LIGGGHTS® - Developer Forum

Topics related to developing with LIGGGHTS® can be discussed here: discussion about implementation details, C++, MPI and debugging tools

Access to all particle coordinates at a specific timestep

Submitted by aishk on Fri, 10/11/2019 - 23:10

Hello,

Is there a way to access the coordinates of all particles in the simulation at a particular timestep in the code? Also, is this done by the neighbor list or the atomVec class? Are they all stored in some array that I can access easily within the source-code?

I tried looking into how dump_custom_vtk works to display the particle coordinates to the .vtk file but I did not find anything.

If this coordinate storage is done by the neighbor list, is there some calculation needed to get global coordinates instead of local coordinates?

Thanks!

easiest way to access young's modulus values

Submitted by mschramm on Thu, 10/03/2019 - 21:57

Hello.
I would like to use the young's modulus that is imported by the user when describing a particle.
What would be the easiest way to access this value?
From looking at fix_check_timestep, I know I want to use
Y = static_cast(modify->find_fix_property("youngsModulus","property/global","peratomtype",max_type,0,style));
but I do not know what style is in this sense.

Invoke the voronoi cell volume using as parameters in other compute

Submitted by Jeff_DONG on Wed, 09/18/2019 - 03:55

Hi,

I am trying to invoke the Voronoi Cell volume from the "compute_voronoi"(voro++ works well) to use as a parameter in the calculation of "easo_capillary.h".

I just found a similar example when doing the "compute_pressure" which also needs to invoke the value of "compute_temperature" for calculation.

My question is that I can only invoke it this way(compute_pressure and compute_temperature)? or does anyone have other suggestions on the invoking?

Very appreciate!

Best Regards,

Jeffrey

RDEM: Radiative heat transfer model in Discrete Element Method

Submitted by wuhaoinet on Mon, 09/09/2019 - 08:43

I have published several paper on radiative heat transfer in packed bed. Thus, I plan to open the code in the numerical simulation.
The numerical code includes radiation with short-range cutoff and the conduction ...
If anyone interested this topic, please discuss on https://github.com/wuhao1938/RDEM

Later, the related code will be in https://github.com/wuhao1938/RDEM

Modification of source code not working

Submitted by aishk on Wed, 08/28/2019 - 20:28

Hi,
I need to extend the source code on LIGGGHTS so I decided to start with looking at adding a few fprintf statements for debugging but any fprintf statements I add to the source code doesn't seem to be showing up on the terminal screen or on the logfile.
This is how I am remaking LIGGGHTS after any modifications:
make clean-auto
make stubs
make auto

Is there something I am doing wrong or not doing?
Any help is appreciated!
Thanks

Multicontact Model Issues

Submitted by amgroen on Wed, 07/31/2019 - 18:48

Is the multicontact model supposed to work with MPI?
I have not had any luck getting the example program in.hydrogel_multicontact to run on Windows.

Due to this error:
Fatal error in MPI_Wait: Message truncated, error stack:
MPI_Wait(request=0x000000EDE92FF220, status0x000000EDE92FF228) failed
Message from rank 1 and tag 0 truncated; 209176 bytes received but buffer size is 117120

CHANGING COHESION_MODEL_SJKR EQUATION

Submitted by amohammad42 on Wed, 07/17/2019 - 18:26

Hello Everyone,

I wanted to change the equation of cohesion model in liggghts, so I changed the equation of Force = Cohesive Energy Density * Area in cohesion_model_sjkr.h (lines 98 & 99). Do I need to change any other file associated with it before compiling ? What file in addition to cohesion_model_sjkr.h I need to change in order to implement this new equation of cohesion?

I am new to liggghts and C++, so please let me know what should I do.

Thanks&Regards,
Abdul.

Pages

Subscribe to RSS - LIGGGHTS® - Developer Forum