📄 cfd recipes chapter 18.mht
字号:
From: <由 Microsoft Internet Explorer 5 保存>
Subject: CFD Recipes: Chapter 18
Date: Fri, 18 Jul 2003 17:28:50 +0800
MIME-Version: 1.0
Content-Type: multipart/related;
boundary="----=_NextPart_000_003F_01C34D52.0DB5A150";
type="text/html"
X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2600.0000
This is a multi-part message in MIME format.
------=_NextPart_000_003F_01C34D52.0DB5A150
Content-Type: text/html;
charset="gb2312"
Content-Transfer-Encoding: quoted-printable
Content-Location: http://capella.colorado.edu/~laney/ch18code.htm
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD><TITLE>CFD Recipes: Chapter 18</TITLE>
<META http-equiv=3DContent-Type content=3D"text/html; charset=3Dgb2312">
<META content=3D"Free CFD Software" name=3DDESCRIPTION>
<META content=3D"MSHTML 6.00.2600.0" name=3DGENERATOR></HEAD>
<BODY text=3D#ffffff vLink=3D#cc9999 aLink=3D#ffcccc link=3D#999999 =
bgColor=3D#006666>
<DIV align=3Dcenter><IMG height=3D42 alt=3D"Cambridge University Press"=20
src=3D"http://capella.colorado.edu/~laney/graphics/cup.gif" width=3D189> =
<BR><IMG=20
height=3D79 alt=3D"Computational Gasdynamics"=20
src=3D"http://capella.colorado.edu/~laney/graphics/title.gif" =
width=3D537>=20
<BR></DIV>
<TABLE cellSpacing=3D0 cellPadding=3D10 width=3D"100%" border=3D0>
<TBODY>
<TR>
<TD vAlign=3Dtop align=3Dright><FONT size=3D+0><A=20
href=3D"http://capella.colorado.edu/~laney/home.htm">Home</A>=20
<P><A=20
=
href=3D"http://capella.colorado.edu/~laney/description.htm">Description</=
A>=20
<P><A href=3D"http://capella.colorado.edu/~laney/faq.htm">FAQ</A>=20
<P><A href=3D"http://capella.colorado.edu/~laney/toc.htm">Table of =
Contents</A>=20
<P><A =
href=3D"http://capella.colorado.edu/~laney/index.htm">Index</A>=20
<P><A =
href=3D"http://capella.colorado.edu/~laney/reviews.htm">Reviews</A>=20
<P><A =
href=3D"http://capella.colorado.edu/~laney/software.htm">Software</A>=20
<P><A =
href=3D"http://capella.colorado.edu/~laney/hwsolns.htm">Solutions to=20
Problems</A>=20
<P><A =
href=3D"http://capella.colorado.edu/~laney/additional.htm">Additional=20
Material</A>=20
<P><A =
href=3D"http://capella.colorado.edu/~laney/errata.htm">Errata </A>
<P><A =
href=3D"http://capella.colorado.edu/~laney/cover.htm">Cover</A>=20
<P><A =
href=3D"http://capella.colorado.edu/~laney/otherbooks.htm">Other Books=20
of Interest</A>=20
<P><A href=3D"http://capella.colorado.edu/~laney/links.htm">Other =
Sites of=20
Interest</A>=20
<P><A href=3D"http://capella.colorado.edu/~laney/author.htm">About =
the=20
Author</A>=20
<P><A =
href=3D"http://capella.colorado.edu/~laney/purchase.htm">Purchase the=20
Book</A> </FONT></P></TD>
<TD vAlign=3Dtop>
<BLOCKQUOTE>
<H1>Software for <I>Computational Gasdynamics</I></H1><BR><IMG=20
=
src=3D"http://capella.colorado.edu/~laney/graphics/cfdrecipes.jpg"=20
border=3D0> <BR><BR>
<HR>
<H2>Chapter 18. Basic Numerical Methods for the Euler =
Equations</H2>
<HR>
<BR>
<TABLE cellSpacing=3D0 cellPadding=3D5 border=3D0>
<TBODY>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/lf.f">Lax-Friedrichs =
Method (18.1)</A> </TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/lw.f">Lax-Wendroff=20
Methods (MacCormack and Richtmyer) (18.1)</A> </TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/sw.f">Steger-Warming =
Flux Split First-Order Upwind Method (18.2)</A> </TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/vanleer.f">Van=20
Leer Flux Split First-Order Upwind Method (18.2)</A> =
</TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/ausm.f">Liou-Steffen =
Flux Split First-Order Upwind Method (AUSM) (18.2)</A> =
</TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/zb.f">Zha-Bilgen=20
Flux Split First-Order Upwind Method (18.2)</A> </TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/bw.f">Beam-Warming=20
Second-Order Upwind Method w/Three Options for Flux Vector =
Splitting(18.2)</A> </TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/godunov.f">Godunov's =
First-Order Upwind Method (18.3)</A> </TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/roefou.f">Roe's=20
First-Order Upwind Method (18.3)</A> </TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/onewave.f">First-Orde=
r=20
Upwind Method Based on One-Wave Solver (18.3)</A> =
</TD></TR>
<TR>
<TD><A=20
=
href=3D"http://capella.colorado.edu/~laney/ch18soft/input.dat">Sample=20
Input File</A> </TD></TR></TBODY></TABLE>
<P>These methods approximate the solution to the Riemann =
problem. See <A=20
href=3D"http://capella.colorado.edu/~laney/ch5code.htm">Chapter =
5</A> for=20
the exact solution. These codes require an ASCII input file =
called=20
input.dat of the following form:
<P>
<TABLE cellSpacing=3D0 cellPadding=3D10 border=3D0>
<TBODY>
<TR>
<TD><I>Lower Limit on x </I></TD>
<TD><I>Upper Limit on x </I></TD>
<TD><I>Final Time </I></TD></TR>
<TR>
<TD><I>Pressure for x<0 when t=3D0</I> </TD>
<TD><I>Density for x<0 when t=3D0</I> </TD>
<TD><I>Speed for x<0 when t=3D0</I> </TD></TR>
<TR>
<TD><I>Pressure for x>0 when t=3D0</I> </TD>
<TD><I>Density for x>0 when t=3D0</I> </TD>
<TD><I>Speed for x>0 when t=3D0</I> =
</TD></TR></TBODY></TABLE>
<P>These codes produce 8 ASCII output files:=20
<UL>
<LI>density.out. Density vs. <I>x</I>=20
<LI>entropy.out. Entropy vs. <I>x</I>=20
<LI>mach.out. Mach number vs. <I>x</I>=20
<LI>massflux.out. Mass flux vs. <I>x</I>=20
<LI>pressure.out. Pressure vs. <I>x</I>=20
<LI>sound.out. Speed-of-sound vs. <I>x</I>=20
<LI>velocity.out. Velocity vs. <I>x</I>=20
<LI>spectral.out. Spectral radius vs. <I>x</I> (the spectral =
radius is=20
the largest wave speed in absolute value) </LI></UL>Using the =
sample=20
input file, these codes will generate most of the figures in =
Chapter 18=20
related to Test Case 1. With a slightly different input file, =
and=20
appropriate changes in the settings for each code, these codes =
will also=20
generate most of the figures in Chapter 18 related to Test Case =
2. The=20
book describes the conditions used to generate each figure in =
precise=20
detail. If you want to duplicate the pictures in Chapter 18 =
exactly, you=20
must make sure to set all parameters correctly including the =
number of=20
appoximation points, the initial CFL number, the coefficient of=20
artificial viscosity (if any), the correct variation on the =
method=20
(where variations exist), the correct smoothing of the initial=20
conditions (where this is an option), and so forth. Most of =
these are=20
controlled by FORTRAN parameters. Make sure to study the source =
code in=20
detail so that you understand the various options that it =
offers.=20
</BLOCKQUOTE></TD></TR></TBODY></TABLE><BR>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -