Skip to content

Commit

Permalink
New pattern: KortewegDeVries1895/kdv.vti
Browse files Browse the repository at this point in the history
  • Loading branch information
timhutton committed Jan 7, 2021
1 parent 151195e commit 176fca7
Show file tree
Hide file tree
Showing 3 changed files with 81 additions and 0 deletions.
1 change: 1 addition & 0 deletions CMakeLists.txt
Original file line number Diff line number Diff line change
Expand Up @@ -267,6 +267,7 @@ set( PATTERN_FILES
Patterns/KuramotoSivashinsky1978/Kuramoto-Sivashinsky_multistable.vti
Patterns/KuramotoSivashinsky1978/Kuramoto-Sivashinsky_travelling_waves.vti
Patterns/KuramotoSivashinsky1978/Kuramoto-Sivashinsky_travelling_waves2.vti
Patterns/KortewegDeVries1895/kdv.vti
)

set( HELP_FILES
Expand Down
1 change: 1 addition & 0 deletions Help/changes.html
Original file line number Diff line number Diff line change
Expand Up @@ -39,6 +39,7 @@ <h3>Changes</h3>
<li><a href="open:Patterns/wave_soliton.vti">wave_soliton.vti</a>
<li><a href="open:Patterns/LotkaVolterra1926/Lotka-Volterra_1D.vti">Lotka-Volterra_1D.vti</a>
<li><a href="open:Patterns/Ginzburg-Landau/complex_Ginzburg-Landau_1D.vti">complex_Ginzburg-Landau_1D.vti</a>
<li><a href="open:Patterns/KortewegDeVries1895/kdv.vti">KortewegDeVries1895/kdv.vti</a>
</ul>
</ul>

Expand Down
79 changes: 79 additions & 0 deletions Patterns/KortewegDeVries1895/kdv.vti
Original file line number Diff line number Diff line change
@@ -0,0 +1,79 @@
<?xml version="1.0"?>
<VTKFile type="ImageData" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
<RD format_version="6">
<description>
The &lt;a href=&quot;https://en.wikipedia.org/wiki/Korteweg%E2%80%93de_Vries_equation&quot;&gt;Korteweg-de Vries (KdV) equation&lt;/a&gt; is a model of waves on shallow water surfaces:

da/dt = k * a * da/dx - d^3a/dx^3

The initial pattern splits into solitons that move at different speeds. The solitons can pass through each other
without being disturbed.

(Eventually instabilities appear because we are using only simple numerical methods.)

</description>
<rule name="advection" type="formula" wrap="1" neighborhood_type="vertex">
<param name="timestep">
0.001000
</param>
<param name="k">
-0.500000
</param>
<formula number_of_chemicals="1">
delta_a = k * a * x_gradient_a - x_deriv3_a;
</formula>
</rule>
<initial_pattern_generator apply_when_loading="true" zero_first="true">
<overlay chemical="a">
<overwrite/>
<gaussian height="1" sigma="0.05">
<point3D x="0.5" y="0.5" z="0.5"/>
</gaussian>
<everywhere/>
</overlay>
</initial_pattern_generator>
<render_settings>
<surface_color r="1" g="1" b="1"/>
<colormap value="spectral"/>
<color_low r="0" g="0" b="1"/>
<color_high r="1" g="0" b="0"/>
<show_color_scale value="false"/>
<show_multiple_chemicals value="true"/>
<active_chemical value="a"/>
<low value="0"/>
<high value="1"/>
<vertical_scale_1D value="30"/>
<vertical_scale_2D value="15"/>
<contour_level value="0.25"/>
<cap_contour value="true"/>
<invert_contour_cap value="false"/>
<use_wireframe value="false"/>
<show_cell_edges value="false"/>
<show_bounding_box value="true"/>
<show_chemical_label value="true"/>
<slice_3D value="true"/>
<slice_3D_axis value="z"/>
<slice_3D_position value="0.5"/>
<show_displacement_mapped_surface value="true"/>
<color_displacement_mapped_surface value="false"/>
<use_image_interpolation value="true"/>
<timesteps_per_render value="1024"/>
<show_phase_plot value="false"/>
<phase_plot_x_axis value="a"/>
<phase_plot_y_axis value="b"/>
<phase_plot_z_axis value="c"/>
<plot_ab_orthogonally value="false"/>
</render_settings>
</RD>
<ImageData WholeExtent="0 255 0 0 0 0" Origin="0 0 0" Spacing="1 1 1" Direction="1 0 0 0 1 0 0 0 1">
<Piece Extent="0 255 0 0 0 0">
<PointData>
<DataArray type="Float64" Name="a" format="binary" RangeMin="0" RangeMax="0">
AQAAAACAAAAACAAAFwAAAA==eJxjYBgFo2AUjIJRMApGwUgDAAgAAAE=
</DataArray>
</PointData>
<CellData>
</CellData>
</Piece>
</ImageData>
</VTKFile>

0 comments on commit 176fca7

Please sign in to comment.