1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
|
##############################################################################
#
# Copyright (c) 2003-2020 by The University of Queensland
# http://www.uq.edu.au
#
# Primary Business: Queensland, Australia
# Licensed under the Apache License, version 2.0
# http://www.apache.org/licenses/LICENSE-2.0
#
# Development until 2012 by Earth Systems Science Computational Center (ESSCC)
# Development 2012-2013 by School of Earth Sciences
# Development from 2014 by Centre for Geoscience Computing (GeoComp)
# Development from 2019 by School of Earth and Environmental Sciences
#
##############################################################################
from __future__ import print_function, division
__copyright__="""Copyright (c) 2003-2020 by The University of Queensland
http://www.uq.edu.au
Primary Business: Queensland, Australia"""
__license__="""Licensed under the Apache License, version 2.0
http://www.apache.org/licenses/LICENSE-2.0"""
__url__="https://launchpad.net/escript-finley"
from esys.escript.modelframe import Model,ParameterSet
from esys.escript.util import exp
import numpy
class GravityForce(ParameterSet):
"""
Sets a gravity force of given direction in given domain:
:note: Instance variable domain - domain of interest (in).
:note: Instance variable domain - `esys.escript.Domain`
:note: Instance variable density - density, default 1 (in).
:note: Instance variable gravity - the gravity constant, default 9.81 (in).
:note: Instance variable direction - the direction of gravity, default [1.,0.,0.] (in).
"""
def __init__(self,**kwargs):
"""
initializes the set
"""
super(GravityForce, self).__init__(**kwargs)
self.declareParameter(domain=None,
gravity=9.81, \
density=1., \
direction=[1.,0.,0.])
def gravity_force(self):
"""
return the gravity force as `density` * `gravity` * `direction`
"""
if isinstance(self.direction,list):
dir=numpy.array(self.direction[:self.domain.getDim()])
else:
dir=self.direction[:self.domain.getDim()]
return self.gravity*self.density*dir
class MaterialTable(ParameterSet):
"""
A simple matrial table which allows setting physical ivar of a model
:ivar density: density(in/out)
:ivar heat_capacity: heat_capacity(in/out)
:ivar thermal_permabilty: permabilty(in/out)
:ivar viscosity: viscosity(in/out)
:ivar radiation_coefficient: (in/out)
"""
def __init__(self,**kwargs):
super(MaterialTable, self).__init__(**kwargs)
self.declareParameter(gravity=9.81, \
density=1., \
heat_capacity=1., \
thermal_permabilty=1., \
radiation_coefficient=0.)
class SimpleEarthModel(ParameterSet):
"""
A simple matrial table run convection models::
density=density0*(1-expansion_coefficient*(temperature-reference_temperature))
viscocity=viscocity0*(exp(alpha*(1/reference_temperature - 1/temperature))
:ivar gravity: gravity constants (9.81) (in)
:ivar reference_temperature: reference temperature (in)
:ivar density0: density at reference temperature (in)
:ivar viscosity0: viscosity0 at reference temperature (in)
:ivar alpha: viscosity contrast (in)
:ivar expansion_coefficient: Raleigh number (in)
:ivar heat_capacity: heat capacity (in)
:ivar thermal_permabilty: permabilty (in)
:ivar temperature: temperature (in)
:ivar viscosity: viscosity (out)
:ivar density: density (out)
"""
def __init__(self,**kwargs):
super(SimpleEarthModel, self).__init__(**kwargs)
self.declareParameter(reference_temperature=1.,
gravity=9.81, \
density0=1., \
expansion_coefficient=0., \
viscosity0=1., \
alpha=0., \
temperature=0.,\
heat_capacity=1., \
thermal_permabilty=1.)
def density(self):
return self.density0*(1-self.expansion_coefficient*(self.temperature-self.reference_temperature))
def viscosity(self):
return self.viscosity0*exp(self.alpha*(1/self.reference_temperature - 1/(self.temperature+1.e-15)))
class SimpleSolidMaterial(MaterialTable):
"""
A simple matrial table which allows setting physical parameters of
a model.
"""
def __init__(self,**kwargs):
super(MaterialTable, self).__init__(**kwargs)
self.declareParameter(lame_lambda=1.,\
lame_my=1.)
class SimpleFluidMaterial(MaterialTable):
"""
A simple matrial table which allows setting physical ivar of a model.
"""
def __init__(self,**kwargs):
super(MaterialTable, self).__init__(**kwargs)
self.declareParameter(viscosity=1., \
hydraulic_conductivity=1.e-4)
# vim: expandtab shiftwidth=4:
|
