Skip to content

d2q9_fracture_dissolution implementation #6

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
Rymarz wants to merge 40 commits into
base: develop
Choose a base branch
from
Open

d2q9_fracture_dissolution implementation #6

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
40 commits
Select commit Hold shift + click to select a range
ac9e1f4
Add d2q9_fracture_dissolution (d2q9_thin_film copy)
Rymarz Sep 30, 2021
4180190
Add concentration with BKG collision
Rymarz Oct 5, 2021
3ed52ef
Add first order/linear reaction
Rymarz Oct 5, 2021
e8ace53
Add h_z evolution equation
Rymarz Oct 5, 2021
dcde969
Correct Dynamics.R
Rymarz Oct 5, 2021
87a43de
Rrun initialization of fields
Rymarz Oct 7, 2021
0713291
Correction of initialization from image
Rymarz Oct 12, 2021
a246f1b
Correction of h_Z in momentum equation
Rymarz Oct 13, 2021
16e1bdc
InitFromFields repair
Nov 8, 2021
3748dc9
Const saturation concentration at the outlet
Nov 17, 2021
a954d96
Change h_Z from field to density
Rymarz Nov 17, 2021
1d283dc
Haight_parameter added
Nov 19, 2021
67555bb
AddSetting C_saturation corrected
Rymarz Nov 19, 2021
77db5bf
Add pressure quantity
Dec 3, 2021
9bd63fc
d2q5 MRT for solute with reaction
Rymarz Dec 16, 2021
1dedf63
Names refactor
Rymarz Dec 20, 2021
39dfd6d
merge conflict resolution
Rymarz Dec 20, 2021
5447531
Names refactor
Rymarz Dec 21, 2021
80a8a61
Concentration d2q9 and TRT
Rymarz Dec 21, 2021
19f60d9
merged and resolved the conflict in fr.xml
Rymarz Jan 29, 2022
623705f
Incompressible equilibrium distirbution implementation
Rymarz Jan 31, 2022
c9908e2
Incompressible equlibrium for flow
Rymarz Feb 2, 2022
bfff3b0
IncompressibleMRT implementation
Rymarz May 10, 2022
cf487f2
Concentration in kg/m3
Rymarz May 16, 2022
bf421e8
Initial copy from d2q9
Rymarz May 27, 2022
12ae3a2
Add Brinkman force
Rymarz May 27, 2022
096d9ea
Brinkman force correction
Rymarz May 27, 2022
5d1030d
d2q9 concentration without source
Rymarz May 28, 2022
4aef790
Adding concentration BCs
Rymarz May 28, 2022
eb6523b
d2q5 concentration
May 28, 2022
c749a45
Reaction added
May 28, 2022
fa8eb37
Incompressible MRT for fluid
Rymarz May 29, 2022
b636b34
Relaxation C changes
Rymarz May 30, 2022
51a7a18
Concentration SRT
Rymarz May 30, 2022
72a8787
NodeTypes refactor
Rymarz May 31, 2022
9462580
Uniform ReactionConstant
Jun 6, 2022
ea2423a
Height_parameter added
Jun 6, 2022
e59c96a
TRT and magicNumber added
Jun 7, 2022
4f8a791
Adding h evolution (naive)
Jun 7, 2022
946498e
Add hEvolution step
Jun 8, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Binary file added example/data/1024x256_initial_hz.png
View file
Open in desktop
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file added example/data/128x32.png
View file
Open in desktop
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file added example/data/32x32.png
View file
Open in desktop
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
45 changes: 45 additions & 0 deletions example/flow/2d/cave.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,45 @@
<?xml version="1.0"?>
<CLBConfig version="2.0" output="output/" permissive="true">
<Units>
<Param name="x" value="0.009765625m" gauge="1"/>
<!-- <Param name="nu" value="1e-6m2/s" gauge="0.1666666666"/> -->
<Param name="Rho" value="1000kg/m3" gauge="1"/>
</Units>
<Geometry nx="1024" ny="512">
<MRT>
<Box />
</MRT>
<WVelocity name="Inlet">
<Inlet />
</WVelocity>
<EPressure name="Outlet">
<Outlet />
</EPressure>
<Inlet nx="1" dx="1">
<Box />
</Inlet>
<Outlet nx="1" dx="-1">
<Box />
</Outlet>
</Geometry>
<RunR>
Solver$Fields$Height[] = rnorm(64*32, mean = 0.02048, sd=0.000002048)

