SEARCH · Engineering Papers

Results for “BALLISTIC TRAJECTORY”

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.

SourceResearch connectionFrom yearThrough year

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

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Along-Trajectory Acoustic Signal Variations Observed During the Hypersonic Re-Entry of the OSIRIS-REx Sample Return Capsule

The re-entry of the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) sample return capsule (SRC) on 24 September 2023 presented a rare opportunity to study atmospheric entry dynamics through a dense network of ground-based infrasound sensors. As the first interplanetary capsule to re-enter over the United States since Stardust in 2006, this event allowed for unprecedented observations of infrasound signals generated during hypersonic descent. We deployed 39 single-sensor stations across Nevada and Utah, strategically distributed to capture signals from distinct trajectory points. Infrasound data were analyzed to examine how signal amplitude and period vary with altitude and propagation path for a nonablating hypersonic object with well-defined physical and aerodynamic properties. Raytracing simulations incorporated atmospheric specifications from the ground-2-space model to estimate source altitudes for observed signals. Results confirmed ballistic arrivals at all stations, with source altitudes ranging from 44 to 62 km along the trajectory. Signal period and amplitude exhibited strong dependence on source altitude, with higher altitudes corresponding to lower amplitudes, longer periods, and reduced high-frequency content. Regression analysis demonstrated strong correlations between signal characteristics and both altitude and propagation geometry. Our results suggest, when attenuation is considered, the amplitude is primarily determined by the source, with the propagation path playing a secondary role over the distances examined. These findings emphasize the utility of controlled SRC re-entries for advancing our understanding of natural meteoroid dynamics, refining atmospheric entry models, and improving methodologies for planetary defense. The OSIRIS-REx SRC campaign represents the most comprehensive infrasound study of a hypersonic re-entry to date, showcasing the potential of coordinated geophysical observational networks for high-energy atmospheric phenomena, including space debris re-entries.

58 GEOSCIENCES↗

NASA NTRS1973

Linear filtering of ballistic-entry-probe data for atmospheric reconstruction.

A Kalman-Schmidt filter is used to estimate atmospheric and trajectory parameters for entry into the Venusian atmosphere. A significantly improved version of the Landing Trajectory Reconstruction (LTR) computer program, used to obtain the estimates is described. Major improvements involve precision and linearity control of numerical differencing, corrected perturbation modeling, and incorporation of a refractivity model. Important results show that gyroscopic data are not usable with LTR, that atmospheric properties are generally well estimated for altitudes less than 100 km, and that the two LTR modes of operation are complementary in performance.

Sabin, M. L.↗

NASA NTRS2018

Supersonic Retropropulsion on Robotic Mars Landers: Selected Design Trades

Many concepts for future robotic Mars lander missions require landing heavier payloads than those landed to date. Mars lander architectures to date have relied on a parachute to help slow the lander; however the effectiveness of a parachute in the thin Martian atmosphere is diminished with heavier payloads unless the diameter of the parachute is increased or it is deployed at a higher Mach number, both of which are significant technical challenges. In addition the parachute can be successfully deployed only within a specific Mach number and dynamic pressure range. Targeting the entry trajectory to hit this “Mach-Q box” imposes constraints on the entry ballistic coefficient, limiting it to ~ 150-200 kg/m^2. Eliminating the parachute from the design requires descent engine ignition at supersonic speeds (Supersonic Retropropulsion, or SRP). SRP increases the propellant requirement, but also allows entry ballistic coefficients of ~600 kg/m^2 or more, with the consequence of significantly increased entry mass and landed payload mass.

Lorenz, Christopher G.↗

NASA NTRSDate not supplied

Keplerian Analysis for Versatile Evaluation of Arbitrary Trajectories

