P.P. Sundaramoorthy
Please Note
20 records found
1
Coupled DSMC Simulation of an Air-Breathing Electric Propulsion Intake and Thruster Interface
Effects of thruster interface, and Chamber Shape on Neutral Gas Delivery
The framework is first validated against the non-gridded conical intake and applied to four passive baseline configurations combining two intake profiles with two discharge-chamber geometries. The intake profile produces sub-percent differences at the throat and discharge-exit planes, as expected because diffuse re-emission removes the directional memory the profile would otherwise carry. The chamber geometry produces a measurable effect on the throat measurement, indicating that the throat responds to the downstream boundary at the level of a passive chamber wall alone.
The framework is then applied to the SITAEL split-ring architecture, scaled to a 16U CubeSat platform and coupled to the RIT-µ3 ion optics at 215 km. The interface is modelled as an absorbing-diffuse boundary at the neutral transparency τn = 0.25, with the central blank resolved as a physical surface. Resolving the interface changes the reported performance. The chamber density is 1.67× the fully absorbing-boundary prediction, the delivered flux is 0.43× it, and the residence time is 4.05× it. These corrections exceed the differences between the intake geometries reported in the literature.
The study then examines three geometric parameters at the exit face against the resolved interface: the blank surface profile, the intake capture area, and the discharge-chamber cross-section. The blank surface profile redirects the returned population within the chamber without altering the flux delivered to the thruster. The intake capture area enters the reported collection efficiency through its normalization alone and leaves the delivered flux unchanged at fixed exit geometry. Only the discharge-chamber cross-section alters the delivered flux through the exit area it presents to the thruster. The beam extraction requirement of the ion optics fixes both parameters that govern that flux. ...
The framework is first validated against the non-gridded conical intake and applied to four passive baseline configurations combining two intake profiles with two discharge-chamber geometries. The intake profile produces sub-percent differences at the throat and discharge-exit planes, as expected because diffuse re-emission removes the directional memory the profile would otherwise carry. The chamber geometry produces a measurable effect on the throat measurement, indicating that the throat responds to the downstream boundary at the level of a passive chamber wall alone.
The framework is then applied to the SITAEL split-ring architecture, scaled to a 16U CubeSat platform and coupled to the RIT-µ3 ion optics at 215 km. The interface is modelled as an absorbing-diffuse boundary at the neutral transparency τn = 0.25, with the central blank resolved as a physical surface. Resolving the interface changes the reported performance. The chamber density is 1.67× the fully absorbing-boundary prediction, the delivered flux is 0.43× it, and the residence time is 4.05× it. These corrections exceed the differences between the intake geometries reported in the literature.
The study then examines three geometric parameters at the exit face against the resolved interface: the blank surface profile, the intake capture area, and the discharge-chamber cross-section. The blank surface profile redirects the returned population within the chamber without altering the flux delivered to the thruster. The intake capture area enters the reported collection efficiency through its normalization alone and leaves the delivered flux unchanged at fixed exit geometry. Only the discharge-chamber cross-section alters the delivered flux through the exit area it presents to the thruster. The beam extraction requirement of the ion optics fixes both parameters that govern that flux.
A building-block approach was adopted. At the coupon level, laminate quality assessments and mechanical tests were conducted on specimens extracted from the curved structure to establish baseline material quality and determine its in-plane material properties. These results showed good agreement with literature values and were used to build a validated LAW25/CRASURV material card with plasticity modelling to capture the matrix-dominated non-linear response.
At the element level, three-point bending tests performed on the curved omega structures showed excellent repeatability, both quantitatively and qualitatively, with consistent cross-sectional deformation, strain field, and load–displacement response across independent instrumentation. A shell-based finite element model, developed using the validated material card in Altair HyperMesh, closely captured this non-linear response at both global and local levels, while identifying limitations in modelling delamination. Sensitivity studies were also carried out to identify the governing parameters. These findings provide an experimental-numerical framework to build robust, validated models that predict the structural response at low computational cost during composite design iterations, and pave the way for efficient modelling using advanced elements to capture the complete structural behaviour. ...
A building-block approach was adopted. At the coupon level, laminate quality assessments and mechanical tests were conducted on specimens extracted from the curved structure to establish baseline material quality and determine its in-plane material properties. These results showed good agreement with literature values and were used to build a validated LAW25/CRASURV material card with plasticity modelling to capture the matrix-dominated non-linear response.
At the element level, three-point bending tests performed on the curved omega structures showed excellent repeatability, both quantitatively and qualitatively, with consistent cross-sectional deformation, strain field, and load–displacement response across independent instrumentation. A shell-based finite element model, developed using the validated material card in Altair HyperMesh, closely captured this non-linear response at both global and local levels, while identifying limitations in modelling delamination. Sensitivity studies were also carried out to identify the governing parameters. These findings provide an experimental-numerical framework to build robust, validated models that predict the structural response at low computational cost during composite design iterations, and pave the way for efficient modelling using advanced elements to capture the complete structural behaviour.
To address this, this work develops an integrated approach that combines experimental determination of thermo-optical and thermal properties with reduced-order thermal modeling. A database of material properties for commonly used CubeSat components is generated through laboratory testing. Further, simplified thermal models are created and validated against experimental results to capture the dominant heat transfer behavior with reduced complexity.
The combined framework allows faster and more reliable preliminary thermal analysis, helping identify thermal issues early in the design process and improving overall confidence in thermal design. ...
To address this, this work develops an integrated approach that combines experimental determination of thermo-optical and thermal properties with reduced-order thermal modeling. A database of material properties for commonly used CubeSat components is generated through laboratory testing. Further, simplified thermal models are created and validated against experimental results to capture the dominant heat transfer behavior with reduced complexity.
The combined framework allows faster and more reliable preliminary thermal analysis, helping identify thermal issues early in the design process and improving overall confidence in thermal design.
Performance and flow characteristics of an active ABEP intake
Consisting of a coaxially coupled Knudsen Pump in continuous operation
Knudsen Pumps are usually operated in the transition flow regimes for its superior performance compared to the free-molecular regime. This study initially investigates the operation of Knudsen Pumps in free-molecular flow regime. This is done because in ABEP applications, the gases are usually in the free-molecular flow regimes. The channel lengths and the ratio of the size of the narrow channel to that of the wide channel is varied in the presented work. The variation of macroscopic properties along the axis of the Knudsen Pump is obtained from the results of the DSMC simulations. The results are compared to 1D analytical solution for some limit cases. It was concluded that increasing the channel lengths and maintaining a transition flow regime inside the wide channel of a Knudsen Pump brings up its compression characteristics while operating in a free-molecular flow regime.
In the second stage of this work, a hybrid intake consisting of a passive converging diffuse wall connected to an active Knudsen pump is investigated. A numerical DSMC model for such a setup is realised and verified using convergence studies. Two different configurations of single stage and five stage Knudsen Pump is studied. In the former case, the mass flow rate through the intake was varied using a given transmissivity for the intake exit and the temperature difference within each channels is varied. It was found that higher compressibility was observed when the total mass flow rate through the thruster was lower than the maximum thermal creep mass flow rate. Whereas the compressibility increased when the temperature difference was increased from Δ𝑇 = 500 to Δ𝑇 = 1000. However, further increase in temperature difference did not result in increased compressibility. In the five stage Knudsen Pump case, two different outlet transmissivities were studied. It was observed that the pressure and mass density drops down after consecutive stages of the Knudsen pump when the mass flow rate through the Knudsen Pump was higher than the maximum thermal creep mass flow rate. This resulted in no significant pressure or density compressibility. However, when the mass flow rate through the thruster was less than the maximum thermal creep mass flow rate, there was progressive increase in the mass density while the pressure continued to drop after consecutive stages. It was concluded that in-order to obtain significant compressibility the mass flow rate through the hybrid intake needs to be very low which in-turn might not yield sufficient thrust. A closed configuration where there is only mass flow during intermittent thruster firing could be ideal for obtaining higher compressibility using the Knudsen Pump ...
Knudsen Pumps are usually operated in the transition flow regimes for its superior performance compared to the free-molecular regime. This study initially investigates the operation of Knudsen Pumps in free-molecular flow regime. This is done because in ABEP applications, the gases are usually in the free-molecular flow regimes. The channel lengths and the ratio of the size of the narrow channel to that of the wide channel is varied in the presented work. The variation of macroscopic properties along the axis of the Knudsen Pump is obtained from the results of the DSMC simulations. The results are compared to 1D analytical solution for some limit cases. It was concluded that increasing the channel lengths and maintaining a transition flow regime inside the wide channel of a Knudsen Pump brings up its compression characteristics while operating in a free-molecular flow regime.
In the second stage of this work, a hybrid intake consisting of a passive converging diffuse wall connected to an active Knudsen pump is investigated. A numerical DSMC model for such a setup is realised and verified using convergence studies. Two different configurations of single stage and five stage Knudsen Pump is studied. In the former case, the mass flow rate through the intake was varied using a given transmissivity for the intake exit and the temperature difference within each channels is varied. It was found that higher compressibility was observed when the total mass flow rate through the thruster was lower than the maximum thermal creep mass flow rate. Whereas the compressibility increased when the temperature difference was increased from Δ𝑇 = 500 to Δ𝑇 = 1000. However, further increase in temperature difference did not result in increased compressibility. In the five stage Knudsen Pump case, two different outlet transmissivities were studied. It was observed that the pressure and mass density drops down after consecutive stages of the Knudsen pump when the mass flow rate through the Knudsen Pump was higher than the maximum thermal creep mass flow rate. This resulted in no significant pressure or density compressibility. However, when the mass flow rate through the thruster was less than the maximum thermal creep mass flow rate, there was progressive increase in the mass density while the pressure continued to drop after consecutive stages. It was concluded that in-order to obtain significant compressibility the mass flow rate through the hybrid intake needs to be very low which in-turn might not yield sufficient thrust. A closed configuration where there is only mass flow during intermittent thruster firing could be ideal for obtaining higher compressibility using the Knudsen Pump
This thesis investigates the feasibility of a semi-automated handheld robotic system for in situ material deposition, leveraging human guidance alongside robotic precision. To this end, a prototype handheld extrusion device was developed, integrating onboard optical flow sensors for real-time positional correction. Its control architecture enables user-guided motion while autonomously compensating for tracking errors, combining human adaptability with the repeatability of robotic control to ensure consistent material deposition.
To evaluate system performance, a surface-level defect modelled as a 175 mm long crack was introduced into a wooden substrate. Material was deposited into the defect under controlled conditions using the prototype, completing the process in 35 s. The device successfully performed real-time correction and deposition, demonstrating the feasibility of augmented handheld additive repair. Average accuracies of 1.05 mm in x and 2.45 mm in y were achieved after movements of about 300 mm in global x and 70 mm in global y. This level of precision demonstrates the basic capabilities of the system, but targeted improvements will be required to meet the demands of industrial deployment.
While developed for composite repair, the underlying technology is versatile and applicable across various domains. Its ability to deposit functional materials with spatial precision suggests potential use cases in structural health monitoring and field-deployable additive manufacturing. Future extensions, such as support for curved surfaces, vision-based localisation, and integrated non-destructive testing methods, could significantly enhance system capability. These developments may enable intelligent, semi-autonomous platforms that combine human intuition with robotic accuracy for material deposition in complex field environments. ...
This thesis investigates the feasibility of a semi-automated handheld robotic system for in situ material deposition, leveraging human guidance alongside robotic precision. To this end, a prototype handheld extrusion device was developed, integrating onboard optical flow sensors for real-time positional correction. Its control architecture enables user-guided motion while autonomously compensating for tracking errors, combining human adaptability with the repeatability of robotic control to ensure consistent material deposition.
To evaluate system performance, a surface-level defect modelled as a 175 mm long crack was introduced into a wooden substrate. Material was deposited into the defect under controlled conditions using the prototype, completing the process in 35 s. The device successfully performed real-time correction and deposition, demonstrating the feasibility of augmented handheld additive repair. Average accuracies of 1.05 mm in x and 2.45 mm in y were achieved after movements of about 300 mm in global x and 70 mm in global y. This level of precision demonstrates the basic capabilities of the system, but targeted improvements will be required to meet the demands of industrial deployment.
While developed for composite repair, the underlying technology is versatile and applicable across various domains. Its ability to deposit functional materials with spatial precision suggests potential use cases in structural health monitoring and field-deployable additive manufacturing. Future extensions, such as support for curved surfaces, vision-based localisation, and integrated non-destructive testing methods, could significantly enhance system capability. These developments may enable intelligent, semi-autonomous platforms that combine human intuition with robotic accuracy for material deposition in complex field environments.
Dataset: https://doi.org/10.4121/uuid:6bd1d545-e071-4f7b-a457-e213909b878f
...
Dataset: https://doi.org/10.4121/uuid:6bd1d545-e071-4f7b-a457-e213909b878f
The selected removal method employs plume impingement: a novel, contactless technique in which high-momentum gas jets are directed at target satellites to induce controlled trajectory changes. This avoids complex capture mechanisms and enables re-entry without requiring physical contact. The mission architecture ensures each target is individually approached and guided toward a controlled re-entry trajectory, before the "shepherd" spacecraft returns to the initial orbit for the next operation. The complete mission concept involves performing numerous controlled approach manoeuvres with each target debris gradually changing its trajectory using a custom-developed momentum transfer and transfer efficiency simulations until the debris reaches an elliptical orbit with a 381km perigee. At its final orbit the debris passively de-orbits within 7.5 months. The developed strategy is proven to be viable for debris ranging from 250- 500 kg and altitudes of 550-630 km. However, the detailed subsystem design is performed for 10 Starlink v1 satellites at 600 km circular orbit as this scenario was deemed most common for an ADR mission considering the abundance of Starlink. The mission is scheduled for deployment in January 2030 and is designed to stay within a 100 million total budget, covering all subsystem development, testing, and launch operations.
Subsystems have been developed in detail: the propulsion subsystem uses bi-propellant thrusters on a two-axis gimbal for precise plume control; the GNC and ADCS subsystems combine LIDAR, IR cameras, reaction wheels, and IMUs for autonomous attitude control, tracking, and detumbling. Power is provided by solar arrays and lithium-ion batteries to operate in the eclipse. Communication is handled via the ESA Estrack network, and all systems are designed for modularity and sustainability. During the final weeks of the DSE exercise, the team will finalise all subsystem designs, refine the plume-based and orbital targeting control algorithms, and integrate all components into a fully verified and validated spacecraft configuration.
...
The selected removal method employs plume impingement: a novel, contactless technique in which high-momentum gas jets are directed at target satellites to induce controlled trajectory changes. This avoids complex capture mechanisms and enables re-entry without requiring physical contact. The mission architecture ensures each target is individually approached and guided toward a controlled re-entry trajectory, before the "shepherd" spacecraft returns to the initial orbit for the next operation. The complete mission concept involves performing numerous controlled approach manoeuvres with each target debris gradually changing its trajectory using a custom-developed momentum transfer and transfer efficiency simulations until the debris reaches an elliptical orbit with a 381km perigee. At its final orbit the debris passively de-orbits within 7.5 months. The developed strategy is proven to be viable for debris ranging from 250- 500 kg and altitudes of 550-630 km. However, the detailed subsystem design is performed for 10 Starlink v1 satellites at 600 km circular orbit as this scenario was deemed most common for an ADR mission considering the abundance of Starlink. The mission is scheduled for deployment in January 2030 and is designed to stay within a 100 million total budget, covering all subsystem development, testing, and launch operations.
Subsystems have been developed in detail: the propulsion subsystem uses bi-propellant thrusters on a two-axis gimbal for precise plume control; the GNC and ADCS subsystems combine LIDAR, IR cameras, reaction wheels, and IMUs for autonomous attitude control, tracking, and detumbling. Power is provided by solar arrays and lithium-ion batteries to operate in the eclipse. Communication is handled via the ESA Estrack network, and all systems are designed for modularity and sustainability. During the final weeks of the DSE exercise, the team will finalise all subsystem designs, refine the plume-based and orbital targeting control algorithms, and integrate all components into a fully verified and validated spacecraft configuration.
Physics-Informed Neural Networks for Aerospace Applications
Advancing Simulation and Surrogate Modelling
Fracture and Impact Dynamics of Beam Lattices
Finite Element Modeling with Penalty Contact for Energy Dissipation in Architected Materials
This research focuses on adapting solar array designs available in literature to 3P PocketQubes, addressing their unique size, power and deployment constraints imposed by the deployer. A 4-panel array design in a wing configuration was designed, facilitating a peak power of 19.7 W. This was supported by a burn-wire mechanism and a torsion spring hinge. Structural analysis was conducted using LS-DYNA and ANSYS Mechanical to simulate deployment impact and launch vibrations, respectively. The analysis evaluated the design’s response to dynamic launch loads and mechanical impacts during deployment. The final assembly weight was 204.8 g, contributing to a specific power of 96.2 W/kg, comparable to many COTS deployable solutions for CubeSats. The procurement cost was found to be €1480 (excluding solar cell assemblies), and can be greatly lowered with higher order quantities for formation flying/distributed missions with multiple PocketQubes.
Additionally, this research examined the impact of solar array deployment configurations on power generation and orbital lifetime. Power generation was assessed using a simplified Python model, later verified through AGI STK simulations. Results indicated that the 𝛽 angle significantly influences power output across different configurations. Moreover, for peak 𝛽 angles representing noon-midnight and dusk-dawn orbits, seasonal variations were studied. This yielded an approximate 30% reduction during summer solstice for dusk-dawn orbits, but no visible change for noon-midnight orbits. For velocity aligned PocketQubes, configurations with panels mounted on the 5 × 5 cm face at a 135° deployment angle were found to be optimal. In contrast, for PocketQubes with pointing capabilities and high peak power requirements, the previously designed wing configuration was recommended.
The study of orbital lifetimes was facilitated by ESA’s DRAMA software and CROC was used to identify the minimum, average (random tumbling scenario), and maximum cross sections, contributing to a wide range of expected orbital lifetimes. In the velocity-aligned case, configurations featuring panels attached long edges (parallel to the drag force) yielded lifetimes suitable for long Earth observation missions. The latter configurations were found to be suited for shorter, technology demonstration missions. The launch date significantly impacted the results of the study, especially for configurations with higher ballistic coefficients. The results from this study however are highly idealistic, as the PocketQube angle of attack is expected to vary, largely influencing the drag area experienced. Assumptions regarding the drag coefficient and the use of NRLMSISE-00 drag model within DRAMA also limit the accuracy of the results, as seen when comparing the simulation and real observations from Delfi-PQ. ...
This research focuses on adapting solar array designs available in literature to 3P PocketQubes, addressing their unique size, power and deployment constraints imposed by the deployer. A 4-panel array design in a wing configuration was designed, facilitating a peak power of 19.7 W. This was supported by a burn-wire mechanism and a torsion spring hinge. Structural analysis was conducted using LS-DYNA and ANSYS Mechanical to simulate deployment impact and launch vibrations, respectively. The analysis evaluated the design’s response to dynamic launch loads and mechanical impacts during deployment. The final assembly weight was 204.8 g, contributing to a specific power of 96.2 W/kg, comparable to many COTS deployable solutions for CubeSats. The procurement cost was found to be €1480 (excluding solar cell assemblies), and can be greatly lowered with higher order quantities for formation flying/distributed missions with multiple PocketQubes.
Additionally, this research examined the impact of solar array deployment configurations on power generation and orbital lifetime. Power generation was assessed using a simplified Python model, later verified through AGI STK simulations. Results indicated that the 𝛽 angle significantly influences power output across different configurations. Moreover, for peak 𝛽 angles representing noon-midnight and dusk-dawn orbits, seasonal variations were studied. This yielded an approximate 30% reduction during summer solstice for dusk-dawn orbits, but no visible change for noon-midnight orbits. For velocity aligned PocketQubes, configurations with panels mounted on the 5 × 5 cm face at a 135° deployment angle were found to be optimal. In contrast, for PocketQubes with pointing capabilities and high peak power requirements, the previously designed wing configuration was recommended.
The study of orbital lifetimes was facilitated by ESA’s DRAMA software and CROC was used to identify the minimum, average (random tumbling scenario), and maximum cross sections, contributing to a wide range of expected orbital lifetimes. In the velocity-aligned case, configurations featuring panels attached long edges (parallel to the drag force) yielded lifetimes suitable for long Earth observation missions. The latter configurations were found to be suited for shorter, technology demonstration missions. The launch date significantly impacted the results of the study, especially for configurations with higher ballistic coefficients. The results from this study however are highly idealistic, as the PocketQube angle of attack is expected to vary, largely influencing the drag area experienced. Assumptions regarding the drag coefficient and the use of NRLMSISE-00 drag model within DRAMA also limit the accuracy of the results, as seen when comparing the simulation and real observations from Delfi-PQ.
For each operating point, the combustor geometry is kept fixed while the primary-zone equivalence ratios and the mixing parameter are optimised. This same methodology is applied at cruise using cruise-specific inlet conditions. The CRN predicts NOx within the measured in-flight range; BFFM2 also falls within this band, while P3T3 remains close to it. Only the CRN resolves the internal mixture structure and reaction pathways.
Water-to-fuel sweeps show that water injection consistently reduces NOx, with the strongest effect at high thrust where baseline temperatures are highest. CO increases mainly at idle and approach due to lower burnout temperatures. Reaction-pathway analysis confirms that thermal NO remains the dominant mechanism and that water suppresses existing pathways by reducing temperature and radicals.
Overall, the CRN provides an accurate and computationally efficient framework for analysing water injection and predicting emissions at both LTO and cruise, while resolving the internal combustor processes that are inaccessible to simpler correlation-based methods. ...
For each operating point, the combustor geometry is kept fixed while the primary-zone equivalence ratios and the mixing parameter are optimised. This same methodology is applied at cruise using cruise-specific inlet conditions. The CRN predicts NOx within the measured in-flight range; BFFM2 also falls within this band, while P3T3 remains close to it. Only the CRN resolves the internal mixture structure and reaction pathways.
Water-to-fuel sweeps show that water injection consistently reduces NOx, with the strongest effect at high thrust where baseline temperatures are highest. CO increases mainly at idle and approach due to lower burnout temperatures. Reaction-pathway analysis confirms that thermal NO remains the dominant mechanism and that water suppresses existing pathways by reducing temperature and radicals.
Overall, the CRN provides an accurate and computationally efficient framework for analysing water injection and predicting emissions at both LTO and cruise, while resolving the internal combustor processes that are inaccessible to simpler correlation-based methods.
A Mission-Agnostic Spacecraft System Simulator Toolkit
Development and Implementation for the ESA EnVision Mission
This thesis, conducted within the ESA EnVision project team at ESTEC, presents the development, implementation and evaluation of PAS3: a Python toolkit for Accessible SPICE-based Spacecraft Simulations. PAS3 streamlines the creation of ad-hoc, simulation-backed models for discipline-specific analysis, for any planetary mission in its implementation phase. It integrates Andrew Annex’ SpiceyPy, a procedural wrapper of NASA/JPL’s NAIF SPICE toolkit, inside an object-oriented framework extending its accessibility to non-developers. The toolkit was used to develop models for EnVision’s attitude, instrument and platform thermals, and power generation, as well as an interactive web-based 3D mission visualization.
PAS3 is an open-source package under the Permissive European Space Agency – ESA Software Community Licence and is available on the European Space Software Repository. ...
This thesis, conducted within the ESA EnVision project team at ESTEC, presents the development, implementation and evaluation of PAS3: a Python toolkit for Accessible SPICE-based Spacecraft Simulations. PAS3 streamlines the creation of ad-hoc, simulation-backed models for discipline-specific analysis, for any planetary mission in its implementation phase. It integrates Andrew Annex’ SpiceyPy, a procedural wrapper of NASA/JPL’s NAIF SPICE toolkit, inside an object-oriented framework extending its accessibility to non-developers. The toolkit was used to develop models for EnVision’s attitude, instrument and platform thermals, and power generation, as well as an interactive web-based 3D mission visualization.
PAS3 is an open-source package under the Permissive European Space Agency – ESA Software Community Licence and is available on the European Space Software Repository.
Investigating Thermal Testing Methods in the Frequency Domain
An Alternate Approach to Thermal Model Correlation
ICARUS
In-Air Capturing Apparatus for Recovering Unpropelled Stages
First, the behaviour of the ACCD as part of a larger towing system is analysed. For this, a two-dimensional (2D), quasi-steady-state Towing Model is developed, enabling a swift exploration of the ACCD’s design space, as well as a preliminary estimation of extreme tow-back loads. This Towing Model shows that the ACCD’s towing position relative to the TA is highly sensitivity to the ACCD’s Centre of Gravity (CoG); in order to stay clear of the TA’s disturbing wake zone, a maximum allowable CoG position of 1 m behind the vehicle’s nose is defined. Additionally, a significant pitch control authority is demonstrated, indicating substantial manoeuvrability. Finally, a study of extreme tow-back conditions results in the definition of a 300 kN axial towing design load - representing a 72% increase compared to previous estimates.
Next, the ACCD’s relative navigation system is studied, used to estimate the position and attitude of the RLV. A comparison between state-of-the-art systems is made, and the so-called VisNav solution is proposed - with a sensor attached to the RLV, and active beacons mounted on the ACCD. In order to identify a representative configuration for this system, a geometric VisNav Model is developed, analysing the visibility of beacons from the sensor’s Point of View (PoV). Based on comparative studies performed with this model, a design with three asymmetric rings of twelve beacons is proposed. Combined with Global Navigation Satellite System (GNSS) and Inertial Navigation System (INS) data, this configuration enables relative position and attitude estimation with an accuracy of 20 cm and 1°.
A representative electromechanical design for the ACCD is then proposed, based on results from the foregoing design space exploration. On the one hand, the selection of Commercial-off-the-Shelf (COTS) avionics is covered, as well as a preliminary dimensioning of the vehicle’s power system. Here, the ACCD’s actuators are identified to be the main drivers of its electronic design in terms of Size, Weight and Power (SWaP) footprint. Additionally, the feasibility of a battery-powered design is demonstrated. On the other hand, a Computer-Aided Design (CAD) study of mechanical subassemblies is performed, which are further analysed in terms of extreme operational loads. Most noticeably, a separation of the docking and release systems is proposed, as combining both into a single mechanism is deemed infeasible, due to the continuous presence of significant towing loads. Based on the presented design, a Bill of Materials (BoM) for the ACCD is established, yielding a total functional mass of 159.14 kg. Because of the substantial impact of the vehicle’s wings and docking system on its overall CoG, an additional 16.3 kg trim mass is required to reach the desired CoG position. Compared to previous studies, the final estimate for the ACCD’s total mass has increased by 33%.
Concluding the study, compliance with the vast majority of the postulated design requirements is demonstrated, while a number of future improvements and complementary studies are identified. These include an aerodynamic redesign and shortening of the ACCD’s shell, to eliminate the need for a trim mass. Additionally, prototyping of this unique vehicle is highly recommended, in order to further study its behaviour, and validate the theoretical design proposed in this work. ...
First, the behaviour of the ACCD as part of a larger towing system is analysed. For this, a two-dimensional (2D), quasi-steady-state Towing Model is developed, enabling a swift exploration of the ACCD’s design space, as well as a preliminary estimation of extreme tow-back loads. This Towing Model shows that the ACCD’s towing position relative to the TA is highly sensitivity to the ACCD’s Centre of Gravity (CoG); in order to stay clear of the TA’s disturbing wake zone, a maximum allowable CoG position of 1 m behind the vehicle’s nose is defined. Additionally, a significant pitch control authority is demonstrated, indicating substantial manoeuvrability. Finally, a study of extreme tow-back conditions results in the definition of a 300 kN axial towing design load - representing a 72% increase compared to previous estimates.
Next, the ACCD’s relative navigation system is studied, used to estimate the position and attitude of the RLV. A comparison between state-of-the-art systems is made, and the so-called VisNav solution is proposed - with a sensor attached to the RLV, and active beacons mounted on the ACCD. In order to identify a representative configuration for this system, a geometric VisNav Model is developed, analysing the visibility of beacons from the sensor’s Point of View (PoV). Based on comparative studies performed with this model, a design with three asymmetric rings of twelve beacons is proposed. Combined with Global Navigation Satellite System (GNSS) and Inertial Navigation System (INS) data, this configuration enables relative position and attitude estimation with an accuracy of 20 cm and 1°.
A representative electromechanical design for the ACCD is then proposed, based on results from the foregoing design space exploration. On the one hand, the selection of Commercial-off-the-Shelf (COTS) avionics is covered, as well as a preliminary dimensioning of the vehicle’s power system. Here, the ACCD’s actuators are identified to be the main drivers of its electronic design in terms of Size, Weight and Power (SWaP) footprint. Additionally, the feasibility of a battery-powered design is demonstrated. On the other hand, a Computer-Aided Design (CAD) study of mechanical subassemblies is performed, which are further analysed in terms of extreme operational loads. Most noticeably, a separation of the docking and release systems is proposed, as combining both into a single mechanism is deemed infeasible, due to the continuous presence of significant towing loads. Based on the presented design, a Bill of Materials (BoM) for the ACCD is established, yielding a total functional mass of 159.14 kg. Because of the substantial impact of the vehicle’s wings and docking system on its overall CoG, an additional 16.3 kg trim mass is required to reach the desired CoG position. Compared to previous studies, the final estimate for the ACCD’s total mass has increased by 33%.
Concluding the study, compliance with the vast majority of the postulated design requirements is demonstrated, while a number of future improvements and complementary studies are identified. These include an aerodynamic redesign and shortening of the ACCD’s shell, to eliminate the need for a trim mass. Additionally, prototyping of this unique vehicle is highly recommended, in order to further study its behaviour, and validate the theoretical design proposed in this work.
This thesis explores two main aspects: (1) optimizing TGFRP transparency through manufacturing parameters and (2) leveraging TGFRP transparency for damage detection and stress analysis. In the first phase, findings demonstrate that achieving high surface smoothness and maximizing light transmission within the green spectrum (520–600 nm) are crucial for enhancing TGFRP transparency. This is best controlled via post-curing, a more stable approach than adding methyl methacrylate (MMA) as suggested in existing literature. Sizing was also found to play a critical role in transparency, ensuring optimal bonding between glass fibers and the matrix. Additionally, minimizing the number of glass fiber fabric layers improved transparency by reducing transmission losses across visible wavelengths.
The second phase investigated the transparency of TGFRP for visualizing internal damage and stress distributions. Confocal microscopy was used to observe cracks at various depths within TGFRP samples, while image processing techniques, such as thresholding, were employed to calculate crack density. The results indicate that a higher crack density corresponds to reduced light transmittance across the visible range. However, some limitations were encountered, as cracks and scratches on the upper layers cast shadows on the lower layers, complicating the detection of deeper damage.
During tensile testing of open-hole specimens, significant opacity changes were observed in high-strain regions, suggesting that opacity correlates with localized stress. Digital image correlation (DIC) and grayscale histogram analysis provided insights into the strain thresholds at which opacity changes begin. Future work could build on this by examining the extent to which opacity changes are due to refractive index variations alone, independent of crack formation, using these strain levels as a baseline.
Furthermore, verifying the reversibility of opacity changes through post-test spectrophotometry could position TGFRPs as highly effective materials for developing innovative and reliable non-destructive testing, damage detection, and stress visualization methods in glass fiber-reinforced polymer composites.
...
This thesis explores two main aspects: (1) optimizing TGFRP transparency through manufacturing parameters and (2) leveraging TGFRP transparency for damage detection and stress analysis. In the first phase, findings demonstrate that achieving high surface smoothness and maximizing light transmission within the green spectrum (520–600 nm) are crucial for enhancing TGFRP transparency. This is best controlled via post-curing, a more stable approach than adding methyl methacrylate (MMA) as suggested in existing literature. Sizing was also found to play a critical role in transparency, ensuring optimal bonding between glass fibers and the matrix. Additionally, minimizing the number of glass fiber fabric layers improved transparency by reducing transmission losses across visible wavelengths.
The second phase investigated the transparency of TGFRP for visualizing internal damage and stress distributions. Confocal microscopy was used to observe cracks at various depths within TGFRP samples, while image processing techniques, such as thresholding, were employed to calculate crack density. The results indicate that a higher crack density corresponds to reduced light transmittance across the visible range. However, some limitations were encountered, as cracks and scratches on the upper layers cast shadows on the lower layers, complicating the detection of deeper damage.
During tensile testing of open-hole specimens, significant opacity changes were observed in high-strain regions, suggesting that opacity correlates with localized stress. Digital image correlation (DIC) and grayscale histogram analysis provided insights into the strain thresholds at which opacity changes begin. Future work could build on this by examining the extent to which opacity changes are due to refractive index variations alone, independent of crack formation, using these strain levels as a baseline.
Furthermore, verifying the reversibility of opacity changes through post-test spectrophotometry could position TGFRPs as highly effective materials for developing innovative and reliable non-destructive testing, damage detection, and stress visualization methods in glass fiber-reinforced polymer composites.
In this thesis, a method was developed to generalize the TSI-USM propagator. The TSI-USM was generalized to allow the propagation of a spacecraft trajectory in a central gravitational field subject to any pre-set thrust profile, without the need to adapt the internal working of the propagator. The three-dimensional thrust profile should be specified a priori, and with respect to time, by the user or by an evolutionary algorithm. The method used to generalize the TSI-USM requires the coefficients of the cubic spline interpolation of the pre-set thrust profile to compute the recurrence relations of the TSI-USM.
Results show that, although the method indeed generalizes the TSI-USM, it is less accurate and requires more CPU time for a given accuracy than existing RK8(7)13M-based propagators. From this result, it can be expected that future attempts to generalize the TSI-USM propagator, without sacrificing its accuracy and computational speed, will be challenging. ...
In this thesis, a method was developed to generalize the TSI-USM propagator. The TSI-USM was generalized to allow the propagation of a spacecraft trajectory in a central gravitational field subject to any pre-set thrust profile, without the need to adapt the internal working of the propagator. The three-dimensional thrust profile should be specified a priori, and with respect to time, by the user or by an evolutionary algorithm. The method used to generalize the TSI-USM requires the coefficients of the cubic spline interpolation of the pre-set thrust profile to compute the recurrence relations of the TSI-USM.
Results show that, although the method indeed generalizes the TSI-USM, it is less accurate and requires more CPU time for a given accuracy than existing RK8(7)13M-based propagators. From this result, it can be expected that future attempts to generalize the TSI-USM propagator, without sacrificing its accuracy and computational speed, will be challenging.
Mass optimization of PocketQube structures and deployers
A method demonstrated on the DelfiPQ PocketQube mission