Neko  0.9.99
A portable framework for high-order spectral element flow simulations
advection.f90
Go to the documentation of this file.
1 ! Copyright (c) 2021-2024, The Neko Authors
2 ! All rights reserved.
3 !
4 ! Redistribution and use in source and binary forms, with or without
5 ! modification, are permitted provided that the following conditions
6 ! are met:
7 !
8 ! * Redistributions of source code must retain the above copyright
9 ! notice, this list of conditions and the following disclaimer.
10 !
11 ! * Redistributions in binary form must reproduce the above
12 ! copyright notice, this list of conditions and the following
13 ! disclaimer in the documentation and/or other materials provided
14 ! with the distribution.
15 !
16 ! * Neither the name of the authors nor the names of its
17 ! contributors may be used to endorse or promote products derived
18 ! from this software without specific prior written permission.
19 !
20 ! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
21 ! "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
22 ! LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
23 ! FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
24 ! COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
25 ! INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
26 ! BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
27 ! LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
28 ! CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
29 ! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
30 ! ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
31 ! POSSIBILITY OF SUCH DAMAGE.
32 !
34 module advection
35  use num_types, only : rp
36  use space, only : space_t
37  use field, only : field_t
38  use coefs, only : coef_t
39  use json_module, only : json_file
40  use field_series, only: field_series_t
42  implicit none
43  private
44 
46  type, public, abstract :: advection_t
47  contains
48  procedure(compute_adv), pass(this), deferred :: compute
49  procedure(compute_scalar_adv), pass(this), deferred :: compute_scalar
50  procedure(advection_free), pass(this), deferred :: free
51  end type advection_t
52 
53  interface
54 
67  module subroutine advection_factory(object, json, coef, &
68  ulag, vlag, wlag, dtlag, tlag, time_scheme, use_dummy, slag)
69  class(advection_t), allocatable, intent(inout) :: object
70  type(json_file), intent(inout) :: json
71  type(coef_t), intent(inout), target :: coef
72  type(field_series_t), intent(in), target :: ulag, vlag, wlag
73  real(kind=rp), intent(in), target :: dtlag(10)
74  real(kind=rp), intent(in), target :: tlag(10)
75  type(time_scheme_controller_t), intent(in), target :: time_scheme
76  logical, optional, intent(in) :: use_dummy
77  type(field_series_t), target, optional, intent(in) :: slag
78  end subroutine advection_factory
79  end interface
80 
81  public :: advection_factory
82 
83  abstract interface
84 
96  subroutine compute_adv(this, vx, vy, vz, fx, fy, fz, Xh, coef, n, dt)
97  import :: advection_t
98  import :: coef_t
99  import :: space_t
100  import :: field_t
101  import :: rp
102  class(advection_t), intent(inout) :: this
103  type(space_t), intent(inout) :: Xh
104  type(coef_t), intent(inout) :: coef
105  type(field_t), intent(inout) :: vx, vy, vz
106  type(field_t), intent(inout) :: fx, fy, fz
107  integer, intent(in) :: n
108  real(kind=rp), intent(in), optional :: dt
109  end subroutine compute_adv
110  end interface
111 
112  abstract interface
113 
124  subroutine compute_scalar_adv(this, vx, vy, vz, s, fs, Xh, coef, n, dt)
125  import :: advection_t
126  import :: coef_t
127  import :: space_t
128  import :: field_t
129  import :: rp
130  class(advection_t), intent(inout) :: this
131  type(field_t), intent(inout) :: vx, vy, vz
132  type(field_t), intent(inout) :: s
133  type(field_t), intent(inout) :: fs
134  type(space_t), intent(inout) :: Xh
135  type(coef_t), intent(inout) :: coef
136  integer, intent(in) :: n
137  real(kind=rp), intent(in), optional :: dt
138  end subroutine compute_scalar_adv
139  end interface
140 
141  abstract interface
142 
143  subroutine advection_free(this)
144  import :: advection_t
145  class(advection_t), intent(inout) :: this
146  end subroutine advection_free
147  end interface
148 
149 end module advection
A factory for advection_t decendants. Both creates and initializes the object.
Definition: advection.f90:96
Add advection operator to the right-hand-side for a scalar.
Definition: advection.f90:124
Subroutines to add advection terms to the RHS of a transport equation.
Definition: advection.f90:34
Coefficients.
Definition: coef.f90:34
Stores a series fields.
Defines a field.
Definition: field.f90:34
integer, parameter, public rp
Global precision used in computations.
Definition: num_types.f90:12
Defines a function space.
Definition: space.f90:34
Compound scheme for the advection and diffusion operators in a transport equation.
Base class for time integration schemes.
Definition: time_scheme.f90:61
Base abstract type for computing the advection operator.
Definition: advection.f90:46
Coefficients defined on a given (mesh, ) tuple. Arrays use indices (i,j,k,e): element e,...
Definition: coef.f90:55
The function space for the SEM solution fields.
Definition: space.f90:62
Implements the logic to compute the time coefficients for the advection and diffusion operators in a ...