Solver$Actions$InitFromFields()
Solver$Actions$InitFromFields()
</RunR>
<Model>
<Param name="VelocityX" value="1e-4m/s" />
<Param name="VelocityX" value="1e-4m/s" zone="Inlet" />
<Param name="Pressure" value="0.0" />
<Param name="C_saturation" value="2.0/m3" />
<Param name="Concentration" value="2.0/m3" />
<Param name="Inlet_concentration" value="0.0" />
<Param name="k" value="2.5e-6m/s" />
<Param name="D" value="10e-9m2/s" />
</Model>
<HDF5 />
<Solve Iterations="1e1">
<HDF5 Iterations="1"/>
</Solve>
</CLBConfig>
76 changes: 76 additions & 0 deletions example/flow/2d/fr.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,76 @@
<?xml version="1.0"?>
<CLBConfig version="2.0" output="output/" permissive="true">
<Units>
<Param name="Viscosity" value="1e-6m2/s" gauge="0.000166666666667" />
<Param name="C_l" value="0.01m" gauge="1" />
<Param name="rho" value="1000kg/m3" gauge="1"/>
</Units>
<Geometry nx="10m" ny="6">
<MRT>
<Box />
</MRT>
<WPressureWConcentration name="Inlet">
<Inlet />
</WPressureWConcentration>
<EPressureEConcentration name="Outlet">
<Outlet />
</EPressureEConcentration>
<Inlet nx="1" dx="1">
<Box />
</Inlet>
<Outlet nx="-1" dx="-1">
<Box />
</Outlet>
</Geometry>
<Model>
<!-- <Param name="VelocityX" value="5e-4m/s"/> -->
<!-- 98.1Pa for each meter -->
<<<<<<< HEAD
<Param name="Pressure" value="0.001Pa" zone="Inlet" />
<Param name="Pressure" value="0" zone="Outlet" />
<Param name="C_saturation" value="2/m3" />
<Param name="Concentration" value="0" />
<Param name="k" value="2.5e-7m/s" />
<Param name="C_solid" value="2.7e4/m3" />
<Param name="Height" value="2e-4m" />
<Param name="Diffusivity" value="1e-6m2/s" />
<!-- <Param name="MagicNumber_D" value="0.25" /> -->
</Model>


<Solve Iterations="100">
<HDF5 Iterations="1" />
</Solve>

<!-- <Param name="k" value="2.5e-7m/s" />
<Param name="Concentration" value="0" zone="Inlet" /> -->

<Solve Iterations="100">
<HDF5 Iterations="1" />
=======
<Param name="Pressure" value="9.81Pa" zone="Inlet"/>
<Param name="Pressure" value="0" zone="Outlet"/>
<Param name="Saturation" value="2.0/m3" />
<Param name="Concentration" value="2.0/m3" />
<Param name="k" value="0" />
<Param name="SolidConcentration" value="2.7e4/m3" />
<Param name="Height" value="2e-4m" />
<Param name="Diffusivity" value="1e-15m2/s" />
<Param name="MagicNumber_D" value="0.25" />
</Model>

<HDF5 />

<Solve Iterations="5000000">
<HDF5 Iterations="100000"/>
</Solve>

<Param name="k" value="2.5e-7m/s" />
<Param name="Concentration" value="0.01/m3" zone="Inlet" />

<Solve Iterations="5000000">
<HDF5 Iterations="100000"/>
>>>>>>> soluteMRT
</Solve>

</CLBConfig>
73 changes: 73 additions & 0 deletions example/flow/2d/fracture.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,73 @@
<?xml version="1.0"?>
<CLBConfig version="2.0" output="output/Pe=100_Da=5e-2_lam=2e-4" permissive="true">
<!-- C = physical_value / gauge -->
<Units>
<Param name="Viscosity" value="1e-6m2/s" gauge="1" />
<Param name="C_l" value="0.02m" gauge="1" />
<Param name="rho" value="1000kg/m3" gauge="1"/>
</Units>
<Geometry nx="1024" ny="128">
<IncompressibleSRT>
<Box />
</IncompressibleSRT>
<WPressureWConcentration name="Inlet">
<Inlet />
</WPressureWConcentration>
<EPressureEConcentration name="Outlet">
<Outlet />
</EPressureEConcentration>
<Inlet nx="1" dx="1">
<Box />
</Inlet>
<Outlet nx="1" dx="-1">
<Box />
</Outlet>
</Geometry>

<!-- Pe = q_0/D = 100 -->
<!-- Da = kh_0 / q_0 = 5e-2 -->
<!-- Lambda = h_0 / L = 1e-5 -->
<!-- Dimensionless -->
<Model>
<Param name="Pressure" value=".1" zone="Inlet" />
<Param name="Pressure" value="-.1" zone="Outlet" />
<Param name="Saturation" value="3/m3" />
<Param name="Concentration" value="3/m3" />
<Param name="k" value="0" />
<Param name="SolidConcentration" value="2.7e4/m3" />
<Param name="Height" value="0.2048" />
<Param name="Diffusivity" value="1" />
<Param name="MagicNumber_D" value="0.25" />
</Model>

<!-- Dimensioned -->
<!-- <Param name="VelocityX" value="5e-4m/s"/> -->
<!-- 98.1Pa for each meter -->
<!-- <Model>
<Param name="Pressure" value="0.01" zone="Inlet"/>
<Param name="Pressure" value="0" zone="Outlet"/>
<Param name="Saturation" value="3.0/m3" />
<Param name="Concentration" value="3.0/m3" />
<Param name="k" value="0" />
<Param name="SolidConcentration" value="2.7e4/m3" />
<Param name="Height" value="2e-4m" />
<Param name="Diffusivity" value="6.73e-10m2/s" />
<Param name="MagicNumber_D" value="0.25" />
</Model> -->
<RunR>
<!-- Solver$Fields$Height_parameter[] = rnorm(100*16, mean = 0.2048, sd=0.0002);
Solver$Actions$InitFromFields(); -->
</RunR>

<HDF5 />

<Solve Iterations="100"></Solve>

<Param name="k" value="2.44140625E-04" />
<Param name="Concentration" value="1/m3" zone="Inlet" />

<Solve Iterations="100">
<HDF5 Iterations="10" />
</Solve>

</CLBConfig>
92 changes: 92 additions & 0 deletions example/flow/2d/fracture_dissolution.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,92 @@
<!-- <?xml version="1.0"?>
<CLBConfig version="2.0" output="output/" permissive="true">
<Geometry nx="1024" ny="512">
<Reaction>
<Box />
</Reaction>
<MRT>
<Box />
</MRT>
<WVelocity name="Inlet">
<Inlet />
</WVelocity>
<EPressure name="Outlet">
<Outlet />
</EPressure>
<Inlet nx="1" dx="5">
<Box />
</Inlet>
<Outlet nx="1" dx="-5">
<Box />
</Outlet>
</Geometry>
<RunR>
Solver$Fields$Height[] = rnorm(64*32, mean = 1, sd=0.001)

Solver$Actions$InitFromFields()
</RunR>
<Model>
<Param name="VelocityX" value="0.1" />
<Param name="Viscosity" value="0.02" />
<Param name="D" value="0.01" />
<Param name="C_saturation" value="0.0" />
<Param name="C_inlet" value="0.0" />
<Param name="k" value="0.0001" />
</Model>
<HDF5 />
<Solve Iterations="1e6">
<HDF5 Iterations="1"/>
</Solve>
</CLBConfig> -->

<?xml version="1.0"?>
<CLBConfig version="2.0" output="output/" permissive="true">
<Units>
<Param name="size" value="0.01m" gauge="1"/>
<Param name="nu" value="1e-5m2/s" gauge="0.1666666666"/>
</Units>
<Geometry nx="20m" ny="10m">
<MRT>
<Box />
</MRT>
<WVelocity name="Inlet">
<Inlet />
</WVelocity>
<EPressure name="Outlet">
<Outlet />
</EPressure>
<Inlet nx="1" dx="5">
<Box />
</Inlet>
<Outlet nx="1" dx="-5">
<Box />
</Outlet>
</Geometry>
<RunR>
Solver$Fields$Height[] = rnorm(100, mean = 0.02, sd=0.01)

Solver$Actions$InitFromFields()
Solver$Actions$InitFromFields()
Solver$Actions$InitFromFields()

</RunR>
<Model>
<Param name="VelocityX" value="0" />
<Param name="VelocityX" value="0.0001m/s" zone="Inlet" />
<Param name="Pressure" value="0.0" />
<Param name="C_saturation" value="2/m3" />
<Param name="C_init" value="0" />
<Param name="C_inlet" value="0.0" />
<Param name="C_outlet" value="2/m3" />
<Param name="D" value="1e-5m2/s" />
<Param name="k" value="1" />
<Param name="C_solid" value="2.7e4/m3" />

</Model>
<HDF5 />
<Solve Iterations="1e9">
<HDF5 Iterations="1e5"/>
</Solve>
</CLBConfig>


8 changes: 8 additions & 0 deletions models/flow/d2q9_fracture/Description.md
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,8 @@
The `d2q9` model is a basic 2D Lattice Boltzmann Method model for flow
simulation.

It is a implementation of a simple Multiple-Relaxation-Time collision
operator.

The model has one [option](/basics/options/): BC. This option turns on
custom fields for setting non-standard profiles on inlets and other things.
Loading