Designing interplanetary missions is an iterative process, with a tight coupling between mission analysis and vehicle design. The ability to simultaneously process both analyses early on can provide significant benefits in both the concept formulation and in overall mission feasibility. To implement such a simultaneous analysis process, there is a need to model interplanetary trajectories rapidly while maintaining an acceptable level of fidelity. Many simple patched conics tools make too many assumptions to satisfactorily address all relevant constraints and objectives, while full-fidelity trajectory design tools often require more time and expertise than is available in early conceptual studies. Hence, some form of middle-ground modeling tool would provide value. This paper introduces the development of such a tool using Python. This tool, Keplerian Analysis for Versatile Evaluation of Arbitrary Trajectories (KAVEAT), implements a patched conics method that enforces continuity of state (time, mass, position, velocity) through all transitions between various spheres of influence on the trajectory. The trajectory is constructed in an object-oriented manner, dynamically assembling unique maneuvers such as departure burns, heliocentric transfers (ballistic or with electric propulsion), and planetary flybys (powered or unpowered). Each maneuver has its own specific set of basic inputs needed to define it. The tool then takes these defined maneuvers, strings them together creating all the interdependencies required, and solves for the ideal trajectory using a gradient-based optimizer.

Katherine T McBrayer↗

NASA NTRS2000

Earth Return Aerocapture for the TransHab/Ellipsled Vehicle

The current architecture being considered by NASA for a human Mars mission involves the use of an aerocapture procedure at Mars arrival and possibly upon Earth return. This technique would be used to decelerate the vehicles and insert them into their desired target orbits, thereby eliminating the need for propulsive orbital insertions. The crew may make the interplanetary journey in a large, inflatable habitat known as the TransHab. It has been proposed that upon Earth return, this habitat be captured into orbit for use on subsequent missions. In this case, the TransHab would be complimented with an aeroshell, which would protect it from heating during the atmospheric entry and provide the vehicle with aerodynamic lift. The aeroshell has been dubbed the "Ellipsled" because of its characteristic shape. This paper reports the results of a preliminary study of the aerocapture of the TransHab/Ellipsled vehicle upon Earth return. Undershoot and overshoot boundaries have been determined for a range of entry velocities, and the effects of variations in the atmospheric density profile, the vehicle deceleration limit, the maximum vehicle roll rate, the target orbit, and the vehicle ballistic coefficient have been examined. A simple, 180 degree roll maneuver was implemented in the undershoot trajectories to target the desired 407 km circular Earth orbit. A three-roll sequence was developed to target not only a specific orbital energy, but also a particular inclination, thereby decreasing propulsive inclination changes and post-aerocapture delta-V requirements. Results show that the TransHab/Ellipsled vehicle has a nominal corridor width of at least 0.7 degrees for entry speeds up to 14.0 km/s. Most trajectories were simulated using continuum flow aerodynamics, but the impact of high-altitude viscous effects was evaluated and found to be minimal. In addition, entry corridor comparisons have been made between the TransHab/Ellipsled and a modified Apollo capsule which is also being considered as the crew return vehicle; because of its slightly higher lift-to-drag ratio, the TransHab has a modest advantage with regard to corridor width. Stagnation-point heating rates and integrated heat loads were determined for a range of vehicle ballistic coefficients and entry velocities.

Muth, W. D.↗

NASA NTRSDate not supplied

Aerocapture Trajectory Design for Uranus Orbiter

