[FLASH-USERS] Continuing to evolve particles that leave the domain

Anshu adubey at lbl.gov
Sat Mar 28 12:38:27 EDT 2015


Try adding another category such as OUTSIDE, with that particle should not get deleted.
Give it a number > maxblocks so that sort puts those particles at the end. When you pass the data structure to particle routines make sure the local count is such that these particles are not
processed. And then add a routine to treat them separately. It has been a couple of years since
I worked with this code, so I might have a few things wrong and a few details missing, but I think 
in general this should work. 

Best,
Anshu

> On Mar 28, 2015, at 9:16 AM, James Guillochon <jfg at ucolick.org> wrote:
> 
> Hi all,
> 
> I'm having a bit of trouble getting this to work. I need to leave a particle assigned to a block at all times, and not marked as NONEXISTENT or LOST, resulting in its deletion. Unfortunately it seems like there's a number of places where particles can be mark as NONEXISTENT in the code, and changing them all seems difficult (there's a only a few places where the particle is marked as "LOST", which I've disabled).
> 
> Anyone have any additional insight here? At the moment, particles are deleted as soon as they leave the domain, even with the changes I've made.
> 
> Thanks,
> - James
> 
>> On Sat, Mar 28, 2015 at 12:15 PM, James Guillochon <guillochon at gmail.com> wrote:
>> Hi all,
>> 
>> I'm having a bit of trouble getting this to work. I need to leave a particle assigned to a block at all times, and not marked as NONEXISTENT or LOST, resulting in its deletion. Unfortunately it seems like there's a number of places where particles can be mark as NONEXISTENT in the code, and changing them all seems difficult (there's a only a few places where the particle is marked as "LOST", which I've disabled).
>> 
>> Anyone have any additional insight here? At the moment, particles are deleted as soon as they leave the domain, even with the changes I've made.
>> 
>> Thanks,
>> - James
>> 
>>> On Fri, Mar 27, 2015 at 4:58 PM, Anshu Dubey <adubey at lbl.gov> wrote:
>>> If you don't need them to interact with the mesh in anyway I don't see
>>> any conceptual problem with advancing them. You will still have to make 
>>> sure that they are not treated the same as particles within domain
>>> in any function. 
>>> 
>>> Anshu
>>> 
>>>> On Fri, Mar 27, 2015 at 1:48 PM, James Guillochon <jfg at ucolick.org> wrote:
>>>> Hi all, I'm using active particles in my calculation, and I would like to continue to evolve the particle trajectories even once the particles have left the domain. I see that particles that leave the domain can be saved after they leave by enabling a flag, but it appears that their evolution at this point is frozen. It looks like a few relatively simple changes in Particles_advance may be able to continue their evolution, but I wasn't sure if evolving these particles would cause issues.
>>>> 
>>>> Has anyone tried to do something like this before?
>>>> 
>>>> Cheers,
>>>> - James
>>>> 
>>>> -- 
>>>> James Guillochon
>>>> Einstein Fellow at the Harvard-Smithsonian CfA
>>>> jguillochon at cfa.harvard.edu
>>> 
>>> 
>>> 
>>> -- 
>>> Anshu Dubey <adubey at lbl.gov>
>>> Mailstop: 50A1148
>>> Building-Room 050A-2154 
>>> 510-486-5242
>> 
>> 
>> 
>> -- 
>> James Guillochon
>> Einstein Fellow at the Harvard-Smithsonian CfA
>> jguillochon at cfa.harvard.edu
> 
> 
> 
> -- 
> James Guillochon
> Einstein Fellow at the Harvard-Smithsonian CfA
> jguillochon at cfa.harvard.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20150328/eb636f16/attachment.htm>


More information about the flash-users mailing list