# SPH Symmetry Conditions

## Multiple SPH Symmetry Conditions

An axi-symmetry condition can be modelized through the use of two conditions with respect to two planes intersecting at the axis of symmetry. A spheric symmetry condition can be modelized through the use of three conditions with respect to three planes intersecting at the center of symmetry.

Nevertheless, these kinds of symmetries are not treated the same way.

Therefore, some characteristics of axi-symmetry (respectively, spheric symmetry) conditions can be closed to the axis of symmetry (respectively, the center of symmetry).

## Kinematic Boundary Condition

With respect to the previous discussion: adding the kinematic boundary condition an explicit way allows to enforce it.

- If "Slide" type, the velocity of the node in direction
"
`Dir`" is set to zero - If "Tied" type, the velocity of the node in all directions is set to zero

`P`

_{1}) and (

`P`

_{2}), the two boundary conditions are modified so that the velocity in the plane normal to common axis of (

`P`

_{1}) and (

`P`

_{2}) will remain zero.

`N`in all directions is set to zero.

It also allows application to the same node, a kinematic boundary condition through a SPH symmetry condition (/SPHBCS) and a standard boundary condition (/BCS) at the same time, as long as the standard boundary condition is not given in a moving skew system, but a fix skew system or the global skew system. The two conditions are then composed the same way.

## Part Mass

You must be advised that when a particle lies on a symmetry plane at time t=0, the mass and the initial volume considered for the particles are respectively:

Where, `m`_{p} is the mass specified into property
set.

When a particle lies on $n$ symmetry planes at time t=0,

Ghost particles built from this particle will get the same initial volume and mass.

When $n>2$ , the previous equation may provide an error on mass and energies output for the part the particles belong to, with respect to the physical model.

## Formulation Level

When a symmetry plane is defined, and even if a kinematic condition is set for all particles lying on the symmetry plane, particles lying at time zero inside the domain are theoretically able to cross the symmetry plane.

If `Ilev`=0, particles crossing the symmetry plane
will not be (progressively) taken into account anymore in the computation, neither
than their symmetric particles which then lie inside the domain.

`Ilev`=1, particles which have crossed the symmetry plane rebound an elastic way upon the symmetry plane: their velocity in the normal direction to the plane is set the opposite.

`Ilev`=1, it is strongly recommended to associate kinematic condition to all particles lying on symmetry plane at time zero, for computational time reasons.

## Maximum Created Number of Ghost Particles

Ghost particles are created at each search for neighbors time within the security distance, and then destroyed when a new search occurs (a new set of ghost particles is then created).

At any search time, all ghost particles which are inside the security distance of any real particle are created.

In practice, some more particles, strictly necessary, are created: a symmetric
particle `G`_{i} to particle `N`_{i} is created, with respect to symmetry plane
$P$
, if
$\exists j$
neighbor of
$i$
:

Where, ${d}_{i}$ and ${d}_{j}$ are the smoothing lengths related to particle $i$ and $j$ .

As long as no real particles cross the symmetry plane (all real particles lie on the same side of the symmetry plane), this criteria is sufficient to get all ghost neighbors of all real particles inside the security distance, since:

$d\left({N}_{i},\left(P\right)\right)\le d\left({G}_{i},{N}_{i}\right)$ for $\forall \left(i,j\right)$

And, $d\left({G}_{i},{N}_{i}\right)\le \sqrt{1+{\alpha}_{sort}}\cdot \mathrm{max}\left({d}_{i}+{d}_{j}\right)$

Particles, which one can expect to remain far from the symmetry plane all along the simulation, will never be symmetrized. This provides a way to over-estimate the number of particles which will be symmetrized at one time.

When a particle `N`_{i} has to be symmetrized with respect to
n conditions, the particle `N`_{i} creates n ghost particles. The following quantity must
remain less than `Maxsph` (since v14.0.220,
`Maxsph` is ignored and the memory is dynamically
allocated).

- $n$
- Number of conditions
- ${n}_{particles}$
- Number of particles to be symmetrized with respect to condition $n$

The default value which is the number of SPH symmetry conditions multiplied by the number of particles will be enough to treat any problem.

## Solid to SPH Options (Sol2SPH)

The solid to SPH option (Sol2SPH) enables you to turn a solid element into particles either in order to increase the time step/robustness of a Lagrangian calculation, while not changing the physics.

## Time Step

- Particle time step (/DT/SPHCEL)
- Nodal time step (/DT/NODA)

In particle time step, stable time step is computed as:

- ${d}_{i}$
- Smoothing length related to particle $i$
- ${c}_{i}$
- Sound speed at location $i$

It is recommended to set time step scale factor $\text{\Delta}{T}_{sca}$ to 0.3.

In nodal time step, stable time step is computed as:

- $m$
- Mass for particles
- ${K}^{*}$
- Stiffness based on SPH interaction

For time step scale factor $\text{\Delta}{T}_{sca}$ , it is recommended to set it to 0.67.

## Thermal Analysis

Heat transfer is now available between SPH particles and finite elements with `I`_{the}=1 in
/INTER/TYPE7 and /INTER/TYPE21; and with
/THERM_STRESS/MAT, thermal expansion in SPH is also
possible.