# 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

Figure 3 indicates that if two kinematic boundary conditions are
applied to N through two symmetry conditions with respect to planes (`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. Note that if one of the two symmetry
conditions is a type "Tied" condition, the velocity of 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

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

`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 be (progressively) not 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:

$d{\left({N}_{i},\left(P\right)\right)}^{2}\le \left(1+{\alpha}_{sort}\right)\cdot \mathrm{max}{\left({d}_{i}+{d}_{j}\right)}^{2}$

With ${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 gives a way to over-estimate the number of particles which will be symmetrized at one time.

`N`

_{i}has to be symmetrized with respect to n conditions, the particle

`N`

_{i}gives birth to 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$

Anyway, 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)

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

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

- $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.