2D Collimator Analysis. Indirect coupled field analysis. Technique. Previous analysis has involved direct coupled fields ANSYS solves the thermal and structural components for each time step as 1 analysis.
Indirect coupled field analysis
Luis gave me several files giving heat deposition values for beam strikes at positions 2mm and 10 mm from the top of a graphite-titanium collimator for 250 GeV and 500 GeV beams.
I rearranged the data into a form acceptable for ANSYS and used data from the Y-Z plane, at x=0.
This should give the worst case scenario being directly in the path of the beam.
Initially I used data from a 250GeV beam 2mm from the top.
Approximate line of 250GeV beam
This is a matlab produced figure showing the heat deposition.
There is a lot of distortion in the figure. The aspect ratio is 300:1.
Unfortunately I couldn’t create a key explaining what heat load the temperatures relate to, other than red is the maximum and blue is the minimum.
After 1ps of heating ANSYS predicted the maximum temperature to be 283.43°C
The distortion in Y-Z is slightly different to that shown on slide 4.
The figure also appears to be smoother. This is because ANSYS has interpolated values between the nodes.
The model was fixed in the z direction.
Theoretically with a TEC of 9.2x10-6m°-1 and ΔT= 283.43°C there should be an expansion of .0026m per metre if the model was free to move.
But as the material is constrained by boundary condition this leads to a strain of -0.26%.
The bottom edge of the area was fixed so it couldn’t move in either direction, set to plane stress so the model was free to move into and out of plane
By multiplying the strain by the Young’s Modulus the stress can be calculated. This works out to be -295.88MPa. This agrees with the result from ANSYS.
I took the temperature figures and applied them to a structural model in two stages.
The first stage ramped linearly from 0°C to the maximum temperature in 1ps with 10 sub steps.
The second stage applied the maximum temperature constantly for 10ns.
A subsequent test showed that the maximum temperature would decrease by 0.045°C after 10ns.
The model is constrained in z by the boundary condition, but it is free to move in x because its set to plane stress.
The inertia of the unheated area causes the model to effectively be constrained in y.
The stress in y and z is predicted by this equation:
This gives a value of -423MPa. This agrees with the ANSYS result.
A mesh size of 50microns was used. This is 5 times larger than the ideal size, but ANSYS wouldn’t allow a 10micron mesh.
This caused some problems, one of which being a non-uniform distribution of nodes along the beam path. ANSYS interpolates the heat load to the nearest node which led to a strange heating pattern through the material. This is an artefact of the meshing.
Graph of stress along reflected wave
The peaks are at the interaction point