Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Testing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60

Numerical simulation of transient hypervelocity flow in an expansion tube

Several numerical simulations of the transient flow of helium in an expansion tube are presented in an effort to identify some of the basic mechanisms which cause the noisy test flows seen in experiments. The calculations were performed with an axisymmetric Navier-Stokes code based on a finite volume formulation and upwinding techniques. Although laminar flow and ideal bursting of the diaphragms was assumed, the simulations showed some of the important features seen in experiments. In particular, the discontinuity in tube diameter of the primary diaphragm station introduced a transverse perturbation to the expanding driver gas and this perturbation was seen to propagate into the test gas under some flow conditions. The disturbances seen in the test flow can be characterized as either small amplitude, low frequency noise possibly introduced during shock compression or large amplitude, high frequency noise associated with the passage of the reflected head of the unsteady expansion.

Jacobs, P. A.↗

Altitude Starting Characteristics of an Afterburner with Autoignition and Hot-streak Ignition

An investigation was conducted in an altitude test chamber at the NACA Lewis Laboratory to determine the altitude starting characteristics of an afterburner with autoignition and with hot-streak ignition. Transient afterburner ignition data were obtained over a range of altitudes from 30,000 to 50,000 feet at a flight Mach number of 0.60. Afterburner ignition with a torch igniter located axially at approximately the midpoint of the combustion chamber was possible over the entire range, but ignition ignition with a torch igniter located in the transition section 1 5/8 inches upstream of the turbine stators proved unsatisfactory at an altitude of 50,000 feet due to the inability to obtain flame through the turbine. Increasing the afterburner-inlet total pressure at a constant afterburner fuel-air ratio decreased the afterburner ignition time. Hot-streak ignition was possible within 2 seconds after the time required to obtain the preset, normal afterburner fuel pressure, whereas autoignition required 4 to 7 seconds for the range of altitudes investigated. Following the ignition there was a period of oscillatory operation existing in the engine-afterburner before steady-state operation was attained. The time required for steady-state stable operation decreased as afterburner inlet total pressure increased. The duration of oscillations also decreased with hot-streak ignition because the fuel-air mixture was ignited before a large volume of combustible mixture was accumulated in the afterburner.

Engines, Control - Turbojet↗

Cold Start of a Radiator Equipped with Titanium-Water Heat Pipes

Radiator panels utilizing titanium-water heat pipes are being considered for lunar applications. A traditional sandwich structure is envisioned where heat pipes are embedded between two high thermal conductivity face sheets. The heat pipe evaporators are to be thermally connected to the heat source through one or more manifolds containing coolant. Initial radiator operation on the lunar surface would likely follow a cold soak where the water in the heat pipes is purposely frozen. To achieve heat pipe operation, it will be necessary to thaw the heat pipes. One option is to allow the sunlight impinging on the surface at sunrise to achieve this goal. Testing was conducted in a thermal vacuum chamber to simulate the lunar sunrise and additional modeling was conducted to identify steady-state and transient response. It was found that sunlight impinging on the radiator surface at sunrise was insufficient to solely achieve the goal of thawing the water in the heat pipes. However, starting from a frozen condition was accomplished successfully by applying power to the evaporators. Start up in this fashion was demonstrated without evaporator dryout. Concern is raised over thawing thermosyphons, vertical heat pipes operating in a gravity field, with no wick in the condenser section. This paper presents the results of the simulated cold start study and identifies future work to support radiator panels equipped with titanium-water heat pipes.

Jaworske, Donald A.↗

Heat Transfer Modeling and Validation for Optically Thick Alumina Fibrous Insulation

Combined radiation/conduction heat transfer through unbonded alumina fibrous insulation was modeled using the diffusion approximation for modeling the radiation component of heat transfer in the optically thick insulation. The validity of the heat transfer model was investigated by comparison to previously reported experimental effective thermal conductivity data over the insulation density range of 24 to 96 kg/cu m, with a pressure range of 0.001 to 750 torr (0.1 to 101.3 x 10(exp 3) Pa), and test sample hot side temperature range of 530 to 1360 K. The model was further validated by comparison to thermal conductivity measurements using the transient step heating technique on an insulation sample at a density of 144 kg/cu m over a pressure range of 0.001 to 760 torr, and temperature range of 290 to 1090 K.

