RADIAL_DEFORMATION

Solution control and techniques

*RADIAL_DEFORMATION
"Optional title"
coid
entype, enid, $t_{beg}$, $t_{end}$

Parameter definition

Variable
Description
coid
Command ID
entype
Entity type
options: P, PS
enid
Entity identification number
$t_{beg}$
Start time
quantity: Time
default: 0
$t_{end}$
End time
quantity: Time
default: 1.0e20

Description

With this command active, the outer surface of an object is only allowed to undergo a uniform radial deformation. It is typically used to stabilize the initial spherical expansion of explosive charges in UNDEX simulations.

Example

Enforcing a radial expansion of an UNDEX charge

In UNDEX simulations, it is recommended to model the detonating explosive with Finite Elements rather than with CFD. The early interaction between explosive and and water (Surface SPH) is more accurate with a Lagrangian FE-SPH contact than with the CFD-SPH coupling.

However, small Finite Elements might distort when interacting with large and heavy SPH nodes. By enforcing a homogeneous radial deformation of the expanding detonation products, these effects can be avoided.

Beyond the initial rapid radial expansion phase, the bubble dynamics is better handled with a CFD approach. In this example, the explosive elements are deactivated after $1ms$ and their state is mapped to the CFD domain (E2C).

*UNIT_SYSTEM SI *PARAMETER %R0 = 0.1, "TNT charge radius" %R1 = 3.0, "Water domain radius" %L = 1.0, "Slab size" %H = 0.2, "Slab thickness" %X = 1.0, "Stand-off distance" %dx = 0.04, "SPH node spacing" %dc = 0.04, "CFD cell size" # # --- TIME --- # *TIME 0.002 # # --- FE MESH --- # *COMPONENT_BOX "Slab" 1, 1, 10, 10, 1 [-%L/2], [-%L/2], [-%X-%H], [%L/2], [%L/2], [-%X] *COMPONENT_SPHERE "TNT charge" 2, 2, 4 0, 0, 0, [%R0] *CHANGE_P-ORDER P, 2, 2 *SMOOTH_MESH P, 2, 60.0 # # --- MATERIAL --- # *MAT_RIGID "Slab" 1, 2500.0 *MAT_OBJECT "TNT" "{c8da5215-7469-4a4e-a3f2-63528907eff7}", 2 *MAT_FLUID "Water" 3, 1000.0, 2.25e9, 0, -1.0e5, 3 *EOS_TAIT 3, 7.0 # # --- PART --- # *PART "Plate" 1, 1 "TNT" 2, 2 # # --- DETONATION OF FE CHARGE --- # *DETONATION 1 P, 2, 0, 0, 0, 0 # # --- RADIAL DEFORMATION OF FE CHARGE --- # *RADIAL_DEFORMATION 1 P, 2 # # --- DEACTIVATE FE CHARGE --- # *ACTIVATE_ELEMENTS 1, P, 2, 0, 1.0e-3 # # --- SPH --- # *PARTICLE_DOMAIN ALL, 0 [-1.5*%R1], [-1.5*%R1], [-1.5*%R1], [1.5*%R1], [1.5*%R1], [1.5*%R1] *GENERATE_PARTICLE_DISTRIBUTION 3 10, 0, [%dx] *GEOMETRY_COMPOSITE 10 11, -12 *GEOMETRY_SPHERE 11 0, 0, 0, [%R0], [%R1] *GEOMETRY_BOX 12 [-%L/2], [-%L/2], [-%X-%H], [%L/2], [%L/2], [-%X] *PARTICLE_SPH "Water" 3 3, 0, 2 # # --- MERGE --- # *MERGE "SPH to slab" SPH, 3, P, 1 "SPH to charge" SPH, 3, P, 2 # # --- CFD --- # *CFD_DOMAIN 1 ALL, 0, [%dc], 0, 1 [-6*%R0], [-6*%R0], [-6*%R0], [6*%R0], [6*%R0], [6*%R0] *CFD_HE "TNT" 2 2, 2, 1 *GEOMETRY_PART "TNT (FE)" 2 2 *END