Introduction: The recently released National Academies Planetary Science and Astrobiology Decadal Survey 2023-2032 [1] identified the Ice Giants as the top priority science destination. While the survey acknowledged the potential for either a Uranus Orbiter and Probe (UOP) mission or a Neptune-Triton Odyssey mission, it ultimately identified the former as the highest priority new flagship mission. UOP missions calls for a launch window of opportunity between 2031-2038 with 12-15 year interplanetary cruise time along with a fully-propulsive Uranus Orbit Insertion burn on the order of a few km/s. However, a mission to Uranus with the same science payload could utilize aerocapture for orbit insertion to achieve both a significant reduction in the interplanetary cruise time and reduction in propulsive burn costs. Why Aerocapture: Aerocapture is a promising propellant and time-saving orbital insertion technique for planetary destinations with an atmosphere. Although not flight-proven, previous aerocapture systems studies in the literature have demonstrated both the validity and robustness of the technique at various planetary destinations. With respect to the Ice Giant planets, Neptune has seen more of the analysis in the literature. For science missions at Neptune, aerocapture can enable 1.4 times more delivered mass to orbit than an all-propulsive mission for the same launch vehicle while reducing interplanetary cruise times by more than 3 years [2] Additionally with modern guidance and control, Neptune aerocapture with blunt-body aeroshells is realizable [3][4]. There are limited papers in the literature investigating Uranus aerocapture with those available providing a preliminary feasibility assessment [5]. Consequently, the two-year funded NASA Space Technology Mission Directorate (STMD)-funded project, titled Aerocapture System as an Enabling Technology for Ice Giants Missions, aims to mature the analysis and technology state of Uranus aerocapture. Trajectory Design: This paper presents the current state of the trajectory design in support of the new aerocapture project. The project design philosophy is inspired from recent Neptune aerocapture studies, which employed modern guidance and control, in the sense that blunt-body aeroshells are analyzed. An assessment of the theoretical flight path angle corridor width is conducted for a range of ballistic coefficients and lift-to-drag ratios for both Space Launch System and Falcon Heavy Launch Vehicle interplanetary trajectory solutions. The results from the corridor width assessment provide an assessment of the aerocapture design-space and qualitative metrics on trajectory design considerations. The Program to Optimize Simulated Trajectories II (POST2) is utilized to run Monte Carlo simulations of Uranus aerocapture three-degree-of-freedom bank angle modulated trajectories using a closed-loop numerical-predictor corrector guidance algorithm. UranusGRAM 2021 is utilized as the atmospheric model [6]. A Uranus-developed aerodatabase, originally derived from Mars Science Laboratory (MSL), is utilized to provide vehicle aerodynamics over a wide range of hypersonic flow regimes. A MSL-derived 70 deg 4.5m diameter sphere-cone aeroshell that houses the UOP payload mass is assumed. Robustness testing and performance analysis is conducted, including the assessment of entry state errors, atmosphere density variations, and aerodynamic dispersions. Post-aerocapture Delta-V and aerothermal statistics are formulated into propellant mass and TPS requirements. The results presented in the paper will demonstrate the trajectory viability of Uranus aerocapture. Preliminary Results: Preliminary trajectory design results indicates successful Uranus aerocapture with a blunt-body aeroshell housing the same payload mass as the UOP mission from an Earth-to-Uranus interplanetary trajectory arriving in less than 7 years. From this interplanetary trajectory, aerocapture provides an orbit insertion Delta-V capability of 6.9 km/s requiring less than 300 m/s for post-aerocapture correction burns (15% of wet mass allocated to propellant). To put this into perspective, the UOP study utilized an Earth-to-Uranus trajectory that arrives in 13 years and requires more than 1000 m/s for fully-propulsive orbit insertion (40% of wet mass allocated to propellant). Achieving the same 6.9 km/s Delta-V capability fully-propulsively is mass prohibitive (97% of wet mass allocated to propellant). Nevertheless, aerocapture has the potential to reduce interplanetary transit times to Uranus by half while delivering the same payload mass to orbit in a reduced propellant mass footprint.

Rohan Deshmukh↗

DOE OSTI2025

Spin-Vibronic Effect in Photoinduced Electron Transfer

The spin-vibronic effect (SVE) accelerates quantum-mechanically forbidden electronic transitions, but experimental manifestations of this phenomenon remain limited. In this contribution, the role of SVE in photoinduced electron transfer (PET) dynamics is probed by coherent vibrational wavepacket (CVWP) motions in Pt(II) dimer-naphthalene diimide donor–acceptor complexes. Metal–metal-to-ligand charge-transfer (MMLCT) excitation of the donor triggers ballistic PET, driving the molecule away from equilibrium. The SVE then directs the PET trajectory from the 1 MMLCT excited state to an intermediate ligand-centered triplet state, modulating the CVWP dynamics along the reaction coordinates. Herzberg–Teller type oscillations indicate that the Pt–Pt stretching vibration, arising from the formation of the excited triplet intermediate, serves as the reaction coordinate for charge separation. Finally, our experimental findings advance the understanding of the interplay between electronic, vibrational, and spin degrees of freedom in ultrafast photoactivated processes.

