Homec4science
Diffusion lammps (master) Browse

lammps/src/SHOCK/fix_msst.h
bc446bb8b01bmaster

Last Change
Raw File

fix_msst.h
View Options

/* -*- c++ -*- ----------------------------------------------------------
LAMMPS - Large-scale Atomic/Molecular Massively Parallel Simulator
http://lammps.sandia.gov, Sandia National Laboratories
Steve Plimpton, sjplimp@sandia.gov
Copyright (2003) Sandia Corporation. Under the terms of Contract
DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government retains
certain rights in this software. This software is distributed under
the GNU General Public License.
See the README file in the top-level LAMMPS directory.
------------------------------------------------------------------------- */
#ifdef FIX_CLASS
FixStyle(msst,FixMSST)
#else
#ifndef FIX_MSST_H
#define FIX_MSST_H
#include "fix.h"
namespace LAMMPS_NS {
class FixMSST : public Fix {
public:
FixMSST(class LAMMPS *, int, char **);
~FixMSST();
int setmask();
void init();
void setup(int);
void initial_integrate(int);
void final_integrate();
double compute_scalar();
double compute_vector(int);
void write_restart(FILE *);
void restart(char *);
int modify_param(int, char **);
double memory_usage();
private:
double dtv,dtf,dthalf; // Full and half step sizes
double boltz,nktv2p, mvv2e; // Boltzmann factor and unit conversions
double total_mass; // Mass of the computational cell
double omega[3]; // Time derivative of the volume
double p_current[3],dilation[3];
double qmass; // Effective cell mass
double mu; // Effective cell viscosity
double tscale; // Converts thermal energy to compressive
// strain ke at simulation start
int dftb; // flag for use with DFTB+
double velocity_sum; // Sum of the velocities squared
double damping; // Damping function for TS force term at
// small volume difference (v0 - vol)
double T0S0; // Initial TS term for DFTB+ simulations
double S_elec,S_elec_1,S_elec_2; // time history of electron entropy
// for DFTB+ simulaitons
double TS_dot; // time derivative of TS term for
// DFTB+ simulations
double **old_velocity; // Saved velocities
int kspace_flag; // 1 if KSpace invoked, 0 if not
int nrigid; // number of rigid fixes
int *rfix; // indices of rigid fixes
char *id_temp,*id_press; // Strings with identifiers of
char *id_pe; // created computes
class Compute *temperature; // Computes created to evaluate
class Compute *pressure; // thermodynamic quantities
class Compute *pe;
int tflag,pflag,vsflag,peflag; // Flags to keep track of computes that
// were created
// shock initial conditions
double e0; // Initial energy
double v0; // Initial volume
double p0; // Initial pressure
double velocity; // Velocity of the shock
double lagrangian_position; // Lagrangian location of computational cell
int direction; // Direction of shock
int p0_set; // Is pressure set
int v0_set; // Is volume set
int e0_set; // Is energy set
double TS_int; // Needed for conserved quantity
// with thermal electronic excitations
double beta; // Energy conservation scaling factor
int maxold; // allocated size of old_velocity
class FixExternal *fix_external; // ptr to fix external
// functions
void couple();
void remap(int);
double compute_etotal();
double compute_vol();
double compute_hugoniot();
double compute_rayleigh();
double compute_lagrangian_speed();
double compute_lagrangian_position();
double compute_vsum();
};
}
#endif
#endif
/* ERROR/WARNING messages:
E: Illegal ... command
Self-explanatory. Check the input script syntax and compare to the
documentation for the command. You can use -echo screen as a
command-line option when running LAMMPS to see the offending line.
E: Fix msst tscale must satisfy 0 <= tscale < 1
Self-explanatory.
E: Fix msst requires a periodic box
Self-explanatory.
E: Cannot use fix msst without per-type mass defined
Self-explanatory.
E: Could not find fix msst compute ID
Self-explanatory.
E: Fix msst compute ID does not compute temperature
Self-explanatory.
E: Fix msst compute ID does not compute pressure
Self-explanatory.
E: Fix msst compute ID does not compute potential energy
Self-explanatory.
E: Could not find fix_modify temperature ID
The compute ID for computing temperature does not exist.
E: Fix_modify temperature ID does not compute temperature
The compute ID assigned to the fix must compute temperature.
W: Temperature for MSST is not for group all
User-assigned temperature to MSST fix does not compute temperature for
all atoms. Since MSST computes a global pressure, the kinetic energy
contribution from the temperature is assumed to also be for all atoms.
Thus the pressure used by MSST could be inaccurate.
E: Could not find fix_modify pressure ID
The compute ID for computing pressure does not exist.
E: Fix_modify pressure ID does not compute pressure
The compute ID assigned to the fix must compute pressure.
*/
c4science · Help