Osney Thermofluids Laboratory Department of Engineering Science Attitude-Dependent Thermal Flux Estimation in the Slip Regime for Geometric Primitives Towards the Improvement of Destructive Re-entry Simulations Nathan Donaldson, MEng DPhil student, University of Oxford
22
Embed
Osney Thermofluids Laboratory Department of …old.esaconferencebureau.com/Custom/15A01/Presentations/Room 1.1... · Osney Thermofluids Laboratory Department of Engineering Science
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Osney Thermofluids Laboratory Department of Engineering Science
Attitude-Dependent Thermal Flux Estimation in the Slip Regime for Geometric Primitives Towards the Improvement of Destructive Re-entry Simulations Nathan Donaldson, MEng DPhil student, University of Oxford
Osney Thermofluids Laboratory Department of Engineering Science
• Rapid debris demise simulation codes make use of simplified aerothermal models
• Heat flux assessments are often tumble-averaged
• Inaccuracies may manifest in demise estimation as a result
• Also, high fidelity, full trajectory simulations are extremely computationally intensive
Motivation
Osney Thermofluids Laboratory Department of Engineering Science
• Heat flux correlations for common primitives are sometimes used to assess survivability of simplified components during re-entry
• However, a set of correlations for edge heating during stable-attitude descent does not currently exist — random tumbling is favoured
• The application of such correlations would allow attitude dependent temperature history calculations for worst-case heating/stabilisation problems to be performed extremely inexpensively
Motivation
Osney Thermofluids Laboratory Department of Engineering Science
• Rapid assessment codes could hence be updated to produce high resolution results without an appreciable change in simulation duration
• The implementation of dynamic tumbling models could consequently be simplified, such that object stabilisation and tumbling rate could be directly assessed
• A comprehensive series of high fidelity analyses was seen as the first step to developing these correlations.
Motivation
Osney Thermofluids Laboratory Department of Engineering Science
•Geometry definition, meshing and DSMC
•Post-processing and data extraction
•Curve fitting and figure generation
Numerical analysis framework
Osney Thermofluids Laboratory Department of Engineering Science
DSMC
Osney Thermofluids Laboratory Department of Engineering Science
DSMC
Osney Thermofluids Laboratory Department of Engineering Science
DSMC data extraction
Pictured are the data extraction line locations for cuboid and cylinder cases. Side elements are pictured in red, front elements in green.
Osney Thermofluids Laboratory Department of Engineering Science
• Heat flux data were extracted from cells at the four locations shown in the previous slide
• A cubic polynomial surface was fit to the data points using pyeq2, and
coefficients computed • Hence, a relation of the form shown below was found between the
pitch, yaw and incident heat flux of each point analysed
Osney Thermofluids Laboratory Department of Engineering Science
Cylinder side face correlation
Pictured are the edge data and fitted surfaces for cylinder cases.
Osney Thermofluids Laboratory Department of Engineering Science
Error relative to local values
Pictured are the relative errors for cylinder correlations. Front and side element correlations are shown from left to right, respectively. The mean
errors for the cuboid and two cylinder correlations were 9.35%, 11.32% and 852.59% respectively (due to small values of normalised heat flux).
Osney Thermofluids Laboratory Department of Engineering Science
Error relative to local maxima
Pictured are the relative errors for cylinder correlations. Front and side element correlations are shown from left to right, respectively. When
normalised to local maxima, mean errors for the cuboid and two cylinder correlations are reduced to 2.08%, 8.58% and 10.80% respectively.
Osney Thermofluids Laboratory Department of Engineering Science
• The correlations may be used trivially as hand calculations in order to make quick, preliminary assessments during a design process
• They are also directly applicable as validation data for aerothermal models
Applying the correlations
Osney Thermofluids Laboratory Department of Engineering Science
• Object-oriented demise codes like DAS, SESAM, and ORSAT can benefit directly from these correlations
• Trends may be applied to primitives, and dimensionalised by multiplying through by the reference heat flux (stagnation at zero incidence)
Applying the correlations
Osney Thermofluids Laboratory Department of Engineering Science
• Spacecraft-oriented codes such as SCARAB can also benefit
• Computationally inexpensive curvature and shape assessment routines may be utilised to evaluate flat faces, rotational symmetry, and sharp corners
• Correlations may then be applied automatically, and dimensionalised by multiplying through by the reference heat flux (stagnation at zero incidence)
Applying the correlations
Osney Thermofluids Laboratory Department of Engineering Science
• The reference heat flux may be calculated from freestream conditions using stagnation heat flux correlation such as Detra-Hidalgo, Tauber, or Brandis-Johnston (or indeed the Fay-Riddell expression, which is derived from first principles)
• High resolution heat flux distributions are determined with minimal computational expense – ultimately, only algebraic expressions are employed
Applying the correlations
Osney Thermofluids Laboratory Department of Engineering Science
• Validation of correlations at different Knudsen numbers
• Modification of correlation equations to reduce peak errors
• Experimental validation using Oxford’s low density wind tunnel (capable of reproducing slip regime flow, Kn ≈ 0.03)
Further work
Osney Thermofluids Laboratory Department of Engineering Science