commit 2b6e5f7 Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 8 21:31:47 2023 -0500 Riser update Calculate fluid volume fraction using IB VF commit 9a4ccf4 Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 8 17:48:24 2023 -0500 Riser case Added the riser case commit 9217d19 Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 8 17:37:46 2023 -0500 Sync w/ Olivier Sync with Olivier's develop branch commit 9d83c8e Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 8 17:36:10 2023 -0500 Squashed commit of the following: commit 0bca3e8 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Feb 7 23:19:41 2023 -0500 Added IB-particle collision from Jesse. Switched all solvers (except MAST!) to new linsol class. ALMOST ALL EXAMPLES WILL NEED TO BE UPDATED! commit 26b4fe2 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Feb 7 22:29:02 2023 -0500 Switched fourier3d to c2c... commit 19ae451 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Feb 7 22:09:34 2023 -0500 Commit what I think is the cleanest implementation of r2r fftw solver - but that requires another strategy for building the operator commit a9becb0 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Feb 7 16:22:59 2023 -0500 p3dfft is also removed as decomposition strategy commit 31219ce Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Feb 7 16:16:07 2023 -0500 May not be working yet - this commit removes p3dfft altogether, replaces it with our in-house fourier solver. commit ff795e1 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Feb 7 09:49:15 2023 -0500 Started adding a multiphase_pipe case. Started adding back in-house fourier solver. Fixed bug with hit. commit ec57c3a Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Mon Feb 6 13:41:15 2023 -0500 Proper non-dimensionalization of the swirl injection case commit 49d5325 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Mon Feb 6 13:17:20 2023 -0500 Moved library directory specification to my .zshrc to remove all reference to 'desjardi'... commit 7937d14 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Mon Feb 6 10:59:56 2023 -0500 Adding a swirl atomizer test case. commit c70dd8b Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Sun Feb 5 21:51:08 2023 -0500 Particles with id=0 do not modify the flow and are not advanced. They're just kinda... there. commit 305694c Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Sun Feb 5 21:35:21 2023 -0500 Removed shock-particle case, it's not ready. commit 4bf8f61 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Sun Feb 5 21:34:04 2023 -0500 Updated shit_p routines to make use of newly defined integration tools and fluid_vol variable. commit af17f44 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Sun Feb 5 20:01:42 2023 -0500 Absorded range of modifications from Jesse's repository. Most important is the use of p3dfft for fourier-based linear solvers. commit f962bc2 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Fri Feb 3 17:17:33 2023 -0500 Clean-up ib_pipe to make it a proper pipe flow LES. Cleaned up Vreman's sgs model. commit 6c4766c Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Fri Feb 3 14:19:32 2023 -0500 Slight clean-up of the turbulent ib_pipe case. commit bf6f539 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Feb 2 22:57:49 2023 -0500 ib_pipe modified to absorb Jesse's changes, and pushed to turbulent conditions. commit d789220 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Feb 2 14:56:07 2023 -0500 Small changes to falling_drop, added new ib_drop case. commit 1ceb7d2 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Wed Feb 1 10:41:13 2023 -0500 Reindent df_class. Need to think about role of rho here for two-phase flows. commit 8461ee6 Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 8 15:03:37 2023 -0500 remove p3dfft decomp Thought I removed this commit 3be151f Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 8 11:29:26 2023 -0500 Update GNUmakefile Remove p3dfft commit 5b56d8f Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 8 08:52:48 2023 -0500 LPT updates Remove solver specific flags. Added optional argument to avoid overlap when injection routine is called. Also fixed a bug with LPT cfl, before the LPT cfl was overwriting flow solver CFL. Now we have option to use collisional CFL and then check max. commit ee68f90 Author: jessecaps <jcaps@umich.edu> Date: Tue Feb 7 14:09:04 2023 -0500 IB collision + linsol integration Replace pipe collision detection with more general (but expensive) IB detection. Also brought linsol into the variable density solvers. commit cac0ade Author: jessecaps <jcaps@umich.edu> Date: Tue Feb 7 14:01:04 2023 -0500 nothing important Back to Olivier's version commit be5f1f9 Author: jessecaps <jcaps@umich.edu> Date: Tue Feb 7 13:47:39 2023 -0500 In-house P3DFFT NGA2 now runs without external dependency for parallel fft library. Slightly slower than using a library, but can likely be optimized commit 5c4d232 Author: jessecaps <jcaps@umich.edu> Date: Sat Feb 4 11:09:40 2023 -0500 Update simulation.f90 monitor the correct mfr each time step commit b7dbc2b Author: jessecaps <jcaps@umich.edu> Date: Sat Feb 4 11:05:08 2023 -0500 P3DFFT allows for general decomposition Implemented P3DFFT compatible parallel decomposition, can now run FFT pressure solver w/ nprocs that don't evenly divide into the number of grid points! For IB Pipe need to add epsilon to cfg%VF, pressure solver doesn't like zero volumes commit 9e56801 Author: jessecaps <jcaps@umich.edu> Date: Fri Feb 3 17:44:31 2023 -0500 Merge branch 'develop' of https://github.com/desjardi/NGA2 into develop commit 3278fd0 Author: jessecaps <jcaps@umich.edu> Date: Thu Feb 2 20:36:28 2023 -0500 Update geometry.f90 I should compile being pushing commit 7aefa45 Author: jessecaps <jcaps@umich.edu> Date: Thu Feb 2 17:16:10 2023 -0500 Update geometry.f90 How did my compiler miss this? commit ef24ad8 Author: jessecaps <jcaps@umich.edu> Date: Thu Feb 2 10:04:56 2023 -0500 Update pfft3d_class.f90 Because Jack insists commit 5842d2e Author: jessecaps <jcaps@umich.edu> Date: Thu Feb 2 09:58:44 2023 -0500 Update pfft3d_class.f90 Less restrictive check, the old check failed when domain length had decimal value commit 2aa1478 Author: jessecaps <jcaps@umich.edu> Date: Thu Feb 2 07:37:32 2023 -0500 Particle-laden pipe flow is working We now have clustering! Added ability to detect collisions with cylinder walls. commit 917ece4 Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 1 14:18:07 2023 -0500 Update pfft3d_class.f90 This version works in 2D + 3D, still issue with 2D when npx=1 commit ec09e9f Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 1 13:08:21 2023 -0500 Update pfft3d_class.f90 3D FFT works in many cases, in parallel num procs must be evenly divisible into nx. Also seems to be a bug in the p3dfft library when nz=1 and npx=1. commit 2d5aba0 Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 1 13:07:14 2023 -0500 Update poissontest Add consistency with NGA solver commit 0750caf Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 1 09:05:18 2023 -0500 Fix conflicts with merge Synced with Olivier's develop branch commit 1247ef2 Merge: 73720a6 e21b0b1 Author: jessecaps <jcaps@umich.edu> Date: Wed Feb 1 08:57:43 2023 -0500 Merge branch 'develop' of https://github.com/desjardi/NGA2 into develop commit 73720a6 Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 31 17:00:46 2023 -0500 Update pfft3d_class.f90 3D FFT is now working in serial and in parallel (tested on the HIT case)! commit e21b0b1 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Jan 31 16:17:16 2023 -0500 Replaced itpi by itpr in ib_tester. commit 48beaf0 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Jan 31 16:14:27 2023 -0500 Cleaned up gradU, strainrate, vorticity to make all solvers more similar. Also removed itpi from incomp and replaced by itpr. commit e9773ce Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 31 14:55:18 2023 -0500 Update pfft3d_class.f90 Serial FFTW implemented and verified it works. commit 3b1639b Author: jessecaps <jcaps@umich.edu> Date: Mon Jan 30 13:56:56 2023 -0500 p3dfft merge Merge Jack's P3DFFT updates commit a8f65c1 Merge: 766c357 ccae70d Author: jessecaps <jcaps@umich.edu> Date: Mon Jan 30 13:56:32 2023 -0500 Merge branch 'p3dfft' into develop commit 766c357 Author: jessecaps <jcaps@umich.edu> Date: Sat Jan 28 10:16:42 2023 -0500 Update lpt_class.f90 Do not include particles w/ id=0 when computing VF commit ccae70d Author: John Wakefield <jwake@umich.edu> Date: Fri Jan 27 21:30:23 2023 -0500 fixed scaling commit 3afc9c9 Author: jessecaps <jcaps@umich.edu> Date: Fri Jan 27 18:48:44 2023 -0500 Update to pipe cases Geometry now computes volume fraction field for IBM, useful for integration to take into account correct volume commit 5781f18 Author: John Wakefield <jwake@umich.edu> Date: Fri Jan 27 13:09:51 2023 -0500 added Poisson equation with point sources test case for p3dfft solver commit 475fade Author: John Wakefield <jwake@umich.edu> Date: Fri Jan 27 13:09:15 2023 -0500 fixed off by one error commit 2fa5f3a Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Fri Jan 27 12:13:08 2023 -0500 Removed extra line from mast commit eebe39d Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Jan 26 17:29:07 2023 -0500 Removed redundant subroutine in low Mach solver due to incorrect conflict resolution. commit 0dc01ce Author: jessecaps <jcaps@umich.edu> Date: Thu Jan 26 14:07:58 2023 -0500 Update Make.p3dfft Not needed commit 3060d19 Author: jessecaps <jcaps@umich.edu> Date: Thu Jan 26 13:05:51 2023 -0500 Update Make.p3dfft Need to point to correct directories commit 0a0a86e Author: jessecaps <jcaps@umich.edu> Date: Thu Jan 26 12:04:22 2023 -0500 fixing merge issues commit 82f803f Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 25 22:26:45 2023 -0500 Update incomp_class.f90 Added itpi_x operator needed by IBM commit 90ca7e4 Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 25 21:43:07 2023 -0500 Remove old example case commit 30016ba Merge: f805f3d 53a06cf Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 25 21:36:34 2023 -0500 Merge branch 'p3dfft' into develop commit f805f3d Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 25 18:16:31 2023 -0500 Resolve issues from previous merge Resolve issues from previous merge commit dd5fb99 Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 25 18:07:45 2023 -0500 Update lowmach_class.f90 cosmetic commit 092cee5 Merge: 8a3d5a9 397c154 Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 25 18:07:13 2023 -0500 Merge branch 'develop' of https://github.com/desjardi/NGA2 into develop commit 8a3d5a9 Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 25 17:56:14 2023 -0500 New riser case New riser case that solves particles in a periodic pipe. Uses a simple IBM for pipe and body forcing. Includes restart capabilities and wallclock timing. TODO: LPT collisions with pipe walls commit 53a06cf Author: John Wakefield <jwake@umich.edu> Date: Wed Jan 25 16:33:23 2023 -0500 fixed processor-local indices in fourier space commit 397c154 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Wed Jan 25 16:32:08 2023 -0500 Squashed commit of the following: commit 49fc8e8 Author: jessecaps <jcaps@umich.edu> Date: Mon Jan 9 15:46:23 2023 -0500 Update simulation.f90 Avoid division by 0 commit 0e80a30 Author: jessecaps <jcaps@umich.edu> Date: Sun Jan 8 11:56:35 2023 -0500 Part inject Particle injection example commit 8e5df3b Author: jessecaps <jcaps@umich.edu> Date: Sun Jan 8 11:53:24 2023 -0500 Update sgsmodel_class.f90 Separate out different SGS models. Currently includes constant coeff Smagornisky, Vreman, and dynamic smag. Dynamic smag still has a bug. commit 4d62407 Author: jessecaps <jcaps@umich.edu> Date: Sun Jan 8 11:51:15 2023 -0500 Direct forcing solver added Direct forcing IBM solver with example. Still needs updates to allow for moving IBs but that is not difficult. commit ec6fad5 Author: jessecaps <jcaps@umich.edu> Date: Sun Jan 8 11:45:05 2023 -0500 Channel case added Olivier's channel case (without non-Newtonian visc part) commit b533e73 Author: jessecaps <jcaps@umich.edu> Date: Sun Jan 8 11:43:23 2023 -0500 Simple IB pipe case Simple IB pipe case. Does not rely on an IBM solver. commit 42c6afb Author: jessecaps <jcaps@umich.edu> Date: Sun Jan 8 11:39:51 2023 -0500 Port new LPT routines to Master Fluidized bed example contains up-to-date changes to LPT. Implicit filtering works, particles feel fluid stress, LPT CFL bug fix, low Mach class has useful routines commit 19cc1c1 Author: jessecaps <jcaps@umich.edu> Date: Sun Jan 8 11:27:29 2023 -0500 Delete matm_class.f90 commit 4849175 Author: jessecaps <jcaps@umich.edu> Date: Sun Dec 18 13:00:37 2022 -0500 LPT updates More progress porting old LPT stuff into lpt_class. Particles now have angular velocity and tangential collisions. By default friction coeff is set to 0, if it is given a value larger than 0 then torque due to collisions is accounted for. Also started making progress to account for lift and fluid torque due to drag. This requires interpolating fluid vorticity. commit 248792c Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 17 10:44:56 2022 -0500 new LPT CFL Added CFLs to LPT, if collisions are used this prevents particles to move more than 1/10 their diameter, otherwise they cannot move more than dx commit bb645e1 Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 17 00:58:08 2022 -0500 Update lpt_class.f90 I think there was a bug in how Olivier implemented the spring stiffness. Updated version is consistent with old NGA and (I think) is the correct version commit 6260927 Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 17 00:25:51 2022 -0500 Update simulation.f90 Particles leave the domain if they are too soft. Lowering Tcol helps commit dc0922e Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 17 00:17:30 2022 -0500 LPT updates Keep track of number of particles added and removed each timestep commit 677da16 Author: jessecaps <jcaps@umich.edu> Date: Fri Dec 16 23:01:25 2022 -0500 New LPT capabilities Added injection routine. This allows for particles to be injected randomly within a circular region. Special care is taken to avoid overlap if collisions are used. Also added separate coeff of restitution for walls. The example case injects particles with diameters sampled from a random lognormal distribution. commit b0a321e Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 10 12:34:36 2022 -0500 diag solver updates The diagonal solver compiles and computes things but it is not working yet. haven't yet switched to the abstract linear solver commit 8f9a8cf Author: jessecaps <jcaps@umich.edu> Date: Fri Dec 9 11:27:45 2022 -0500 New diagonal solver class Added a diagonal solver class, taken from old NGA. Haven't yet tested. Will test on implicit filter in lpt class commit ab391b2 Author: jessecaps <jcaps@umich.edu> Date: Wed Dec 7 15:31:42 2022 -0500 Output part vector to EnSight Added vectors to particle mesh following Max's suggestions commit c2aebbf Author: jessecaps <jcaps@umich.edu> Date: Wed Dec 7 08:50:43 2022 -0500 Updates to fluidized bed Still working on it commit 55eeef5 Author: jessecaps <jcaps@umich.edu> Date: Wed Dec 7 08:49:45 2022 -0500 Update .gitignore We should ignore makefiles. Doesn't seem to be working though.... commit be802f8 Author: jessecaps <jcaps@umich.edu> Date: Wed Dec 7 08:44:19 2022 -0500 fbed input Updates to fluidized bed input commit 7305337 Author: jessecaps <jcaps@umich.edu> Date: Wed Dec 7 08:42:42 2022 -0500 Update lpt_class.f90 Working on fluidized bed example. Started adding ability to do implicit filtering of projected terms. Still need to interpolate fluid stress to particle location. Case does not run using 3 or 4 procs in x, but runs with 1, 3, and 5... commit a311987 Author: jessecaps <jcaps@umich.edu> Date: Mon Dec 5 14:45:42 2022 -0500 Delete .#lpt_class.f90 remove garbage file commit a769c4e Merge: b674f28 2126f7a Author: jessecaps <80769669+jessecaps@users.noreply.github.com> Date: Mon Dec 5 20:31:32 2022 +0100 Merge pull request #1 from rgrawe/develop Develop commit 2126f7a Merge: c6d27d3 b674f28 Author: rgrawe <rgrawe@umich.edu> Date: Mon Dec 5 12:13:06 2022 -0500 Merge branch 'jcaps' of https://github.com/jessecaps/NGA2 into jcaps commit c6d27d3 Author: rgrawe <rgrawe@umich.edu> Date: Mon Dec 5 12:02:17 2022 -0500 Fluidized Bed Initialization With Grid Problems Attempting to initialize fluidized bed but there's something wrong with the grid. commit b674f28 Author: jessecaps <jcaps@umich.edu> Date: Wed Nov 30 15:06:21 2022 -0500 Update simulation.f90 Need to specify drag model in HIT LPT initialization commit d9b56f6 Author: jessecaps <jcaps@umich.edu> Date: Wed Nov 30 14:17:23 2022 -0500 Update lpt_class.f90 Drag correction was being canceled out commit c7e8178 Author: jessecaps <jcaps@umich.edu> Date: Tue Nov 29 14:00:10 2022 -0500 Update lpt_class.f90 Added option to interpolate temperature to particle. Also made drag model public. commit 09a1e4f Author: jessecaps <jcaps@umich.edu> Date: Thu Nov 17 22:21:55 2022 -0600 Added HIT from new_linsol branch commit adaf88f Author: jessecaps <jcaps@umich.edu> Date: Mon Nov 7 15:59:50 2022 +0100 Update to MAST Separate where/when viz arrays are computed in MAST commit d72f61d Author: jessecaps <jcaps@umich.edu> Date: Mon Nov 7 15:59:07 2022 +0100 Fixed Tenneti correlation Terms in Tenneti should all be multiplied by extra fVF commit c334b71 Author: jessecaps <jcaps@umich.edu> Date: Mon Nov 7 15:17:44 2022 +0100 Particle mod Starting to create ability to choose different drag models commit 714298c Author: jessecaps <jcaps@umich.edu> Date: Mon Nov 7 15:17:08 2022 +0100 Update to MAST Compute various terms that are useful for visualization commit 23fd21b Author: jessecaps <jcaps@umich.edu> Date: Mon Nov 7 14:10:18 2022 +0100 Example MAST cases from Olivier's develop branch Example MAST cases from Olivier's develop branch commit ecbfe62 Author: Jesse Capecelatro <> Date: Mon Nov 7 13:14:47 2022 +0100 Some updates to Makefiles commit 6aed7db Author: Jesse Capecelatro <> Date: Mon Nov 7 12:59:20 2022 +0100 Shock-particle example case commit f808d59 Author: Jesse Capecelatro <> Date: Mon Nov 7 12:59:02 2022 +0100 Olivier's changes from develop. Also updated LPT class to include energy source, temperature, etc. Still need to add heat transfer between phases. MAST now gets source from particles. Need to check energy equation. commit b0d6e7b Author: John Wakefield <jwake@umich.edu> Date: Wed Jan 25 14:30:58 2023 -0500 fixed ref_stencil for non-root processes commit 30c26b7 Author: John Wakefield <jwake@umich.edu> Date: Wed Jan 25 12:07:52 2023 -0500 current (not working) p3dfft class; same error as library example when run with valgrind commit dda8479 Author: Fabien Evrard <fa.evrard@gmail.com> Date: Wed Jan 25 11:21:51 2023 -0500 Patch that fixes compilation using GNU compilers commit c8761db Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 24 14:18:38 2023 -0500 Update lpt_class.f90 Make stress and src terms in LPT optional commit 350d467 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Mon Jan 23 10:16:23 2023 -0500 Reset falling_drop case to original commit 20e7136 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Mon Jan 23 10:08:15 2023 -0500 Renamed some dev files commit 642c411 Merge: de6f407 8ec29c4 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Jan 17 15:08:59 2023 -0500 Merge pull request desjardi#27 from mbkuhn/levitating_drop commit 8ec29c4 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Jan 16 21:32:04 2023 -0700 working w/ stretched mesh, added print statements commit 514de66 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Jan 16 21:17:38 2023 -0700 working version with uniform mesh commit a6a05af Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Jan 16 12:22:57 2023 -0700 compiling with forcing commit 4a5c64d Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Jan 16 11:37:47 2023 -0700 compiling version - without forcing commit 2224e25 Author: jessecaps <jcaps@umich.edu> Date: Sat Jan 14 14:09:52 2023 -0500 Remove case More work needs to be done before we can simulate shock-particle interactions in NGA2 commit 8420943 Author: jessecaps <jcaps@umich.edu> Date: Fri Jan 13 12:57:56 2023 -0500 Update diag_class.f90 Added buffers to avoid warning messages printed to the screen commit aad15ec Merge: 0616fb5 49fc8e8 Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 10 10:20:07 2023 -0500 Merge branch 'master' into develop commit 0616fb5 Author: jessecaps <jcaps@umich.edu> Date: Sat Jan 7 19:08:56 2023 -0500 Create GNUmakefile Makefile for ib_pipe commit 9ec1db7 Author: jessecaps <jcaps@umich.edu> Date: Sat Jan 7 18:43:41 2023 -0500 Channel update Channel compiles using new lin_sol stuff. commit 4a52436 Author: jessecaps <jcaps@umich.edu> Date: Sat Jan 7 18:30:31 2023 -0500 channel case Added Olivier's channel case (without non-Newtonian visc). Modified forcing in IB_pipe to be consistent. commit ced0d79 Author: jessecaps <jcaps@umich.edu> Date: Fri Jan 6 16:40:49 2023 -0500 Update lowmach_class.f90 Added routine to calculate fluid vorticity. This will be needed when using lift in LPT commit 644f8b8 Author: jessecaps <jcaps@umich.edu> Date: Thu Jan 5 11:24:40 2023 -0500 Updates to SGS Max's changes separating out different SGS models and including Vreman model. commit 7a2b63e Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 4 20:56:25 2023 -0500 Delete matm_class.f90 commit 5780cd5 Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 4 16:49:53 2023 -0500 LPT updates I prefer to have src terms as optional arguments in lpt_advance, this avoids us having to compute these terms if one-way coupled, and also allows us to eventually include sources for energy. Also reduced the inlet velocity in the fluidized bed test case to get proper bubbling. commit cd0cbd4 Author: jessecaps <jcaps@umich.edu> Date: Wed Jan 4 11:00:11 2023 -0500 Update lpt_class.f90 Nothing important. commit 2f95259 Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 3 22:30:45 2023 -0500 LPT implicit filter now works and doesn't print warnings Store recbuf in buf to avoid annoying warnings. Also by default the implicit filter should be true, it is cheaper. commit 319148b Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 3 21:17:12 2023 -0500 Update lpt_class.f90 Fixed a bug with the implicit filter. It is now working! commit 65dd294 Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 3 18:34:23 2023 -0500 Update simulation.f90 Time velocity and pressure separately. This can probably be written cleaner commit 4b3f51c Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 3 18:05:24 2023 -0500 Update simulation.f90 Slightly cleaner way to handle timing commit 042b250 Author: jessecaps <jcaps@umich.edu> Date: Tue Jan 3 17:41:23 2023 -0500 Update simulation.f90 Added timers to fluidized bed case, might be handy later on commit b78a1f3 Author: jessecaps <jcaps@umich.edu> Date: Mon Jan 2 17:25:11 2023 -0500 Update GNUmakefile part_inject does not require IRL commit 0328a00 Author: jessecaps <jcaps@umich.edu> Date: Mon Jan 2 17:03:51 2023 -0500 Part_inject makefil Makefile for part.inject requires updating .gitignore commit 12523eb Author: jessecaps <jcaps@umich.edu> Date: Sun Jan 1 15:32:40 2023 -0500 Update README proper reference commit 85bd226 Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 31 18:41:44 2022 -0500 ib_pipe updates Small edits to ib pipe case commit c4645da Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 31 14:36:46 2022 -0500 Add fluctuations to pipe Basic idea, this can be improved commit 55478fd Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 31 14:27:07 2022 -0500 Simple IB pipe case A very simple implementation of pipe flow using a direct forcing-type scheme. commit 5705bd5 Author: jessecaps <jcaps@umich.edu> Date: Sat Dec 31 14:25:27 2022 -0500 Rename variable To be consistent with other classes commit 057fa6c Author: jessecaps <jcaps@umich.edu> Date: Thu Dec 29 22:58:47 2022 -0500 Update df_class.f90 Calculate center of mass and volume of each IBM object based on marker location and size. Does not work when object crosses a periodic boundary! commit 90577ff Author: jessecaps <jcaps@umich.edu> Date: Thu Dec 29 16:27:58 2022 -0500 Update simulation.f90 Important to sync U after adding IBM correction commit ab9668a Author: jessecaps <jcaps@umich.edu> Date: Sun Dec 25 21:58:05 2022 -0500 LPT and DF update CFLs for particles should be based on dx not dxm commit 9d3b0be Author: jessecaps <jcaps@umich.edu> Date: Sun Dec 25 12:35:37 2022 -0500 IBM update Added force computation acting on each 'object' and changes interpolation subroutine to be a function commit 620e9f7 Author: jessecaps <jcaps@umich.edu> Date: Sun Dec 25 09:46:26 2022 -0500 Update df_class.f90 cleanup commit ea4b12b Author: jessecaps <jcaps@umich.edu> Date: Fri Dec 23 16:37:35 2022 -0500 Update df_class.f90 Added levelset array to direct forcing class. This is only used for visualization and eventually collision detection with LPT commit bdba18b Author: jessecaps <jcaps@umich.edu> Date: Fri Dec 23 16:36:59 2022 -0500 Update mathtools.f90 Added safe arctan function commit 563b1d0 Author: jessecaps <jcaps@umich.edu> Date: Fri Dec 23 16:03:35 2022 -0500 direct forcing update Make ID integer not integer8, unlikely we'll have THAT many IBM objects commit d6d3e64 Author: jessecaps <jcaps@umich.edu> Date: Fri Dec 23 12:46:03 2022 -0500 New direct forcing IBM solver Created a simple immersed boundary solver based on direct forcing. The marker particle ID determines the number of objects. At the moment objects cannot move, but will eventually add an advance routine to update objects. commit e3d440f Author: jessecaps <jcaps@umich.edu> Date: Tue Dec 20 20:16:42 2022 -0500 Fluidize bed works! Collisions need to be properly resolved otherwise everything goes crazy... Added LPT CFL and some other small things. This case now bubbles commit 4cf4506 Author: jessecaps <jcaps@umich.edu> Date: Tue Dec 20 16:51:00 2022 -0500 Remove print statement Remove print statement commit cf9d4a1 Author: jessecaps <jcaps@umich.edu> Date: Tue Dec 20 13:34:29 2022 -0500 More LPT changes LPT now slips over particles with id<=0 to allow for "glued" particles. Renamed p%col to p%Acol (acceleration due to collisions). Also LPT now reads in fluid stress and interpolates to particle position. commit 9180874 Author: jessecaps <jcaps@umich.edu> Date: Tue Dec 20 13:07:58 2022 -0500 Update lowmach_class.f90 Added routine that computes divergence of fluid stress. This will be needed by LPT to interpolate stress to particle location. commit de6f407 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Fri Dec 16 12:10:58 2022 +0100 Volume, mass, momentum conserving version of tpns. Should be KE conserving spatially too. No stabilization here. commit 802a47a Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Dec 8 17:34:35 2022 -0500 1- started work on variable density jet case 2- rewrote strain rate calculation for compactness 3- falling drop case implements a new momentum conservation strategy - still debugging this! commit 70a97a2 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Wed Nov 16 17:07:53 2022 -0500 Added some basic postprocessing to filmKH case commit c21b48a Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Nov 10 12:59:44 2022 -0500 New filmKH case, modifications to falling drop case commit 452ae98 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Oct 11 16:15:23 2022 -0400 Minor cosmetic changes. commit 91a5f82 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Oct 11 10:50:11 2022 -0400 Some basic directory changes - more importantly, replaced lstdc++ (deprecated) by lc++ commit 8d45c16 Merge: 08751be 5a80cba Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Oct 4 10:36:19 2022 -0400 Merge pull request desjardi#23 from mbkuhn/LJSC commit 08751be Merge: c89f3b8 90be61b Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Oct 4 10:27:34 2022 -0400 Merge branch 'pr/mbkuhn/22' into develop commit 5a80cba Author: Michael Kuhn <mbk225@cornell.edu> Date: Sat Sep 10 20:14:13 2022 -0600 correcting the specification of the jet BC locator. fixes the incorrect pressure solution near the base of the jet. commit c89f3b8 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Sep 8 18:04:37 2022 -0400 In the process of adding bquick commit 8af422c Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Wed Aug 31 16:45:51 2022 -0400 New grad(U) implementation based on face velocities. commit f56e6c7 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Aug 29 23:00:16 2022 -0600 Setup for liquid jet(s) in supersonic crossflow commit 90be61b Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Aug 29 22:47:42 2022 -0600 simplifying how material models are referred, allowing fewer parameters commit 02c8156 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Aug 29 22:45:09 2022 -0600 removing some left-behind code in the shock_water_cylinder case and adjusting to new BC system commit dfcaadd Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Aug 29 22:43:50 2022 -0600 lots of changes to how boundary conditions work - heavy reliance on the shift array, alongside integers that specify the correct part of the array to use, based on the BC - ad-hoc BCs for VF and face pressures - BCs can be specified by face or by cell - fixing some oversights in the treatment of wall cells commit 7f45d85 Author: Michael Kuhn <mbk225@cornell.edu> Date: Sat Aug 20 23:48:26 2022 -0600 additional changes to fully set up the cell-centered 'c' option for iterator class commit 783742b Author: Michael Kuhn <mbk225@cornell.edu> Date: Fri Aug 12 21:22:07 2022 -0600 neumann for volume fraction at outlet. eliminates NaN issue at boundary when droplet is leaving domain. commit 04c0a12 Author: Michael Kuhn <mbk225@cornell.edu> Date: Fri Aug 12 21:21:13 2022 -0600 setting up cell-centered option face='c' for iterator class. added to cell-centered solvers (volume fraction and scalars) when adding boundary conditions. commit 7f046c8 Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Aug 10 22:04:00 2022 -0600 in defining SL flag, +1/-1 in loop over domain doesn't work well for non-periodic boundaries. commit 47082c2 Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Aug 10 22:02:26 2022 -0600 changing solver for momentum equations to get rid of asymmetry commit d8a7671 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Aug 8 22:58:25 2022 -0600 correcting shock sensor definition commit 110cbe2 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Aug 8 22:02:59 2022 -0600 goes with previous commit commit 083b963 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Aug 8 22:00:33 2022 -0600 better monitors and specifying pressure relaxation model in examples commit b6b7a5c Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Aug 4 13:42:00 2022 -0400 Increased time step size for channel case. commit fad6ee3 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Aug 4 13:37:35 2022 -0400 Fixed enormous strain rate bug... Added non-Newtonian model to channel case. commit bf7ee20 Author: Michael Kuhn <mbk225@cornell.edu> Date: Thu Jul 21 23:08:20 2022 -0600 now, relaxation models are actually required commit 8cf295e Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Jul 18 23:49:47 2022 -0600 removind unused declarations commit 26aeb75 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Jul 18 23:43:04 2022 -0600 relaxation models can be specified and are now required commit a67f234 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Jul 18 23:24:16 2022 -0600 thermo-mechanical relaxation is working commit b9c86e4 Author: Michael Kuhn <mbk225@cornell.edu> Date: Sun Jul 17 23:06:58 2022 -0600 changes to locally update Temperature where it should be updated: in pressure relax and energy fix commit 32ccda1 Author: Michael Kuhn <mbk225@cornell.edu> Date: Sun Jul 17 18:50:11 2022 -0600 output non-dimensional numbers upon initialization commit cddbc21 Author: Michael Kuhn <mbk225@cornell.edu> Date: Sun Jul 17 17:32:12 2022 -0600 spurious currents input file - results look good commit 0267bad Author: Michael Kuhn <mbk225@cornell.edu> Date: Sat Jul 16 17:05:14 2022 -0600 gravity implementation improved - better energy terms - picking implicit solver that works for this input file - fixing level_set_drop_center - adding prompt to read gravity vector commit d20fa20 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Jul 14 16:50:21 2022 -0400 Turbulent channel works okay commit 409c9bb Merge: 06b2cf8 900a96e Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Thu Jul 14 09:23:55 2022 -0400 Merge pull request desjardi#21 from mbkuhn/viscous_explicit commit 06b2cf8 Merge: c614fc9 96b8e02 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Wed Jul 13 11:08:06 2022 -0400 Merge pull request desjardi#20 from austinhan/develop commit c614fc9 Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Wed Jul 13 10:46:24 2022 -0400 Implementation of a channel case. Could use more refined forcing and better initial conditions. commit 900a96e Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Jul 13 08:41:45 2022 -0600 removing duplicate of elvira reconstruction commit fdd3129 Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Jul 13 07:52:20 2022 -0600 changing tolerance for inertial case and ensight interval for thermal cylinder commit a4438d2 Author: Michael Kuhn <mbk225@cornell.edu> Date: Tue Jul 12 21:59:03 2022 -0600 modifying ensight output to be more useful commit 802e0fc Author: Michael Kuhn <mbk225@cornell.edu> Date: Tue Jul 12 21:57:57 2022 -0600 adding convergence monitor and tweaking simulation monitor to include max temperature commit 669a372 Author: Michael Kuhn <mbk225@cornell.edu> Date: Tue Jul 12 19:37:14 2022 -0600 small fixes to setup of case and communicating bulk modulus for the sake of sl_face access. commit 50e3fe6 Author: Michael Kuhn <mbk225@cornell.edu> Date: Sun Jul 10 23:05:10 2022 -0600 important update to the semi-Lagrangian flux flag (formerly sl_face). Now broken into 3 arrays, which allow it to properly be communicated using the sync function. commit ad1e3e8 Author: Michael Kuhn <mbk225@cornell.edu> Date: Sun Jul 10 22:02:20 2022 -0600 adding monitors and fixing some loop indices commit fdffbb5 Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Jul 6 22:13:26 2022 -0600 removing some redundant code commit c830294 Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Jul 6 22:12:53 2022 -0600 output functionality for temperature and density profiles, and an important fix for apply_bcond so things still work for fully periodic simulations. still in the middle of debugging and validating. commit d05b9c2 Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Jun 29 23:49:46 2022 -0600 setting up MAST_drop case - intended for testing a variety of things. commit fa9a46f Author: Michael Kuhn <mbk225@cornell.edu> Date: Thu Jun 16 23:17:18 2022 -0600 calculate the viscous CFL, and use the material model to get the largest local kinematic viscosity if it is included as an argument. commit 0f36883 Author: Michael Kuhn <mbk225@cornell.edu> Date: Thu Jun 16 22:17:55 2022 -0600 changing solvers for better convergence and incorporating a kinetic_energy monitor file for the taylor vortex commit c2fe36e Author: Michael Kuhn <mbk225@cornell.edu> Date: Thu Jun 16 21:28:09 2022 -0600 correcting indices in the viscous calculation commit fb61d49 Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Jun 15 23:23:43 2022 -0600 Taylor-vortex test case commit 6bad5dd Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Jun 15 23:23:11 2022 -0600 fixed a big typo! this was really affecting all 2D simulations that involved y. commit 6b7174e Author: Michael Kuhn <mbk225@cornell.edu> Date: Tue Jun 14 23:33:50 2022 -0600 further changes for clarity commit f9542fa Author: Michael Kuhn <mbk225@cornell.edu> Date: Tue Jun 14 23:22:57 2022 -0600 modifying water_air_rarefaction to compile properly with recent changes. commit 0ea59da Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon Jun 13 09:11:15 2022 -0600 explicit diffusion step implemented, with viscous and heat diffusion included. Temperature array and calculation added, along with method to select from available models for material parameters. Needs to be validated. commit fb183aa Author: Michael Kuhn <mbk225@cornell.edu> Date: Fri Jun 3 08:57:08 2022 -0600 added ELVIRA reconstruction method. May not be desired in the long run, but should be useful in development. commit 9d653ab Author: Michael Kuhn <mbk225@cornell.edu> Date: Wed Jul 13 07:47:40 2022 -0600 First-order moments need to be reset after the interface reconstruction and VOF have changed. Hugely important. commit a731028 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon May 30 14:22:21 2022 -0600 little changes commit 1f6f30a Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon May 30 13:48:41 2022 -0600 correcting flaw in calculation - barycenters need to be primitively advected commit ac43b70 Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon May 30 13:03:10 2022 -0600 minor changes to shock_water_cylinder case commit ea4271f Author: Michael Kuhn <mbk225@cornell.edu> Date: Mon May 30 13:01:43 2022 -0600 adding the advection of barycenters. no more segfault in debug for LVIRA reconstruction. commit 6e4d2ce Author: Olivier Desjardins <olivier.desjardins@cornell.edu> Date: Tue Jul 12 09:28:20 2022 -0400 Output of PLIC mesh for film_RT case

Uh oh!

Branch selector

User selector

Datepicker

Commit History

Commits on Feb 9, 2023

Commits on Jan 9, 2023

Commits on Jan 8, 2023

Commits on Dec 18, 2022

Commits on Dec 17, 2022

Commits on Dec 10, 2022

Commits on Dec 9, 2022

Commits on Dec 7, 2022

Commits on Dec 5, 2022

Commits on Nov 30, 2022

Commits on Nov 29, 2022

Commits on Nov 18, 2022

Commits on Nov 7, 2022

Pagination