Added documentation as-is

.gitattributes

@@ -0,0 +1 @@
+*.png binary

README.md

@@ -0,0 +1,35 @@
+# TaSMET
+Thermoacoustic System Modeling Environment Twente
+
+Welcome to the ThermoAcoustic System Modeling Environment Twente, or
+TASMET. TASMET is a computer code to model thermoacoustic (TA)
+engines, refrigerators and combined systems by providing nonlinear
+models for laminar/turbulent oscillating flow in ducts, heat
+exchangers and stacks/regenerators. A coupling to the mechanical and
+electrical domain is provided with a piston model.
+
+- For a quick glance on the possibilities of TaSMET, please look into
+the User's guide provided with the code.
+- For the installation guide, please look into the [Wiki][Wiki]
+- Examples of working models will be provided. We are currently
+  working on that. Once a regenerator model is implemented, we will
+  provide examples of Stirling models.
+
+
+[Wiki]: https://github.com/amdj/tasmet/wiki
+
+## Example models
+
+Example models for TaSMET will be provided in the
+[tasmet_examples][ta_ex] repository.
+
+
+
+[ta_ex]: https://github.com/amdj/tasmet_examples
+
+
+## Release notes version 0.1
+- Stacks cannot yet be connected to ConnectorVolumes
+- Stacks cannot yet be connected to Pistons
+- VelocityBc probably still has a bug
+

doc/bib/TaSMET.bib

