Compare commits

18 Commits
p1 ... p3

Author SHA1 Message Date
562a5a713c fix tables 2026-01-28 16:21:23 +08:00
47e5f68759 update sa 2026-01-28 15:50:08 +08:00
a1d41b049c rename 2026-01-21 14:26:00 +08:00
6a178db069 clean 2026-01-21 12:20:17 +08:00
1498209b0e latex 2026-01-20 03:45:15 +08:00
101e8b7e8e update task3 2026-01-20 01:55:46 +08:00
d49a10a042 update readmes 2026-01-19 22:57:22 +08:00
b4a307e5a8 update readme 2026-01-19 22:18:35 +08:00
a14695cea8 update readme 2026-01-19 20:48:58 +08:00
a01123b5b5 fix colors 2026-01-19 19:43:57 +08:00
b98ba91fc1 fix color 2026-01-19 19:38:38 +08:00
22003ea65e P1: fixs 2026-01-19 15:35:56 +08:00
86daefbea0 P3: visual 2026-01-19 15:05:55 +08:00
30a57d6012 P3: readme fix1 2026-01-19 14:32:53 +08:00
73ae99885a P3: sens figure 2026-01-19 14:15:51 +08:00
d4d7f7d94b P3: sens and readme 2026-01-19 12:59:03 +08:00
ffc51b2927 P3: main 2026-01-19 11:57:19 +08:00
c6e22729c7 P3: readme 2026-01-19 11:31:46 +08:00
171 changed files with 17063 additions and 127 deletions

2
.gitignore vendored Normal file
View File

@@ -0,0 +1,2 @@
__pycache__
.DS_Store

BIN
last.pdf Normal file

Binary file not shown.

View File

@@ -0,0 +1,7 @@
auto_roxygenize_for_build_and_reload="0"
auto_roxygenize_for_build_package="1"
auto_roxygenize_for_check="1"
live_preview_website="1"
makefile_args=""
preview_website="1"
website_output_format="all"

View File

@@ -0,0 +1,9 @@
{
"sortOrder": [
{
"columnIndex": 2,
"ascending": true
}
],
"path": "C:/Users/Zhang/Desktop/21美赛/美赛模板/MCM_Latex模板(2021最新修改)"
}

View File

@@ -0,0 +1,3 @@
{
"activeTab": 1
}

View File

@@ -0,0 +1,14 @@
{
"left": {
"splitterpos": 409,
"topwindowstate": "NORMAL",
"panelheight": 848,
"windowheight": 886
},
"right": {
"splitterpos": 531,
"topwindowstate": "NORMAL",
"panelheight": 848,
"windowheight": 886
}
}

View File

@@ -0,0 +1,5 @@
{
"TabSet1": 0,
"TabSet2": 0,
"TabZoom": {}
}

View File

@@ -0,0 +1,5 @@

View File

@@ -0,0 +1 @@
{"active_set":"","sets":[]}

View File

@@ -0,0 +1,26 @@
{
"id": "461C3D72",
"path": "C:/Users/Zhang/Desktop/21美赛/美赛模板/MCM_Latex模板(2021最新修改)/mcmthesis-demo.tex",
"project_path": "mcmthesis-demo.tex",
"type": "tex",
"hash": "657170375",
"contents": "",
"dirty": false,
"created": 1611800986918.0,
"source_on_save": false,
"relative_order": 1,
"properties": {
"source_window_id": "",
"Source": "Source",
"cursorPosition": "80,4",
"scrollLine": "0"
},
"folds": "",
"lastKnownWriteTime": 1612445849,
"encoding": "UTF-8",
"collab_server": "",
"source_window": "",
"last_content_update": 1612445849571,
"read_only": false,
"read_only_alternatives": []
}

View File