Daryabeigi, Kamran↗

Output Feedback M-MRAC Backstepping With Aerospace Applications

The paper presents a certainty equivalence output feedback backstepping adaptive control design method for the systems of any relative degree with unmatched uncertainties without over-parametrization. It uses a fast prediction model to estimate the unknown parameters, which is independent of the control design. It is shown that the system's input and output tracking errors can be systematically decreased by the proper choice of the design parameters. The approach is applied to aerospace control problems and tested in numerical simulations.

Output feedback adaptive control↗

Development of Multi-Physics Dynamics Models for High-Frequency Large-Amplitude Structural Response Simulation

An analytic approach is demonstrated to reveal potential pyroshock-driven dynamic effects causing power losses in the Thermo-Electric (TE) module bars of the Mars Science Laboratory (MSL) Multi-Mission Radioisotope Thermoelectric Generator (MMRTG). This study utilizes high-fidelity finite element analysis with SIERRA/PRESTO codes to estimate wave propagation effects due to large-amplitude suddenly-applied pyro shock loads in the MMRTG. A high fidelity model of the TE module bar was created with approximately 30 million degrees-of-freedom (DOF). First, a quasi-static preload was applied on top of the TE module bar, then transient tri-axial acceleration inputs were simultaneously applied on the preloaded module. The applied input acceleration signals were measured during MMRTG shock qualification tests performed at the Jet Propulsion Laboratory. An explicit finite element solver in the SIERRA/PRESTO computational environment, along with a 3000 processor parallel super -computing framework at NASA AMES, was used for the simulation. The simulation results were investigated both qualitatively and quantitatively. The predicted shock wave propagation results provide detailed structural responses throughout the TE module bar, and key insights into the dynamic response (i.e., loads, displacements, accelerations) of critical internal spring/piston compression systems, TE materials, and internal component interfaces in the MMRTG TE module bar. They also provide confidence on the viability of this high-fidelity modeling scheme to accurately predict shock wave propagation patterns within complex structures. This analytic approach is envisioned for modeling shock sensitive hardware susceptible to intense shock environments positioned near shock separation devices in modern space vehicles and systems.

Shock Waves↗

Thermal Stability Testing of Two-Phase Thermal Control Hardware for the Surface Water Ocean Topography Mission

The thermal architecture for the Surface Water Ocean Topography mission utilizes loop heat pipes and constant conductance heat pipes to transport waste heat (> 1000 Watts) from the instrument electronics to the radiator. The main thermal design risk is the ability to maintain a temporal stability of <0.05 Celsius per minute in a low earth orbit environment. The stringent thermal requirements are part of the overall error budget needed to meet the primary mission science objectives. A testbed was developed to simulate flight-like loads and environments in order to validate the thermal subsystem can meet the temporal stability requirements. Preliminary testing showed that the thermal control system can meet the stability requirement and that loop heat pipes are actually helping to maintain stability during transient sink temperature changes for this specific flight application and boundary conditions.

Somawardhana, Ruwan P.↗

Some aspects of algorithm performance and modeling in transient thermal analysis of structures

The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms with variable time steps, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, have been selected and finite-element models of the configurations are described. These problems include a Space Shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system, and a model of the Space Shuttle Orbiter wing. Results generally indicate a preference for implicit over explicit algorithms for solution of transient structural heat transfer problems when the governing equations are 'stiff' (typical of many practical problems such as insulated metal structures).

Adelman, H. M.↗

Design and Testing of an Active Heat Rejection Radiator with Digital Turn-Down Capability