@@ -0,0 +1,25 @@
+
+@article{de_jong_modeling_2015,
+	title = {Modeling of thermoacoustic systems using the nonlinear frequency domain method},
+	journal = {To be published},
+	author = {de Jong, J. A. and Wijnant, Y. H. and Wilcox, D. and de Boer, A},
+	year = {2015}
+}
+
+@article{de_jong_frequency_2015,
+	title = {A frequency domain method for solving periodic nonlinear acoustic problems},
+	journal = {To be published},
+	author = {de Jong, J. A. and Wijnant, Y. H. and Wilcox, D. and de Boer, A},
+	year = {2015}
+}
+
+@phdthesis{jong_numerical_2015,
+	address = {Enschede},
+	title = {Numerical modeling of thermoacoustic systems},
+	url = {http://doc.utwente.nl/96275/},
+	abstract = {The subject of this thesis is a relatively new class of heat engines and refrigerators, called thermoacoustic (TA) systems. TA systems have gained commercial interest due to their low number of moving parts and potentially high efficiency. In the case of a TA engine, heat is converted to acoustic power. This power can subsequently be converted to electricity using a ?reversed? loudspeaker, called a linear alternator. In a TA refrigerator, a speaker or linear alternator is used to generate a strong acoustic wave, which is used to pump heat. To achieve competitive power densities, thermoacoustic systems are generally run at such high amplitudes, that performance deteriorating nonlinear effects can no longer be neglected. To accu- rately predict performance in the nonlinear regime, nonlinear models are required. This thesis describes two contributions to the field of thermoacoustic system modeling. Firstly, a one-dimensional heat transfer model has been developed. This model can be used to estimate the performance of often used parallel-plate heat exchangers for thermoacoustic systems. These heat exchangers are located close to the stack or regenerator of a TA system and are responsible for the heat in/output required to let the system execute its thermodynamic cycle. The results of the model show a good match with a different heat transfer model from the literature, and the model provides guidelines for future heat exchanger design. Secondly, a nonlinear frequency domain method is developed with which the initial transient start-up process can be skipped in the simulations. The method can be used to directly simulate a TA system in its periodic steady-state. This significantly reduces computational cost, since the initial transient regime often involves several hundred oscillation cycles. The method is applied to a one-dimensional nonlinear model of TA systems. The model is used to simulate an experimental standing wave thermoacoustic engine from the literature. The obtained results are in agreement with literature results.},
+	school = {Universiteit Twente},
+	author = {Jong, Jan Albertus de},
+	year = {2015},
+	file = {Jong - 2015 - Numerical modeling of thermoacoustic systems.pdf:/home/anne/.literature/storage/GQTWDG7B/Jong - 2015 - Numerical modeling of thermoacoustic systems.pdf:application/pdf}
+}

doc/fig/duct_vars.eps

@@ -0,0 +1,219 @@
+%!PS-Adobe-3.0 EPSF-3.0
+%%Creator: Ipelib 70108 (Ipe 7.1.8)
+%%CreationDate: D:20150928213252
+%%LanguageLevel: 2
+%%BoundingBox: 126 703 387 769
+%%HiResBoundingBox: 126 703.8 386 768.2
+%%EndComments
+%%BeginProlog
+%%BeginResource: procset ipe 7.0 70108
+/ipe 40 dict def ipe begin
+/np { newpath } def
+/m { moveto } def
+/l { lineto } def
+/c { curveto } def
+/h { closepath } def
+/re { 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto
+      neg 0 rlineto closepath } def
+/d { setdash } def
+/w { setlinewidth } def
+/J { setlinecap } def
+/j { setlinejoin } def
+/cm { [ 7 1 roll ] concat } def
+/q { gsave } def
+/Q { grestore } def
+/g { setgray } def
+/G { setgray } def
+/rg { setrgbcolor } def
+/RG { setrgbcolor } def
+/S { stroke } def
+/f* { eofill } def
+/f { fill } def
+/ipeMakeFont {
+  exch findfont
+  dup length dict begin
+    { 1 index /FID ne { def } { pop pop } ifelse } forall
+    /Encoding exch def
+    currentdict
+  end
+  definefont pop
+} def
+/ipeFontSize 0 def
+/Tf { dup /ipeFontSize exch store selectfont } def
+/Td { translate } def
+/BT { gsave } def
+/ET { grestore } def
+/TJ { 0 0 moveto { dup type /stringtype eq
+ { show } { ipeFontSize mul -0.001 mul 0 rmoveto } ifelse
+} forall } def
+<<
+/PatternType 1
+/PaintType 2
+/TilingType 2
+/BBox [ 0 0 100 4 ]
+/XStep 99
+/YStep 4
+/PaintProc { pop 0 0 100 1 re fill} bind
+>>
+[ 0.5 -0.866025 0.866025 0.5 0 0 ]
+makepattern
+/Pat103 exch def
+<<
+/PatternType 1
+/PaintType 2
+/TilingType 2
+/BBox [ 0 0 100 4 ]
+/XStep 99
+/YStep 4
+/PaintProc { pop 0 0 100 1 re fill} bind
+>>
+[ 0.866025 0.5 -0.5 0.866025 0 0 ]
+makepattern
+/Pat104 exch def
+/patg { [/Pattern /DeviceGray ] setcolorspace setcolor } def
+/patrg { [/Pattern /DeviceRGB ] setcolorspace setcolor } def
+end
+%%EndResource
+%%EndProlog
+%%BeginSetup
+ipe begin
+%%EndSetup
+0 J 1 j
+q 128 768 m
+384 768 l
+384 704 l
+128 704 l
+128 768 l
+0.4 w S
+Q
+q q 161.8 736 m
+161.8 736.994 160.994 737.8 160 737.8 c
+159.006 737.8 158.2 736.994 158.2 736 c
+158.2 735.006 159.006 734.2 160 734.2 c
+160.994 734.2 161.8 735.006 161.8 736 c
+f*
+Q
+Q
+q q 225.8 736 m
+225.8 736.994 224.994 737.8 224 737.8 c
+223.006 737.8 222.2 736.994 222.2 736 c
+222.2 735.006 223.006 734.2 224 734.2 c
+224.994 734.2 225.8 735.006 225.8 736 c
+f*
+Q
+Q
+q q 289.8 736 m
+289.8 736.994 288.994 737.8 288 737.8 c
+287.006 737.8 286.2 736.994 286.2 736 c
+286.2 735.006 287.006 734.2 288 734.2 c
+288.994 734.2 289.8 735.006 289.8 736 c
+f*
+Q
+Q
+q q 353.8 736 m
+353.8 736.994 352.994 737.8 352 737.8 c
+351.006 737.8 350.2 736.994 350.2 736 c
+350.2 735.006 351.006 734.2 352 734.2 c
+352.994 734.2 353.8 735.006 353.8 736 c
+f*
+Q
+Q
+q q 190.2 734.2 m
+193.8 734.2 l
+193.8 737.8 l
+190.2 737.8 l
+h f*
+Q
+Q
+q q 254.2 734.2 m
+257.8 734.2 l
+257.8 737.8 l
+254.2 737.8 l
+h f*
+Q
+Q
+q q 318.2 734.2 m
+321.8 734.2 l
+321.8 737.8 l
+318.2 737.8 l
+h f*
+Q
+Q
+q q 126.5 734.5 m
+129.5 734.5 l
+129.5 737.5 l
+126.5 737.5 l
+h 1 g
+f*
+Q
+q 126.2 734.2 m
+129.8 734.2 l
+129.8 737.8 l
+126.2 737.8 l
+h 126.8 734.8 m
+129.2 734.8 l
+129.2 737.2 l
+126.8 737.2 l
+h f*
+Q
+Q
+q q 382.5 734.5 m
+385.5 734.5 l
+385.5 737.5 l
+382.5 737.5 l
+h 1 g
+f*
+Q
+q 382.2 734.2 m
+385.8 734.2 l
+385.8 737.8 l
+382.2 737.8 l
+h 382.8 734.8 m
+385.2 734.8 l
+385.2 737.2 l
+382.8 737.2 l
+h f*
+Q
+Q
+q 192 704 m
+192 768 l
+[4] 0 d 0.4 w S
+Q
+q 256 704 m
+256 768 l
+[4] 0 d 0.4 w S
+Q
+q 320 704 m
+320 768 l
+[4] 0 d 0.4 w S
+Q
+showpage
+%%BeginIpeXml: /FlateDecode
+%GhUDE9on$e&A@sB&3;*BgCj0j(YJcl=u6D'lF#@KD::DJPc.FLm^b^]gH-Z,RC&oomHbV!fYc\.
+%j5Y)HNq@m&'Ko>s1K!!4Go26\JC]mH7``4"GNd7`d;0DBS:"`+<gYMN\I3l^$&;@rVVs:8i\6>n
+%</;"<195>M(H@:lMM&GX?I?>EEf*Ied(<[SC:70Mos%U*hKk$Cib;0KqH[Ai^K95+C%eKCUj.E-
+%^eTk0TIo3Zoi/:T%PV:(LdG'7V@a/JE?$4+b\:_Y-QhK.)J5?U"b%oT/$qG+=n>*<WogEm'$_9N
+%].l9Hags;k[j`e_I$h+*>]]bE#!JoIUo_$l(8QN:Y3uGZi25&J*CL;Mbqc(oFoN.<8aArOcUDS,
+%XF>D:/Jl[E2,eEHRX6F$Tc.ZBQ<'*U$IWu!IaQjo%ihrQd?@fKdSA/u<k&dgh\B=eb,LQ%Iu6g9
+%E[+mpe!%=DN<JL<@a^b$/509sPBgU5emAnTIW,'h[u/AWW>"tD#rP_NUjn7nDfV%Sg7"VWr@Q.b
+%s(p^A8\14^XtSSUOeCC^R`%6h\HuN$Zq!cKO`m,Ej_3BA]GSACen1K,Zr;a(8Gc:%g3]Pn$g0ak
+%4XsP?ji\0ZP?oX'D=:AtM6i?Xq6H_N&cM?Uf&P&(@5Kn%YYt5j)s3(D;#U1=O0=VlF[C""<c:eB
+%]WK[H-'+#-!'LsDJ2S[,Y/jCp/SY(<BSOiZ%ic'0)e&eZVOIJ>Q1UTf*UQ.l+&p-D^gCPc=F%7V
+%72i!I4uus`H9pt@W>FU=(2LupE$<tR8N3PM<3\2Dlm.9Zm)+`*%qYun('Hl)aH:KqSufJM49Fp[
+%HeI[.96-XQ_&_@<Kf\E&Yj*>]DCM&d=<7#G<d4$gM*-\B7gkn9l#<a9hpSAq1hDuAhZ:G:O%1QU
+%E8-%j.$^&.cAM0l?![7])aJ12c8K5L/JN&r1rPts2$`s=+BWK9!iJ(s\Ji#7L9iut_.dBg=N:1i
+%<MR>Li&=IH%.ZZAnj:qu@V8]=ILL\60j,9Q^t`ISE%J7A0bh&=4'FP\@au+*+i9G7BsJ_P\"Ban
+%H@Y'm[r:'+bqcUQ3j4j@<7,42q!!3c4HH&QScs%`o.cc)<N/aYG2GH^XU"YCCaDTH2l,<b+%B5a
+%MT_X6?YrI@5g_DH0LG:X'\_ajT-+P/n4R92`n=i@qV7%$bkJ"!_7>?m]kMe3c2i?s>9;:,GLk)M
+%Z[VoLO5GDoE4Y;Mm12+GDL#'N*7Ukue5hC&'TIh54:60K@E!^*:S,/)=*2KU<4LQ[?]FfP!7>fh
+%:AF3^4q]9&$8oX;9C$b7HJ]Rc(s0HH^5d:%,hlFp85RjdI2ikAA>(.0h7lcR(@QV6_?Z(+(7R>:
+%p96^_a1):CfKmR19-&12e@V_T*r=g)KBE^]HbiS-K$TgtiP8_]5`nG$NC_bl"M<mIKER<J:+c5r
+%6lW`';9R=Eb$7p.8&54DUY76X*#6(jZ(DEaREMHI-8ugTRI'dmQ)dr9o_?%+$IYRKVj0k+oY`Z4
+%*l9rX-Wt9sMap:^]C4=^]\f\Q:6)\Ch[-L[6j'.$IR:WHmbEW")trKgoG`BOoq:`$@E$@b!F>pR
+%f^Cq,5e.si&oDVuL$BV1mQF^a\$,cA2>)HLX@003M@elh"uAT5U$C&Zjlc6]?_kLtGnldei\bE/
+%']2;gIqpjf=hYL\L*uIjLgkbdI^u*:bi8*q_00._0B:T.Gr;dS#R]!Np>!A[5r"Ke75l'_MNGOt
+%Jk%%?O$R2*2u1.g;2!h[2UB+$G2=8el?rCeq%"FO;`k~>
+%%EndIpeXml
+%%Trailer
+end
+%%EOF

doc/gas_spring.lyx

@@ -0,0 +1,108 @@
+#LyX 2.1 created this file. For more info see http://www.lyx.org/
+\lyxformat 474
+\begin_document
+\begin_header
+\textclass article
+\use_default_options true
+\maintain_unincluded_children false
+\language english
+\language_package default
+\inputencoding auto
+\fontencoding global
+\font_roman default
+\font_sans default
+\font_typewriter default
+\font_math auto
+\font_default_family default
+\use_non_tex_fonts false
+\font_sc false
+\font_osf false
+\font_sf_scale 100
+\font_tt_scale 100
+\graphics default
+\default_output_format default
+\output_sync 0
+\bibtex_command default
+\index_command default
+\paperfontsize default
+\use_hyperref false
+\papersize default
+\use_geometry false
+\use_package amsmath 1
+\use_package amssymb 1
+\use_package cancel 1
+\use_package esint 1
+\use_package mathdots 1
+\use_package mathtools 1
+\use_package mhchem 1
+\use_package stackrel 1
+\use_package stmaryrd 1
+\use_package undertilde 1
+\cite_engine basic
+\cite_engine_type default
+\biblio_style plain
+\use_bibtopic false
+\use_indices false
+\paperorientation portrait
+\suppress_date false
+\justification true
+\use_refstyle 1
+\index Index
+\shortcut idx
+\color #008000
+\end_index
+\secnumdepth 3
+\tocdepth 3
+\paragraph_separation indent
+\paragraph_indentation default
+\quotes_language english
+\papercolumns 1
+\papersides 1
+\paperpagestyle default
+\tracking_changes false
+\output_changes false
+\html_math_output 0
+\html_css_as_file 0
+\html_be_strict false
+\end_header
+
+\begin_body
+
+\begin_layout Standard
+\begin_inset Formula $K=\frac{\partial F}{\partial x}=S_{l}\frac{\partial p}{\partial x}=S_{l}\frac{\partial p}{\partial\rho}\frac{\partial\rho}{\partial x}=S_{l}c_{0}^{2}$
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\begin_inset Formula $\rho=\rho_{0}\frac{V_{0}}{V}$
+\end_inset
+
+, 
+\begin_inset Formula $V=V_{0}+S_{l}x$
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\begin_inset Formula $\frac{\partial\rho}{\partial x}=-S_{l}\frac{\rho_{0}}{V_{0}}$
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+Such that
+\end_layout
+
+\begin_layout Standard
+\begin_inset Formula $K=-S_{l}^{2}\frac{\rho_{0}c_{0}^{2}}{V_{0}}=-S_{l}^{2}\frac{p_{0}}{\gamma V_{0}}$
+\end_inset
+
+
+\end_layout
+
+\end_body
+\end_document