<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=big5">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
Dear Erik,
<div>Thanks for your reply.It is conceptually helpful.</div>
<div>I think I would use other codes to deal with the exterior after the in<span style="font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols; font-size:16px">terior
is simulated by ET like you said.</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols; font-size:16px">Best regards,</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols; font-size:16px">Chia-Hui</span></div>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>寄件者:</b> Erik Schnetter <schnetter@protonmail.com><br>
<b>寄件日期:</b> 2018年10月11日 下午 09:14:49<br>
<b>收件者:</b> 林家暉<br>
<b>副本:</b> Einstein Toolkit Users<br>
<b>主旨:</b> Re: [Users] separation of simulation domain</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Chia-Hui<br>
<br>
I think you are asking a question about physics, not about the Einstein Toolkit.<br>
<br>
To simulate a system, you need to have a mathematical description to start with. If you want to combine a general relativistic interior with a Newtonian (?) exterior, then you first need to come up with a physically consistent description, for example a consistent
set of PDEs with respective interface conditions. I believe this is very difficult. At the least this is a mathematical question, not one about software engineering.<br>
<br>
What we have been doing in the past is to simple use general relativity everywhere, even in the exterior. If you use adaptive mesh refinement (AMR), then you can have a central region where small length and time scales are resolve, and an exterior region that
is much larger and where only coarser length and time scales are resolved. As Roland mentions, using approximately spherical coordinates in the exterior via Llama allows you to have a constant angular resolution there, which is convenient for tracking gravitational
waves generated by a point source.<br>
<br>
If you need to resolve hours and days, then in general the following approximations hold:<br>
- from the point of view of the interior, the exterior is constant in time<br>
- from the point of view of the exterior, only the temporally averaged interior is relevant<br>
<br>
People would thus first perform a calculation of the interior, and once this has reached a "steady state", use this to motivate certain inner boundary conditions for a second simulation that treats only the exterior.<br>
<br>
-erik<br>
<br>
> On Oct 11, 2018, at 0:29 , 林家暉 <r06222015@ntu.edu.tw> wrote:<br>
> <br>
> Hello,<br>
> I have a question about whether it is possible to separate the domain of the simulation. I want to simulate with a larger scale such as (1 million km )^3 , however it is too computation-consuming to simulate the whole domain with GR and it is not necessary
to consider GR because GR could be ignored in such large scale. But in the region near center with high density , GR is important .
<br>
> Therefore I am wondering if it is possible that GR is included in the center region ,say (1000 km)^3, but not included in the other region. However another problem is that timescale is also very different. ET's time scale is about ~ ms but what I need is
around hours to days. Is there some suggestion for me to deal with it ?<br>
> Thanks,<br>
> Best regards,<br>
> Chia-Hui <br>
<br>
-- <br>
Erik Schnetter <schnetter@protonmail.com><br>
<a href="http://www.perimeterinstitute.ca/personal/eschnetter/">http://www.perimeterinstitute.ca/personal/eschnetter/</a><br>
<br>
<br>
<br>
</div>
</span></font>
</body>
</html>