NASA's proposed lunar lander, Altair, will be exposed to vastly different external environment temperatures. The challenges to the active thermal control system (ATCS) are compounded by unfavorable transients in the internal waste heat dissipation profile: the lowest heat load occurs in the coldest environment while peak loads coincide with the warmest environment. The current baseline for this fluid is a 50/50 inhibited propylene glycol/water mixture with a freeze temperature around -35 C. While the overall size of the radiator's heat rejection area is dictated by the worst case hot scenario, a turn-down feature is necessary to tolerate the worst case cold scenario. A radiator with digital turn-down capability is being designed as a robust means to maintain cabin environment and equipment temperatures while minimizing mass and power consumption. It utilizes active valving to isolate and render ineffective any number of parallel flow tubes which span across the ATCS radiator. Several options were assessed in a trade-study to accommodate flow tube isolation and how to deal with the stagnant fluid that would otherwise remain in the tube. Bread-board environmental tests were conducted for options to drain the fluid from a turned-down leg as well an option to allow a leg to freeze/thaw. Each drain option involved a positive displacement gear pump with different methods of providing a pressure head to feed it. Test results showed that a start-up heater used to generate vapor at the tube inlet held the most promise for tube evacuation. Based on these test results and conclusions drawn from the trade-study, a full-scale radiator design is being worked for the Altair mission profile.

Sunada, Eric↗

GT SUITE Fuel Cell Model Validation with Power Module Test Data

The AES Modular Power Systems (AMPS) Fuel Cell team at the NASA Glenn Research Center (GRC) has created a Microsoft Excel model of a regenerative fuel cell (RFC) system for potential lunar and Martian applications. Due to some limitations of this model, including the ability to predict transient effects in the thermal and fluidic performance of the system, a system-level transient thermal/fluids model was desired. GT-SUITE is a transient multi-physics simulation tool that is popular in the automotive industry and has a built in template for proton exchange membrane (PEM) fuel cells. However, terrestrial fuel cells used in the automotive industry vary significantly from aerospace fuel cells. A validation effort was undertaken to compare GT SUITE model results to actual data collected during the 2015 AMPS Power Module Demonstration in which the AMPS team used a PEM fuel cell system to power a 478 kg rover assembly at the Dunes test site located at GRC. The purpose of the model was to evaluate the efficacy of using GT SUITE to model a non-flow through PEM fuel cell system. The validation effort examined electrochemical performance including fuel cell stack voltage and power supplied, the thermal performance of the stack namely stack temperature and heat rejection rates, and the fluidic performance of the system including pressures, reactant consumption, and water production rates.

Gilligan, Ryan↗

Durability of YSZ Coated Ti2AlC in 1300 °C Mach 0.3 Burner Rig Tests

A thermal barrier coating system survived burner rig testing at 1300 °C for 500 h. A 160 µm thick yttria stabilized zirconia (YSZ) coating was applied to a Ti2AlC MAX phase bar sample by plasma spray physical vapor deposition (PS-PVD) and tested face-on in an atmospheric Mach 0.3 jet fuel burner, using 5-h thermal cycles. No thermal barrier coating (TBC) spallation or recession was observed, only a 2.4 mg/cm2 mass gain. The modest weight gain precluded severe volatility losses under high velocity burner conditions. The coating surface exhibited colonies of (111)flourite fiber-textured columns separated by craze patterns, with no visible moisture attack. The metastable tetragonal t' YSZ phase was obtained initially, transitioning to equilibrium teq and cubic YSZ, but with little detrimental monoclinic. The thickness of the alumina TGO was ~21 to 23 μm under the heated YSZ face and ~13 to 15 μm on the uncoated, cooler backside. The backside exhibited removal of initial transient TiO2 nodules and partial etching of the underlying Al2O3 scale by volatile hydroxides formed in high temperature, high velocity water vapor. Aerodynamic forces produced some bending of the cantilevered sample via creep. The test indicated exceptional stability of YSZ coatings on Ti2AlC under turbine conditions, with thermal expansion matching playing a key role. The purpose of this study was to demonstrate long term durability of YSZ/MAX phase system in aggressive high temperature burner rig testing. MAX phases have been keenly studied because of their unique crystal structure and intriguing properties (Refs. 1 and 2). Having Mn+1(Al,Si)(C,N)n general composition, they are defined as ceramics, but possess unusual desirable attributes such as high conductivity, thermal shock resistance, easy machinability, and deformation tolerance. The mechanical properties derive from weak M-(Al,Si) bonding in the basal plane that leads to sliding and kinking in preference to catastrophic crack growth. Like most ceramics they are phase stable at high temperatures, generally up to 1500 °C. High temperature oxidation resistance is excellent for alumina-forming Ti3AlC2, Ti2AlC, and Cr2AlC, as reviewed by Tallman, et al. (Ref. 3). Compatibility with α-Al2O3 scales is further enhanced in cyclic exposures by a close matching of thermal expansion coefficients, (Ref. 4) i.e., (~9.3, 10.2, 11.3×10–6/K for Al2O3, Ti2AlC, and YSZ, to be discussed). Turbine environments generally contain 10 percent water vapor in the combustion gases, therefore moisture effects can be a concern for some materials (Ref. 5). Furnace tests of MAX phases in high temperature steam generally showed little effect on Al2O3 scale growth (Ref. 6). However, high velocity and high pressure gas can influence scale losses by the formation of volatile reaction products, such as TiO(OH)2 and Al(OH)3 (Refs. 7 to 10). This phenomenon had been discussed for 1100 to 1300 °C high pressure burner rig tests of Ti2AlC (Ref. 11). A single cubic growth rate parameter kcubic was measurably lower than comparable furnace TGA data, but it could be matched reasonably well if corrected for a slight volatility term. In general, a two-parameter cubic-linear growth-volatility law was believed to apply. Corresponding scale volatility loss rates, directly measured at 1300 °C on a pre-oxidized sample, were moderate (0.012 mg/cm2/h) and largely attributed to removal of the initial TiO2 transient scale. A related CH4 burner study of high purity Cr2AlC MAX phase demonstrated 1200 °C durability after 500 rapid (5 min. heat and 2 min. cool) thermal shock cycling (29 h hot time) (Ref. 12). Heating and cooling rates were ~1000 and 500 °C per minute, with a gas velocity of 5 m/s, producing a 75 °C/mm gradient. A 7 μm Al2O3 surface scale and a 13 μm Cr7C3 depletion zone formed with no signs of failure. No evidence of scale volatility was evident, although weight change was not provided, the velocity was moderate, and the total hot time was not extensive. The same high gradient BRT was used to produce 1400 °C surface temperatures for a YSZ/Cr2AlC/IN738 system in the first study of MAX phases used as bond coats for thermal barrier coatings (TBC) (Ref. 25). Here TBC failure was reported after 745 cycles, with only a 1.5 μm Al2O3 scale entrained within a porous, Cr7C3 bondcoat depletion phase. YSZ thermal barrier coatings have been considered to be a compatible complement to Al-MAX phases because of thermal expansion matching and extremely low volatility in water vapor. Initial studies showed superior oxidative stability up to 1300 °C, for long times (at least 500 h) for Ti2AlC substrates and less (268 h) for Cr2AlC, while withstanding large alumina TGO scale thickness (~35 to 40 μm) (Refs. 13 and 14). By comparison, typical superalloy systems can only survive 1150 °C maximum interface temperatures for extended periods, with a maximum sustained TGO below 10 μm (Ref. 15). High temperature SiC based systems are known to form slow-growing SiO2 scales. But these are subject to rate enhancement and volatile Si(OH)4 products in the presence of water vapor, as described comprehensively by Opila, et al. (Refs. 5, 16 to 19). Net weight losses are generally observed in high velocity, high pressure burner rig studies (e.g., 0.084 mg/cm2/h at 1300 °C) (Ref. 20). Furthermore, the loss rates have been shown from chemical physics to scale with v1/2 and pH2O 2 (Ref. 16). Low activity, moisture-resistant environmental barrier coatings (EBC), such as rare earth silicates, are needed to prevent substrate recession under turbine conditions (Refs. 21 to 23).

Smialek, James L.↗

Innovative Comparison of Transient Ignition Temperature at the Booster Interface, New Stainless Steel Pyrovalve Primer Chamber Assembly "V" (PCA) Design Versus the Current Aluminum "Y" PCA Design

An assessment of four spacecraft pyrovalve anomalies that occurred during ground testing was conducted by the NASA Engineering & Safety Center (NESC) in 2008. In all four cases, a common aluminum (Al) primer chamber assembly (PCA) was used with dual NASA Standard Initiators (NSIs) and the nearly simultaneous (separated by less than 80 microseconds) firing of both initiators failed to ignite the booster charge. The results of the assessment and associated test program were reported in AIAA Paper AIAA-2008-4798, NESC Independent Assessment of Pyrovalve Ground Test Anomalies. As a result of the four Al PCA anomalies, and the test results and findings of the NESC assessment, the Mars Science Laboratory (MSL) project team decided to make changes to the PCA. The material for the PCA body was changed from aluminum (Al) to stainless steel (SS) to avoid melting, distortion, and potential leakage of the NSI flow passages when the device functioned. The flow passages, which were interconnected in a Y-shaped configuration (Y-PCA) in the original design, were changed to a V-shaped configuration (V-PCA). The V-shape was used to more efficiently transfer energy from the NSIs to the booster. Development and qualification testing of the new design clearly demonstrated faster booster ignition times compared to the legacy AL Y-PCA design. However, the final NESC assessment report recommended that the SS V-PCA be experimentally characterized and quantitatively compared to the Al Y-PCA design. This data was deemed important for properly evaluating the design options for future NASA projects. This test program has successfully quantified the improvement of the SS V-PCA over the Al Y-PCA. A phase B of the project was also conducted and evaluated the effect of firing command skew and enlargement of flame channels to further assist spacecraft applications.

Saulsberry, Regor L.↗

Development and Operation of an Automatic Rotor Trim Control System for use During the UH-60 Individual Blade Control Wind Tunnel Test

A full-scale wind tunnel test to evaluate the effects of Individual Blade Control (IBC) on the performance, vibration, noise and loads of a UH-60A rotor was recently completed in the National Full-Scale Aerodynamics Complex (NFAC) 40- by 80-Foot Wind Tunnel [1]. A key component of this wind tunnel test was an automatic rotor trim control system that allowed the rotor trim state to be set more precisely, quickly and repeatably than was possible with the rotor operator setting the trim condition manually. The trim control system was also able to maintain the desired trim condition through changes in IBC actuation both in open- and closed-loop IBC modes, and through long-period transients in wind tunnel flow. This ability of the trim control system to automatically set and maintain a steady rotor trim enabled the effects of different IBC inputs to be compared at common trim conditions and to perform these tests quickly without requiring the rotor operator to re-trim the rotor. The trim control system described in this paper was developed specifically for use during the IBC wind tunnel test

Theodore, Colin R.↗

Thermal Considerations for Reducing the Cooldown and Warmup Duration of the James Webb Space Telescope OTIS Cryo-Vacuum Test

