
- 90
- 206 933
Fernando Zigunov
Приєднався 13 бер 2010
LxLaser 2024 - One-Shot Omnidirectional Pressure Integration Through Matrix Inversion
This is a talk I gave at LxLaser 2024 on a work I've done at Los Alamos National Lab. If you have interest or any questions about the work, feel free to comment or send me an e-mail!
Referenced papers:
Zigunov & Charonko (2024, ArXiv) "One-shot omnidirectional pressure integration through matrix inversion"
arxiv.org/pdf/2402.09988
Zigunov & Charonko (2024, Measurement Science and Technology) "A fast, matrix-based method to perform omnidirectional pressure integration"
iopscience.iop.org/article/10.1088/1361-6501/ad2da5/meta
DOI 10.1088/1361-6501/ad2da5
Our Github repo with the files for the open-source solver:
github.com/lanl/pressure-osmosis
Referenced papers:
Zigunov & Charonko (2024, ArXiv) "One-shot omnidirectional pressure integration through matrix inversion"
arxiv.org/pdf/2402.09988
Zigunov & Charonko (2024, Measurement Science and Technology) "A fast, matrix-based method to perform omnidirectional pressure integration"
iopscience.iop.org/article/10.1088/1361-6501/ad2da5/meta
DOI 10.1088/1361-6501/ad2da5
Our Github repo with the files for the open-source solver:
github.com/lanl/pressure-osmosis
Переглядів: 135
Відео
Matlab Tutorial for the OSMODI (Pressure from PIV) solver
Переглядів 2113 місяці тому
I made this video going through what is the kind of pre-processing to successfully obtain pressure from PIV with our new One-Shot Matrix Omnidirectional Integration solver (OS-MODI). This is a slow-paced tutorial, so you can catch all details. You can fast-forward me as much as you can =) Timestamps: 0:00 Where to find this compiled code and 'installation' 1:45 2D PIV, Average Pressure using th...
6 hours of jet noise captured at 290,000 fps
Переглядів 4787 місяців тому
Yes, this is the raw footage from the entire memory of our Phantom T2410 camera, because WHY NOT? Jet velocity is 100m/s (air jet at ambient temperature); and the video captures the shadowgraph of the flow (highlighting the Laplacian of density). The video resolution is reduced because we captured it at 290,000 fps. The total length of the video is 2.3 seconds, which stretch into 6.18 hours whe...
AIAA Scitech 2023 - Power Efficiency Analysis of a Co-Flow Jet Airfoil in CruiseConditions
Переглядів 1,4 тис.Рік тому
This is a presentation I gave at AIAA Aviation 2022 / AIAA Scitech 2023. A conference proceeding paper was generated for this work. Once it is published (by Jan 2023), I'll attach the link here. Thanks for watching and let me know if you have any questions about our work!
Scanning Stereoscopic PIV - An easy 3D flow measurement technique without Tomo-PIV!
Переглядів 4,1 тис.Рік тому
This video describes a technique I developed with my friends at Florida State University to overcome the limitations of Tomo-PIV for volumetric measurements of flow fields in supersonic tunnels. Any aerodynamicist that has worked with Tomo-PIV understands it is very challenging to get Tomo-PIV measurements in supersonic flows. This is a workaround I came up with (it was a shower thought!). The ...
The Spectral Proper Orthogonal Decomposition
Переглядів 13 тис.2 роки тому
I made this video in an attempt to popularize the Spectral POD technique. It is an incredibly powerful analysis tool for understanding the data coming from a multitude of sensors. It elevates the Fourier Transform to a whole new level; hence I call it "The Mother of All Fourier Transforms"! I'm linking below the works of the greatest minds in this topic. I attempted to keep the mathematical ter...
The Theory Behind the Acoustic Intensity Probe
Переглядів 1,4 тис.3 роки тому
Acoustic intensity, or sound intensity, is a very useful quantity that enables measurements of sound power of equipment in the field. Here I'm taking a step back and talking about the basic math behind the probe and demonstrating in practice how these measurements are made. Paper referenced: J. Y. Chung - Cross‐spectral method of measuring acoustic intensity without error caused by instrument p...
Understanding POD: the Proper Orthogonal Decomposition
Переглядів 22 тис.3 роки тому
This was a lot of fun to make! 3blue1brown has inspired me a lot to make a math video with cool animations! This is my take on the Proper Orthogonal Decomposition, an incredibly powerful and useful algorithm for data analysis ! I'm a fluid dynamicist by trade, so my examples come from there. Hope you found this useful! This was produced with Manim, a python library developed by 3blue1brown and ...
Supersonic Impinging Jets - Acoustic Resonance
Переглядів 4533 роки тому
This is a video I submitted to the Gallery of Fluid Motion on 2018. Just posting here to have it in my library!
Pressure Sensitive Paint [PSP] - An amazing tool for aerodynamicists!
Переглядів 3,9 тис.3 роки тому
This is a compilation of my learning process as I used Prof. Keisuke Asai's Pressure Sensitive Paint [PSP] formulation to perform aerodynamic experiments in the Slanted Cylinder model at subsonic conditions. Paper coming up! Check out the papers by Prof. Asai's group for details on the working principle, processing and the technique: A review of pressure-sensitive paint for high-speed and unste...
AIAA Scitech 2021 Talk - Experimental optimization of actuator location for active flow control
Переглядів 2373 роки тому
AIAA Scitech 2021 Talk - Experimental optimization of actuator location for active flow control
PhD Defense Talk - Experimental optimization of active flow control
Переглядів 1,1 тис.3 роки тому
PhD Defense Talk - Experimental optimization of active flow control
Experimental Optimization of Microjet Actuator Location for Active Flow Control
Переглядів 4003 роки тому
Experimental Optimization of Microjet Actuator Location for Active Flow Control
Get rid of ALL small bubbles when casting resin with a bit of Fluid Mechanics!
Переглядів 3244 роки тому
Get rid of ALL small bubbles when casting resin with a bit of Fluid Mechanics!
An EPIC print timelapse of the Voronoi bunny
Переглядів 6874 роки тому
An EPIC print timelapse of the Voronoi bunny
Make EPIC panning timelapses with a DIY tripod slider!
Переглядів 1,8 тис.4 роки тому
Make EPIC panning timelapses with a DIY tripod slider!
PIV movie of a NACA0020 airfoil at Re_c=47000
Переглядів 3344 роки тому
PIV movie of a NACA0020 airfoil at Re_c=47000
Visualizing sound! Shadowgraph of an acoustic levitator.
Переглядів 1,4 тис.4 роки тому
Visualizing sound! Shadowgraph of an acoustic levitator.
Assembling an FPGA driven ultrasonic phased array
Переглядів 6434 роки тому
Assembling an FPGA driven ultrasonic phased array
3d printed Reverse Lens Adapter for Macro Photography
Переглядів 1,7 тис.4 роки тому
3d printed Reverse Lens Adapter for Macro Photography
Installing the Inventor Assembly Tools
Переглядів 1,3 тис.4 роки тому
Installing the Inventor Assembly Tools
A little app to see the branch cuts of functions
Переглядів 1184 роки тому
A little app to see the branch cuts of functions
Branch cut of the sqrt(z^2-A) function
Переглядів 1,5 тис.4 роки тому
Branch cut of the sqrt(z^2-A) function
How to adjust the Scheimpflug adaptor in practice
Переглядів 1 тис.5 років тому
How to adjust the Scheimpflug adaptor in practice
In your setup, you have the cameras below looking "up" into the laser, would it work just as well if they were looking "down" the laser?
@@will6699 yes, it would. You may face some issues because the backscatter signal is weaker than the forward scattering signal (i.e particles will look much dimmer). But I've built setups like that and it works; you just have to deal with having less light available.
@fzigunov Thanks for the quick reply! I just started a new job and they want me to start looking into PIV so I'm just trying to figure it out in my head 🤣
Truly exceptional teaching. If only more lecturers had your skills!
Always been skeptical of the derivation of pressure from the fluid velocity via PIV but tbh it's been above my knowledge. This method makes a lot of sense! I'll share with my colleague doing PIV!
So relaxing, I am falling asleep
Explanation=0👎
Amazing video!
Excellent Video and pretty easy and clear to understand the operations of blocks due to visualization!!
Don't remember, when I hit like on a youtube video the last time. Probably years ago. But this one deserves one.
nice! i used a fully digital implementation of this on a microcontroller. it means i can have a microcontroller be sending an approximation of a floating point value (from a pid controller) to drive a relay, but only switching at the cycle period at most :) in this case switching eventually leads to mechanical wear, so having something that switches less when further from 50% is very nice.
Nice. I wanna add those led infinite mirrors on mine to make the basisc portal man.
Sir, the .mex file is missing in Github, and the file has also been removed from Matlab works...
Hi, nurachand, I apologize for the problem. I think my collaborator has not uploaded the compiled files to the Github page yet. We had some delays due to LANL's open-sourcing process. I can share the file via personal communication until this problem is fixed, if you're interested.
@@fzigunov@fzigunov Sure, I sincerely want to understand the Pressure estimation methodology from PIV (velocity) data. I have 2D PIV files in .vec format. How do you think I could communicate with you?
@nurachand, try the link below with the files in the meantime. The .vec files will work if you follow the procedure in the tutorial. If you face issues, please let me know. I can be contacted at fzigunov (at) syr.edu. drive.google.com/file/d/1uh29GmNYVBLysGIssJu8EmvoPPcqA448/view?usp=drive_link
@@fzigunov Thanks , i will try and let you know.
Link to the Github Open-Source code: github.com/3dfernando/pressure-osmosis If you use our code, I'd appreciate you cite us! Zigunov & Charonko (2024, ArXiv) "One-shot omnidirectional pressure integration through matrix inversion" arxiv.org/pdf/2402.09988 Zigunov & Charonko (2024, Measurement Science and Technology) "A fast, matrix-based method to perform omnidirectional pressure integration" iopscience.iop.org/article/10... DOI 10.1088/1361-6501/ad2da5
This is a talk I gave at LxLaser 2024 on a work I've done at Los Alamos National Lab. If you have interest or any questions about the work, feel free to comment or send me an e-mail! Link to the Github Open-Source code: github.com/3dfernando/pressure-osmosis Referenced papers: Zigunov & Charonko (2024, ArXiv) "One-shot omnidirectional pressure integration through matrix inversion" arxiv.org/pdf/2402.09988 Zigunov & Charonko (2024, Measurement Science and Technology) "A fast, matrix-based method to perform omnidirectional pressure integration" iopscience.iop.org/article/10.1088/1361-6501/ad2da5/meta DOI 10.1088/1361-6501/ad2da5
great Explaination Can you please tell me how to calculate the value of Sampling and feedback capacitor in 2nd order SD-ADC 2MHz samplig frequency and NFFT is 8192.
What does it mean by spatial patterns where the data is most temporally correlated? Like the spatial patterns across consecutive time? Or across like periodically when a time step repeats if something is periodic..?
Very good question. If the correlation between some variable X and Y is high (close to +1), then X and Y move together (i.e., when X is up, Y is also up). If correlation is low, (close to -1), we get the opposite effect and when X is high, Y is low. For POD, it calculates from the data given how they are most correlated (i.e., instead of X and Y, you have many X variables, all contained in the matrix A). So if in a mode a hot spot shows in two different regions of an image (on a high rank mode), it is safe to say they have a high correlation (i.e. , they vary in unison). Hope this helps understand a little better.
@@fzigunov Yes thank you!!
Great work!
So all waves in the universe behave the same. Light and darkness, sound and stillness all the same just have different weight. That section where sound expand but I also contracts to something the steady motion it was in... Nothing does coming from something nothing is something.
Good explanation, but just one note. It will be easy for understanding if the final counter will be signed and bidirectional (up/down counter), so when the analog input is near zero and output of the comparator is "+ - + - + - ...." the output of accumulator is about zero too. I know that it is just a matter of signed and unsigned numbers, but anyway :)
silicon detectors in phones can see light out to about 1100nm. They put filters in to stop it on good cameras. Sometimes the filters can't block all the IR
Yeah I think the laser light is so intense that even though the sensor only responds to 0.1% it still is enough to show as a bright spot. Very interesting!
Record a VR headset (one that uses sensor satellites) or even a modern game controller or TV remote with your phones camera. You will see them flashing all over with their IR LEDs
Well you're explaining what is happening operationally (this specific form) but not the why or how it does it's actually function of "noise shaping". That being said the theory would be beyond most.
Great video Fernando. Is there a website that I can view relevant code to the fluid dynamic example you showed? I would love to explore that more. Thanks!
What an amazing explanation. Thank you very much!
Amazing as always
excellent visualization! Very helpful!💖
Great explanation and visualization!
Impressive!
Impressive! How is the sound related to the shown shadowgraphs?
Bernd, it was a very simple thing, really. The time series of pixel intensity is played back at 3000 samples/s. Since the video is played at 30 samples/s, there is a mismatch between the length of the records (the sound is 100x shorter), so I simply "tacked on" more pixels to the end of the sound. Since it's noise, you can't notice the difference. The idea here was simply to show the data in a different manner. Most of the sound played back is ultrasonic, really.
Yes, this is the raw footage from the entire memory of our Phantom T2410 camera, because WHY NOT? Jet velocity is 100m/s (air jet at ambient temperature); and the video captures the shadowgraph of the flow (highlighting the Laplacian of density). The video resolution is reduced because we captured it at 290,000 fps. The total length of the video is 2.3 seconds, which stretch into 6.18 hours when re-played at 30 fps. Compression artifacts are courtesy of UA-cam. I can't do anything about that unfortunately 😐 The noise in the background is actually coming from the video. The intensity of each pixel is used as a small microphone to produce a sound wave that is slowed down to audible frequencies. Multiple pixels are used to make the length of the audio match the video; but it is all procedurally generated. The pictures of the shear layer turbulence are so mesmerizing that I had to share this with the world. Maybe you can use this as a background for study or sleep, who knows? 😂
Crazy
Whaaat
Damn that's insane
Great video and all but man why the 2 girls 1 cup music
Thanks for checking this short! Links to the referred papers: "A continuously scanning spatiotemporal averaging method for obtaining volumetric mean flow measurements with stereoscopic PIV" (Zigunov et al, ExIF, 2023) - link.springer.com/article/10.1007/s00348-023-03596-w "A fast, matrix-based method to perform omnidirectional pressure integration" (Zigunov and Charonko, MST, 2024) - iopscience.iop.org/article/10.1088/1361-6501/ad2da5 "One-shot omnidirectional pressure integration through matrix inversion" (Zigunov and Charonko, Under Review, 2024) arxiv.org/pdf/2402.09988.pdf
Thanks! Wouldve been nicer is you show how an analog signal was converted to digital bits by this ADC.
Thanks for checking our short! Here's our two papers on the method. First paper demonstrates the math integrals; second paper eliminates the iterative part of the method to solve for pressure in a single matrix inversion. A fast, matrix-based method to perform omnidirectional pressure integration arxiv.org/abs/2311.16935 One-shot omnidirectional pressure integration through matrix inversion arxiv.org/abs/2402.09988
Hello Fernando, Interesting technique described here. I was wondering whether HFSB's combined with counter illuminating lasers can solve the problem of intensity drop that you are describing. Ofc you would have to switch to particle tracking and multi-double-pulse techniques such as Shake The Box, to resolve high speed flows. Again, the time average would be captured though...
Hi, Niko, thanks for the comment. HFSBs are brighter and definitely a solution for low-speed flows. I have not seen HFSBs being used in a supersonic flow; and I think that has to do with a few aspects: (1) Evaporation rate is actually very fast - For normal water in room air it is about ~2mm/day or 23nm/s. This means a bubble with a very thin shell will disappear in a really short time. Though I haven't looked closely into whether a linear scaling applies all the way to the nano-scales. (2) If the bubbles are to last a few seconds, they would need to have shells of ~100's of nm; but that's about the size of an oil droplet used in supersonic flows (~1000 nm). So the gain in brightness for a given frequency response is likely not there; (3) It is just in general not something people have done - Generating bubbles at the single-digit micron diameter scale. There are engineering challenges in doing that with any repeatability/volume and making sure the diameter distribution is not too wide. For a typical supersonic tunnel in a university, you'd have to generate particles at a rate of ~1e12 to 1e13/s (our 4x4" M=3 tunnel requires that many). That's already hard with droplet condensation systems; let alone with repeatable bubbles. I think it is a very rich research direction; there's a lot to learn here. It may not be impossible to do; and definitely would have great benefits to image supersonic flows in 3D. Perhaps as a community we can move towards that direction!
what is the difference between this method and SLS(scanning light sheet)?? Just curious, i am preparing an exam about this stuff
Thanks for the comment! Scanning light sheet is moving only the laser, without moving the camera. In that setup you need to require a depth of field of the thickness of the entire volume (not only the light sheet) such that the particles are in focus, meaning the intensity of the particle images is ~(1/thickness)^2 due to this focusing requirement; and the amount of light required is much larger (potentially hundreds or thousands of times more light). For a supersonic flow you'd need kilojoules of laser energy per pulse to do PIV in a large volume; which is just too expensive to have in a fluids lab.
Thanks a lot....This is the best POD explaination I have seen.... Very helpful... Keep it up.
This was really cool! I was wondering if you could tell me how much of each part I need, because I want to replicate this project.
Hi, Bob, I'm glad you found it neat! See the link in my description for the STL files. The only part that needs a qty is the belt links, which are approximately 45 for this conveyor. The other parts can be visually counted (it's either 1, 2 or 4 of each. I no longer have the original design files, it's been a while I built this!
@@fzigunov Thanks so much!
Details on the paper: link.springer.com/article/10.1007/s00348-023-03596-w
Hey nice video man! you always surprise me!
Thanks!! I saw your talk at APS on the lobed nozzle, very nice work!
Good description of the basics. I would like to add that Nyquist refers to frequency and not phase (in recreating a representation of the input) if the sampler is asynchronous to the process, which is usually the case. You can prove this to yourself by making a small circuit with a sample and hold and slowly bring up the frequency of the input in relation to the sample frequency. At a sample rate of 1/2 the highest frequency (which is hard to get correct and will wonder in relation to a tightly regulated sample clock), asynchronously the samples are almost void of phase information. So, from a practical standpoint the first task is to preclude aliasing which means a low pass filter in relation to the sample period, with good process signal attenuation of around 40dB minimum. Then in designing the controller you would like the sample rate, to maintain phase, of at best 20x the highest process frequency. Many industrial off-the-shelf PID controllers provide this 20x sample rate, many call it oversampling as a marketing point. Now with this over sampling and elimination of aliasing we can design and tune a classic control law in the Z domain. If there are multiple feedback loops we must have synchronous sampling of all required process variables. From here we can use the classical control system stability criteria.
you're right, very cringe humor lmao
Man this animation was incredible, thankyou!
Great video. Do you need time resolved PIV measurements to conduct spectral POD? Thanks
Hi, Honcho, you're correct, you need time-resolved data (PIV, PSP, Shadowgraph, etc.) to perform SPOD
THanks!
This video is gold.
Very good visualization of the interpretation of the V matrix. I just needed it, thank you
THIS IS AMAZING PLEASE KEEP DOING THIS