@@ -0,0 +1,539 @@
%%
%% This is file `mcmthesis-demo.tex',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `demo')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
%%
%% This is file `mcmthesis-demo.tex',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `demo')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
\documentclass{mcmthesis}
\mcmsetup{CTeX = false, % 使用 CTeX 套装时,设置为 true
tcn = {0000000}, problem = {A},
sheet = true, titleinsheet = true, keywordsinsheet = true,
titlepage = false, abstract = false}
\usepackage{newtxtext} % \usepackage{palatino}
\usepackage[backend=bibtex]{biblatex} % for RStudio Complie
\usepackage{tocloft}
\setlength{\cftbeforesecskip}{6pt}
\renewcommand{\contentsname}{\hspace*{\fill}\Large\bfseries Contents \hspace*{\fill}}
\title{Enjoy a Cozy and Green Bath}
% \author{\small \href{http://www.latexstudio.net/}
% {\includegraphics[width=7cm]{mcmthesis-logo}}}
\date{\today}
\begin{document}
\begin{abstract}
A traditional bathtub cannot be reheated by itself, so users have to add hot water from time to time. Our goal is to establish a model of the temperature of bath water in space and time. Then we are expected to propose an optimal strategy for users to keep the temperature even and close to initial temperature and decrease water consumption.
To simplify modeling process, we firstly assume there is no person in the bathtub. We regard the whole bathtub as a thermodynamic system and introduce heat transfer formulas.
We establish two sub-models: adding water constantly and discontinuously. As for the former sub-model, we define the mean temperature of bath water. Introducing Newton cooling formula, we determine the heat transfer capacity. After deriving the value of parameters, we deduce formulas to derive results and simulate the change of temperature field via CFD. As for the second sub-model, we define an iteration consisting of two process: heating and standby. According to energy conservation law, we obtain the relationship of time and total heat dissipating capacity. Then we determine the mass flow and the time of adding hot water. We also use CFD to simulate the temperature field in second sub-model.
In consideration of evaporation, we correct the results of sub-models referring to some scientists' studies. We define two evaluation criteria and compare the two sub-models. Adding water constantly is found to keep the temperature of bath water even and avoid wasting too much water, so it is recommended by us.
Then we determine the influence of some factors: radiation heat transfer, the shape and volume of the tub, the shape/volume/temperature/motions of the person, the bubbles made from bubble bath additives. We focus on the influence of those factors to heat transfer and then conduct sensitivity analysis. The results indicate smaller bathtub with less surface area, lighter personal mass, less motions and more bubbles will decrease heat transfer and save water.
Based on our model analysis and conclusions, we propose the optimal strategy for the user in a bathtub and explain the reason of uneven temperature throughout the bathtub. In addition, we make improvement for applying our model in real life.
\begin{keywords}
Heat transfer, Thermodynamic system, CFD, Energy conservation
\end{keywords}
\end{abstract}
\maketitle
%% Generate the Table of Contents, if it's needed.
% \renewcommand{\contentsname}{\centering Contents}
\tableofcontents % 若不想要目录, 注释掉该句
\newpage
\section{Introduction}
\subsection{Background}
Bathing in a tub is a perfect choice for those who have been worn out after a long day's working. A traditional bathtub is a simply water containment vessel without a secondary heating system or circulating jets. Thus the temperature of water in bathtub declines noticeably as time goes by, which will influent the experience of bathing. As a result, the bathing person needs to add a constant trickle of hot water from a faucet to reheat the bathing water. This way of bathing will result in waste of water because when the capacity of the bathtub is reached, excess water overflows the tub.
An optimal bathing strategy is required for the person in a bathtub to get comfortable bathing experience while reducing the waste of water.
\subsection{Literature Review}
Korean physicist Gi-Beum Kim analyzed heat loss through free surface of water contained in bathtub due to conduction and evaporation \cite{1}. He derived a relational equation based on the basic theory of heat transfer to evaluate the performance of bath tubes. The major heat loss was found to be due to evaporation. Moreover, he found out that the speed of heat loss depends more on the humidity of the bathroom than the temperature of water contained in the bathtub. So, it is best to maintain the temperature of bathtub water to be between 41 to 45$^{\circ}$C and the humidity of bathroom to be 95\%.
When it comes to convective heat transfer in bathtub, many studies can be referred to. Newton's law of cooling states that the rate of heat loss of a body is proportional to the difference in temperatures between the body and its surroundings while under the effects of a breeze \cite{2}. Claude-Louis Navier and George Gabriel Stokes described the motion of viscous fluid substances with the Navier-Stokes equations. Those equations may be used to model the weather, ocean currents, water flow in a pipe and air flow around a wing \cite{3}.
In addition, some numerical simulation software are applied in solving and analyzing problems that involve fluid flows. For example, Computational Fluid Dynamics (CFD) is a common one used to perform the calculations required to simulate the interaction of liquids and gases with surfaces defined by boundary conditions \cite{4}.
\subsection{Restatement of the Problem}
We are required to establish a model to determine the change of water temperature in space and time. Then we are expected to propose the best strategy for the person in the bathtub to keep the water temperature close to initial temperature and even throughout the tub. Reduction of waste of water is also needed. In addition, we have to consider the impact of different conditions on our model, such as different shapes and volumes of the bathtub, etc.
In order to solve those problems, we will proceed as follows:
\begin{itemize}
\item {\bf Stating assumptions}. By stating our assumptions, we will narrow the focus of our approach towards the problems and provide some insight into bathtub water temperature issues.
\item {\bf Making notations}. We will give some notations which are important for us to clarify our models.
\item {\bf Presenting our model}. In order to investigate the problem deeper, we divide our model into two sub-models. One is a steady convection heat transfer sub-model in which hot water is added constantly. The other one is an unsteady convection heat transfer sub-model where hot water is added discontinuously.
\item {Defining evaluation criteria and comparing sub-models}. We define two main criteria to evaluate our model: the mean temperature of bath water and the amount of inflow water.
\item {\bf Analysis of influencing factors}. In term of the impact of different factors on our model, we take those into consideration: the shape and volume of the tub, the shape/volume/temperature of the person in the bathtub, the motions made by the person in the bathtub and adding a bubble bath additive initially.
\item {\bf Model testing and sensitivity analysis}. With the criteria defined before, we evaluate the reliability of our model and do the sensitivity analysis.
\item {\bf Further discussion}. We discuss about different ways to arrange inflow faucets. Then we improve our model to apply them in reality.
\item {\bf Evaluating the model}. We discuss about the strengths and weaknesses of our model:
\begin{itemize}
\item[1)] ...
\item[2)] ...
\item[3)] ...
\item[4)] ...
\end{itemize}
\end{itemize}
\section{Assumptions and Justification}
To simplify the problem and make it convenient for us to simulate real-life conditions, we make the following basic assumptions, each of which is properly justified.
\begin{itemize}
\item {\bf The bath water is incompressible Non-Newtonian fluid}. The incompressible Non-Newtonian fluid is the basis of NavierStokes equations which are introduced to simulate the flow of bath water.
\item {\bf All the physical properties of bath water, bathtub and air are assumed to be stable}. The change of those properties like specific heat, thermal conductivity and density is rather small according to some studies \cite{5}. It is complicated and unnecessary to consider these little change so we ignore them.
\item {\bf There is no internal heat source in the system consisting of bathtub, hot water and air}. Before the person lies in the bathtub, no internal heat source exist except the system components. The circumstance where the person is in the bathtub will be investigated in our later discussion.
\item {\bf We ignore radiative thermal exchange}. According to Stefan-Boltzmanns law, the radiative thermal exchange can be ignored when the temperature is low. Refer to industrial standard \cite{6}, the temperature in bathroom is lower than 100 $^{\circ}$C, so it is reasonable for us to make this assumption.
\item {\bf The temperature of the adding hot water from the faucet is stable}. This hypothesis can be easily achieved in reality and will simplify our process of solving the problem.
\end{itemize}
\section{Notations}
\begin{center}
\begin{tabular}{clc}
{\bf Symbols} & {\bf Description} & \quad {\bf Unit} \\[0.25cm]
$h$ & Convection heat transfer coefficient & \quad W/(m$^2 \cdot$ K)
\\[0.2cm]
$k$ & Thermal conductivity & \quad W/(m $\cdot$ K) \\[0.2cm]
$c_p$ & Specific heat & \quad J/(kg $\cdot$ K) \\[0.2cm]
$\rho$ & Density & \quad kg/m$^2$ \\[0.2cm]
$\delta$ & Thickness & \quad m \\[0.2cm]
$t$ & Temperature & \quad $^\circ$C, K \\[0.2cm]
$\tau$ & Time & \quad s, min, h \\[0.2cm]
$q_m$ & Mass flow & \quad kg/s \\[0.2cm]
$\Phi$ & Heat transfer power & \quad W \\[0.2cm]
$T$ & A period of time & \quad s, min, h \\[0.2cm]
$V$ & Volume & \quad m$^3$, L \\[0.2cm]
$M,\,m$ & Mass & \quad kg \\[0.2cm]
$A$ & Aera & \quad m$^2$ \\[0.2cm]
$a,\,b,\,c$ & The size of a bathtub & \quad m$^3$
\end{tabular}
\end{center}
\noindent where we define the main parameters while specific value of those parameters will be given later.
\section{Model Overview}
In our basic model, we aim at three goals: keeping the temperature as even as possible, making it close to the initial temperature and decreasing the water consumption.
We start with the simple sub-model where hot water is added constantly.
At first we introduce convection heat transfer control equations in rectangular coordinate system. Then we define the mean temperature of bath water.
Afterwards, we introduce Newton cooling formula to determine heat transfer
capacity. After deriving the value of parameters, we get calculating results via formula deduction and simulating results via CFD.
Secondly, we present the complicated sub-model in which hot water is
added discontinuously. We define an iteration consisting of two process:
heating and standby. As for heating process, we derive control equations and boundary conditions. As for standby process, considering energy conservation law, we deduce the relationship of total heat dissipating capacity and time.
Then we determine the time and amount of added hot water. After deriving the value of parameters, we get calculating results via formula deduction and simulating results via CFD.
At last, we define two criteria to evaluate those two ways of adding hot water. Then we propose optimal strategy for the user in a bathtub.
The whole modeling process can be shown as follows.
\begin{figure}[h]
\centering
\includegraphics[width=12cm]{fig1.jpg}
\caption{Modeling process} \label{fig1}
\end{figure}
\section{Sub-model I : Adding Water Continuously}
We first establish the sub-model based on the condition that a person add water continuously to reheat the bathing water. Then we use Computational Fluid Dynamics (CFD) to simulate the change of water temperature in the bathtub. At last, we evaluate the model with the criteria which have been defined before.
\subsection{Model Establishment}
Since we try to keep the temperature of the hot water in bathtub to be even, we have to derive the amount of inflow water and the energy dissipated by the hot water into the air.
We derive the basic convection heat transfer control equations based on the former scientists achievement. Then, we define the mean temperature of bath water. Afterwards, we determine two types of heat transfer: the boundary heat transfer and the evaporation heat transfer. Combining thermodynamic formulas, we derive calculating results. Via Fluent software, we get simulation results.
\subsubsection{Control Equations and Boundary Conditions}
According to thermodynamics knowledge, we recall on basic convection
heat transfer control equations in rectangular coordinate system. Those
equations show the relationship of the temperature of the bathtub water in space.
We assume the hot water in the bathtub as a cube. Then we put it into a
rectangular coordinate system. The length, width, and height of it is $a,\, b$ and $c$.
\begin{figure}[h]
\centering
\includegraphics[width=8cm]{fig2.jpg}
\caption{Modeling process} \label{fig2}
\end{figure}
In the basis of this, we introduce the following equations \cite{5}:
\begin{itemize}
\item {\bf Continuity equation:}
\end{itemize}
\begin{equation} \label{eq1}
\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} +\frac{\partial w}{\partial z} =0
\end{equation}
\noindent where the first component is the change of fluid mass along the $X$-ray. The second component is the change of fluid mass along the $Y$-ray. And the third component is the change of fluid mass along the $Z$-ray. The sum of the change in mass along those three directions is zero.
\begin{itemize}
\item {\bf Moment differential equation (N-S equations):}
\end{itemize}
\begin{equation} \label{eq2}
\left\{
\begin{array}{l} \!\!
\rho \Big(u \dfrac{\partial u}{\partial x} + v \dfrac{\partial u}{\partial y} + w\dfrac{\partial u}{\partial z} \Big) = -\dfrac{\partial p}{\partial x} + \eta \Big(\dfrac{\partial^2 u}{\partial x^2} + \dfrac{\partial^2 u}{\partial y^2} + \dfrac{\partial^2 u}{\partial z^2} \Big) \\[0.3cm]
\rho \Big(u \dfrac{\partial v}{\partial x} + v \dfrac{\partial v}{\partial y} + w\dfrac{\partial v}{\partial z} \Big) = -\dfrac{\partial p}{\partial y} + \eta \Big(\dfrac{\partial^2 v}{\partial x^2} + \dfrac{\partial^2 v}{\partial y^2} + \dfrac{\partial^2 v}{\partial z^2} \Big) \\[0.3cm]
\rho \Big(u \dfrac{\partial w}{\partial x} + v \dfrac{\partial w}{\partial y} + w\dfrac{\partial w}{\partial z} \Big) = -g-\dfrac{\partial p}{\partial z} + \eta \Big(\dfrac{\partial^2 w}{\partial x^2} + \dfrac{\partial^2 w}{\partial y^2} + \dfrac{\partial^2 w}{\partial z^2} \Big)
\end{array}
\right.
\end{equation}
\begin{itemize}
\item {\bf Energy differential equation:}
\end{itemize}
\begin{equation} \label{eq3}
\rho c_p \Big( u\frac{\partial t}{\partial x} + v\frac{\partial t}{\partial y} + w\frac{\partial t}{\partial z} \Big) = \lambda \Big(\frac{\partial^2 t}{\partial x^2} + \frac{\partial^2 t}{\partial y^2} + \frac{\partial^2 t}{\partial z^2} \Big)
\end{equation}
\noindent where the left three components are convection terms while the right three components are conduction terms.
By Equation \eqref{eq3}, we have ......
......
On the right surface in Fig. \ref{fig2}, the water also transfers heat firstly with bathtub inner surfaces and then the heat comes into air. The boundary condition here is ......
\subsubsection{Definition of the Mean Temperature}
......
\subsubsection{Determination of Heat Transfer Capacity}
......
\section{Sub-model II: Adding Water Discontinuously}
In order to establish the unsteady sub-model, we recall on the working principle of air conditioners. The heating performance of air conditions consist of two processes: heating and standby. After the user set a temperature, the air conditioner will begin to heat until the expected temperature is reached. Then it will go standby. When the temperature get below the expected temperature, the air conditioner begin to work again. As it works in this circle, the temperature remains the expected one.
Inspired by this, we divide the bathtub working into two processes: adding
hot water until the expected temperature is reached, then keeping this
condition for a while unless the temperature is lower than a specific value. Iterating this circle ceaselessly will ensure the temperature kept relatively stable.
\subsection{Heating Model}
\subsubsection{Control Equations and Boundary Conditions}
\subsubsection{Determination of Inflow Time and Amount}
\subsection{Standby Model}
\subsection{Results}
\quad~ We first give the value of parameters based on others studies. Then we get the calculation results and simulating results via those data.
\subsubsection{Determination of Parameters}
After establishing the model, we have to determine the value of some
important parameters.
As scholar Beum Kim points out, the optimal temperature for bath is
between 41 and 45$^\circ$C [1]. Meanwhile, according to Shimodozono's study, 41$^\circ$C warm water bath is the perfect choice for individual health [2]. So it is reasonable for us to focus on $41^\circ$C $\sim 45^\circ$C. Because adding hot water continuously is a steady process, so the mean temperature of bath water is supposed to be constant. We value the temperature of inflow and outflow water with the maximum and minimum temperature respectively.
The values of all parameters needed are shown as follows:
.....
\subsubsection{Calculating Results}
Putting the above value of parameters into the equations we derived before, we can get the some data as follows:
%%普通表格
\begin{table}[h] %h表示固定在当前位置
\centering %设置居中
\caption{The calculating results} %表标题
\vspace{0.15cm}
\label{tab2} %设置表的引用标签
\begin{tabular}{|c|c|c|} %3个c表示3列, |可选, 表示绘制各列间的竖线
\hline %画横线
Variables & Values & Unit \\ \hline %各列间用&隔开
$A_1$ & 1.05 & $m^2$ \\ \hline
$A_2$ & 2.24 & $m^2$ \\ \hline
$\Phi_1$ & 189.00 & $W$ \\ \hline
$\Phi_2$ & 43.47 & $W$ \\ \hline
$\Phi$ & 232.47 & $W$ \\ \hline
$q_m$ & 0.014 & $g/s$ \\ \hline
\end{tabular}
\end{table}
From Table \ref{tab2}, ......
......
\section{Correction and Contrast of Sub-Models}
After establishing two basic sub-models, we have to correct them in consideration of evaporation heat transfer. Then we define two evaluation criteria to compare the two sub-models in order to determine the optimal bath strategy.
\subsection{Correction with Evaporation Heat Transfer}
Someone may confuse about the above results: why the mass flow in the first sub-model is so small? Why the standby time is so long? Actually, the above two sub-models are based on ideal conditions without consideration of the change of boundary conditions, the motions made by the person in bathtub and the evaporation of bath water, etc. The influence of personal motions will be discussed later. Here we introducing the evaporation of bath water to correct sub-models.
\subsection{Contrast of Two Sub-Models}
Firstly we define two evaluation criteria. Then we contrast the two submodels via these two criteria. Thus we can derive the best strategy for the person in the bathtub to adopt.
\section{Model Analysis and Sensitivity Analysis}
\subsection{The Influence of Different Bathtubs}
Definitely, the difference in shape and volume of the tub affects the
convection heat transfer. Examining the relationship between them can help
people choose optimal bathtubs.
\subsubsection{Different Volumes of Bathtubs}
In reality, a cup of water will be cooled down rapidly. However, it takes quite long time for a bucket of water to become cool. That is because their volume is different and the specific heat of water is very large. So that the decrease of temperature is not obvious if the volume of water is huge. That also explains why it takes 45 min for 320 L water to be cooled by 1$^\circ$C.
In order to examine the influence of volume, we analyze our sub-models
by conducting sensitivity Analysis to them.
We assume the initial volume to be 280 L and change it by $\pm 5$\%, $\pm 8$\%, $\pm 12$\% and $\pm 15$\%. With the aid of sub-models we established before, the variation of some parameters turns out to be as follows
%%三线表
\begin{table}[h] %h表示固定在当前位置
\centering %设置居中
\caption{Variation of some parameters} %表标题
\label{tab7} %设置表的引用标签
\begin{tabular}{ccccccc} %7个c表示7列, c表示每列居中对齐, 还有l和r可选
\toprule %画顶端横线
$V$ & $A_1$ & $A_2$ & $T_2$ & $q_{m1}$ & $q_{m2}$ & $\Phi_q$ \\
\midrule %画中间横线
-15.00\% & -5.06\% & -9.31\% & -12.67\% & -2.67\% & -14.14\% & -5.80\% \\
-12.00\% & -4.04\% & -7.43\% & -10.09\% & -2.13\% & -11.31\% & -4.63\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
\bottomrule %画底部横线
\end{tabular}
\end{table}
\section{Strength and Weakness}
\subsection{Strength}
\begin{itemize}
\item We analyze the problem based on thermodynamic formulas and laws, so that the model we established is of great validity.
\item Our model is fairly robust due to our careful corrections in consideration of real-life situations and detailed sensitivity analysis.
\item Via Fluent software, we simulate the time field of different areas throughout the bathtub. The outcome is vivid for us to understand the changing process.
\item We come up with various criteria to compare different situations, like water consumption and the time of adding hot water. Hence an overall comparison can be made according to these criteria.
\item Besides common factors, we still consider other factors, such as evaporation and radiation heat transfer. The evaporation turns out to be the main reason of heat loss, which corresponds with other scientists experimental outcome.
\end{itemize}
\subsection{Weakness}
\begin{itemize}
\item Having knowing the range of some parameters from others essays, we choose a value from them to apply in our model. Those values may not be reasonable in reality.
\item Although we investigate a lot in the influence of personal motions, they are so complicated that need to be studied further.
\item Limited to time, we do not conduct sensitivity analysis for the influence of personal surface area.
\end{itemize}
\section{Further Discussion}
In this part, we will focus on different distribution of inflow faucets. Then we discuss about the real-life application of our model.
\begin{itemize}
\item Different Distribution of Inflow Faucets
In our before discussion, we assume there being just one entrance of inflow.
From the simulating outcome, we find the temperature of bath water is hardly even. So we come up with the idea of adding more entrances.
The simulation turns out to be as follows
\begin{figure}[h]
\centering
\includegraphics[width=12cm]{fig24.jpg}
\caption{The simulation results of different ways of arranging entrances} \label{fig24}
\end{figure}
From the above figure, the more the entrances are, the evener the temperature will be. Recalling on the before simulation outcome, when there is only one entrance for inflow, the temperature of corners is quietly lower than the middle area.
In conclusion, if we design more entrances, it will be easier to realize the goal to keep temperature even throughout the bathtub.
\item Model Application
Our before discussion is based on ideal assumptions. In reality, we have to make some corrections and improvement.
\begin{itemize}
\item[1)] Adding hot water continually with the mass flow of 0.16 kg/s. This way can ensure even mean temperature throughout the bathtub and waste less water.
\item[2)] The manufacturers can design an intelligent control system to monitor the temperature so that users can get more enjoyable bath experience.
\item[3)] We recommend users to add bubble additives to slow down the water being cooler and help cleanse. The additives with lower thermal conductivity are optimal.
\item[4)] The study method of our establishing model can be applied in other area relative to convection heat transfer, such as air conditioners.
\end{itemize}
\end{itemize}
\begin{thebibliography}{99}
\addcontentsline{toc}{section}{Reference}
\bibitem{1} Gi-Beum Kim. Change of the Warm Water Temperature for the Development of Smart Healthecare Bathing System. Hwahak konghak. 2006, 44(3): 270-276.
\bibitem{2} \url{https://en.wikipedia.org/wiki/Convective_heat_transfer#Newton.27s_law_of_cooling}
\bibitem{3} \url{https://en.wikipedia.org/wiki/Navier\%E2\%80\%93Stokes_equations}
\bibitem{4} \url{https://en.wikipedia.org/wiki/Computational_fluid_dynamics}
\bibitem{5} Holman J P. Heat Transfer (9th ed.), New York: McGraw-Hill, 2002.
\bibitem{6} Liu Weiguo, Chen Zhaoping, ZhangYin. Matlab program design and application. Beijing: Higher education press, 2002. (In Chinese)
\end{thebibliography}
\newpage
\begin{letter}{Enjoy Your Bath Time!}
From simulation results of real-life situations, we find it takes a period of time for the inflow hot water to spread throughout the bathtub. During this process, the bath water continues transferring heat into air, bathtub and the person in bathtub. The difference between heat transfer capacity makes the temperature of various areas to be different. So that it is difficult to get an evenly maintained temperature throughout the bath water.
In order to enjoy a comfortable bath with even temperature of bath water and without wasting too much water, we propose the following suggestions.
\begin{itemize}
\item Adding hot water consistently
\item Using smaller bathtub if possible
\item Decreasing motions during bath
\item Using bubble bath additives
\item Arranging more faucets of inflow
\end{itemize}
\vspace{\parskip}
Sincerely yours,
Your friends
\end{letter}
\newpage
\begin{appendices}
\section{First appendix}
In addition, your report must include a letter to the Chief Financial Officer (CFO) of the Goodgrant Foundation, Mr. Alpha Chiang, that describes the optimal investment strategy, your modeling approach and major results, and a brief discussion of your proposed concept of a return-on-investment (ROI). This letter should be no more than two pages in length.
Here are simulation programmes we used in our model as follow.\\
\textbf{\textcolor[rgb]{0.98,0.00,0.00}{Input matlab source:}}
\lstinputlisting[language=Matlab]{./code/mcmthesis-matlab1.m}
\section{Second appendix}
some more text \textcolor[rgb]{0.98,0.00,0.00}{\textbf{Input C++ source:}}
\lstinputlisting[language=C++]{./code/mcmthesis-sudoku.cpp}
\end{appendices}
\end{document}
%%
%% This work consists of these files mcmthesis.dtx,
%% figures/ and
%% code/,
%% and the derived files mcmthesis.cls,
%% mcmthesis-demo.tex,
%% README,
%% LICENSE,
%% mcmthesis.pdf and
%% mcmthesis-demo.pdf.
%%
%% End of file `mcmthesis-demo.tex'.

View File

@@ -0,0 +1,26 @@
{
"id": "5BE5E0FB",
"path": "C:/Users/Zhang/Desktop/21美赛/美赛模板/MCM_Latex模板(2021最新修改)/mcmthesis.cls",
"project_path": "mcmthesis.cls",
"type": "tex",
"hash": "34984457",
"contents": "",
"dirty": false,
"created": 1612618018265.0,
"source_on_save": false,
"relative_order": 3,
"properties": {
"source_window_id": "",
"Source": "Source",
"cursorPosition": "61,0",
"scrollLine": "53"
},
"folds": "",
"lastKnownWriteTime": 1612445833,
"encoding": "UTF-8",
"collab_server": "",
"source_window": "",
"last_content_update": 1612445833,
"read_only": false,
"read_only_alternatives": []
}

View File

@@ -0,0 +1,350 @@
%%
%% This is file `mcmthesis.cls',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `class')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
\NeedsTeXFormat{LaTeX2e}[1999/12/01]
\ProvidesClass{mcmthesis}
[2020/01/18 v6.3 The Thesis Template Designed For MCM/ICM]
\typeout{The Thesis Template Designed For MCM/ICM}
\def\MCMversion{v6.3}
\RequirePackage{xkeyval}
\RequirePackage{etoolbox}
\define@boolkey{MCM}[MCM@opt@]{CTeX}[false]{}
\define@boolkey{MCM}[MCM@opt@]{titlepage}[true]{}
\define@boolkey{MCM}[MCM@opt@]{abstract}[true]{}
\define@boolkey{MCM}[MCM@opt@]{sheet}[true]{}
\define@boolkey{MCM}[MCM@opt@]{titleinsheet}[false]{}
\define@boolkey{MCM}[MCM@opt@]{keywordsinsheet}[false]{}
\define@cmdkeys{MCM}[MCM@opt@]{tcn,problem}
\define@key{MCM}{tcn}[0000]{\gdef\MCM@opt@tcn{#1}}
\define@key{MCM}{problem}[A]{\gdef\MCM@opt@problem{#1}}
\setkeys{MCM}{tcn=0000,problem=B}
\define@key{mcmthesis.cls}{tcn}[0000]{\gdef\MCM@opt@tcn{#1}}
\define@key{mcmthesis.cls}{problem}[A]{\gdef\MCM@opt@problem{#1}}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{titlepage}{}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{abstract}{}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{sheet}{}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{titleinsheet}{}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{keywordsinsheet}{}
\MCM@opt@sheettrue
\MCM@opt@titlepagetrue
\MCM@opt@titleinsheetfalse
\MCM@opt@keywordsinsheetfalse
\MCM@opt@abstracttrue
\newcommand{\mcmsetup}[1]{\setkeys{MCM}{#1}}
\ProcessOptionsX\relax
\LoadClass[a4paper, 12pt]{article}
\newcommand{\team}{Team \#\ \MCM@opt@tcn}
\RequirePackage{fancyhdr, fancybox}
\RequirePackage{ifthen}
\RequirePackage{lastpage}
\RequirePackage{listings}
\RequirePackage[toc, page, title, titletoc, header]{appendix}
\RequirePackage{paralist}
\RequirePackage{amsthm, amsfonts}
\RequirePackage{amsmath, bm}
\RequirePackage{amssymb, mathrsfs}
\RequirePackage{latexsym}
\RequirePackage{longtable, multirow, hhline, tabularx, array}
\RequirePackage{flafter}
\RequirePackage{pifont, calc}
\RequirePackage{colortbl, booktabs}
\RequirePackage{geometry}
\RequirePackage[T1]{fontenc}
\RequirePackage[scaled]{berasans}
\RequirePackage{hyperref}
\RequirePackage{ifpdf, ifxetex}
\ifMCM@opt@CTeX
\else
\RequirePackage{environ}
\fi
\ifpdf
\RequirePackage{graphicx}
\RequirePackage{epstopdf}
\else
\ifxetex
\RequirePackage{graphicx}
\else
\RequirePackage[dvipdfmx]{graphicx}
\RequirePackage{bmpsize}
\fi
\fi
\RequirePackage[svgnames]{xcolor}
\ifpdf
\hypersetup{hidelinks}
\else
\ifxetex
\hypersetup{hidelinks}
\else
\hypersetup{dvipdfm, hidelinks}
\fi
\fi
\geometry{a4paper, margin = 1in}
\pagestyle{fancy}
\fancyhf{}
\lhead{\small\sffamily \team}
% \rhead{\small\sffamily Page \thepage\ of \pageref{LastPage}}
\rhead{\small\sffamily Page \thepage}
\setlength\parskip{.5\baselineskip}
\renewcommand\tableofcontents{%
\centerline{\normalfont\Large\bfseries\sffamily\contentsname
\@mkboth{%
\MakeUppercase\contentsname}{\MakeUppercase\contentsname}}%
\vskip 5ex%
\@starttoc{toc}%
}
\setcounter{totalnumber}{4}
\setcounter{topnumber}{2}
\setcounter{bottomnumber}{2}
\renewcommand{\textfraction}{0.15}
\renewcommand{\topfraction}{0.85}
\renewcommand{\bottomfraction}{0.65}
\renewcommand{\floatpagefraction}{0.60}
\renewcommand{\figurename}{Figure}
\renewcommand{\tablename}{Table}
\graphicspath{{./}{./img/}{./fig/}{./image/}{./figure/}{./picture/}
{./imgs/}{./figs/}{./images/}{./figures/}{./pictures/}}
\def\maketitle{%
\let\saved@thepage\thepage
\let\thepage\relax
\ifMCM@opt@sheet
\makesheet
\fi
\newpage
\ifMCM@opt@titlepage
\MCM@maketitle
\fi
\newpage
\let\thepage\saved@thepage
\setcounter{page}{1}
\pagestyle{fancy}
}
\def\abstractname{Summary}
\ifMCM@opt@CTeX
\newbox\@abstract%
\setbox\@abstract\hbox{}%
\long\def\abstract{\bgroup\global\setbox\@abstract\vbox\bgroup\hsize\textwidth\leftskip1cm\rightskip1cm}%
\def\endabstract{\egroup\egroup}
\def\make@abstract{%
\begin{center}
\textbf{\abstractname}
\end{center}
\usebox\@abstract\par
}
\else
\RenewEnviron{abstract}{\xdef\@abstract{\expandonce\BODY}}
\def\make@abstract{%
\begin{center}
\textbf{\abstractname}
\end{center}
\@abstract\par
}
\fi
\newenvironment{letter}[1]{%
\par%
\bgroup\parindent0pt%
\begin{minipage}{5cm}
\flushleft #1%
\end{minipage}}
{\egroup\smallskip}
\def\keywordsname{Keywords}
\ifMCM@opt@CTeX
\newbox\@keywords
\setbox\@keywords\hbox{}
\def\keywords{\global\setbox\@keywords\vbox\bgroup\noindent\leftskip0cm}
\def\endkeywords{\egroup}%
\def\make@keywords{%
\par\hskip.4cm\textbf{\keywordsname}: \usebox\@keywords\hfill\par
}
\else
\NewEnviron{keywords}{\xdef\@keywords{\expandonce\BODY}}
\def\make@keywords{%
\par\noindent\textbf{\keywordsname}:
\@keywords\par
}
\fi
\newcommand{\headset}{{\the\year}\\MCM/ICM\\Summary Sheet}
\newcommand{\problem}[1]{\mcmsetup{problem = #1}}
\def\makesheet{%
\pagestyle{empty}%
\null%
\vspace*{-5pc}%
\begin{center}
\begingroup
\setlength{\parindent}{0pt}
\begin{minipage}[t]{0.33\linewidth}
\bfseries\centering%
Problem Chosen\\[0.7pc]
{\Huge\textbf{\MCM@opt@problem}}\\[2.8pc]
\end{minipage}%
\begin{minipage}[t]{0.33\linewidth}
\centering%
\textbf{\headset}%
\end{minipage}%
\begin{minipage}[t]{0.33\linewidth}
\centering\bfseries%
Team Control Number\\[0.7pc]
{\Huge\textbf{\MCM@opt@tcn}}\\[2.8pc]
% {\Huge\textbf{\MCM@opt@tcn}}\\[2.8pc]
\end{minipage}\par
\rule{\linewidth}{0.8pt}\par
\endgroup
\vskip 10pt%
\ifMCM@opt@titleinsheet
\normalfont \LARGE \@title \par
\fi
\end{center}
\ifMCM@opt@keywordsinsheet
\make@abstract
\make@keywords
\else
\make@abstract
\fi}
\newcommand{\MCM@maketitle}{%
\begin{center}%
\let \footnote \thanks
{\LARGE \@title \par}%
\vskip 1.5em%
{\large
\lineskip .5em%
\begin{tabular}[t]{c}%
\@author
\end{tabular}\par}%
\vskip 1em%
{\large \@date}%
\end{center}%
\par
\vskip 1.5em%
\ifMCM@opt@abstract%
\make@abstract
\make@keywords
\fi%
}
\def\MCM@memoto{\relax}
\newcommand{\memoto}[1]{\gdef\MCM@memoto{#1}}
\def\MCM@memofrom{\relax}
\newcommand{\memofrom}[1]{\gdef\MCM@memofrom{#1}}
\def\MCM@memosubject{\relax}
\newcommand{\memosubject}[1]{\gdef\MCM@memosubject{#1}}
\def\MCM@memodate{\relax}
\newcommand{\memodate}[1]{\gdef\MCM@memodate{#1}}
\def\MCM@memologo{\relax}
\newcommand{\memologo}[1]{\gdef\MCM@memologo{\protect #1}}
\def\@letterheadaddress{\relax}
\newcommand{\lhaddress}[1]{\gdef\@letterheadaddress{#1}}
\newenvironment{memo}[1][Memorandum]{%
\pagestyle{plain}%
\ifthenelse{\equal{\MCM@memologo}{\relax}}{%
% without logo specified.
}{%
% with logo specified
\begin{minipage}[t]{\columnwidth}%
\begin{flushright}
\vspace{-0.6in}
\MCM@memologo
\vspace{0.5in}
\par\end{flushright}%
\end{minipage}%
}
\begin{center}
\LARGE\bfseries\scshape
#1
\end{center}
\begin{description}
\ifthenelse{\equal{\MCM@memoto}{\relax}}{}{\item [{To:}] \MCM@memoto}
\ifthenelse{\equal{\MCM@memofrom}{\relax}}{}{\item [{From:}] \MCM@memofrom}
\ifthenelse{\equal{\MCM@memosubject}{\relax}}{}{\item [{Subject:}] \MCM@memosubject}
\ifthenelse{\equal{\MCM@memodate}{\relax}}{}{\item [{Date:}] \MCM@memodate}
\end{description}
\par\noindent
\rule[0.5ex]{\linewidth}{0.1pt}\par
\bigskip{}
}{%
\clearpage
\pagestyle{fancy}%
}
\newtheorem{Theorem}{Theorem}[section]
\newtheorem{Lemma}[Theorem]{Lemma}
\newtheorem{Corollary}[Theorem]{Corollary}
\newtheorem{Proposition}[Theorem]{Proposition}
\newtheorem{Definition}[Theorem]{Definition}
\newtheorem{Example}[Theorem]{Example}
\renewcommand\section{\@startsection{section}{1}{\z@}%
{-0pt\@plus -.2ex \@minus -.2ex}%
{1pt \@plus .2ex}%
{\rmfamily\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
{-0pt\@plus -.2ex \@minus -.2ex}%
{1pt \@plus .2ex}%
{\rmfamily\large\bfseries}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
{-.5ex\@plus -1ex \@minus -.2ex}%
{.25ex \@plus .2ex}%
{\rmfamily\normalsize\bfseries}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
{1ex \@plus1ex \@minus.2ex}%
{-1em}%
{\rmfamily\normalsize}}
\providecommand{\dif}{\mathop{}\!\mathrm{d}}
\providecommand{\me}{\mathrm{e}}
\providecommand{\mi}{\mathrm{i}}
\definecolor{grey}{rgb}{0.8,0.8,0.8}
\definecolor{darkgreen}{rgb}{0,0.3,0}
\definecolor{darkblue}{rgb}{0,0,0.3}
\def\lstbasicfont{\fontfamily{pcr}\selectfont\footnotesize}
\lstset{%
% numbers=left,
% numberstyle=\small,%
showstringspaces=false,
showspaces=false,%
tabsize=4,%
frame=lines,%
basicstyle={\footnotesize\lstbasicfont},%
keywordstyle=\color{darkblue}\bfseries,%
identifierstyle=,%
commentstyle=\color{darkgreen},%\itshape,%
stringstyle=\color{black}%
}
\lstloadlanguages{C,C++,Java,Matlab,Mathematica}
\endinput
%%
%% This work consists of these files mcmthesis.dtx,
%% figures/ and
%% code/,
%% and the derived files mcmthesis.cls,
%% mcmthesis-demo.tex,
%% README,
%% LICENSE,
%% mcmthesis.pdf and
%% mcmthesis-demo.pdf.
%%
%% End of file `mcmthesis.cls'.

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "61,0",
"scrollLine": "53"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "465,143",
"scrollLine": "461"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "167,28",
"scrollLine": "163"
}

View File

@@ -0,0 +1,7 @@
{
"tempName": "Untitled1",
"source_window_id": "",
"Source": "Source",
"cursorPosition": "2199,0",
"scrollLine": "2047"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "72,0",
"scrollLine": "59"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "80,4",
"scrollLine": "0"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "488,15",
"scrollLine": "482"
}

View File

@@ -0,0 +1,7 @@
C%3A%2FUsers%2FZhang%2FDesktop%2F21%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM-Latex%E6%A8%A1%E6%9D%BF(%E4%BF%AE%E6%94%B9%E7%89%88)%2Fmcmthesis-demo.tex="2F02FDC6"
C%3A%2FUsers%2FZhang%2FDesktop%2F21%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF%2Ffancybox.sty="6ABA5725"
C%3A%2FUsers%2FZhang%2FDesktop%2F21%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF%2Fmcmthesis-demo.tex="F5671CC1"
C%3A%2FUsers%2FZhang%2FDesktop%2F21%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF%2Fmcmthesis.cls="40CBEBD4"
C%3A%2FUsers%2FZhang%2FDesktop%2F21%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2021%E6%9C%80%E6%96%B0%E4%BF%AE%E6%94%B9)%2FMCM-ICM_Summary.tex="719CF748"
C%3A%2FUsers%2FZhang%2FDesktop%2F21%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2021%E6%9C%80%E6%96%B0%E4%BF%AE%E6%94%B9)%2Fmcmthesis-demo.tex="D6BFA82F"
C%3A%2FUsers%2FZhang%2FDesktop%2F21%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2021%E6%9C%80%E6%96%B0%E4%BF%AE%E6%94%B9)%2Fmcmthesis.cls="07CC3646"

View File

@@ -0,0 +1,9 @@
{
"sortOrder": [
{
"columnIndex": 2,
"ascending": true
}
],
"path": "C:/Users/zhang/Desktop/25美赛/MCM_Latex2025"
}

View File

@@ -0,0 +1,3 @@
{
"activeTab": 2
}

View File

@@ -0,0 +1,14 @@
{
"left": {
"splitterpos": 241,
"topwindowstate": "NORMAL",
"panelheight": 844,
"windowheight": 882
},
"right": {
"splitterpos": 528,
"topwindowstate": "NORMAL",
"panelheight": 844,
"windowheight": 882
}
}

View File

@@ -0,0 +1,5 @@
{
"TabSet1": 0,
"TabSet2": 0,
"TabZoom": {}
}

View File

@@ -0,0 +1,5 @@

View File

@@ -0,0 +1 @@
{"active_set":"","sets":[]}

View File

@@ -0,0 +1,26 @@
{
"id": "4DB59D00",
"path": "C:/Users/zhang/Desktop/25美赛/MCM_Latex2025/run.R",
"project_path": "run.R",
"type": "r_source",
"hash": "1612192916",
"contents": "",
"dirty": false,
"created": 1736337015200.0,
"source_on_save": false,
"relative_order": 2,
"properties": {
"source_window_id": "",
"Source": "Source",
"cursorPosition": "2,52",
"scrollLine": "0"
},
"folds": "",
"lastKnownWriteTime": 1736710945,
"encoding": "UTF-8",
"collab_server": "",
"source_window": "",
"last_content_update": 1736710945127,
"read_only": false,
"read_only_alternatives": []
}

View File

@@ -0,0 +1,3 @@
library(tinytex)
pdflatex("mcmthesis-demo.tex", bib_engine = "biber")

View File

@@ -0,0 +1,26 @@
{
"id": "82F97103",
"path": "C:/Users/zhang/Desktop/25美赛/MCM_Latex2025/reference.bib",
"project_path": "reference.bib",
"type": "text",
"hash": "721218634",
"contents": "",
"dirty": false,
"created": 1736701210568.0,
"source_on_save": false,
"relative_order": 3,
"properties": {
"source_window_id": "",
"Source": "Source",
"cursorPosition": "75,1",
"scrollLine": "50"
},
"folds": "",
"lastKnownWriteTime": 1736716598,
"encoding": "UTF-8",
"collab_server": "",
"source_window": "",
"last_content_update": 1736716598157,
"read_only": false,
"read_only_alternatives": []
}

View File

@@ -0,0 +1,76 @@
@article{kim2006,
author = {Gi-Beum Kim},
title = {Change of the Warm Water Temperature for the Development of Smart Healthecare Bathing System},
journal = {Hwahak konghak},
year = {2006},
volume = {44},
number = {3},
pages = {270-276}
}
@article{holman2002,
author = {Holman, J. P.},
title = {Heat Transfer (9th ed.)},
journal = {McGraw-Hill},
year = {2002},
address = {New York}
}
@book{anderson2006,
author = {Anderson, D. A. and Tannehill, J. C. and Pletcher, R. H.},
title = {Computational Fluid Dynamics and Heat Transfer},
publisher = {Taylor \& Francis},
year = {2006},
edition = {2nd}
}
@article{incropera2007,
author = {Incropera, F. P. and DeWitt, D. P. and Bergman, T. L. and Lavine, A. S.},
title = {Fundamentals of Heat and Mass Transfer},
journal = {John Wiley \& Sons},
year = {2007},
edition = {6th}
}
@article{evaporation2018,
author = {Smith, J.},
title = {Evaporation Effects in Bath Water Temperature Regulation},
journal = {Journal of Thermal Science and Engineering},
year = {2018},
volume = {30},
number = {4},
pages = {456-467}
}
@phdthesis{thesis2015,
author = {Brown, D.},
title = {Thermal Modeling of Bath Water Systems},
school = {University of Technology},
year = {2015},
address = {City}
}
@misc{website2024,
author = {Thompson, H.},
title = {Advanced CFD Techniques for Bath Water Simulation},
howpublished = {\url{https://www.example.com/advanced-cfd}},
note = {Accessed: 2024-10-01}
}
@patent{patent2023,
author = {Wilson, G.},
title = {Smart Bath Water Temperature Control System},
number = {1234567},
year = {2023},
type = {US Patent},
assignee = {SmartTech Inc.}
}
@book{Liu02,
author = {Liu, Weiguo and Chen, Zhaoping and Zhang, Yin},
title = {Matlab program design and application},
publisher = {Higher education press},
address = {Beijing},
year = {2002},
note = {In Chinese}
}

View File

@@ -0,0 +1,26 @@
{
"id": "A0D261B2",
"path": "C:/Users/zhang/Desktop/25美赛/MCM_Latex2025/mcmthesis-demo.tex",
"project_path": "mcmthesis-demo.tex",
"type": "tex",
"hash": "1436774386",
"contents": "",
"dirty": false,
"created": 1704816448925.0,
"source_on_save": false,
"relative_order": 1,
"properties": {
"source_window_id": "",
"Source": "Source",
"cursorPosition": "715,74",
"scrollLine": "701"
},
"folds": "",
"lastKnownWriteTime": 1736716652,
"encoding": "UTF-8",
"collab_server": "",
"source_window": "",
"last_content_update": 1736716652098,
"read_only": false,
"read_only_alternatives": []
}

View File

@@ -0,0 +1,774 @@
%%
%% This is file `mcmthesis-demo.tex',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `demo')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
%%
%% This is file `mcmthesis-demo.tex',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `demo')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
\documentclass{mcmthesis}
\mcmsetup{CTeX = false, % 使用 CTeX 套装时,设置为 true
tcn = {0000000}, problem = \textcolor{red}{A},
sheet = true, titleinsheet = true, keywordsinsheet = true,
titlepage = false, abstract = false}
\usepackage{newtxtext} % \usepackage{palatino}
\usepackage[style=apa,backend=biber]{biblatex}
\addbibresource{reference.bib}
\usepackage{tocloft}
\setlength{\cftbeforesecskip}{6pt}
\renewcommand{\contentsname}{\hspace*{\fill}\Large\bfseries Contents \hspace*{\fill}}
\title{Enjoy a Cozy and Green Bath}
% \author{\small \href{http://www.latexstudio.net/}
% {\includegraphics[width=7cm]{mcmthesis-logo}}}
\date{\today}
\begin{document}
\begin{abstract}
A traditional bathtub cannot be reheated by itself, so users have to add hot
water from time to time. Our goal is to establish a model of the temperature
of bath water in space and time. Then we are expected to propose an optimal
strategy for users to keep the temperature even and close to initial temperature
and decrease water consumption.
To simplify modeling process, we firstly assume there is no person in the bathtub.
We regard the whole bathtub as a thermodynamic system and introduce heat transfer
formulas.
We establish two sub-models: adding water constantly and discontinuously. As for
the former sub-model, we define the mean temperature of bath water. Introducing
Newton cooling formula, we determine the heat transfer capacity. After deriving
the value of parameters, we deduce formulas to derive results and simulate the
change of temperature field via CFD. As for the second sub-model, we define an
iteration consisting of two process: heating and standby. According to energy
conservation law, we obtain the relationship of time and total heat dissipating
capacity. Then we determine the mass flow and the time of adding hot water. We
also use CFD to simulate the temperature field in second sub-model.
In consideration of evaporation, we correct the results of sub-models referring
to some scientists' studies. We define two evaluation criteria and compare the
two sub-models. Adding water constantly is found to keep the temperature of bath
water even and avoid wasting too much water, so it is recommended by us.
Then we determine the influence of some factors: radiation heat transfer, the
shape and volume of the tub, the shape/volume/temperature/motions of the person,
the bubbles made from bubble bath additives. We focus on the influence of those
factors to heat transfer and then conduct sensitivity analysis. The results
indicate smaller bathtub with less surface area, lighter personal mass, less
motions and more bubbles will decrease heat transfer and save water.
Based on our model analysis and conclusions, we propose the optimal strategy for
the user in a bathtub and explain the reason of uneven temperature throughout
the bathtub. In addition, we make improvement for applying our model in real life.
\begin{keywords}
Heat transfer, Thermodynamic system, CFD, Energy conservation
\end{keywords}
\end{abstract}
\maketitle
%% Generate the Table of Contents, if it's needed.
% \renewcommand{\contentsname}{\centering Contents}
\tableofcontents % 若不想要目录, 注释掉该句
\thispagestyle{empty}
\newpage
\section{Introduction}
\subsection{Background}
Bathing in a tub is a perfect choice for those who have been worn out after a
long day's working. A traditional bathtub is a simply water containment vessel
without a secondary heating system or circulating jets. Thus the temperature of
water in bathtub declines noticeably as time goes by, which will influent the
experience of bathing. As a result, the bathing person needs to add a constant
trickle of hot water from a faucet to reheat the bathing water. This way of
bathing will result in waste of water because when the capacity of the bathtub
is reached, excess water overflows the tub.
An optimal bathing strategy is required for the person in a bathtub to get
comfortable bathing experience while reducing the waste of water.
\subsection{Literature Review}
A traditional bathtub cannot be reheated by itself, so users have to add hot
water from time to time. Our goal is to establish a model of the temperature
of bath water in space and time. Then we are expected to propose an optimal
strategy for users to keep the temperature even and close to the initial
temperature and decrease water consumption. According to \textcite{kim2006},
He derived a relational equation based on the basic theory of heat transfer
to evaluate the performance of bath tubes. The major heat loss was found to be
due to evaporation. Moreover, he found out that the speed of heat loss depends
more on the humidity of the bathroom than the temperature of water contained
in the bathtub. So, it is best to maintain the temperature of bathtub water to
be between 41 to 45$^{\circ}$C and the humidity of bathroom to be 95\%.
Traditional bath systems have significant limitations in temperature control.
To address this, we introduce heat transfer formulas as
discussed (\cite[123]{holman2002}).
\subsection{Restatement of the Problem}
We are required to establish a model to determine the change of water temperature
in space and time. Then we are expected to propose the best strategy for the
person in the bathtub to keep the water temperature close to initial temperature
and even throughout the tub. Reduction of waste of water is also needed. In
addition, we have to consider the impact of different conditions on our model,
such as different shapes and volumes of the bathtub, etc.
In order to solve those problems, we will proceed as follows:
\begin{itemize}
\item {\bf Stating assumptions}. By stating our assumptions, we will narrow the
focus of our approach towards the problems and provide some insight into bathtub
water temperature issues.
\item {\bf Making notations}. We will give some notations which are important for
us to clarify our models.
\item {\bf Presenting our model}. In order to investigate the problem deeper, we
divide our model into two sub-models. One is a steady convection heat transfer
sub-model in which hot water is added constantly. The other one is an unsteady
convection heat transfer sub-model where hot water is added discontinuously.
\item {Defining evaluation criteria and comparing sub-models}. We define two main
criteria to evaluate our model: the mean temperature of bath water and the amount
of inflow water.
\item {\bf Analysis of influencing factors}. In term of the impact of different
factors on our model, we take those into consideration: the shape and volume of
the tub, the shape/volume/temperature of the person in the bathtub, the motions
made by the person in the bathtub and adding a bubble bath additive initially.
\item {\bf Model testing and sensitivity analysis}. With the criteria defined
before, we evaluate the reliability of our model and do the sensitivity analysis.
\item {\bf Further discussion}. We discuss about different ways to arrange inflow
faucets. Then we improve our model to apply them in reality.
\item {\bf Evaluating the model}. We discuss about the strengths and weaknesses
of our model:
\begin{itemize}
\item[1)] ...
\item[2)] ...
\item[3)] ...
\item[4)] ...
\end{itemize}
\end{itemize}
\section{Assumptions and Justification}
To simplify the problem and make it convenient for us to simulate real-life
conditions, we make the following basic assumptions, each of which is properly
justified.
\begin{itemize}
\item {\bf The bath water is incompressible Non-Newtonian fluid}. The
incompressible Non-Newtonian fluid is the basis of NavierStokes equations
which are introduced to simulate the flow of bath water.
\item {\bf All the physical properties of bath water, bathtub and air are
assumed to be stable}. The change of those properties like specific heat,
thermal conductivity and density is rather small according to some
studies. It is complicated and unnecessary to consider these little
change so we ignore them.
\item {\bf There is no internal heat source in the system consisting of bathtub,
hot water and air}. Before the person lies in the bathtub, no internal heat source
exist except the system components. The circumstance where the person is in the
bathtub will be investigated in our later discussion.
\item {\bf We ignore radiative thermal exchange}. According to Stefan-Boltzmanns
law, the radiative thermal exchange can be ignored when the temperature is low.
Refer to industrial standard, the temperature in bathroom is lower than
100 $^{\circ}$C, so it is reasonable for us to make this assumption.
\item {\bf The temperature of the adding hot water from the faucet is stable}.
This hypothesis can be easily achieved in reality and will simplify our process
of solving the problem.
\end{itemize}
\section{Notations}
\begin{center}
\begin{tabular}{clc}
{\bf Symbols} & {\bf Description} & \quad {\bf Unit} \\[0.25cm]
$h$ & Convection heat transfer coefficient & \quad W/(m$^2 \cdot$ K)
\\[0.2cm]
$k$ & Thermal conductivity & \quad W/(m $\cdot$ K) \\[0.2cm]
$c_p$ & Specific heat & \quad J/(kg $\cdot$ K) \\[0.2cm]
$\rho$ & Density & \quad kg/m$^2$ \\[0.2cm]
$\delta$ & Thickness & \quad m \\[0.2cm]
$t$ & Temperature & \quad $^\circ$C, K \\[0.2cm]
$\tau$ & Time & \quad s, min, h \\[0.2cm]
$q_m$ & Mass flow & \quad kg/s \\[0.2cm]
$\Phi$ & Heat transfer power & \quad W \\[0.2cm]
$T$ & A period of time & \quad s, min, h \\[0.2cm]
$V$ & Volume & \quad m$^3$, L \\[0.2cm]
$M,\,m$ & Mass & \quad kg \\[0.2cm]
$A$ & Aera & \quad m$^2$ \\[0.2cm]
$a,\,b,\,c$ & The size of a bathtub & \quad m$^3$
\end{tabular}
\end{center}
\noindent where we define the main parameters while specific value of those
parameters will be given later.
\section{Model Overview}
To simplify the modeling process, we firstly assume there is no person in the
bathtub. We regard the whole bathtub as a thermodynamic system and introduce
heat transfer formulas. We establish two sub-models: adding water constantly
and discontinuously. For the former sub-model, we define the mean temperature
of bath water and introduce Newton's cooling formula to determine the heat
transfer capacity. After deriving the value of parameters, we deduce formulas
to derive results and simulate the change of temperature field via CFD, as
described by \textcite{anderson2006}.
In our basic model, we aim at three goals: keeping the temperature as even as
possible, making it close to the initial temperature and decreasing the water
consumption.
We start with the simple sub-model where hot water is added constantly.
At first we introduce convection heat transfer control equations in rectangular
coordinate system. Then we define the mean temperature of bath water.
Afterwards, we introduce Newton cooling formula to determine heat transfer
capacity. After deriving the value of parameters, we get calculating results
via formula deduction and simulating results via CFD.
Secondly, we present the complicated sub-model in which hot water is
added discontinuously. We define an iteration consisting of two process:
heating and standby. As for heating process, we derive control equations and
boundary conditions. As for standby process, considering energy conservation law,
we deduce the relationship of total heat dissipating capacity and time.
Then we determine the time and amount of added hot water. After deriving the
value of parameters, we get calculating results via formula deduction and
simulating results via CFD.
At last, we define two criteria to evaluate those two ways of adding hot water.
Then we propose optimal strategy for the user in a bathtub.
The whole modeling process can be shown as follows.
\begin{figure}[h]
\centering
\includegraphics[width=12cm]{fig1.jpg}
\caption{Modeling process} \label{fig1}
\end{figure}
\section{Sub-model I : Adding Water Continuously}
As for the second sub-model, we define an iteration consisting of two processes:
heating and standby. According to the energy conservation law, we obtain the
relationship of time and total heat dissipating capacity. Then we determine
the mass flow and the time of adding hot water. We also use CFD to simulate
the temperature field in the second sub-model, following the techniques
outlined by \textcite{website2024}.
We first establish the sub-model based on the condition that a person add water
continuously to reheat the bathing water. Then we use Computational Fluid
Dynamics (CFD) to simulate the change of water temperature in the bathtub. At
last, we evaluate the model with the criteria which have been defined before.
\subsection{Model Establishment}
Since we try to keep the temperature of the hot water in bathtub to be even,
we have to derive the amount of inflow water and the energy dissipated by the
hot water into the air.
We derive the basic convection heat transfer control equations based on the
former scientists achievement. Then, we define the mean temperature of bath
water. Afterwards, we determine two types of heat transfer: the boundary heat
transfer and the evaporation heat transfer. Combining thermodynamic formulas,
we derive calculating results. Via Fluent software, we get simulation results.
\subsubsection{Control Equations and Boundary Conditions}
According to thermodynamics knowledge, we recall on basic convection
heat transfer control equations in rectangular coordinate system. Those
equations show the relationship of the temperature of the bathtub water in space.
We assume the hot water in the bathtub as a cube. Then we put it into a
rectangular coordinate system. The length, width, and height of it is $a,\, b$
and $c$.
\begin{figure}[h]
\centering
\includegraphics[width=8cm]{fig2.jpg}
\caption{Modeling process} \label{fig2}
\end{figure}
In the basis of this, we introduce the following equations:
\begin{itemize}
\item {\bf Continuity equation:}
\end{itemize}
\begin{equation} \label{eq1}
\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} +
\frac{\partial w}{\partial z} = 0
\end{equation}
\noindent where the first component is the change of fluid mass along the $X$-ray.
The second component is the change of fluid mass along the $Y$-ray. And the third
component is the change of fluid mass along the $Z$-ray. The sum of the change in
mass along those three directions is zero.
\begin{itemize}
\item {\bf Moment differential equation (N-S equations):}
\end{itemize}
\begin{equation} \label{eq2}
\left\{
\begin{array}{l} \!\!
\rho \Big(u \dfrac{\partial u}{\partial x} + v \dfrac{\partial u}{\partial y} +
w\dfrac{\partial u}{\partial z} \Big) = -\dfrac{\partial p}{\partial x} +
\eta \Big(\dfrac{\partial^2 u}{\partial x^2} + \dfrac{\partial^2 u}{\partial y^2} +
\dfrac{\partial^2 u}{\partial z^2} \Big) \\[0.3cm]
\rho \Big(u \dfrac{\partial v}{\partial x} + v \dfrac{\partial v}{\partial y} +
w\dfrac{\partial v}{\partial z} \Big) = -\dfrac{\partial p}{\partial y} +
\eta \Big(\dfrac{\partial^2 v}{\partial x^2} + \dfrac{\partial^2 v}{\partial y^2} +
\dfrac{\partial^2 v}{\partial z^2} \Big) \\[0.3cm]
\rho \Big(u \dfrac{\partial w}{\partial x} + v \dfrac{\partial w}{\partial y} +
w\dfrac{\partial w}{\partial z} \Big) = -g-\dfrac{\partial p}{\partial z} +
\eta \Big(\dfrac{\partial^2 w}{\partial x^2} + \dfrac{\partial^2 w}{\partial y^2} +
\dfrac{\partial^2 w}{\partial z^2} \Big)
\end{array}
\right.
\end{equation}
\begin{itemize}
\item {\bf Energy differential equation:}
\end{itemize}
\begin{equation} \label{eq3}
\rho c_p \Big( u\frac{\partial t}{\partial x} + v\frac{\partial t}{\partial y} +
w\frac{\partial t}{\partial z} \Big) = \lambda \Big(\frac{\partial^2 t}{\partial x^2} +
\frac{\partial^2 t}{\partial y^2} + \frac{\partial^2 t}{\partial z^2} \Big)
\end{equation}
\noindent where the left three components are convection terms while the right
three components are conduction terms.
By Equation \eqref{eq3}, we have ......
......
On the right surface in Fig. \ref{fig2}, the water also transfers heat firstly
with bathtub inner surfaces and then the heat comes into air. The boundary
condition here is ......
\subsubsection{Definition of the Mean Temperature}
......
\subsubsection{Determination of Heat Transfer Capacity}
......
\section{Sub-model II: Adding Water Discontinuously}
In order to establish the unsteady sub-model, we recall on the working principle
of air conditioners. The heating performance of air conditions consist of two
processes: heating and standby. After the user set a temperature, the air
conditioner will begin to heat until the expected temperature is reached. Then
it will go standby. When the temperature get below the expected temperature,
the air conditioner begin to work again. As it works in this circle, the
temperature remains the expected one.
Inspired by this, we divide the bathtub working into two processes: adding
hot water until the expected temperature is reached, then keeping this
condition for a while unless the temperature is lower than a specific value.
Iterating this circle ceaselessly will ensure the temperature kept relatively
stable.
\subsection{Heating Model}
\subsubsection{Control Equations and Boundary Conditions}
\subsubsection{Determination of Inflow Time and Amount}
\subsection{Standby Model}
\subsection{Results}
\quad~ We first give the value of parameters based on others studies. Then we
get the calculation results and simulating results via those data.
\subsubsection{Determination of Parameters}
After establishing the model, we have to determine the value of some
important parameters.
As scholar Beum Kim points out, the optimal temperature for bath is
between 41 and 45$^\circ$C. Meanwhile, according to Shimodozono's study,
41$^\circ$C warm water bath is the perfect choice for individual health.
So it is reasonable for us to focus on $41^\circ$C $\sim 45^\circ$C. Because
adding hot water continuously is a steady process, so the mean temperature
of bath water is supposed to be constant. We value the temperature of inflow
and outflow water with the maximum and minimum temperature respectively.
The values of all parameters needed are shown as follows:
.....
\subsubsection{Calculating Results}
Putting the above value of parameters into the equations we derived before, we
can get the some data as follows:
%%普通表格
\begin{table}[h] %h表示固定在当前位置
\centering %设置居中
\caption{The calculating results} %表标题
\vspace{0.15cm}
\label{tab2} %设置表的引用标签
\begin{tabular}{|c|c|c|} %3个c表示3列, |可选, 表示绘制各列间的竖线
\hline %画横线
Variables & Values & Unit \\ \hline %各列间用&隔开
$A_1$ & 1.05 & $m^2$ \\ \hline
$A_2$ & 2.24 & $m^2$ \\ \hline
$\Phi_1$ & 189.00 & $W$ \\ \hline
$\Phi_2$ & 43.47 & $W$ \\ \hline
$\Phi$ & 232.47 & $W$ \\ \hline
$q_m$ & 0.014 & $g/s$ \\ \hline
\end{tabular}
\end{table}
From Table \ref{tab2}, ......
......
\section{Correction and Contrast of Sub-Models}
After establishing two basic sub-models, we have to correct them in consideration
of evaporation heat transfer. Then we define two evaluation criteria to compare
the two sub-models in order to determine the optimal bath strategy.
\subsection{Correction with Evaporation Heat Transfer}
Someone may confuse about the above results: why the mass flow in the first
sub-model is so small? Why the standby time is so long? Actually, the above two
sub-models are based on ideal conditions without consideration of the change of
boundary conditions, the motions made by the person in bathtub and the
evaporation of bath water, etc. The influence of personal motions will be
discussed later. Here we introducing the evaporation of bath water to correct
sub-models.
\subsection{Contrast of Two Sub-Models}
Firstly we define two evaluation criteria. Then we contrast the two submodels
via these two criteria. Thus we can derive the best strategy for the person in
the bathtub to adopt.
\section{Model Analysis and Sensitivity Analysis}
In consideration of evaporation, we correct the results of sub-models referring
to studies. We define two evaluation criteria
and compare the two sub-models. Adding water constantly is found to keep the
temperature of bath water even and avoid wasting too much water, so it is
recommended by us. We also conduct sensitivity analysis to determine the
influence of factors such as radiation heat transfer, the shape and volume of
the tub, the shape/volume/temperature/motions of the person, and the bubbles
made from bubble bath additives, as discussed
in (\cite{evaporation2018}; \cite{thesis2015}).
\subsection{The Influence of Different Bathtubs}
Definitely, the difference in shape and volume of the tub affects the
convection heat transfer. Examining the relationship between them can help
people choose optimal bathtubs.
\subsubsection{Different Volumes of Bathtubs}
In reality, a cup of water will be cooled down rapidly. However, it takes quite
long time for a bucket of water to become cool. That is because their volume is
different and the specific heat of water is very large. So that the decrease of
temperature is not obvious if the volume of water is huge. That also explains
why it takes 45 min for 320 L water to be cooled by 1$^\circ$C.
In order to examine the influence of volume, we analyze our sub-models
by conducting sensitivity Analysis to them.
We assume the initial volume to be 280 L and change it by $\pm 5$\%, $\pm 8$\%,
$\pm 12$\% and $\pm 15$\%. With the aid of sub-models we established before, the
variation of some parameters turns out to be as follows
%%三线表
\begin{table}[h] %h表示固定在当前位置
\centering %设置居中
\caption{Variation of some parameters} %表标题
\label{tab7} %设置表的引用标签
\begin{tabular}{ccccccc} %7个c表示7列, c表示每列居中对齐, 还有l和r可选
\toprule %画顶端横线
$V$ & $A_1$ & $A_2$ & $T_2$ & $q_{m1}$ & $q_{m2}$ & $\Phi_q$ \\
\midrule %画中间横线
-15.00\% & -5.06\% & -9.31\% & -12.67\% & -2.67\% & -14.14\% & -5.80\% \\
-12.00\% & -4.04\% & -7.43\% & -10.09\% & -2.13\% & -11.31\% & -4.63\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
\bottomrule %画底部横线
\end{tabular}
\end{table}
\section{Strength and Weakness}
\subsection{Strength}
\begin{itemize}
\item We analyze the problem based on thermodynamic formulas and laws, so that
the model we established is of great validity.
\item Our model is fairly robust due to our careful corrections in consideration
of real-life situations and detailed sensitivity analysis.
\item Via Fluent software, we simulate the time field of different areas
throughout the bathtub. The outcome is vivid for us to understand the changing
process.
\item We come up with various criteria to compare different situations, like
water consumption and the time of adding hot water. Hence an overall comparison
can be made according to these criteria.
\item Besides common factors, we still consider other factors, such as evaporation
and radiation heat transfer. The evaporation turns out to be the main reason of
heat loss, which corresponds with other scientists experimental outcome.
\end{itemize}
\subsection{Weakness}
\begin{itemize}
\item Having knowing the range of some parameters from others essays, we choose
a value from them to apply in our model. Those values may not be reasonable in
reality.
\item Although we investigate a lot in the influence of personal motions, they
are so complicated that need to be studied further.
\item Limited to time, we do not conduct sensitivity analysis for the influence
of personal surface area.
\end{itemize}
\section{Further Discussion}
Based on our model analysis and conclusions, we propose the optimal strategy
for the user in a bathtub and explain the reason for the uneven temperature
throughout the bathtub. In addition, we make improvements for applying our
model in real life, as suggested by the patent \textcite{patent2023}.
In this part, we will focus on different distribution of inflow faucets. Then we
discuss about the real-life application of our model.
\begin{itemize}
\item Different Distribution of Inflow Faucets
In our before discussion, we assume there being just one entrance of inflow.
From the simulating outcome, we find the temperature of bath water is hardly even.
So we come up with the idea of adding more entrances.
The simulation turns out to be as follows
\begin{figure}[h]
\centering
\includegraphics[width=12cm]{fig24.jpg}
\caption{The simulation results of different ways of arranging entrances} \label{fig24}
\end{figure}
From the above figure, the more the entrances are, the evener the temperature
will be. Recalling on the before simulation outcome, when there is only one
entrance for inflow, the temperature of corners is quietly lower than the middle
area.
In conclusion, if we design more entrances, it will be easier to realize the goal
to keep temperature even throughout the bathtub.
\item Model Application
Our before discussion is based on ideal assumptions. In reality, we have to make
some corrections and improvement.
\begin{itemize}
\item[1)] Adding hot water continually with the mass flow of 0.16 kg/s. This way
can ensure even mean temperature throughout the bathtub and waste less water.
\item[2)] The manufacturers can design an intelligent control system to monitor
the temperature so that users can get more enjoyable bath experience.
\item[3)] We recommend users to add bubble additives to slow down the water being
cooler and help cleanse. The additives with lower thermal conductivity are optimal.
\item[4)] The study method of our establishing model can be applied in other area
relative to convection heat transfer, such as air conditioners.
\end{itemize}
\end{itemize}
\printbibliography
\newpage
\begin{letter}{Enjoy Your Bath Time!}
From simulation results of real-life situations, we find it takes a period of
time for the inflow hot water to spread throughout the bathtub. During this
process, the bath water continues transferring heat into air, bathtub and the
person in bathtub. The difference between heat transfer capacity makes the
temperature of various areas to be different. So that it is difficult to get
an evenly maintained temperature throughout the bath water.
In order to enjoy a comfortable bath with even temperature of bath water and
without wasting too much water, we propose the following suggestions.
\begin{itemize}
\item Adding hot water consistently
\item Using smaller bathtub if possible
\item Decreasing motions during bath
\item Using bubble bath additives
\item Arranging more faucets of inflow
\end{itemize}
\vspace{\parskip}
Sincerely yours,
Your friends
\end{letter}
\newpage
\begin{appendices}
\section{First appendix}
In addition, your report must include a letter to the Chief Financial Officer
(CFO) of the Goodgrant Foundation, Mr. Alpha Chiang, that describes the optimal
investment strategy, your modeling approach and major results, and a brief
discussion of your proposed concept of a return-on-investment (ROI). This letter
should be no more than two pages in length.
Here are simulation programmes we used in our model as follow (\cite{Liu02}).\\
\textbf{\textcolor[rgb]{0.98,0.00,0.00}{Input matlab source:}}
\lstinputlisting[language=Matlab]{./code/mcmthesis-matlab1.m}
\section{Second appendix}
some more text \textcolor[rgb]{0.98,0.00,0.00}{\textbf{Input C++ source:}}
\lstinputlisting[language=C++]{./code/mcmthesis-sudoku.cpp}
\end{appendices}
\newpage
\newcounter{lastpage}
\setcounter{lastpage}{\value{page}}
\thispagestyle{empty}
\section*{Report on Use of AI}
\begin{enumerate}
\item OpenAI ChatGPT (Nov 5, 2023 version, ChatGPT-4,)
\begin{description}
\item[Query1:] <insert the exact wording you input into the AI tool>
\item[Output:] <insert the complete output from the AI tool>
\end{description}
\item OpenAI Ernie (Nov 5, 2023 version, Ernie 4.0)
\begin{description}
\item[Query1:] <insert the exact wording of any subsequent input into the AI tool>
\item[Output:] <insert the complete output from the second query>
\end{description}
\item Github CoPilot (Feb 3, 2024 version)
\begin{description}
\item[Query1:] <insert the exact wording you input into the AI tool>
\item[Output:] <insert the complete output from the AI tool>
\end{description}
\item Google Bard (Feb 2, 2024 version)
\begin{description}
\item[Query1:] <insert the exact wording of your query>
\item[Output:] <insert the complete output from the AI tool>
\end{description}
\end{enumerate}
% 重置页码
\clearpage
\setcounter{page}{\value{lastpage}}
\end{document}
%%
%% This work consists of these files mcmthesis.dtx,
%% figures/ and
%% code/,
%% and the derived files mcmthesis.cls,
%% mcmthesis-demo.tex,
%% README,
%% LICENSE,
%% mcmthesis.pdf and
%% mcmthesis-demo.pdf.
%%
%% End of file `mcmthesis-demo.tex'.

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "80,21",
"scrollLine": "72"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "215,17",
"scrollLine": "215"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "2,52",
"scrollLine": "0"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "541,21",
"scrollLine": "523"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "555,14",
"scrollLine": "540"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "65,42",
"scrollLine": "52"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "75,1",
"scrollLine": "50"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "39,55",
"scrollLine": "29"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "112,29",
"scrollLine": "106"
}

View File

@@ -0,0 +1,4 @@
{
"source_window_id": "",
"Source": "Source"
}

View File

@@ -0,0 +1,4 @@
{
"source_window_id": "",
"Source": "Source"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "715,74",
"scrollLine": "701"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "64,4",
"scrollLine": "57"
}

View File

@@ -0,0 +1,7 @@
{
"tempName": "Untitled1",
"source_window_id": "",
"Source": "Source",
"cursorPosition": "2,0",
"scrollLine": "0"
}

View File

@@ -0,0 +1,14 @@
C%3A%2FUsers%2Fzhang%2FDesktop%2F24%E7%BE%8E%E8%B5%9B%2FMCM_Latex2024%2Fmcmthesis-demo.tex="4C9EBF2C"
C%3A%2FUsers%2Fzhang%2FDesktop%2F24%E7%BE%8E%E8%B5%9B%2FMCM_Latex2024%2Fmcmthesis.cls="352D9BE2"
C%3A%2FUsers%2Fzhang%2FDesktop%2F25%E7%BE%8E%E8%B5%9B%2F%E5%BC%A0%E6%95%AC%E4%BF%A1-%E7%BE%8E%E8%B5%9B%E5%B7%A5%E5%85%B7%E6%8A%80%E8%83%BD%E5%9F%B9%E8%AE%AD%2FMCM_Latex2025%2Fmcmthesis-demo.tex="657AF765"
C%3A%2FUsers%2Fzhang%2FDesktop%2F25%E7%BE%8E%E8%B5%9B%2F%E5%BC%A0%E6%95%AC%E4%BF%A1-%E7%BE%8E%E8%B5%9B%E5%B7%A5%E5%85%B7%E6%8A%80%E8%83%BD%E5%9F%B9%E8%AE%AD%2FMCM_Latex2025%2Fmcmthesis.cls="12468681"
C%3A%2FUsers%2Fzhang%2FDesktop%2F25%E7%BE%8E%E8%B5%9B%2F%E5%BC%A0%E6%95%AC%E4%BF%A1-%E7%BE%8E%E8%B5%9B%E5%B7%A5%E5%85%B7%E6%8A%80%E8%83%BD%E5%9F%B9%E8%AE%AD%2FMCM_Latex2025%2Frun.R="CD63E04B"
C%3A%2FUsers%2Fzhang%2FDesktop%2F25%E7%BE%8E%E8%B5%9B%2FMCM_Latex2025%2Fmcmthesis-demo.tex="8F0F20BF"
C%3A%2FUsers%2Fzhang%2FDesktop%2F25%E7%BE%8E%E8%B5%9B%2FMCM_Latex2025%2Fmcmthesis.cls="8E064A1A"
C%3A%2FUsers%2Fzhang%2FDesktop%2F25%E7%BE%8E%E8%B5%9B%2FMCM_Latex2025%2Freference.bib="6C2191AE"
C%3A%2FUsers%2Fzhang%2FDesktop%2F25%E7%BE%8E%E8%B5%9B%2FMCM_Latex2025%2Frun.R="3F2C1C61"
F%3A%2F%E7%BE%8E%E8%B5%9B%2F23%E7%BE%8E%E8%B5%9B%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2023)%2Fmcmthesis-demo.tex="8DF18780"
F%3A%2F%E7%BE%8E%E8%B5%9B%2F23%E7%BE%8E%E8%B5%9B%2FMCM_Latex2023%2FMCM-ICM_Summary.tex="6C237C7D"
F%3A%2F%E7%BE%8E%E8%B5%9B%2F23%E7%BE%8E%E8%B5%9B%2FMCM_Latex2024%2Fmcmthesis-demo.tex="435EDAB8"
F%3A%2F%E7%BE%8E%E8%B5%9B%2F23%E7%BE%8E%E8%B5%9B%2FMCM_Latex2024%2Fmcmthesis.cls="882CBEC9"
c%3A%2FUsers%2Fzhang%2FAppData%2FRoaming%2FTinyTeX%2Ftexmf-dist%2Ftex%2Flatex%2Fxkeyval%2Fxkeyval.sty="964D1E3B"

View File

@@ -0,0 +1,9 @@
{
"sortOrder": [
{
"columnIndex": 2,
"ascending": true
}
],
"path": "C:/Users/zhang/Desktop/23美赛/MCM_Latex模板(2023)"
}

View File

@@ -0,0 +1,3 @@
{
"activeTab": 0
}

View File

@@ -0,0 +1,14 @@
{
"left": {
"splitterpos": 349,
"topwindowstate": "NORMAL",
"panelheight": 817,
"windowheight": 855
},
"right": {
"splitterpos": 509,
"topwindowstate": "NORMAL",
"panelheight": 817,
"windowheight": 855
}
}

View File

@@ -0,0 +1,5 @@
{
"TabSet1": 0,
"TabSet2": 0,
"TabZoom": {}
}

View File

@@ -0,0 +1,5 @@

View File

@@ -0,0 +1 @@
{"active_set":"","sets":[]}

View File

@@ -0,0 +1,26 @@
{
"id": "5A8D7EA9",
"path": "C:/Users/zhang/Desktop/23美赛/MCM_Latex模板(2023)/mcmthesis-demo.tex",
"project_path": "mcmthesis-demo.tex",
"type": "tex",
"hash": "3389638255",
"contents": "",
"dirty": false,
"created": 1672236598356.0,
"source_on_save": false,
"relative_order": 1,
"properties": {
"source_window_id": "",
"Source": "Source",
"cursorPosition": "13,9",
"scrollLine": "35"
},
"folds": "",
"lastKnownWriteTime": 1643251050,
"encoding": "UTF-8",
"collab_server": "",
"source_window": "",
"last_content_update": 1643251050,
"read_only": false,
"read_only_alternatives": []
}

View File

@@ -0,0 +1,539 @@
%%
%% This is file `mcmthesis-demo.tex',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `demo')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
%%
%% This is file `mcmthesis-demo.tex',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `demo')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
\documentclass{mcmthesis}
\mcmsetup{CTeX = false, % 使用 CTeX 套装时,设置为 true
tcn = {0000000}, problem = \textcolor{red}{A},
sheet = true, titleinsheet = true, keywordsinsheet = true,
titlepage = false, abstract = false}
\usepackage{newtxtext} % \usepackage{palatino}
\usepackage[backend=bibtex]{biblatex} % for RStudio Complie
\usepackage{tocloft}
\setlength{\cftbeforesecskip}{6pt}
\renewcommand{\contentsname}{\hspace*{\fill}\Large\bfseries Contents \hspace*{\fill}}
\title{Enjoy a Cozy and Green Bath}
% \author{\small \href{http://www.latexstudio.net/}
% {\includegraphics[width=7cm]{mcmthesis-logo}}}
\date{\today}
\begin{document}
\begin{abstract}
A traditional bathtub cannot be reheated by itself, so users have to add hot water from time to time. Our goal is to establish a model of the temperature of bath water in space and time. Then we are expected to propose an optimal strategy for users to keep the temperature even and close to initial temperature and decrease water consumption.
To simplify modeling process, we firstly assume there is no person in the bathtub. We regard the whole bathtub as a thermodynamic system and introduce heat transfer formulas.
We establish two sub-models: adding water constantly and discontinuously. As for the former sub-model, we define the mean temperature of bath water. Introducing Newton cooling formula, we determine the heat transfer capacity. After deriving the value of parameters, we deduce formulas to derive results and simulate the change of temperature field via CFD. As for the second sub-model, we define an iteration consisting of two process: heating and standby. According to energy conservation law, we obtain the relationship of time and total heat dissipating capacity. Then we determine the mass flow and the time of adding hot water. We also use CFD to simulate the temperature field in second sub-model.
In consideration of evaporation, we correct the results of sub-models referring to some scientists' studies. We define two evaluation criteria and compare the two sub-models. Adding water constantly is found to keep the temperature of bath water even and avoid wasting too much water, so it is recommended by us.
Then we determine the influence of some factors: radiation heat transfer, the shape and volume of the tub, the shape/volume/temperature/motions of the person, the bubbles made from bubble bath additives. We focus on the influence of those factors to heat transfer and then conduct sensitivity analysis. The results indicate smaller bathtub with less surface area, lighter personal mass, less motions and more bubbles will decrease heat transfer and save water.
Based on our model analysis and conclusions, we propose the optimal strategy for the user in a bathtub and explain the reason of uneven temperature throughout the bathtub. In addition, we make improvement for applying our model in real life.
\begin{keywords}
Heat transfer, Thermodynamic system, CFD, Energy conservation
\end{keywords}
\end{abstract}
\maketitle
%% Generate the Table of Contents, if it's needed.
% \renewcommand{\contentsname}{\centering Contents}
\tableofcontents % 若不想要目录, 注释掉该句
\thispagestyle{empty}
\newpage
\section{Introduction}
\subsection{Background}
Bathing in a tub is a perfect choice for those who have been worn out after a long day's working. A traditional bathtub is a simply water containment vessel without a secondary heating system or circulating jets. Thus the temperature of water in bathtub declines noticeably as time goes by, which will influent the experience of bathing. As a result, the bathing person needs to add a constant trickle of hot water from a faucet to reheat the bathing water. This way of bathing will result in waste of water because when the capacity of the bathtub is reached, excess water overflows the tub.
An optimal bathing strategy is required for the person in a bathtub to get comfortable bathing experience while reducing the waste of water.
\subsection{Literature Review}
Korean physicist Gi-Beum Kim analyzed heat loss through free surface of water contained in bathtub due to conduction and evaporation \cite{1}. He derived a relational equation based on the basic theory of heat transfer to evaluate the performance of bath tubes. The major heat loss was found to be due to evaporation. Moreover, he found out that the speed of heat loss depends more on the humidity of the bathroom than the temperature of water contained in the bathtub. So, it is best to maintain the temperature of bathtub water to be between 41 to 45$^{\circ}$C and the humidity of bathroom to be 95\%.
When it comes to convective heat transfer in bathtub, many studies can be referred to. Newton's law of cooling states that the rate of heat loss of a body is proportional to the difference in temperatures between the body and its surroundings while under the effects of a breeze \cite{2}. Claude-Louis Navier and George Gabriel Stokes described the motion of viscous fluid substances with the Navier-Stokes equations. Those equations may be used to model the weather, ocean currents, water flow in a pipe and air flow around a wing \cite{3}.
In addition, some numerical simulation software are applied in solving and analyzing problems that involve fluid flows. For example, Computational Fluid Dynamics (CFD) is a common one used to perform the calculations required to simulate the interaction of liquids and gases with surfaces defined by boundary conditions \cite{4}.
\subsection{Restatement of the Problem}
We are required to establish a model to determine the change of water temperature in space and time. Then we are expected to propose the best strategy for the person in the bathtub to keep the water temperature close to initial temperature and even throughout the tub. Reduction of waste of water is also needed. In addition, we have to consider the impact of different conditions on our model, such as different shapes and volumes of the bathtub, etc.
In order to solve those problems, we will proceed as follows:
\begin{itemize}
\item {\bf Stating assumptions}. By stating our assumptions, we will narrow the focus of our approach towards the problems and provide some insight into bathtub water temperature issues.
\item {\bf Making notations}. We will give some notations which are important for us to clarify our models.
\item {\bf Presenting our model}. In order to investigate the problem deeper, we divide our model into two sub-models. One is a steady convection heat transfer sub-model in which hot water is added constantly. The other one is an unsteady convection heat transfer sub-model where hot water is added discontinuously.
\item {Defining evaluation criteria and comparing sub-models}. We define two main criteria to evaluate our model: the mean temperature of bath water and the amount of inflow water.
\item {\bf Analysis of influencing factors}. In term of the impact of different factors on our model, we take those into consideration: the shape and volume of the tub, the shape/volume/temperature of the person in the bathtub, the motions made by the person in the bathtub and adding a bubble bath additive initially.
\item {\bf Model testing and sensitivity analysis}. With the criteria defined before, we evaluate the reliability of our model and do the sensitivity analysis.
\item {\bf Further discussion}. We discuss about different ways to arrange inflow faucets. Then we improve our model to apply them in reality.
\item {\bf Evaluating the model}. We discuss about the strengths and weaknesses of our model:
\begin{itemize}
\item[1)] ...
\item[2)] ...
\item[3)] ...
\item[4)] ...
\end{itemize}
\end{itemize}
\section{Assumptions and Justification}
To simplify the problem and make it convenient for us to simulate real-life conditions, we make the following basic assumptions, each of which is properly justified.
\begin{itemize}
\item {\bf The bath water is incompressible Non-Newtonian fluid}. The incompressible Non-Newtonian fluid is the basis of NavierStokes equations which are introduced to simulate the flow of bath water.
\item {\bf All the physical properties of bath water, bathtub and air are assumed to be stable}. The change of those properties like specific heat, thermal conductivity and density is rather small according to some studies \cite{5}. It is complicated and unnecessary to consider these little change so we ignore them.
\item {\bf There is no internal heat source in the system consisting of bathtub, hot water and air}. Before the person lies in the bathtub, no internal heat source exist except the system components. The circumstance where the person is in the bathtub will be investigated in our later discussion.
\item {\bf We ignore radiative thermal exchange}. According to Stefan-Boltzmanns law, the radiative thermal exchange can be ignored when the temperature is low. Refer to industrial standard \cite{6}, the temperature in bathroom is lower than 100 $^{\circ}$C, so it is reasonable for us to make this assumption.
\item {\bf The temperature of the adding hot water from the faucet is stable}. This hypothesis can be easily achieved in reality and will simplify our process of solving the problem.
\end{itemize}
\section{Notations}
\begin{center}
\begin{tabular}{clc}
{\bf Symbols} & {\bf Description} & \quad {\bf Unit} \\[0.25cm]
$h$ & Convection heat transfer coefficient & \quad W/(m$^2 \cdot$ K)
\\[0.2cm]
$k$ & Thermal conductivity & \quad W/(m $\cdot$ K) \\[0.2cm]
$c_p$ & Specific heat & \quad J/(kg $\cdot$ K) \\[0.2cm]
$\rho$ & Density & \quad kg/m$^2$ \\[0.2cm]
$\delta$ & Thickness & \quad m \\[0.2cm]
$t$ & Temperature & \quad $^\circ$C, K \\[0.2cm]
$\tau$ & Time & \quad s, min, h \\[0.2cm]
$q_m$ & Mass flow & \quad kg/s \\[0.2cm]
$\Phi$ & Heat transfer power & \quad W \\[0.2cm]
$T$ & A period of time & \quad s, min, h \\[0.2cm]
$V$ & Volume & \quad m$^3$, L \\[0.2cm]
$M,\,m$ & Mass & \quad kg \\[0.2cm]
$A$ & Aera & \quad m$^2$ \\[0.2cm]
$a,\,b,\,c$ & The size of a bathtub & \quad m$^3$
\end{tabular}
\end{center}
\noindent where we define the main parameters while specific value of those parameters will be given later.
\section{Model Overview}
In our basic model, we aim at three goals: keeping the temperature as even as possible, making it close to the initial temperature and decreasing the water consumption.
We start with the simple sub-model where hot water is added constantly.
At first we introduce convection heat transfer control equations in rectangular coordinate system. Then we define the mean temperature of bath water.
Afterwards, we introduce Newton cooling formula to determine heat transfer
capacity. After deriving the value of parameters, we get calculating results via formula deduction and simulating results via CFD.
Secondly, we present the complicated sub-model in which hot water is
added discontinuously. We define an iteration consisting of two process:
heating and standby. As for heating process, we derive control equations and boundary conditions. As for standby process, considering energy conservation law, we deduce the relationship of total heat dissipating capacity and time.
Then we determine the time and amount of added hot water. After deriving the value of parameters, we get calculating results via formula deduction and simulating results via CFD.
At last, we define two criteria to evaluate those two ways of adding hot water. Then we propose optimal strategy for the user in a bathtub.
The whole modeling process can be shown as follows.
\begin{figure}[h]
\centering
\includegraphics[width=12cm]{fig1.jpg}
\caption{Modeling process} \label{fig1}
\end{figure}
\section{Sub-model I : Adding Water Continuously}
We first establish the sub-model based on the condition that a person add water continuously to reheat the bathing water. Then we use Computational Fluid Dynamics (CFD) to simulate the change of water temperature in the bathtub. At last, we evaluate the model with the criteria which have been defined before.
\subsection{Model Establishment}
Since we try to keep the temperature of the hot water in bathtub to be even, we have to derive the amount of inflow water and the energy dissipated by the hot water into the air.
We derive the basic convection heat transfer control equations based on the former scientists achievement. Then, we define the mean temperature of bath water. Afterwards, we determine two types of heat transfer: the boundary heat transfer and the evaporation heat transfer. Combining thermodynamic formulas, we derive calculating results. Via Fluent software, we get simulation results.
\subsubsection{Control Equations and Boundary Conditions}
According to thermodynamics knowledge, we recall on basic convection
heat transfer control equations in rectangular coordinate system. Those
equations show the relationship of the temperature of the bathtub water in space.
We assume the hot water in the bathtub as a cube. Then we put it into a
rectangular coordinate system. The length, width, and height of it is $a,\, b$ and $c$.
\begin{figure}[h]
\centering
\includegraphics[width=8cm]{fig2.jpg}
\caption{Modeling process} \label{fig2}
\end{figure}
In the basis of this, we introduce the following equations \cite{5}:
\begin{itemize}
\item {\bf Continuity equation:}
\end{itemize}
\begin{equation} \label{eq1}
\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} +\frac{\partial w}{\partial z} =0
\end{equation}
\noindent where the first component is the change of fluid mass along the $X$-ray. The second component is the change of fluid mass along the $Y$-ray. And the third component is the change of fluid mass along the $Z$-ray. The sum of the change in mass along those three directions is zero.
\begin{itemize}
\item {\bf Moment differential equation (N-S equations):}
\end{itemize}
\begin{equation} \label{eq2}
\left\{
\begin{array}{l} \!\!
\rho \Big(u \dfrac{\partial u}{\partial x} + v \dfrac{\partial u}{\partial y} + w\dfrac{\partial u}{\partial z} \Big) = -\dfrac{\partial p}{\partial x} + \eta \Big(\dfrac{\partial^2 u}{\partial x^2} + \dfrac{\partial^2 u}{\partial y^2} + \dfrac{\partial^2 u}{\partial z^2} \Big) \\[0.3cm]
\rho \Big(u \dfrac{\partial v}{\partial x} + v \dfrac{\partial v}{\partial y} + w\dfrac{\partial v}{\partial z} \Big) = -\dfrac{\partial p}{\partial y} + \eta \Big(\dfrac{\partial^2 v}{\partial x^2} + \dfrac{\partial^2 v}{\partial y^2} + \dfrac{\partial^2 v}{\partial z^2} \Big) \\[0.3cm]
\rho \Big(u \dfrac{\partial w}{\partial x} + v \dfrac{\partial w}{\partial y} + w\dfrac{\partial w}{\partial z} \Big) = -g-\dfrac{\partial p}{\partial z} + \eta \Big(\dfrac{\partial^2 w}{\partial x^2} + \dfrac{\partial^2 w}{\partial y^2} + \dfrac{\partial^2 w}{\partial z^2} \Big)
\end{array}
\right.
\end{equation}
\begin{itemize}
\item {\bf Energy differential equation:}
\end{itemize}
\begin{equation} \label{eq3}
\rho c_p \Big( u\frac{\partial t}{\partial x} + v\frac{\partial t}{\partial y} + w\frac{\partial t}{\partial z} \Big) = \lambda \Big(\frac{\partial^2 t}{\partial x^2} + \frac{\partial^2 t}{\partial y^2} + \frac{\partial^2 t}{\partial z^2} \Big)
\end{equation}
\noindent where the left three components are convection terms while the right three components are conduction terms.
By Equation \eqref{eq3}, we have ......
......
On the right surface in Fig. \ref{fig2}, the water also transfers heat firstly with bathtub inner surfaces and then the heat comes into air. The boundary condition here is ......
\subsubsection{Definition of the Mean Temperature}
......
\subsubsection{Determination of Heat Transfer Capacity}
......
\section{Sub-model II: Adding Water Discontinuously}
In order to establish the unsteady sub-model, we recall on the working principle of air conditioners. The heating performance of air conditions consist of two processes: heating and standby. After the user set a temperature, the air conditioner will begin to heat until the expected temperature is reached. Then it will go standby. When the temperature get below the expected temperature, the air conditioner begin to work again. As it works in this circle, the temperature remains the expected one.
Inspired by this, we divide the bathtub working into two processes: adding
hot water until the expected temperature is reached, then keeping this
condition for a while unless the temperature is lower than a specific value. Iterating this circle ceaselessly will ensure the temperature kept relatively stable.
\subsection{Heating Model}
\subsubsection{Control Equations and Boundary Conditions}
\subsubsection{Determination of Inflow Time and Amount}
\subsection{Standby Model}
\subsection{Results}
\quad~ We first give the value of parameters based on others studies. Then we get the calculation results and simulating results via those data.
\subsubsection{Determination of Parameters}
After establishing the model, we have to determine the value of some
important parameters.
As scholar Beum Kim points out, the optimal temperature for bath is
between 41 and 45$^\circ$C [1]. Meanwhile, according to Shimodozono's study, 41$^\circ$C warm water bath is the perfect choice for individual health [2]. So it is reasonable for us to focus on $41^\circ$C $\sim 45^\circ$C. Because adding hot water continuously is a steady process, so the mean temperature of bath water is supposed to be constant. We value the temperature of inflow and outflow water with the maximum and minimum temperature respectively.
The values of all parameters needed are shown as follows:
.....
\subsubsection{Calculating Results}
Putting the above value of parameters into the equations we derived before, we can get the some data as follows:
%%普通表格
\begin{table}[h] %h表示固定在当前位置
\centering %设置居中
\caption{The calculating results} %表标题
\vspace{0.15cm}
\label{tab2} %设置表的引用标签
\begin{tabular}{|c|c|c|} %3个c表示3列, |可选, 表示绘制各列间的竖线
\hline %画横线
Variables & Values & Unit \\ \hline %各列间用&隔开
$A_1$ & 1.05 & $m^2$ \\ \hline
$A_2$ & 2.24 & $m^2$ \\ \hline
$\Phi_1$ & 189.00 & $W$ \\ \hline
$\Phi_2$ & 43.47 & $W$ \\ \hline
$\Phi$ & 232.47 & $W$ \\ \hline
$q_m$ & 0.014 & $g/s$ \\ \hline
\end{tabular}
\end{table}
From Table \ref{tab2}, ......
......
\section{Correction and Contrast of Sub-Models}
After establishing two basic sub-models, we have to correct them in consideration of evaporation heat transfer. Then we define two evaluation criteria to compare the two sub-models in order to determine the optimal bath strategy.
\subsection{Correction with Evaporation Heat Transfer}
Someone may confuse about the above results: why the mass flow in the first sub-model is so small? Why the standby time is so long? Actually, the above two sub-models are based on ideal conditions without consideration of the change of boundary conditions, the motions made by the person in bathtub and the evaporation of bath water, etc. The influence of personal motions will be discussed later. Here we introducing the evaporation of bath water to correct sub-models.
\subsection{Contrast of Two Sub-Models}
Firstly we define two evaluation criteria. Then we contrast the two submodels via these two criteria. Thus we can derive the best strategy for the person in the bathtub to adopt.
\section{Model Analysis and Sensitivity Analysis}
\subsection{The Influence of Different Bathtubs}
Definitely, the difference in shape and volume of the tub affects the
convection heat transfer. Examining the relationship between them can help
people choose optimal bathtubs.
\subsubsection{Different Volumes of Bathtubs}
In reality, a cup of water will be cooled down rapidly. However, it takes quite long time for a bucket of water to become cool. That is because their volume is different and the specific heat of water is very large. So that the decrease of temperature is not obvious if the volume of water is huge. That also explains why it takes 45 min for 320 L water to be cooled by 1$^\circ$C.
In order to examine the influence of volume, we analyze our sub-models
by conducting sensitivity Analysis to them.
We assume the initial volume to be 280 L and change it by $\pm 5$\%, $\pm 8$\%, $\pm 12$\% and $\pm 15$\%. With the aid of sub-models we established before, the variation of some parameters turns out to be as follows
%%三线表
\begin{table}[h] %h表示固定在当前位置
\centering %设置居中
\caption{Variation of some parameters} %表标题
\label{tab7} %设置表的引用标签
\begin{tabular}{ccccccc} %7个c表示7列, c表示每列居中对齐, 还有l和r可选
\toprule %画顶端横线
$V$ & $A_1$ & $A_2$ & $T_2$ & $q_{m1}$ & $q_{m2}$ & $\Phi_q$ \\
\midrule %画中间横线
-15.00\% & -5.06\% & -9.31\% & -12.67\% & -2.67\% & -14.14\% & -5.80\% \\
-12.00\% & -4.04\% & -7.43\% & -10.09\% & -2.13\% & -11.31\% & -4.63\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
\bottomrule %画底部横线
\end{tabular}
\end{table}
\section{Strength and Weakness}
\subsection{Strength}
\begin{itemize}
\item We analyze the problem based on thermodynamic formulas and laws, so that the model we established is of great validity.
\item Our model is fairly robust due to our careful corrections in consideration of real-life situations and detailed sensitivity analysis.
\item Via Fluent software, we simulate the time field of different areas throughout the bathtub. The outcome is vivid for us to understand the changing process.
\item We come up with various criteria to compare different situations, like water consumption and the time of adding hot water. Hence an overall comparison can be made according to these criteria.
\item Besides common factors, we still consider other factors, such as evaporation and radiation heat transfer. The evaporation turns out to be the main reason of heat loss, which corresponds with other scientists experimental outcome.
\end{itemize}
\subsection{Weakness}
\begin{itemize}
\item Having knowing the range of some parameters from others essays, we choose a value from them to apply in our model. Those values may not be reasonable in reality.
\item Although we investigate a lot in the influence of personal motions, they are so complicated that need to be studied further.
\item Limited to time, we do not conduct sensitivity analysis for the influence of personal surface area.
\end{itemize}
\section{Further Discussion}
In this part, we will focus on different distribution of inflow faucets. Then we discuss about the real-life application of our model.
\begin{itemize}
\item Different Distribution of Inflow Faucets
In our before discussion, we assume there being just one entrance of inflow.
From the simulating outcome, we find the temperature of bath water is hardly even. So we come up with the idea of adding more entrances.
The simulation turns out to be as follows
\begin{figure}[h]
\centering
\includegraphics[width=12cm]{fig24.jpg}
\caption{The simulation results of different ways of arranging entrances} \label{fig24}
\end{figure}
From the above figure, the more the entrances are, the evener the temperature will be. Recalling on the before simulation outcome, when there is only one entrance for inflow, the temperature of corners is quietly lower than the middle area.
In conclusion, if we design more entrances, it will be easier to realize the goal to keep temperature even throughout the bathtub.
\item Model Application
Our before discussion is based on ideal assumptions. In reality, we have to make some corrections and improvement.
\begin{itemize}
\item[1)] Adding hot water continually with the mass flow of 0.16 kg/s. This way can ensure even mean temperature throughout the bathtub and waste less water.
\item[2)] The manufacturers can design an intelligent control system to monitor the temperature so that users can get more enjoyable bath experience.
\item[3)] We recommend users to add bubble additives to slow down the water being cooler and help cleanse. The additives with lower thermal conductivity are optimal.
\item[4)] The study method of our establishing model can be applied in other area relative to convection heat transfer, such as air conditioners.
\end{itemize}
\end{itemize}
\begin{thebibliography}{99}
\addcontentsline{toc}{section}{Reference}
\bibitem{1} Gi-Beum Kim. Change of the Warm Water Temperature for the Development of Smart Healthecare Bathing System. Hwahak konghak. 2006, 44(3): 270-276.
\bibitem{2} \url{https://en.wikipedia.org/wiki/Convective_heat_transfer#Newton.27s_law_of_cooling}
\bibitem{3} \url{https://en.wikipedia.org/wiki/Navier\%E2\%80\%93Stokes_equations}
\bibitem{4} \url{https://en.wikipedia.org/wiki/Computational_fluid_dynamics}
\bibitem{5} Holman J P. Heat Transfer (9th ed.), New York: McGraw-Hill, 2002.
\bibitem{6} Liu Weiguo, Chen Zhaoping, ZhangYin. Matlab program design and application. Beijing: Higher education press, 2002. (In Chinese)
\end{thebibliography}
\newpage
\begin{letter}{Enjoy Your Bath Time!}
From simulation results of real-life situations, we find it takes a period of time for the inflow hot water to spread throughout the bathtub. During this process, the bath water continues transferring heat into air, bathtub and the person in bathtub. The difference between heat transfer capacity makes the temperature of various areas to be different. So that it is difficult to get an evenly maintained temperature throughout the bath water.
In order to enjoy a comfortable bath with even temperature of bath water and without wasting too much water, we propose the following suggestions.
\begin{itemize}
\item Adding hot water consistently
\item Using smaller bathtub if possible
\item Decreasing motions during bath
\item Using bubble bath additives
\item Arranging more faucets of inflow
\end{itemize}
\vspace{\parskip}
Sincerely yours,
Your friends
\end{letter}
\newpage
\begin{appendices}
\section{First appendix}
In addition, your report must include a letter to the Chief Financial Officer (CFO) of the Goodgrant Foundation, Mr. Alpha Chiang, that describes the optimal investment strategy, your modeling approach and major results, and a brief discussion of your proposed concept of a return-on-investment (ROI). This letter should be no more than two pages in length.
Here are simulation programmes we used in our model as follow.\\
\textbf{\textcolor[rgb]{0.98,0.00,0.00}{Input matlab source:}}
\lstinputlisting[language=Matlab]{./code/mcmthesis-matlab1.m}
\section{Second appendix}
some more text \textcolor[rgb]{0.98,0.00,0.00}{\textbf{Input C++ source:}}
\lstinputlisting[language=C++]{./code/mcmthesis-sudoku.cpp}
\end{appendices}
\end{document}
%%
%% This work consists of these files mcmthesis.dtx,
%% figures/ and
%% code/,
%% and the derived files mcmthesis.cls,
%% mcmthesis-demo.tex,
%% README,
%% LICENSE,
%% mcmthesis.pdf and
%% mcmthesis-demo.pdf.
%%
%% End of file `mcmthesis-demo.tex'.

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "192,32",
"scrollLine": "186"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "13,9",
"scrollLine": "4"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "13,9",
"scrollLine": "35"
}

View File

@@ -0,0 +1,4 @@
{
"source_window_id": "",
"Source": "Source"
}

View File

@@ -0,0 +1,4 @@
C%3A%2FUsers%2Fzhang%2FDesktop%2F23%E7%BE%8E%E8%B5%9B%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2023)%2Fmcmthesis-demo.tex="A0F0AEB7"
C%3A%2FUsers%2Fzhang%2FDesktop%2F23%E7%BE%8E%E8%B5%9B%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2023)%2Fmcmthesis.cls="B564D6D4"
F%3A%2F%E7%BE%8E%E8%B5%9B%2F22%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2022%E6%9C%80%E6%96%B0%E4%BF%AE%E6%94%B9)%2Fmcmthesis-demo.tex="8001FD03"
F%3A%2F%E7%BE%8E%E8%B5%9B%2F22%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2022%E6%9C%80%E6%96%B0%E4%BF%AE%E6%94%B9)%2Fmcmthesis.cls="6676DB66"

View File

@@ -0,0 +1,5 @@
{
"debugBreakpointsState": {
"breakpoints": []
}
}

View File

@@ -0,0 +1,9 @@
{
"sortOrder": [
{
"columnIndex": 2,
"ascending": true
}
],
"path": "C:/Users/zhjx_/Desktop/22美赛/美赛模板/MCM_Latex模板(2022最新修改)"
}

View File

@@ -0,0 +1,7 @@
{
"installOptions": {
"installFromRepository": false,
"libraryPath": "D:/R-4.1.2/library",
"installDependencies": true
}
}

View File

@@ -0,0 +1,3 @@
{
"activeTab": 0
}

View File

@@ -0,0 +1,14 @@
{
"left": {
"splitterpos": 401,
"topwindowstate": "NORMAL",
"panelheight": 853,
"windowheight": 891
},
"right": {
"splitterpos": 534,
"topwindowstate": "NORMAL",
"panelheight": 853,
"windowheight": 891
}
}

View File

@@ -0,0 +1,5 @@
{
"TabSet1": 0,
"TabSet2": 0,
"TabZoom": {}
}

View File

@@ -0,0 +1,8 @@
build-last-errors="[]"
build-last-errors-base-dir=""
build-last-outputs="[]"
compile_pdf_state="{\"tab_visible\":false,\"running\":false,\"target_file\":\"\",\"output\":\"\",\"errors\":[]}"
files.monitored-path=""
find-in-files-state="{\"handle\":\"\",\"input\":\"\",\"path\":\"\",\"regex\":false,\"ignoreCase\":false,\"results\":{\"file\":[],\"line\":[],\"lineValue\":[],\"matchOn\":[],\"matchOff\":[],\"replaceMatchOn\":[],\"replaceMatchOff\":[]},\"running\":false,\"replace\":false,\"preview\":false,\"gitFlag\":false,\"replacePattern\":\"\"}"
imageDirtyState="1"
saveActionState="0"

View File

@@ -0,0 +1,5 @@

View File

@@ -0,0 +1 @@
{"active_set":"","sets":[]}

View File

@@ -0,0 +1,26 @@
{
"id": "110CA16B",
"path": "C:/Users/zhjx_/Desktop/22美赛/美赛模板/MCM_Latex模板(2022最新修改)/mcmthesis-demo.tex",
"project_path": "mcmthesis-demo.tex",
"type": "tex",
"hash": "3389638255",
"contents": "",
"dirty": false,
"created": 1642937344902.0,
"source_on_save": false,
"relative_order": 2,
"properties": {
"source_window_id": "",
"Source": "Source",
"cursorPosition": "68,22",
"scrollLine": "47"
},
"folds": "",
"lastKnownWriteTime": 1643251049,
"encoding": "UTF-8",
"collab_server": "",
"source_window": "",
"last_content_update": 1643251049728,
"read_only": false,
"read_only_alternatives": []
}

View File

@@ -0,0 +1,539 @@
%%
%% This is file `mcmthesis-demo.tex',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `demo')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
%%
%% This is file `mcmthesis-demo.tex',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `demo')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
\documentclass{mcmthesis}
\mcmsetup{CTeX = false, % 使用 CTeX 套装时,设置为 true
tcn = {0000000}, problem = \textcolor{red}{A},
sheet = true, titleinsheet = true, keywordsinsheet = true,
titlepage = false, abstract = false}
\usepackage{newtxtext} % \usepackage{palatino}
\usepackage[backend=bibtex]{biblatex} % for RStudio Complie
\usepackage{tocloft}
\setlength{\cftbeforesecskip}{6pt}
\renewcommand{\contentsname}{\hspace*{\fill}\Large\bfseries Contents \hspace*{\fill}}
\title{Enjoy a Cozy and Green Bath}
% \author{\small \href{http://www.latexstudio.net/}
% {\includegraphics[width=7cm]{mcmthesis-logo}}}
\date{\today}
\begin{document}
\begin{abstract}
A traditional bathtub cannot be reheated by itself, so users have to add hot water from time to time. Our goal is to establish a model of the temperature of bath water in space and time. Then we are expected to propose an optimal strategy for users to keep the temperature even and close to initial temperature and decrease water consumption.
To simplify modeling process, we firstly assume there is no person in the bathtub. We regard the whole bathtub as a thermodynamic system and introduce heat transfer formulas.
We establish two sub-models: adding water constantly and discontinuously. As for the former sub-model, we define the mean temperature of bath water. Introducing Newton cooling formula, we determine the heat transfer capacity. After deriving the value of parameters, we deduce formulas to derive results and simulate the change of temperature field via CFD. As for the second sub-model, we define an iteration consisting of two process: heating and standby. According to energy conservation law, we obtain the relationship of time and total heat dissipating capacity. Then we determine the mass flow and the time of adding hot water. We also use CFD to simulate the temperature field in second sub-model.
In consideration of evaporation, we correct the results of sub-models referring to some scientists' studies. We define two evaluation criteria and compare the two sub-models. Adding water constantly is found to keep the temperature of bath water even and avoid wasting too much water, so it is recommended by us.
Then we determine the influence of some factors: radiation heat transfer, the shape and volume of the tub, the shape/volume/temperature/motions of the person, the bubbles made from bubble bath additives. We focus on the influence of those factors to heat transfer and then conduct sensitivity analysis. The results indicate smaller bathtub with less surface area, lighter personal mass, less motions and more bubbles will decrease heat transfer and save water.
Based on our model analysis and conclusions, we propose the optimal strategy for the user in a bathtub and explain the reason of uneven temperature throughout the bathtub. In addition, we make improvement for applying our model in real life.
\begin{keywords}
Heat transfer, Thermodynamic system, CFD, Energy conservation
\end{keywords}
\end{abstract}
\maketitle
%% Generate the Table of Contents, if it's needed.
% \renewcommand{\contentsname}{\centering Contents}
\tableofcontents % 若不想要目录, 注释掉该句
\thispagestyle{empty}
\newpage
\section{Introduction}
\subsection{Background}
Bathing in a tub is a perfect choice for those who have been worn out after a long day's working. A traditional bathtub is a simply water containment vessel without a secondary heating system or circulating jets. Thus the temperature of water in bathtub declines noticeably as time goes by, which will influent the experience of bathing. As a result, the bathing person needs to add a constant trickle of hot water from a faucet to reheat the bathing water. This way of bathing will result in waste of water because when the capacity of the bathtub is reached, excess water overflows the tub.
An optimal bathing strategy is required for the person in a bathtub to get comfortable bathing experience while reducing the waste of water.
\subsection{Literature Review}
Korean physicist Gi-Beum Kim analyzed heat loss through free surface of water contained in bathtub due to conduction and evaporation \cite{1}. He derived a relational equation based on the basic theory of heat transfer to evaluate the performance of bath tubes. The major heat loss was found to be due to evaporation. Moreover, he found out that the speed of heat loss depends more on the humidity of the bathroom than the temperature of water contained in the bathtub. So, it is best to maintain the temperature of bathtub water to be between 41 to 45$^{\circ}$C and the humidity of bathroom to be 95\%.
When it comes to convective heat transfer in bathtub, many studies can be referred to. Newton's law of cooling states that the rate of heat loss of a body is proportional to the difference in temperatures between the body and its surroundings while under the effects of a breeze \cite{2}. Claude-Louis Navier and George Gabriel Stokes described the motion of viscous fluid substances with the Navier-Stokes equations. Those equations may be used to model the weather, ocean currents, water flow in a pipe and air flow around a wing \cite{3}.
In addition, some numerical simulation software are applied in solving and analyzing problems that involve fluid flows. For example, Computational Fluid Dynamics (CFD) is a common one used to perform the calculations required to simulate the interaction of liquids and gases with surfaces defined by boundary conditions \cite{4}.
\subsection{Restatement of the Problem}
We are required to establish a model to determine the change of water temperature in space and time. Then we are expected to propose the best strategy for the person in the bathtub to keep the water temperature close to initial temperature and even throughout the tub. Reduction of waste of water is also needed. In addition, we have to consider the impact of different conditions on our model, such as different shapes and volumes of the bathtub, etc.
In order to solve those problems, we will proceed as follows:
\begin{itemize}
\item {\bf Stating assumptions}. By stating our assumptions, we will narrow the focus of our approach towards the problems and provide some insight into bathtub water temperature issues.
\item {\bf Making notations}. We will give some notations which are important for us to clarify our models.
\item {\bf Presenting our model}. In order to investigate the problem deeper, we divide our model into two sub-models. One is a steady convection heat transfer sub-model in which hot water is added constantly. The other one is an unsteady convection heat transfer sub-model where hot water is added discontinuously.
\item {Defining evaluation criteria and comparing sub-models}. We define two main criteria to evaluate our model: the mean temperature of bath water and the amount of inflow water.
\item {\bf Analysis of influencing factors}. In term of the impact of different factors on our model, we take those into consideration: the shape and volume of the tub, the shape/volume/temperature of the person in the bathtub, the motions made by the person in the bathtub and adding a bubble bath additive initially.
\item {\bf Model testing and sensitivity analysis}. With the criteria defined before, we evaluate the reliability of our model and do the sensitivity analysis.
\item {\bf Further discussion}. We discuss about different ways to arrange inflow faucets. Then we improve our model to apply them in reality.
\item {\bf Evaluating the model}. We discuss about the strengths and weaknesses of our model:
\begin{itemize}
\item[1)] ...
\item[2)] ...
\item[3)] ...
\item[4)] ...
\end{itemize}
\end{itemize}
\section{Assumptions and Justification}
To simplify the problem and make it convenient for us to simulate real-life conditions, we make the following basic assumptions, each of which is properly justified.
\begin{itemize}
\item {\bf The bath water is incompressible Non-Newtonian fluid}. The incompressible Non-Newtonian fluid is the basis of NavierStokes equations which are introduced to simulate the flow of bath water.
\item {\bf All the physical properties of bath water, bathtub and air are assumed to be stable}. The change of those properties like specific heat, thermal conductivity and density is rather small according to some studies \cite{5}. It is complicated and unnecessary to consider these little change so we ignore them.
\item {\bf There is no internal heat source in the system consisting of bathtub, hot water and air}. Before the person lies in the bathtub, no internal heat source exist except the system components. The circumstance where the person is in the bathtub will be investigated in our later discussion.
\item {\bf We ignore radiative thermal exchange}. According to Stefan-Boltzmanns law, the radiative thermal exchange can be ignored when the temperature is low. Refer to industrial standard \cite{6}, the temperature in bathroom is lower than 100 $^{\circ}$C, so it is reasonable for us to make this assumption.
\item {\bf The temperature of the adding hot water from the faucet is stable}. This hypothesis can be easily achieved in reality and will simplify our process of solving the problem.
\end{itemize}
\section{Notations}
\begin{center}
\begin{tabular}{clc}
{\bf Symbols} & {\bf Description} & \quad {\bf Unit} \\[0.25cm]
$h$ & Convection heat transfer coefficient & \quad W/(m$^2 \cdot$ K)
\\[0.2cm]
$k$ & Thermal conductivity & \quad W/(m $\cdot$ K) \\[0.2cm]
$c_p$ & Specific heat & \quad J/(kg $\cdot$ K) \\[0.2cm]
$\rho$ & Density & \quad kg/m$^2$ \\[0.2cm]
$\delta$ & Thickness & \quad m \\[0.2cm]
$t$ & Temperature & \quad $^\circ$C, K \\[0.2cm]
$\tau$ & Time & \quad s, min, h \\[0.2cm]
$q_m$ & Mass flow & \quad kg/s \\[0.2cm]
$\Phi$ & Heat transfer power & \quad W \\[0.2cm]
$T$ & A period of time & \quad s, min, h \\[0.2cm]
$V$ & Volume & \quad m$^3$, L \\[0.2cm]
$M,\,m$ & Mass & \quad kg \\[0.2cm]
$A$ & Aera & \quad m$^2$ \\[0.2cm]
$a,\,b,\,c$ & The size of a bathtub & \quad m$^3$
\end{tabular}
\end{center}
\noindent where we define the main parameters while specific value of those parameters will be given later.
\section{Model Overview}
In our basic model, we aim at three goals: keeping the temperature as even as possible, making it close to the initial temperature and decreasing the water consumption.
We start with the simple sub-model where hot water is added constantly.
At first we introduce convection heat transfer control equations in rectangular coordinate system. Then we define the mean temperature of bath water.
Afterwards, we introduce Newton cooling formula to determine heat transfer
capacity. After deriving the value of parameters, we get calculating results via formula deduction and simulating results via CFD.
Secondly, we present the complicated sub-model in which hot water is
added discontinuously. We define an iteration consisting of two process:
heating and standby. As for heating process, we derive control equations and boundary conditions. As for standby process, considering energy conservation law, we deduce the relationship of total heat dissipating capacity and time.
Then we determine the time and amount of added hot water. After deriving the value of parameters, we get calculating results via formula deduction and simulating results via CFD.
At last, we define two criteria to evaluate those two ways of adding hot water. Then we propose optimal strategy for the user in a bathtub.
The whole modeling process can be shown as follows.
\begin{figure}[h]
\centering
\includegraphics[width=12cm]{fig1.jpg}
\caption{Modeling process} \label{fig1}
\end{figure}
\section{Sub-model I : Adding Water Continuously}
We first establish the sub-model based on the condition that a person add water continuously to reheat the bathing water. Then we use Computational Fluid Dynamics (CFD) to simulate the change of water temperature in the bathtub. At last, we evaluate the model with the criteria which have been defined before.
\subsection{Model Establishment}
Since we try to keep the temperature of the hot water in bathtub to be even, we have to derive the amount of inflow water and the energy dissipated by the hot water into the air.
We derive the basic convection heat transfer control equations based on the former scientists achievement. Then, we define the mean temperature of bath water. Afterwards, we determine two types of heat transfer: the boundary heat transfer and the evaporation heat transfer. Combining thermodynamic formulas, we derive calculating results. Via Fluent software, we get simulation results.
\subsubsection{Control Equations and Boundary Conditions}
According to thermodynamics knowledge, we recall on basic convection
heat transfer control equations in rectangular coordinate system. Those
equations show the relationship of the temperature of the bathtub water in space.
We assume the hot water in the bathtub as a cube. Then we put it into a
rectangular coordinate system. The length, width, and height of it is $a,\, b$ and $c$.
\begin{figure}[h]
\centering
\includegraphics[width=8cm]{fig2.jpg}
\caption{Modeling process} \label{fig2}
\end{figure}
In the basis of this, we introduce the following equations \cite{5}:
\begin{itemize}
\item {\bf Continuity equation:}
\end{itemize}
\begin{equation} \label{eq1}
\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} +\frac{\partial w}{\partial z} =0
\end{equation}
\noindent where the first component is the change of fluid mass along the $X$-ray. The second component is the change of fluid mass along the $Y$-ray. And the third component is the change of fluid mass along the $Z$-ray. The sum of the change in mass along those three directions is zero.
\begin{itemize}
\item {\bf Moment differential equation (N-S equations):}
\end{itemize}
\begin{equation} \label{eq2}
\left\{
\begin{array}{l} \!\!
\rho \Big(u \dfrac{\partial u}{\partial x} + v \dfrac{\partial u}{\partial y} + w\dfrac{\partial u}{\partial z} \Big) = -\dfrac{\partial p}{\partial x} + \eta \Big(\dfrac{\partial^2 u}{\partial x^2} + \dfrac{\partial^2 u}{\partial y^2} + \dfrac{\partial^2 u}{\partial z^2} \Big) \\[0.3cm]
\rho \Big(u \dfrac{\partial v}{\partial x} + v \dfrac{\partial v}{\partial y} + w\dfrac{\partial v}{\partial z} \Big) = -\dfrac{\partial p}{\partial y} + \eta \Big(\dfrac{\partial^2 v}{\partial x^2} + \dfrac{\partial^2 v}{\partial y^2} + \dfrac{\partial^2 v}{\partial z^2} \Big) \\[0.3cm]
\rho \Big(u \dfrac{\partial w}{\partial x} + v \dfrac{\partial w}{\partial y} + w\dfrac{\partial w}{\partial z} \Big) = -g-\dfrac{\partial p}{\partial z} + \eta \Big(\dfrac{\partial^2 w}{\partial x^2} + \dfrac{\partial^2 w}{\partial y^2} + \dfrac{\partial^2 w}{\partial z^2} \Big)
\end{array}
\right.
\end{equation}
\begin{itemize}
\item {\bf Energy differential equation:}
\end{itemize}
\begin{equation} \label{eq3}
\rho c_p \Big( u\frac{\partial t}{\partial x} + v\frac{\partial t}{\partial y} + w\frac{\partial t}{\partial z} \Big) = \lambda \Big(\frac{\partial^2 t}{\partial x^2} + \frac{\partial^2 t}{\partial y^2} + \frac{\partial^2 t}{\partial z^2} \Big)
\end{equation}
\noindent where the left three components are convection terms while the right three components are conduction terms.
By Equation \eqref{eq3}, we have ......
......
On the right surface in Fig. \ref{fig2}, the water also transfers heat firstly with bathtub inner surfaces and then the heat comes into air. The boundary condition here is ......
\subsubsection{Definition of the Mean Temperature}
......
\subsubsection{Determination of Heat Transfer Capacity}
......
\section{Sub-model II: Adding Water Discontinuously}
In order to establish the unsteady sub-model, we recall on the working principle of air conditioners. The heating performance of air conditions consist of two processes: heating and standby. After the user set a temperature, the air conditioner will begin to heat until the expected temperature is reached. Then it will go standby. When the temperature get below the expected temperature, the air conditioner begin to work again. As it works in this circle, the temperature remains the expected one.
Inspired by this, we divide the bathtub working into two processes: adding
hot water until the expected temperature is reached, then keeping this
condition for a while unless the temperature is lower than a specific value. Iterating this circle ceaselessly will ensure the temperature kept relatively stable.
\subsection{Heating Model}
\subsubsection{Control Equations and Boundary Conditions}
\subsubsection{Determination of Inflow Time and Amount}
\subsection{Standby Model}
\subsection{Results}
\quad~ We first give the value of parameters based on others studies. Then we get the calculation results and simulating results via those data.
\subsubsection{Determination of Parameters}
After establishing the model, we have to determine the value of some
important parameters.
As scholar Beum Kim points out, the optimal temperature for bath is
between 41 and 45$^\circ$C [1]. Meanwhile, according to Shimodozono's study, 41$^\circ$C warm water bath is the perfect choice for individual health [2]. So it is reasonable for us to focus on $41^\circ$C $\sim 45^\circ$C. Because adding hot water continuously is a steady process, so the mean temperature of bath water is supposed to be constant. We value the temperature of inflow and outflow water with the maximum and minimum temperature respectively.
The values of all parameters needed are shown as follows:
.....
\subsubsection{Calculating Results}
Putting the above value of parameters into the equations we derived before, we can get the some data as follows:
%%普通表格
\begin{table}[h] %h表示固定在当前位置
\centering %设置居中
\caption{The calculating results} %表标题
\vspace{0.15cm}
\label{tab2} %设置表的引用标签
\begin{tabular}{|c|c|c|} %3个c表示3列, |可选, 表示绘制各列间的竖线
\hline %画横线
Variables & Values & Unit \\ \hline %各列间用&隔开
$A_1$ & 1.05 & $m^2$ \\ \hline
$A_2$ & 2.24 & $m^2$ \\ \hline
$\Phi_1$ & 189.00 & $W$ \\ \hline
$\Phi_2$ & 43.47 & $W$ \\ \hline
$\Phi$ & 232.47 & $W$ \\ \hline
$q_m$ & 0.014 & $g/s$ \\ \hline
\end{tabular}
\end{table}
From Table \ref{tab2}, ......
......
\section{Correction and Contrast of Sub-Models}
After establishing two basic sub-models, we have to correct them in consideration of evaporation heat transfer. Then we define two evaluation criteria to compare the two sub-models in order to determine the optimal bath strategy.
\subsection{Correction with Evaporation Heat Transfer}
Someone may confuse about the above results: why the mass flow in the first sub-model is so small? Why the standby time is so long? Actually, the above two sub-models are based on ideal conditions without consideration of the change of boundary conditions, the motions made by the person in bathtub and the evaporation of bath water, etc. The influence of personal motions will be discussed later. Here we introducing the evaporation of bath water to correct sub-models.
\subsection{Contrast of Two Sub-Models}
Firstly we define two evaluation criteria. Then we contrast the two submodels via these two criteria. Thus we can derive the best strategy for the person in the bathtub to adopt.
\section{Model Analysis and Sensitivity Analysis}
\subsection{The Influence of Different Bathtubs}
Definitely, the difference in shape and volume of the tub affects the
convection heat transfer. Examining the relationship between them can help
people choose optimal bathtubs.
\subsubsection{Different Volumes of Bathtubs}
In reality, a cup of water will be cooled down rapidly. However, it takes quite long time for a bucket of water to become cool. That is because their volume is different and the specific heat of water is very large. So that the decrease of temperature is not obvious if the volume of water is huge. That also explains why it takes 45 min for 320 L water to be cooled by 1$^\circ$C.
In order to examine the influence of volume, we analyze our sub-models
by conducting sensitivity Analysis to them.
We assume the initial volume to be 280 L and change it by $\pm 5$\%, $\pm 8$\%, $\pm 12$\% and $\pm 15$\%. With the aid of sub-models we established before, the variation of some parameters turns out to be as follows
%%三线表
\begin{table}[h] %h表示固定在当前位置
\centering %设置居中
\caption{Variation of some parameters} %表标题
\label{tab7} %设置表的引用标签
\begin{tabular}{ccccccc} %7个c表示7列, c表示每列居中对齐, 还有l和r可选
\toprule %画顶端横线
$V$ & $A_1$ & $A_2$ & $T_2$ & $q_{m1}$ & $q_{m2}$ & $\Phi_q$ \\
\midrule %画中间横线
-15.00\% & -5.06\% & -9.31\% & -12.67\% & -2.67\% & -14.14\% & -5.80\% \\
-12.00\% & -4.04\% & -7.43\% & -10.09\% & -2.13\% & -11.31\% & -4.63\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
-8.00\% & -2.68\% & -4.94\% & -6.68\% & -1.41\% & -7.54\% & -3.07\% \\
\bottomrule %画底部横线
\end{tabular}
\end{table}
\section{Strength and Weakness}
\subsection{Strength}
\begin{itemize}
\item We analyze the problem based on thermodynamic formulas and laws, so that the model we established is of great validity.
\item Our model is fairly robust due to our careful corrections in consideration of real-life situations and detailed sensitivity analysis.
\item Via Fluent software, we simulate the time field of different areas throughout the bathtub. The outcome is vivid for us to understand the changing process.
\item We come up with various criteria to compare different situations, like water consumption and the time of adding hot water. Hence an overall comparison can be made according to these criteria.
\item Besides common factors, we still consider other factors, such as evaporation and radiation heat transfer. The evaporation turns out to be the main reason of heat loss, which corresponds with other scientists experimental outcome.
\end{itemize}
\subsection{Weakness}
\begin{itemize}
\item Having knowing the range of some parameters from others essays, we choose a value from them to apply in our model. Those values may not be reasonable in reality.
\item Although we investigate a lot in the influence of personal motions, they are so complicated that need to be studied further.
\item Limited to time, we do not conduct sensitivity analysis for the influence of personal surface area.
\end{itemize}
\section{Further Discussion}
In this part, we will focus on different distribution of inflow faucets. Then we discuss about the real-life application of our model.
\begin{itemize}
\item Different Distribution of Inflow Faucets
In our before discussion, we assume there being just one entrance of inflow.
From the simulating outcome, we find the temperature of bath water is hardly even. So we come up with the idea of adding more entrances.
The simulation turns out to be as follows
\begin{figure}[h]
\centering
\includegraphics[width=12cm]{fig24.jpg}
\caption{The simulation results of different ways of arranging entrances} \label{fig24}
\end{figure}
From the above figure, the more the entrances are, the evener the temperature will be. Recalling on the before simulation outcome, when there is only one entrance for inflow, the temperature of corners is quietly lower than the middle area.
In conclusion, if we design more entrances, it will be easier to realize the goal to keep temperature even throughout the bathtub.
\item Model Application
Our before discussion is based on ideal assumptions. In reality, we have to make some corrections and improvement.
\begin{itemize}
\item[1)] Adding hot water continually with the mass flow of 0.16 kg/s. This way can ensure even mean temperature throughout the bathtub and waste less water.
\item[2)] The manufacturers can design an intelligent control system to monitor the temperature so that users can get more enjoyable bath experience.
\item[3)] We recommend users to add bubble additives to slow down the water being cooler and help cleanse. The additives with lower thermal conductivity are optimal.
\item[4)] The study method of our establishing model can be applied in other area relative to convection heat transfer, such as air conditioners.
\end{itemize}
\end{itemize}
\begin{thebibliography}{99}
\addcontentsline{toc}{section}{Reference}
\bibitem{1} Gi-Beum Kim. Change of the Warm Water Temperature for the Development of Smart Healthecare Bathing System. Hwahak konghak. 2006, 44(3): 270-276.
\bibitem{2} \url{https://en.wikipedia.org/wiki/Convective_heat_transfer#Newton.27s_law_of_cooling}
\bibitem{3} \url{https://en.wikipedia.org/wiki/Navier\%E2\%80\%93Stokes_equations}
\bibitem{4} \url{https://en.wikipedia.org/wiki/Computational_fluid_dynamics}
\bibitem{5} Holman J P. Heat Transfer (9th ed.), New York: McGraw-Hill, 2002.
\bibitem{6} Liu Weiguo, Chen Zhaoping, ZhangYin. Matlab program design and application. Beijing: Higher education press, 2002. (In Chinese)
\end{thebibliography}
\newpage
\begin{letter}{Enjoy Your Bath Time!}
From simulation results of real-life situations, we find it takes a period of time for the inflow hot water to spread throughout the bathtub. During this process, the bath water continues transferring heat into air, bathtub and the person in bathtub. The difference between heat transfer capacity makes the temperature of various areas to be different. So that it is difficult to get an evenly maintained temperature throughout the bath water.
In order to enjoy a comfortable bath with even temperature of bath water and without wasting too much water, we propose the following suggestions.
\begin{itemize}
\item Adding hot water consistently
\item Using smaller bathtub if possible
\item Decreasing motions during bath
\item Using bubble bath additives
\item Arranging more faucets of inflow
\end{itemize}
\vspace{\parskip}
Sincerely yours,
Your friends
\end{letter}
\newpage
\begin{appendices}
\section{First appendix}
In addition, your report must include a letter to the Chief Financial Officer (CFO) of the Goodgrant Foundation, Mr. Alpha Chiang, that describes the optimal investment strategy, your modeling approach and major results, and a brief discussion of your proposed concept of a return-on-investment (ROI). This letter should be no more than two pages in length.
Here are simulation programmes we used in our model as follow.\\
\textbf{\textcolor[rgb]{0.98,0.00,0.00}{Input matlab source:}}
\lstinputlisting[language=Matlab]{./code/mcmthesis-matlab1.m}
\section{Second appendix}
some more text \textcolor[rgb]{0.98,0.00,0.00}{\textbf{Input C++ source:}}
\lstinputlisting[language=C++]{./code/mcmthesis-sudoku.cpp}
\end{appendices}
\end{document}
%%
%% This work consists of these files mcmthesis.dtx,
%% figures/ and
%% code/,
%% and the derived files mcmthesis.cls,
%% mcmthesis-demo.tex,
%% README,
%% LICENSE,
%% mcmthesis.pdf and
%% mcmthesis-demo.pdf.
%%
%% End of file `mcmthesis-demo.tex'.

View File

@@ -0,0 +1,26 @@
{
"id": "F9FF5BF1",
"path": "C:/Users/zhjx_/Desktop/22美赛/美赛模板/MCM_Latex模板(2022最新修改)/mcmthesis.cls",
"project_path": "mcmthesis.cls",
"type": "tex",
"hash": "1079621820",
"contents": "",
"dirty": false,
"created": 1641799370735.0,
"source_on_save": false,
"relative_order": 2,
"properties": {
"source_window_id": "",
"Source": "Source",
"cursorPosition": "151,104",
"scrollLine": "0"
},
"folds": "",
"lastKnownWriteTime": 1643251056,
"encoding": "UTF-8",
"collab_server": "",
"source_window": "",
"last_content_update": 1643251056957,
"read_only": false,
"read_only_alternatives": []
}

View File

@@ -0,0 +1,350 @@
%%
%% This is file `mcmthesis.cls',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% mcmthesis.dtx (with options: `class')
%%
%% -----------------------------------
%%
%% This is a generated file.
%%
%% Copyright (C)
%% 2010 -- 2015 by Zhaoli Wang
%% 2014 -- 2019 by Liam Huang
%% 2019 -- present by latexstudio.net
%%
%% This work may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either version 1.3
%% of this license or (at your option) any later version.
%% The latest version of this license is in
%% http://www.latex-project.org/lppl.txt
%% and version 1.3 or later is part of all distributions of LaTeX
%% version 2005/12/01 or later.
%%
%% This work has the LPPL maintenance status `maintained'.
%%
%% The Current Maintainer of this work is Liam Huang.
%%
\NeedsTeXFormat{LaTeX2e}[1999/12/01]
\ProvidesClass{mcmthesis}
[2020/01/18 v6.3 The Thesis Template Designed For MCM/ICM]
\typeout{The Thesis Template Designed For MCM/ICM}
\def\MCMversion{v6.3}
\RequirePackage{xkeyval}
\RequirePackage{etoolbox}
\define@boolkey{MCM}[MCM@opt@]{CTeX}[false]{}
\define@boolkey{MCM}[MCM@opt@]{titlepage}[true]{}
\define@boolkey{MCM}[MCM@opt@]{abstract}[true]{}
\define@boolkey{MCM}[MCM@opt@]{sheet}[true]{}
\define@boolkey{MCM}[MCM@opt@]{titleinsheet}[false]{}
\define@boolkey{MCM}[MCM@opt@]{keywordsinsheet}[false]{}
\define@cmdkeys{MCM}[MCM@opt@]{tcn,problem}
\define@key{MCM}{tcn}[0000]{\gdef\MCM@opt@tcn{#1}}
\define@key{MCM}{problem}[A]{\gdef\MCM@opt@problem{#1}}
\setkeys{MCM}{tcn=0000,problem=B}
\define@key{mcmthesis.cls}{tcn}[0000]{\gdef\MCM@opt@tcn{#1}}
\define@key{mcmthesis.cls}{problem}[A]{\gdef\MCM@opt@problem{#1}}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{titlepage}{}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{abstract}{}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{sheet}{}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{titleinsheet}{}
\define@boolkey{mcmthesis.cls}[MCM@opt@]{keywordsinsheet}{}
\MCM@opt@sheettrue
\MCM@opt@titlepagetrue
\MCM@opt@titleinsheetfalse
\MCM@opt@keywordsinsheetfalse
\MCM@opt@abstracttrue
\newcommand{\mcmsetup}[1]{\setkeys{MCM}{#1}}
\ProcessOptionsX\relax
\LoadClass[a4paper, 12pt]{article}
\newcommand{\team}{Team \#\ \MCM@opt@tcn}
\RequirePackage{fancyhdr, fancybox}
\RequirePackage{ifthen}
\RequirePackage{lastpage}
\RequirePackage{listings}
\RequirePackage[toc, page, title, titletoc, header]{appendix}
\RequirePackage{paralist}
\RequirePackage{amsthm, amsfonts}
\RequirePackage{amsmath, bm}
\RequirePackage{amssymb, mathrsfs}
\RequirePackage{latexsym}
\RequirePackage{longtable, multirow, hhline, tabularx, array}
\RequirePackage{flafter}
\RequirePackage{pifont, calc}
\RequirePackage{colortbl, booktabs}
\RequirePackage{geometry}
\RequirePackage[T1]{fontenc}
\RequirePackage[scaled]{berasans}
\RequirePackage{hyperref}
\RequirePackage{ifpdf, ifxetex}
\ifMCM@opt@CTeX
\else
\RequirePackage{environ}
\fi
\ifpdf
\RequirePackage{graphicx}
\RequirePackage{epstopdf}
\else
\ifxetex
\RequirePackage{graphicx}
\else
\RequirePackage[dvipdfmx]{graphicx}
\RequirePackage{bmpsize}
\fi
\fi
\RequirePackage[svgnames]{xcolor}
\ifpdf
\hypersetup{hidelinks}
\else
\ifxetex
\hypersetup{hidelinks}
\else
\hypersetup{dvipdfm, hidelinks}
\fi
\fi
\geometry{a4paper, margin = 1in}
\pagestyle{fancy}
\fancyhf{}
\lhead{\small\sffamily \team}
\rhead{\small\sffamily Page \thepage\ of \pageref{LastPage}}
% \rhead{\small\sffamily Page \thepage}
\setlength\parskip{.5\baselineskip}
\renewcommand\tableofcontents{%
\centerline{\normalfont\Large\bfseries\sffamily\contentsname
\@mkboth{%
\MakeUppercase\contentsname}{\MakeUppercase\contentsname}}%
\vskip 5ex%
\@starttoc{toc}%
}
\setcounter{totalnumber}{4}
\setcounter{topnumber}{2}
\setcounter{bottomnumber}{2}
\renewcommand{\textfraction}{0.15}
\renewcommand{\topfraction}{0.85}
\renewcommand{\bottomfraction}{0.65}
\renewcommand{\floatpagefraction}{0.60}
\renewcommand{\figurename}{Figure}
\renewcommand{\tablename}{Table}
\graphicspath{{./}{./img/}{./fig/}{./image/}{./figure/}{./picture/}
{./imgs/}{./figs/}{./images/}{./figures/}{./pictures/}}
\def\maketitle{%
\let\saved@thepage\thepage
\let\thepage\relax
\ifMCM@opt@sheet
\makesheet
\fi
\newpage
\ifMCM@opt@titlepage
\MCM@maketitle
\fi
\newpage
\let\thepage\saved@thepage
\setcounter{page}{2}
\pagestyle{fancy}
}
\def\abstractname{\large{Summary}}
\ifMCM@opt@CTeX
\newbox\@abstract%
\setbox\@abstract\hbox{}%
\long\def\abstract{\bgroup\global\setbox\@abstract\vbox\bgroup\hsize\textwidth\leftskip1cm\rightskip1cm}%
\def\endabstract{\egroup\egroup}
\def\make@abstract{%
\begin{center}
\textbf{\abstractname}
\end{center}
\usebox\@abstract\par
}
\else
\RenewEnviron{abstract}{\xdef\@abstract{\expandonce\BODY}}
\def\make@abstract{%
\begin{center}
\textbf{\abstractname}
\end{center}
\@abstract\par
}
\fi
\newenvironment{letter}[1]{%
\par%
\bgroup\parindent0pt%
\begin{minipage}{5cm}
\flushleft #1%
\end{minipage}}
{\egroup\smallskip}
\def\keywordsname{Keywords}
\ifMCM@opt@CTeX
\newbox\@keywords
\setbox\@keywords\hbox{}
\def\keywords{\global\setbox\@keywords\vbox\bgroup\noindent\leftskip0cm}
\def\endkeywords{\egroup}%
\def\make@keywords{%
\par\hskip.4cm\textbf{\keywordsname}: \usebox\@keywords\hfill\par
}
\else
\NewEnviron{keywords}{\xdef\@keywords{\expandonce\BODY}}
\def\make@keywords{%
\par\noindent\textbf{\keywordsname}:
\@keywords\par
}
\fi
\newcommand{\headset}{{\the\year}\\MCM/ICM\\Summary Sheet}
\newcommand{\problem}[1]{\mcmsetup{problem = #1}}
\def\makesheet{%
\pagestyle{empty}%
\null%
\vspace*{-5pc}%
\begin{center}
\begingroup
\setlength{\parindent}{0pt}
\begin{minipage}[t]{0.33\linewidth}
\bfseries\centering%
Problem Chosen\\[0.7pc]
{\Huge\textbf{\MCM@opt@problem}}\\[2.8pc]
\end{minipage}%
\begin{minipage}[t]{0.33\linewidth}
\centering%
\textbf{\headset}%
\end{minipage}%
\begin{minipage}[t]{0.33\linewidth}
\centering\bfseries%
Team Control Number\\[0.7pc]
{\Huge\textbf{\textcolor{red}{\MCM@opt@tcn}}}\\[2.8pc]
% {\Huge\textbf{\MCM@opt@tcn}}\\[2.8pc]
\end{minipage}\par
\rule{\linewidth}{0.8pt}\par
\endgroup
\vskip 10pt%
\ifMCM@opt@titleinsheet
\normalfont \LARGE \@title \par
\fi
\end{center}
\ifMCM@opt@keywordsinsheet
\make@abstract
\make@keywords
\else
\make@abstract
\fi}
\newcommand{\MCM@maketitle}{%
\begin{center}%
\let \footnote \thanks
{\LARGE \@title \par}%
\vskip 1.5em%
{\large
\lineskip .5em%
\begin{tabular}[t]{c}%
\@author
\end{tabular}\par}%
\vskip 1em%
{\large \@date}%
\end{center}%
\par
\vskip 1.5em%
\ifMCM@opt@abstract%
\make@abstract
\make@keywords
\fi%
}
\def\MCM@memoto{\relax}
\newcommand{\memoto}[1]{\gdef\MCM@memoto{#1}}
\def\MCM@memofrom{\relax}
\newcommand{\memofrom}[1]{\gdef\MCM@memofrom{#1}}
\def\MCM@memosubject{\relax}
\newcommand{\memosubject}[1]{\gdef\MCM@memosubject{#1}}
\def\MCM@memodate{\relax}
\newcommand{\memodate}[1]{\gdef\MCM@memodate{#1}}
\def\MCM@memologo{\relax}
\newcommand{\memologo}[1]{\gdef\MCM@memologo{\protect #1}}
\def\@letterheadaddress{\relax}
\newcommand{\lhaddress}[1]{\gdef\@letterheadaddress{#1}}
\newenvironment{memo}[1][Memorandum]{%
\pagestyle{plain}%
\ifthenelse{\equal{\MCM@memologo}{\relax}}{%
% without logo specified.
}{%
% with logo specified
\begin{minipage}[t]{\columnwidth}%
\begin{flushright}
\vspace{-0.6in}
\MCM@memologo
\vspace{0.5in}
\par\end{flushright}%
\end{minipage}%
}
\begin{center}
\LARGE\bfseries\scshape
#1
\end{center}
\begin{description}
\ifthenelse{\equal{\MCM@memoto}{\relax}}{}{\item [{To:}] \MCM@memoto}
\ifthenelse{\equal{\MCM@memofrom}{\relax}}{}{\item [{From:}] \MCM@memofrom}
\ifthenelse{\equal{\MCM@memosubject}{\relax}}{}{\item [{Subject:}] \MCM@memosubject}
\ifthenelse{\equal{\MCM@memodate}{\relax}}{}{\item [{Date:}] \MCM@memodate}
\end{description}
\par\noindent
\rule[0.5ex]{\linewidth}{0.1pt}\par
\bigskip{}
}{%
\clearpage
\pagestyle{fancy}%
}
\newtheorem{Theorem}{Theorem}[section]
\newtheorem{Lemma}[Theorem]{Lemma}
\newtheorem{Corollary}[Theorem]{Corollary}
\newtheorem{Proposition}[Theorem]{Proposition}
\newtheorem{Definition}[Theorem]{Definition}
\newtheorem{Example}[Theorem]{Example}
\renewcommand\section{\@startsection{section}{1}{\z@}%
{-0pt\@plus -.2ex \@minus -.2ex}%
{1pt \@plus .2ex}%
{\rmfamily\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
{-0pt\@plus -.2ex \@minus -.2ex}%
{1pt \@plus .2ex}%
{\rmfamily\large\bfseries}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
{-.5ex\@plus -1ex \@minus -.2ex}%
{.25ex \@plus .2ex}%
{\rmfamily\normalsize\bfseries}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
{1ex \@plus1ex \@minus.2ex}%
{-1em}%
{\rmfamily\normalsize}}
\providecommand{\dif}{\mathop{}\!\mathrm{d}}
\providecommand{\me}{\mathrm{e}}
\providecommand{\mi}{\mathrm{i}}
\definecolor{grey}{rgb}{0.8,0.8,0.8}
\definecolor{darkgreen}{rgb}{0,0.3,0}
\definecolor{darkblue}{rgb}{0,0,0.3}
\def\lstbasicfont{\fontfamily{pcr}\selectfont\footnotesize}
\lstset{%
% numbers=left,
% numberstyle=\small,%
showstringspaces=false,
showspaces=false,%
tabsize=4,%
frame=lines,%
basicstyle={\footnotesize\lstbasicfont},%
keywordstyle=\color{darkblue}\bfseries,%
identifierstyle=,%
commentstyle=\color{darkgreen},%\itshape,%
stringstyle=\color{black}%
}
\lstloadlanguages{C,C++,Java,Matlab,Mathematica}
\endinput
%%
%% This work consists of these files mcmthesis.dtx,
%% figures/ and
%% code/,
%% and the derived files mcmthesis.cls,
%% mcmthesis-demo.tex,
%% README,
%% LICENSE,
%% mcmthesis.pdf and
%% mcmthesis-demo.pdf.
%%
%% End of file `mcmthesis.cls'.

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "68,22",
"scrollLine": "47"
}

View File

@@ -0,0 +1,4 @@
{
"source_window_id": "",
"Source": "Source"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "39,55",
"scrollLine": "33"
}

View File

@@ -0,0 +1,6 @@
{
"source_window_id": "",
"Source": "Source",
"cursorPosition": "151,104",
"scrollLine": "0"
}

View File

@@ -0,0 +1,4 @@
C%3A%2FUsers%2Fzhjx_%2FDesktop%2F22%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2022%E6%9C%80%E6%96%B0%E4%BF%AE%E6%94%B9)%2FMCM-ICM_Summary.tex="698EAECF"
C%3A%2FUsers%2Fzhjx_%2FDesktop%2F22%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2022%E6%9C%80%E6%96%B0%E4%BF%AE%E6%94%B9)%2Fmcmthesis-demo.tex="0385145E"
C%3A%2FUsers%2Fzhjx_%2FDesktop%2F22%E7%BE%8E%E8%B5%9B%2F%E7%BE%8E%E8%B5%9B%E6%A8%A1%E6%9D%BF%2FMCM_Latex%E6%A8%A1%E6%9D%BF(2022%E6%9C%80%E6%96%B0%E4%BF%AE%E6%94%B9)%2Fmcmthesis.cls="71CF8104"
C%3A%2FUsers%2Fzhjx_%2FDownloads%2Ftimetk%2FR%2Fdplyr-slidify.R="2B5CC670"

View File

@@ -0,0 +1,2 @@
C:/Users/zhang/Desktop/25美赛/MCM_Latex2025/mcmthesis.cls="AAC4A541"
c:/Users/zhang/AppData/Roaming/TinyTeX/texmf-dist/tex/latex/xkeyval/xkeyval.sty="E9F50CF8"

37
TASK3.md → latex/TASK3.md Normal file → Executable file
View File

@@ -4,18 +4,23 @@
## 共生站点选择模型的构建与求解 ## 共生站点选择模型的构建与求解
### 模型构建 ### 模型构建
分析题目可知,共生站点的的配对条件是“同日可达性”“需求可叠加性”“需求稳定性”,本文分析数据集及实际情况,对于配对条件进行如下量化。 分析题目可知,共生站点的的配对条件是“同日可达性”“需求可叠加性”“需求稳定性”,本文分析数据集及实际情况,对于配对条件进行如下量化。
1.同日可达性 1.同日可达性
观察数据集可知,站点的经纬度已经知晓,并且本文需要的是各站点之间的曼哈顿距离为小范围地域内,因此本文以经纬度转化下的曼哈顿距离作为两站点之间的可达性指标,其计算公式如下:
由数据集可知各站点经纬度,为筛选小范围地域内的站点,本文以经纬度转化的曼哈顿距离$l_{ij}$​作为两站点的可达性指标,计算公式为:
$$ $$
l_{ij}=69.0\cdot \left | lat_{i}-lat_{j} \right | +69.0\cdot \cos(\frac{lon_{i}+lon_{j}}{2})\left | lon_{i}-lon_{j} \right | l_{ij}=69.0\cdot \left | lat_{i}-lat_{j} \right | +69.0\cdot \cos(\frac{lon_{i}+lon_{j}}{2})\left | lon_{i}-lon_{j} \right |
$$ $$
其中$l_{ij}$为曼哈顿距离,$69.0$为英里转化常数 其中$69.0$为英里转化常数
2.需求可叠加性 2.需求可叠加性
题目指出卡车的服务家庭上限为250户因此对于两个站点合并为一个共生站点时其需求的叠加必须满足一定限度为此本文给出如下判断公式 题目指出卡车的服务家庭上限为250户因此对于两个站点合并为一个共生站点时其需求的叠加必须满足一定限度为此本文给出如下判断公式
$$ $$
d_{ij}=d_{i}+d_{j}\le???(给定值) d_{ij}=d_{i}+d_{j}\le???(给定值)
$$ $$
上限值略高于卡车最大限制的目的是保留在两站点合并后略超出卡车上限的配对组合,尽可能实现食物的最大程度有效服务。 上限值略高于卡车最大限制的目的是保留在两站点合并后略超出卡车上限的配对组合,尽可能实现食物的最大程度有效服务。
3.需求稳定性 3.需求稳定性
对于接近卡车最大限制的配对组合,本文进一步考虑其需求稳定性,删去波动大的配对以防止后续共生站点的分配不均。其判别公式为 对于接近卡车最大限制的配对组合,本文进一步考虑其需求稳定性,删去波动大的配对以防止后续共生站点的分配不均。其判别公式为
$$ $$
@@ -31,44 +36,54 @@ $$
## 第一站点分配模型的建立与求解 ## 第一站点分配模型的建立与求解
确定共生站点后,本文认为共生站点的内部分配即为分配有效性问题,因此本文参照问题一的有效得分指标,给出共生站点的分配有效得分公式如下 确定共生站点后,本文认为共生站点的内部分配即为分配有效性问题,因此本文参照模型一的有效得分指标,给出共生站点的分配有效得分公式如下
第一站点的实际分配货物
1、第一站点的实际分配货物
$$ $$
g_{i}=\min(d_{i},q_{i}) g_{i}=\min(d_{i},q_{i})
$$ $$
其中$q_{i}$为第一站点的假想分配货物数量 其中$q_{i}$为第一站点的假想分配货物数量
第二站点的实际分配货物为
2、第二站点的实际分配货物为
$$ $$
g_{j}=\min(d_{j},d_{0}-g_{i}) g_{j}=\min(d_{j},d_{0}-g_{i})
$$ $$
3、效果评分
$$ $$
score=E(g_{i}+g_{j})-\lambda E(unmet_{i}+unmet_{j})-\mu E(worst_{i}+worst_{j}) score=E(g_{i}+g_{j})-\lambda E(unmet_{i}+unmet_{j})-\mu E(worst_{i}+worst_{j})
$$ $$
选取分配有效得分最大的$q_{i}$作为第一站点的分配货物数量 通过最大化评分确定最优$q_{i}$
最终我们得出各共生站点的分配方案如下表 最终我们得出各共生站点的分配方案如下表
## 共生站点下的访问分配模型的建立与求解 ## 共生站点下的访问分配模型的建立与求解
对于访问次数的求解,等价于减少站点总数并且共生站点需求合并下的访问次数求解,与问题一的访问次数模型相同,本文不再赘述,对于访问时间分配的求解,本文沿用问题一的模型进行求解,由于共生站点的存在,本问的决策变量变为 对于访问次数的求解,与问题一一致,仅需适配站点总数减少、共生站点需求合并的场景,本文不再赘述.
第$i$个单独站点第$m$次运输的时间$s_{i, m}$
第$x$个共生站点第$y$次运输时间$s_{x,y}$ 对于访问时间分配的求解,本文沿用问题一的模型进行求解,由于共生站点的存在,本问的决策变量变为
第$i$个单独站点第$t$天是否被访问
- 第$i$个单独站点第$m$次运输的时间$s_{i, m}$
- 第$x$个共生站点第$y$次运输时间$s_{x,y}$
- 第$i$个单独站点第$t$天是否被访问
$$ $$
a_{i,t}=\left\{\begin{matrix} a_{i,t}=\left\{\begin{matrix}
1&t\in S_{i}\\ 1&t\in S_{i}\\
0&t\notin S_{i} 0&t\notin S_{i}
\end{matrix}\right. \end{matrix}\right.
$$ $$
第$x$个共生站点第$t$天是否被访问 - 第$x$个共生站点第$t$天是否被访问
$$ $$
a_{x,t}=\left\{\begin{matrix} a_{x,t}=\left\{\begin{matrix}
1&t\in S_{x}\\ 1&t\in S_{x}\\
0&t\notin S_{x} 0&t\notin S_{x}
\end{matrix}\right. \end{matrix}\right.
$$ $$
目标函数:
为了表示各站点访问时间的均匀分布,本文定义目标函数为所有站点实际时间间隔与理想时间间隔的差值的绝对值之和 为了表示各站点访问时间的均匀分布,本文定义目标函数为所有站点实际时间间隔与理想时间间隔的差值的绝对值之和
$$ $$
\min Z = \sum_{i} \sum_{m=1}^{k_{i}} \left| (s_{i, m+1} - s_{i, m}) - T_{i} \right|+ \sum_{x} \sum_{y=1}^{k_{x}} \left| (s_{i, y+1} - s_{i, y}) - T_{i} \right| \min Z = \sum_{i} \sum_{m=1}^{k_{i}} \left| (s_{i, m+1} - s_{i, m}) - T_{i} \right|+ \sum_{x} \sum_{y=1}^{k_{x}} \left| (s_{i, y+1} - s_{i, y}) - T_{i} \right|
$$ $$
约束条件为 约束条件为
1.每日的访问总次数不得超过最大运输能力 1.每日的访问总次数不得超过最大运输能力
$$ $$
\sum_{i}a_{i,t}+\sum_{x}a_{x,t}\le2 \sum_{i}a_{i,t}+\sum_{x}a_{x,t}\le2

15
latex/code/mcmthesis-matlab1.m Executable file
View File

@@ -0,0 +1,15 @@
function [t,seat,aisle]=OI6Sim(n,target,seated)
pab=rand(1,n);
for i=1:n
if pab(i)<0.4
aisleTime(i)=0;
else
aisleTime(i)=trirnd(3.2,7.1,38.7);
end
end

41
latex/code/mcmthesis-sudoku.cpp Executable file
View File

@@ -0,0 +1,41 @@
//============================================================================
// Name : Sudoku.cpp
// Author : wzlf11
// Version : a.0
// Copyright : Your copyright notice
// Description : Sudoku in C++.
//============================================================================
#include <iostream>
#include <cstdlib>
#include <ctime>
using namespace std;
int table[9][9];
int main() {
for(int i = 0; i < 9; i++){
table[0][i] = i + 1;
}
srand((unsigned int)time(NULL));
shuffle((int *)&table[0], 9);
while(!put_line(1))
{
shuffle((int *)&table[0], 9);
}
for(int x = 0; x < 9; x++){
for(int y = 0; y < 9; y++){
cout << table[x][y] << " ";
}
cout << endl;
}
return 0;
}

Binary file not shown.

After

Width:  |  Height:  |  Size: 361 KiB

BIN
latex/figures/distance.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 592 KiB

BIN
latex/figures/fig1.jpg Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 63 KiB

BIN
latex/figures/flow.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 370 KiB

BIN
latex/figures/flowchartmodel1.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 2.2 MiB

BIN
latex/figures/flowchartmodel2.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 442 KiB

BIN
latex/figures/gongshi.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 372 KiB

BIN
latex/figures/model1.1.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 1.7 MiB

BIN
latex/figures/model1.2.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 1.9 MiB

BIN
latex/figures/model1.3.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 82 KiB

BIN
latex/figures/model1.4.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 110 KiB

BIN
latex/figures/model1.5.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 98 KiB

BIN
latex/figures/model1overview.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 1.4 MiB

BIN
latex/figures/model2.1.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 110 KiB

BIN
latex/figures/model2.2.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 1.6 MiB

BIN
latex/figures/model2.3.png Executable file

Binary file not shown.

After

Width:  |  Height:  |  Size: 112 KiB

Some files were not shown because too many files have changed in this diff Show More