The James Webb Space Telescope (JWST), set to launch in 2018, is NASAs next-generation flagship telescope. The Optical Telescope Element (OTE) and Integrated Science Instrument Module (ISIM) contain all of the optical surfaces and instruments to capture and analyze the telescopes infrared targets. The integrated OTE and ISIM are denoted as OTIS, and will be tested as a single unit in a critical thermal-vacuum test in mid-2017 at NASA Johnson Space Centers Chamber A facility. The payload will be evaluated for workmanship and functionality in a 20K simulated flight environment during this thermal-vacuum test. However, the sheer thermal mass of the OTIS payload as well as the restrictive gradient, rate, and contamination-related constraints placed on test components precludes rapid cooldown or warmup to its steady-state cryo-balance condition. Hardware safety considerations precludes injection of helium gas for free molecular heat transfer. Initial thermal analysis predicted that transient radiative cooldown from ambient temperatures, while meeting all limits and constraints, would take 33.3 days; warmup similarly would take 28.4 days. This paper discusses methods used to reduce transition times from the original predictions through modulation of boundary temperatures and environmental conditions. By optimizing helium shroud transition rates and heater usage, as well as rigorously re-examining previously imposed constraints, savings of up to three days on cooldown and up to a week on warmup can be achieved. The efficiencies gained through these methods allow the JWST thermal test team to create faster cooldown and warmup profiles, thus reducing the overall test duration and cost, while keeping all of the required test operations.

Thermal↗

Thermal Considerations for Reducing the Cooldown and Warmup Duration of the James Webb Space Telescope OTIS Cryo-Vacuum Test

The James Webb Space Telescope (JWST), set to launch in 2018, is NASAs next-generation flagship telescope. The Optical Telescope Element (OTE) and Integrated Science Instrument Module (ISIM) contain all of the optical surfaces and instruments to capture and analyze the telescopes infrared targets. The integrated OTE and ISIM are denoted as OTIS, and will be tested as a single unit in a critical thermal-vacuum test in mid-2017 at NASA Johnson Space Centers Chamber A facility. The payload will be evaluated for workmanship and functionality in a 20K simulated flight environment during this thermal-vacuum test. However, the sheer thermal mass of the OTIS payload as well as the restrictive gradient, rate, and contamination-related constraints placed on test components precludes rapid cooldown or warmup to its steady-state cryo-balance condition. Hardware safety considerations precludes injection of helium gas for free molecular heat transfer. Initial thermal analysis predicted that transient radiative cooldown from ambient temperatures, while meeting all limits and constraints, would take 33.3 days; warmup similarly would take 28.4 days. This paper discusses methods used to reduce transition times from the original predictions through modulation of boundary temperatures and environmental conditions. By optimizing helium shroud transition rates and heater usage, as well as rigorously re-examining previously imposed constraints, savings of up to three days on cooldown and up to a week on warmup can be achieved. The efficiencies gained through these methods allow the JWST thermal test team to create faster cooldown and warmup profiles, thus reducing the overall test duration and cost, while keeping all of the required test operations.

JWST↗

Assessment of tbe Performance of Ablative Insulators Under Realistic Solid Rocket Motor Operating Conditions (a Doctoral Dissertation)

Ablative insulators are used in the interior surfaces of solid rocket motors to prevent the mechanical structure of the rocket from failing due to intense heating by the high-temperature solid-propellant combustion products. The complexity of the ablation process underscores the need for ablative material response data procured from a realistic solid rocket motor environment, where all of the potential contributions to material degradation are present and in their appropriate proportions. For this purpose, the present study examines ablative material behavior in a laboratory-scale solid rocket motor. The test apparatus includes a planar, two-dimensional flow channel in which flat ablative material samples are installed downstream of an aluminized solid propellant grain and imaged via real-time X-ray radiography. In this way, the in-situ transient thermal response of an ablator to all of the thermal, chemical, and mechanical erosion mechanisms present in a solid rocket environment can be observed and recorded. The ablative material is instrumented with multiple micro-thermocouples, so that in-depth temperature histories are known. Both total heat flux and thermal radiation flux gauges have been designed, fabricated, and tested to characterize the thermal environment to which the ablative material samples are exposed. These tests not only allow different ablative materials to be compared in a realistic solid rocket motor environment but also improve the understanding of the mechanisms that influence the erosion behavior of a given ablative material.

Martin, Heath Thomas↗

Molecular flow venting of a volume with an outgassing or desorbing source

A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution.

Michael Woronowicz↗

Molecular flow venting of a volume with an outgassing or desorbing source

"A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution."

Michael Woronowicz↗