Limiter & Rhs

Project.toml

@@ -1,7 +1,7 @@
 name = "SegregatedVMSSolver"
 uuid = "90220dfc-fac3-4908-b895-51c1a1f4cb85"
 authors = ["Carlo Brunelli"]
-version = "3.0.0"
+version = "3.0.1"
 
 [deps]
 CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
@@ -32,27 +32,29 @@ Trapz = "592b5752-818d-11e9-1e9a-2b8ca4a44cd1"
 UnPack = "3a884ed6-31ef-47d7-9d2a-63182c4928ed"
 
 [compat]
-CSV = "0.10.15"
-DataFrames = "1.7.0"
-DataInterpolations = "6.6.0"
-DelimitedFiles = "1.9.1"
-FFTW = "1.8.0"
-FileIO = "1.16.6"
-FillArrays ="=1.13.0"
-Gridap = "=0.17.20"
-GridapDistributed = "=0.3.2"
-GridapGmsh = "=0.7.1"
+CSV = "=0.10.15"
+DataFrames = "=1.7.0"
+DataInterpolations = "=8.0.1"
+DelimitedFiles = "=1.9.1"
+FFTW = "=1.9.0"
+FileIO = "=1.17.0"
+FillArrays = "=1.13.0"
+Gridap = "=0.18.12"
+GridapDistributed = "=0.4.7"
+GridapGmsh = "=0.7.2"
 GridapPETSc = "=0.5.3"
-MPI ="=0.20.22"
-NearestNeighbors ="0.4.21"
-Parameters = "0.12.3"
-PartitionedArrays ="=0.3.4"
-Revise ="3.7.1"
-ScatteredInterpolation = "0.3.6"
 LinearAlgebra = "1.11.0"
+MPI = "=0.20.22"
+NearestNeighbors = "=0.4.21"
+Parameters = "=0.12.3"
+PartitionedArrays = "=0.3.5"
+Pkg = "1.11.0"
+Random = "1.11.0"
+Revise = "=3.8.0"
+ScatteredInterpolation = "=0.3.6"
 SparseArrays = "1.11.0"
+Statistics = "1.11.1"
 SyntheticEddyMethod = "=0.4.4"
-Trapz ="2.0.3"
 Test = "1.11.0"
-UnPack = "1.0.2"
-julia = "1.11.1"
+Trapz = "=2.0.3"
+UnPack = "=1.0.2"

TGV2D.jl

@@ -0,0 +1,61 @@
+using PartitionedArrays
+using SegregatedVMSSolver
+using SegregatedVMSSolver.ParametersDef
+using SegregatedVMSSolver.SolverOptions
+using MPI
+
+
+
+
+
+t0 = 0.0
+dt = 0.01
+tF = 10*dt
+vortex_diameter = 1.0
+N = 16
+Re = 1600
+D = 2
+
+backend = with_debug
+
+rank_partition = (2,2)
+
+
+
+
+
+# solver_options = petsc_options(; vel_ksp="gmres", vel_pc="gamg", pres_ksp="cg", pres_pc="gamg")
+solver_options = " -vel_ksp_type gmres -vel_pc_type gamg -vel_ksp_rtol 1.e-10 -vel_ksp_converged_reason \
+  -pres_ksp_type cg -pres_pc_type gamg -pres_ksp_rtol 1.e-6 -pres_ksp_converged_reason \
+  -ksp_atol 0.0 -vec_type cuda -mat_type aijcusparse -log_view"
+sprob = StabilizedProblem(VMS(3))
+timep = TimeParameters(t0=t0, dt=dt, tF=tF)
+
+physicalp = PhysicalParameters(Re=Re, c=vortex_diameter)
+solverp = SolverParameters(matrix_freq_update=1, Number_Skip_Expansion=10e6, M=40,
+petsc_options=solver_options)
+exportp = ExportParameters(printinitial=true, printmodel=true, 
+vtu_export = ["uh","ph","uh_analytic", "ph_analytic"], extra_export=["VelocityError","PressureError"])
+
+
+
+
+meshp = MeshParameters(rank_partition, D; N=N, L=0.5 * vortex_diameter)
+simparams = SimulationParameters(timep, physicalp, solverp, exportp)
+bc_tgv = Periodic(meshp, physicalp)
+
+
+
+mcase = TaylorGreen(bc_tgv, meshp, simparams, sprob)
+
+
+
+# Create folder and file
+# mkdir("Log")
+# open("Log/PrintSim.txt", "w") do file
+# end
+
+
+
+SegregatedVMSSolver.solve(mcase, backend)
+