Kim, Pyosang [Argonne National Laboratory (ANL), A↗

NASA NTRS1958

A Feasibility Study of the Flare-Cylinder Configuration as a Reentry Body Shape for an Intermediate Range Ballistic Missile

A study has been made of a flare-cylinder configuration to investigate its feasibility as a reentry body of an intermediate range ballistic missile. Factors considered were heating, weight, stability, and impact velocity. A series of trajectories covering the possible range of weight-drag ratios were computed for simple truncated nose shapes of varying pointedness, and hence varying weight-drag ratios. Four trajectories were chosen for detailed temperature computation from among those trajectories estimated to be possible. Temperature calculations were made for both "conventional" (for example, copper, Inconel, and stainless steel) and "unconventional" (for example, beryllium and graphite) materials. Results of the computations showed that an impact Mach number of 0.5 was readily obtainable for a body constructed from conventional materials. A substantial increase in subsonic impact velocity above a Mach number of 0.5 was possible without exceeding material temperature limits. A weight saving of up to 134 pounds out of 822 was possible with unconventional materials. This saving represents 78 percent of the structural weight. Supersonic impact would require construction of the body from unconventional materials but appeared to be well within the range of attainability.

Garland, B. J.↗

NASA NTRS2016

HIAD on ULA (HULA) Orbital Reentry Flight Experiment Concept

This paper describes a proposed orbital velocity reentry flight test of a Hypersonic Inflatable Aerodynamic Decelerator (HIAD). The flight test builds upon ground development activities that continue to advance the materials, design, and manufacturing techniques for the inflatable structure and flexible thermal protection system (F-TPS) that comprise the inflatable heat shield. While certain aspects of material and system performance can be assessed using a variety of ground testing capabilities, only orbital velocity energy on a trajectory through the gradient density of the atmosphere can impart the combined aerodynamic and aeroheating design environments in real time. To achieve this at limited cost, the HIAD would be delivered to a spin-stabilized entry trajectory as a secondary payload on the Centaur stage of a United Launch Alliance (ULA) Atlas V launch vehicle. Initial trajectory studies indicate that the combination of launch vehicle capability and achievable reentry vehicle ballistic numbers make this a strategic opportunity for technology development. This 4 to 6 meter diameter scale aeroshell flight, referred to as HIAD on ULA (HULA), would also contribute to ULA asset recovery development. ULA has proposed that a HIAD be utilized as part of the Sensible, Modular, Autonomous Return Technology (SMART) initiative to enable recovery of the Vulcan launch vehicle booster main engines [1], including a Mid-Air Recovery (MAR) to gently return these assets for reuse. Whereas HULA will attain valuable aerothermal and structural response data toward advancing HIAD technology, it may also provide a largest-to-date scaled flight test of the MAR operation, which in turn would allow the examination of a nearly pristine post-entry aeroshell. By utilizing infrared camera imaging, HULA will also attain aft-side thermal response data, enhancing understanding of the aft side aerothermal environment, an area of high uncertainty. The aeroshell inflation will utilize a heritage design compressed gas system to minimize development costs. The data will be captured to both an onboard recorder and a recorder that is jettisoned and recovered separately from the reentry vehicle to mitigate risk. This paper provides an overview, including the architecture and flight concept of operations, for the proposed HULA flight experiment.

Dinonno, J. M.↗

NASA NTRS1972

Mission strategy for combined comet-asteroid flybys.

Review of the currently consolidating interest in one or more relatively inexpensive small body exploration missions, and discussion of a proposed ballistic flyby of two asteroids and the periodic comet Forbes, to be launched in 1977. Trajectory data and target encounter conditions are presented for a selected set of missions which are compatible with the propulsion capabilities of a modified Pioneer F & G spacecraft.

Brooks, D. R.↗

NASA NTRS2014

Improving FermiI Orbit Determination and Prediction in an Uncertain Atmospheric Drag Environment

Orbit determination and prediction of the Fermi Gamma-ray Space Telescope trajectory is strongly impacted by the unpredictability and variability of atmospheric density and the spacecrafts ballistic coefficient. Operationally, Global Positioning System point solutions are processed with an extended Kalman filter for orbit determination, and predictions are generated for conjunction assessment with secondary objects. When these predictions are compared to Joint Space Operations Center radar-based solutions, the close approach distance between the two predictions can greatly differ ahead of the conjunction. This work explores strategies for improving prediction accuracy and helps to explain the prediction disparities. Namely, a tuning analysis is performed to determine atmospheric drag modeling and filter parameters that can improve orbit determination as well as prediction accuracy. A 45 improvement in three-day prediction accuracy is realized by tuning the ballistic coefficient and atmospheric density stochastic models, measurement frequency, and other modeling and filter parameters.

orbit prediction↗

NASA NTRS2017

High Ballistic Coefficient Mars EDL with Supersonic Retropropulsion

Future Mars missions will require the ability to land increasingly heavy systems with high precision, possibly at high altitude landing sites. A key obstacle in landing systems with high ballistic coefficient is Mars’ very thin atmosphere, approximately one hundred times less dense than Earth’s atmosphere. As a result, hypersonic deceleration occurs slowly throughout atmospheric entry and additional deceleration methods are required before landing. Previous landers up to and including MSL relied on Viking-heritage technologies including a Disk-Gap-Band (DGB) supersonic parachute and a blunt body aeroshell. Safe and effective deployment of a DGB parachute can be accomplished only within prescribed bounds on dynamic pressure and Mach number.2 In addition to the restrictive constraint this imposes on the entry trajectory, there is a concern that for a parachute to effectively slow an entry vehicle with a high ballistic coefficient, its diameter must be large enough to impose certification costs, modeling uncertainty, or mission risk that could be unacceptable. This has prompted investigation into alternative mission concepts based on use of supersonic retropropulsion (SRP) without a parachute. This work investigates SRP as an enabling technology for Mars EDL of high ballistic coefficient vehicles.

Noyes, Connor↗

NASA NTRS2010

Targeting Low-Energy Ballistic Lunar Transfers

Numerous low-energy ballistic transfers exist between the Earth and Moon that require less fuel than conventional transfers, but require three or more months of transfer time. An entirely ballistic lunar transfer departs the Earth from a particular declination at some time in order to arrive at the Moon at a given time along a desirable approach. Maneuvers may be added to the trajectory in order to adjust the Earth departure to meet mission requirements. In this paper, we characterize the (Delta)V cost required to adjust a low-energy ballistic lunar transfer such that a spacecraft may depart the Earth at a desirable declination, e.g., 28.5(white bullet), on a designated date. This study identifies the optimal locations to place one or two maneuvers along a transfer to minimize the (Delta)V cost of the transfer. One practical application of this study is to characterize the launch period for a mission that aims to launch from a particular launch site, such as Cape Canaveral, Florida, and arrive at a particular orbit at the Moon on a given date using a three-month low-energy transfer.

Low-energy transfers↗

DOE OSTI2021

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2001

Aerothermal Analysis of the Project Fire II Afterbody Flow

Computational fluid dynamics (CFD) is used to simulate the wake flow and afterbody heating of the Project Fire II ballistic reentry to Earth at 11.4 km/sec. Laminar results are obtained over a portion of the trajectory between the initial heat pulse and peak afterbody heating. Although non-catalytic forebody convective heating results are in excellent agreement with previous computations, initial predictions of afterbody heating were about a factor of two below the experimental values. Further analysis suggests that significant catalysis may be occurring on the afterbody heat shield. Computations including finite-rate catalysis on the afterbody surface are in good agreement with the data over the early portion of the trajectory, but are conservative near the peak afterbody heating point, especially on the rear portion of the conical frustum. Further analysis of the flight data from Fire II shows that peak afterbody heating occurs before peak forebody heating, a result that contradicts computations and flight data from other entry vehicles. This result suggests that another mechanism, possibly pyrolysis, may be occurring during the later portion of the trajectory, resulting in less total heat transfer than the current predictions.

Wright, Michael J.↗

NASA NTRS2014

Improving Fermi Orbit Determination and Prediction in an Uncertain Atmospheric Drag Environment

Orbit determination and prediction of the Fermi Gamma-ray Space Telescope trajectory is strongly impacted by the unpredictability and variability of atmospheric density and the spacecraft's ballistic coefficient. Operationally, Global Positioning System point solutions are processed with an extended Kalman filter for orbit determination, and predictions are generated for conjunction assessment with secondary objects. When these predictions are compared to Joint Space Operations Center radar-based solutions, the close approach distance between the two predictions can greatly differ ahead of the conjunction. This work explores strategies for improving prediction accuracy and helps to explain the prediction disparities. Namely, a tuning analysis is performed to determine atmospheric drag modeling and filter parameters that can improve orbit determination as well as prediction accuracy. A 45% improvement in three-day prediction accuracy is realized by tuning the ballistic coefficient and atmospheric density stochastic models, measurement frequency, and other modeling and filter parameters.

orbit prediction↗

NASA NTRS2009

Tracking Debris Shed by a Space-Shuttle Launch Vehicle

The DEBRIS software predicts the trajectories of debris particles shed by a space-shuttle launch vehicle during ascent, to aid in assessing potential harm to the space-shuttle orbiter and crew. The user specifies the location of release and other initial conditions for a debris particle. DEBRIS tracks the particle within an overset grid system by means of a computational fluid dynamics (CFD) simulation of the local flow field and a ballistic simulation that takes account of the mass of the particle and its aerodynamic properties in the flow field. The computed particle trajectory is stored in a file to be post-processed by other software for viewing and analyzing the trajectory. DEBRIS supplants a prior debris tracking code that took .15 minutes to calculate a single particle trajectory: DEBRIS can calculate 1,000 trajectories in .20 seconds on a desktop computer. Other improvements over the prior code include adaptive time-stepping to ensure accuracy, forcing at least one step per grid cell to ensure resolution of all CFD-resolved flow features, ability to simulate rebound of debris from surfaces, extensive error checking, a builtin suite of test cases, and dynamic allocation of memory.

Stuart, Phillip C.↗

DOE OSTI2021

(U) A Theoretical Study of Asay Foil Trajectories (V.02)

We consider the trajectory of a generic Asay foil ejecta momentum diagnostic sensor, for scenarios where ejecta are produced at a planar surface and fly ballistically through a perfect vacuum to the sensor. To do so, we build upon a previously established mathematical formalism derived for the analytic study of stationary sensors (i.e., piezopins). First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source areal mass function (“source model”). This defines an integro-differential equation (IDE) for the foil trajectory. When ejecta production is instantaneous - as is generally assumed in momentum diagnostic data analyses - the IDE leads to an implicit and easily calculable closed-form solution for the foil trajectory in a perfect system, as long as the ejecta particle velocity distribution is twice-integrable. General properties of the instant-production trajectory solution indicate the existence of a boundary condition the particle velocity distribution must satisfy in order for the analytically predicted foil trajectories to be compatible with certain features commonly observed in foil data. This boundary condition is identical to one derived previously from a consideration of piezopin data. Armed with the analytic solution for instant production, we also consider various techniques used to extract time-dependent cumulative ejecta masses from foil trajectories, and derive an expression for the error imposed by using an approximated equation of motion. This analytic trajectory solution furthermore makes it possible to examine the common practice of presenting inferred cumulative ejecta masses as a function of a normalized implied velocity; we derive conditions under which this methodology is and is not meaningful. We also propose a strategy for extending the instant-production trajectory solution to time-dependent source functions.

42 ENGINEERING↗