TGV2D_limiter.jl

@@ -0,0 +1,82 @@
+using PartitionedArrays
+using SegregatedVMSSolver
+using SegregatedVMSSolver.ParametersDef
+using SegregatedVMSSolver.SolverOptions
+using MPI
+using Gridap
+
+
+
+    N = 128
+    dt = 0.025
+
+    order = 2
+    t0 =0.0
+    vortex_diameter = 1.0
+    L = 0.5 * vortex_diameter
+
+    u_conv = 0.2
+    Re = 1_000
+
+    D = 2
+
+    CFL_eff = u_conv*N * dt * order / (2*L)
+
+    backend = with_debug
+
+    rank_partition = (4,4)
+    tF = 0.1 #dt*10
+
+
+    println("------------------------")
+    println("CFL = $(CFL_eff)")
+    println("------------------------")
+
+    function petsc_options_cstm()
+        return " -vel_ksp_type gmres -vel_pc_type gamg  -vel_ksp_rtol 1.e-8 -vel_ksp_converged_reason \
+        -pres_ksp_type cg -pres_pc_type gamg  -pres_ksp_rtol 1.e-8 -pres_ksp_converged_reason"
+    end
+
+
+    my_lim = CustomLimiter() do dt, uun, G, GG, ν
+        h = L /(order*N)
+        dt_eff = h/u_conv
+        return max(dt_eff,dt)
+
+    end
+
+
+
+    sprob = StabilizedProblem(VMS(order), TensorFormulation(dt_limiter = my_lim), false)
+
+
+
+
+    # sprob = StabilizedProblem(method=VMS(order), coeff_method=TensorFormulation(r=2, Ci=[4, 36]), skew=false)
+    timep = TimeParameters(t0=t0,dt=dt,tF=tF)
+
+    physicalp = PhysicalParameters(Re=Re,c=vortex_diameter, u_in_mag=u_conv)
+
+    solverp = SolverParameters(matrix_freq_update = 1, Number_Skip_Expansion=10e6, M = 40, a_err_threshold=10_000,
+    petsc_options = petsc_options_cstm())
+    exportp = ExportParameters(printinitial=true,printmodel=true, extra_export=["VelocityError","PressureError"])
+
+
+    meshp= MeshParameters(rank_partition,D;N=N,L=L)
+
+
+    simparams = SimulationParameters(timep,physicalp,solverp,exportp)
+
+    params_tvg = TaylorGreenParameters(Vs=u_conv, Ua=0.3*u_conv, Va = 0.2*u_conv)
+    bc_tgv = Periodic(meshp,physicalp,params_tvg ) 
+
+
+
+    mcase = TaylorGreen(bc_tgv, meshp,simparams,sprob)
+
+
+    SegregatedVMSSolver.solve(mcase,backend)
+
+
+
+

paper/paper.bib

@@ -7,7 +7,7 @@ @article{Badia:2020
   number = {52},
   pages = {2520},
   author = {Santiago Badia and Francesc Verdugo},
-  title = {Gridap: An extensible Finite Element toolbox in Julia},
+  title = {Gridap: An extensible Finite Element toolbox in {J}ulia},
   journal = {Journal of Open Source Software}
 }
 
@@ -20,7 +20,7 @@ @article{Verdugo:2022
   volume = {276},
   pages = {108341},
   author = {Francesc Verdugo and Santiago Badia},
-  title = {The software design of Gridap: A Finite Element package based on the Julia {JIT} compiler},
+  title = {The software design of {Gridap}: A Finite Element package based on the {J}ulia {JIT} compiler},
   journal = {Computer Physics Communications}
 }
 
@@ -33,7 +33,7 @@ @article{BadiaD:2022
   number = {74},
   pages = {4157},
   author = {Santiago Badia and Alberto F. Martín and Francesc Verdugo},
-  title = {GridapDistributed: a massively parallel finite element toolbox in Julia},
+  title = {GridapDistributed: a massively parallel finite element toolbox in {J}ulia},
   journal = {Journal of Open Source Software}
 }
 
@@ -72,14 +72,77 @@ @article{Hughes:2000
   doi = {10.1007/s007910050051},
 }
 
-@article{Banyai:2016,
+@phdthesis{Banyai:2016,
    author = {Tamas Banyai},
    isbn = {9782875161093},
-   institution = {Université Libre de Bruxelles},
-   title = {Development of Stabilized Finite Element Method for Numerical Simulation of Turbulent Incompressible Single and Eulerian-Eulerian Two-Phase Flows},
+   school = {Université Libre de Bruxelles, Ecole polytechnique de Bruxelles},
+   address = {Bruxelles},
+   title = {Development of Stabilized Finite Element Method for Numerical Simulation of Turbulent Incompressible Single and {E}ulerian-{E}ulerian Two-Phase Flows},
    year = {2016},
+   type = {PhD thesis},
+   month = aug,
+   url = {http://hdl.handle.net/2013/ULB-DIPOT:oai:dipot.ulb.ac.be:2013/235110}
 }
 
 
- @article{botella_peyret:1998, 
- title={Benchmark spectral results on the lid-driven cavity flow}, volume={27}, DOI={10.1016/s0045-7930(98)00002-4}, number={4}, journal={Computers and Fluids}, author={Botella, O. and Peyret, R.}, year={1998}, pages={421–433}} 
+@article{botella_peyret:1998, 
+  title={Benchmark spectral results on the lid-driven cavity flow}, 
+  volume={27}, 
+  DOI={10.1016/s0045-7930(98)00002-4}, 
+  number={4}, 
+  journal={Computers and Fluids}, 
+  author={Botella, O. and Peyret, R.}, 
+  year={1998}, 
+  pages={421–433}
+}
+
+
+
+@misc{Trixi:2025,
+  title={{T}rixi.jl: {A}daptive high-order numerical simulations
+         of hyperbolic {PDE}s in {J}ulia},
+  author={Schlottke-Lakemper, Michael and Gassner, Gregor J and
+          Ranocha, Hendrik and Winters, Andrew R and Chan, Jesse
+          and Rueda-Ramírez, Andrés},
+  year={2025},
+  howpublished={\url{https://github.com/trixi-framework/Trixi.jl}},
+  doi={10.5281/zenodo.3996439}
+}
+
+
+
+@misc{WaterLily:2024,
+      title={WaterLily.jl: A differentiable and backend-agnostic Julia solver to simulate incompressible viscous flow and dynamic bodies}, 
+      author={Gabriel D. Weymouth and Bernat Font},
+      year={2024},
+      eprint={2407.16032},
+      archivePrefix={arXiv},
+      primaryClass={physics.flu-dyn},
+      url={https://arxiv.org/abs/2407.16032}, 
+}
+
+
+@article{Oceananigans:2020,
+ doi = {10.21105/joss.02018}, url = {https://doi.org/10.21105/joss.02018}, year = {2020}, publisher = {The Open Journal}, volume = {5}, number = {53}, pages = {2018}, author = {Ali Ramadhan and Gregory LeClaire Wagner and Chris Hill and Jean-Michel Campin and Valentin Churavy and Tim Besard and Andre Souza and Alan Edelman and Raffaele Ferrari and John Marshall}, title = {Oceananigans.jl: Fast and friendly geophysical fluid dynamics on GPUs}, journal = {Journal of Open Source Software} }
+
+
+
+ @misc{Verdugo:2025,
+  author       = {Francesc Verdugo and Alberto F. Martin and HRvanElderen and GeliezaK and Jordi Manyer Fuertes and Fredrik Ekre and Joroks and Dennis Ogiermann and Yung-sheng Tu and Hendrik Ranocha and Oriol Colomés and Mosè Giordano and Oscar Smith and Robin Richardson},
+  title        = {PartitionedArrays.jl: PartitionedArrays-v0.5.10},
+  year         = 2025,
+  publisher    = {Zenodo},
+  version      = {v0.5.10},
+  doi          = {10.5281/zenodo.15266760},
+  url          = {https://github.com/PartitionedArrays/PartitionedArrays.jl}
+}
+
+@misc{GridapPETSc:2025,
+  author       = {Francesc Verdugo and Santiago Badia and Alberto F. Martin and Jordi Manyer Fuertes and Eric Neiva and the Gridap Team},
+  title        = {GridapPETSc.jl: A PETSc backend for Gridap.jl},
+  year         = {2025},
+  howpublished = {\url{https://github.com/gridap/GridapPETSc.jl}},
+}
+
+@misc{@BrunelliGmsh:2025,
+  author = {Carlo Brunelli}, title={AirfoilGmsh.jl: A Julia packgage for generating airfoil meshes using Gmsh}, year={2025}, howpublished={\url{https://github.com/carlodev/AirfoilGmsh.jl}}}