Search NASASearch

SEARCH · Search NASA

Results for “consistent discretization”

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 127 records · Page 7

An empirically derived figure of merit for the quality of overall task performance

The need to develop an operationally relevant figure of merit for the quality of performance of a complex system such as an aircraft cockpit stems from a hypothesized dissociation between measures of performance and those of workload. Performance can be measured in terms of time, errors, or a combination of these. In most tasks performed by expert operators, errors are relatively rare and often corrected in time to avoid consequences. Moreover, perfect performance is seldom necessary to accomplish a particular task. Moreover, how well an expert performs a complex task consisting of a series of discrete cognitive tasks superimposed on a continuous task, such as flying an aircraft, does not depend on how well each discrete task is performed, but on their smooth sequencing. This makes the amount of time spent on each subtask of paramount importance in measuring overall performance, since smooth sequencing requires a minimum amount of time spent on each task. Quality consists in getting tasks done within a crucial time interval while maintaining acceptable continuous task performance. Thus, a figure of merit for overall quality of performance should be primarily a measure of time to perform discrete subtasks combined with a measure of basic vehicle control. Thus, the proposed figure of merit requires doing a task analysis on a series of performance, or runs, of a particular task, listing each discrete task and its associated time, and calculating the mean and standard deviation of these times, along with the mean and standard deviation of tracking error for the whole task. A set of simulator data on 30 runs of a landing task was obtained and a figure of merit will be calculated for each run. The figure of merit will be compared for voice and data link, so that the impact of this technology on total crew performance (not just communication performance) can be assessed. The effect of data link communication on other cockpit tasks will also be considered.

Lemay, Moira

An algebra of discrete event processes

This report deals with an algebraic framework for modeling and control of discrete event processes. The report consists of two parts. The first part is introductory, and consists of a tutorial survey of the theory of concurrency in the spirit of Hoare's CSP, and an examination of the suitability of such an algebraic framework for dealing with various aspects of discrete event control. To this end a new concurrency operator is introduced and it is shown how the resulting framework can be applied. It is further shown that a suitable theory that deals with the new concurrency operator must be developed. In the second part of the report the formal algebra of discrete event control is developed. At the present time the second part of the report is still an incomplete and occasionally tentative working paper.

Heymann, Michael

Analytical study of blockage- and lift-interference corrections for slotted tunnels obtained by the substitution of an equivalent homogeneous boundary for the discrete slots

The solid-blockage interference for a doublet on the tunnel axis and the boundary interference for lifting wings in circular, rectangular, and two-dimensional slotted tunnels have been calculated by substituting an equivalent homogeneous boundary for the physical boundary of discrete slots. In the case of small wings, the interference calculated with the assumption of homogeneity has been found to be consistent with that calculated for the discrete slots for as few as four slots in a circular tunnel. Furthermore, available experimental results for blockage interference are consistent with the results of the present analysis. As a consequence of the assumption of homogeneity it is possible to express the interference of multislotted tunnels as a function of a single parameter which combines the effects of two physical variables: the ratio of open to total slotted wall perimeter and the number of slots. A curve is presented which permits the rapid evaluation of this parameter and numerical results for lift and blockage interference are plotted against the parameter.

Davis, Don D , Jr

Acquisition of Skill Proficiency Over Multiple Sessions of a Novel Rover Simulation

Following long-duration exploration transits, adaptive changes in sensorimotor function may impair the crew's ability to safely perform manual control tasks such as operating pressurized rovers. Postflight performance will also be influenced by the level of preflight skill proficiency they have attained. The purpose of this study was to characterize the acquisition of skills in a motion-based rover simulation over multiple sessions, and to investigate the effects of varying the simulation scenarios. METHODS: Twenty healthy subjects were tested in 5 sessions, with 1-3 days between sessions. Each session consisted of a serial presentation of 8 discrete tasks to be completed as quickly and accurately as possible. Each task consisted of 1) perspective-taking, using a map that defined a docking target, 2) navigation toward the target around a Martian outpost, and 3) docking a side hatch of the rover to a visually guided target. The simulator utilized a Stewart-type motion base (CKAS, Australia), single-seat cabin with triple scene projection covering 150 deg horizontal by 50 deg vertical, and joystick controller. Subjects were randomly assigned to a control group (tasks identical in the first 4 sessions) or a varied-practice group. The dependent variables for each task included accuracy toward the target and time to completion. RESULTS: The greatest improvements in time to completion occurred during the docking phase. The varied-practice group showed more improvement in perspective-taking accuracy. Perspective-taking accuracy was also affected by the relative orientation of the rover to the docking target. Skill acquisition was correlated with self-ratings of previous gaming experience. DISCUSSION: Varying task selection and difficulty will optimize the preflight acquisition of skills when performing novel operational tasks. Simulation of operational manual control will provide functionally relevant evidence regarding the impact of sensorimotor adaptation on early surface operations and what countermeasures are needed. Learning Objective: The use of a motion-based simulation to investigate decrements in the proficiency to operate pressurized rovers during early surface operations of space exploration missions, along with the acquisition of skill proficiency needed during the preflight phase of the mission.

Dean, S. L.

Hyperspectral Fluorescence and Reflectance Imaging Instrument

The system is a single hyperspectral imaging instrument that has the unique capability to acquire both fluorescence and reflectance high-spatial-resolution data that is inherently spatially and spectrally registered. Potential uses of this instrument include plant stress monitoring, counterfeit document detection, biomedical imaging, forensic imaging, and general materials identification. Until now, reflectance and fluorescence spectral imaging have been performed by separate instruments. Neither a reflectance spectral image nor a fluorescence spectral image alone yields as much information about a target surface as does a combination of the two modalities. Before this system was developed, to benefit from this combination, analysts needed to perform time-consuming post-processing efforts to co-register the reflective and fluorescence information. With this instrument, the inherent spatial and spectral registration of the reflectance and fluorescence images minimizes the need for this post-processing step. The main challenge for this technology is to detect the fluorescence signal in the presence of a much stronger reflectance signal. To meet this challenge, the instrument modulates artificial light sources from ultraviolet through the visible to the near-infrared part of the spectrum; in this way, both the reflective and fluorescence signals can be measured through differencing processes to optimize fluorescence and reflectance spectra as needed. The main functional components of the instrument are a hyperspectral imager, an illumination system, and an image-plane scanner. The hyperspectral imager is a one-dimensional (line) imaging spectrometer that includes a spectrally dispersive element and a two-dimensional focal plane detector array. The spectral range of the current imaging spectrometer is between 400 to 1,000 nm, and the wavelength resolution is approximately 3 nm. The illumination system consists of narrowband blue, ultraviolet, and other discrete wavelength light-emitting-diode (LED) sources and white-light LED sources designed to produce consistently spatially stable light. White LEDs provide illumination for the measurement of reflectance spectra, while narrowband blue and UV LEDs are used to excite fluorescence. Each spectral type of LED can be turned on or off depending on the specific remote-sensing process being performed. Uniformity of illumination is achieved by using an array of LEDs and/or an integrating sphere or other diffusing surface. The image plane scanner uses a fore optic with a field of view large enough to provide an entire scan line on the image plane. It builds up a two-dimensional image in pushbroom fashion as the target is scanned across the image plane either by moving the object or moving the fore optic. For fluorescence detection, spectral filtering of a narrowband light illumination source is sometimes necessary to minimize the interference of the source spectrum wings with the fluorescence signal. Spectral filtering is achieved with optical interference filters and absorption glasses. This dual spectral imaging capability will enable the optimization of reflective, fluorescence, and fused datasets as well as a cost-effective design for multispectral imaging solutions. This system has been used in plant stress detection studies and in currency analysis.

Ryan, Robert E.

Numerical methods for the calculation of three-dimensional nozzle exhaust flow fields

Numerical codes developed for the calculation of three-dimensional nozzle exhaust flow fields associated with hypersonic airbreathing aircraft are described. Both codes employ reference plane grid networks with respect to three coordinate systems. Program CHAR3D is a characteristic code utilizing a new wave preserving network within the reference planes, while program BIGMAC is a finite difference code utilizing conservation variables and a one-sided difference algorithm. Secondary waves are numerically captured by both codes, while the underexpansion shock and plume boundary are treated discretely. The exhaust gas properties consist of hydrogen-air combustion product mixtures in local chemical equilibrium. Nozzle contours are treated by a newly developed geometry package based on dual cubic splines. Results are presented for simple configurations demonstrating two- and three-dimensional multiple wave interactions.

Dash, S. M.

The temporal structure of RF radiation from lightning

Radiation from lightning in the RF band from 3-300 MHz has been examined. Simultaneous measurements were made of the RF radiation from lightning together with records of fast and slow field changes. Continuous analogue recordings were made with a system having 300 kHz of bandwidth in the RF channels. The temporal history of RF radiation of these frequencies consists of a sequence of discrete pulses. The data reveal a distinct pattern in the radiation which is independent of frequency and depends on the type of lightning flash: Cloud-to-ground flashes are characterized by an abrupt beginning associated with the stepped leader, whereas cloud-to-cloud flashes begin with a slower train of noise pulses more typical of the end of both types of flash. An exception to this pattern is cloud-to-ground flashes preceded by a breakdown phase, in which case the radiation begins like a cloud-to-cloud flash.

Levine, D. M.

Cross correlation analysis of medium energy gamma rays for the Northern Hemisphere

In the cross correlation method the observed gamma rays are compared with the expected telescope response for a discrete celestial source. The background consists of the atmospheric flux with its maximum near the horizon, the cosmic diffuse flux, and neutron induced gamma rays in the telescope. In sharp contrast to the background, a celestial source produces an asymmetric azimuthal response which varies predictably in time as the source moves through the telescope's aperture. This contrast serves as the basis of the cross correlation technique. Continuous data of 47.5 hr were obtained during a balloon flight from Palestine, TX from 0930 UT on September 30, 1978 to 2300 UT on October 1, 1978. The Crab Nebula-Anticenter region was observed on two consecutive days. A number of other medium energy source candidates also crossed the field-of-view. The obtained results are discussed.

Long, J.

Transient vibration test criteria for spacecraft hardware

Transient vibration test criteria, developed for spacecraft hardware, provide a test rationale to verify the capability of the hardware to withstand the low and mid frequency transient vibration environments induced by launch vehicle events. A test method, consisting of a series of discrete frequency, limited cycle, modulated sine wave pulses, was developed to avoid the slow swept sine drawbacks, yet provide a repeatable test that would excite all frequencies. The shape of the waveform is that of the classic response of the mass of a one degree of freedom system when it is base-excited by an exponentially decayed sine wave transient. Criteria were developed to define pulse amplitudes, shapes, and center frequencies from spacecraft loads analyses. Test tolerance criteria were also developed and specified. The transient vibration test criteria were implemented on spacecraft flight hardware and provided a more realistic test simulation (i.e., less conservative) for qualification of spacecraft hardware without risk of undertest.

Kern, D. L.

Optical scanner

An optical scanner for imaging lines in an object plane onto a linear array in a focal plane either continuously or discretely is described. The scanner consists of a set of four mutually perpendicularly oriented plane corner mirrors which provide a reflecting path that describes a parallelogram. In addition, there is a plane parallel scanning mirror with a front and back reflecting surface located midway between the first and fourth corner mirrors. It is oriented so that in the mid-scan position it is parallel to the first corner mirror, and therefore perpendicular to the fourth corner mirror. As the scan mirror rotates, rays incident from a plurality of lines in the object plane are selectively directed through the optical system arriving at a common intersection on the back surface of the scanning mirror where the rays are colinearly directed toward a lens and then imaged onto the linear array in the focal plane. A set of compensating mirrors may be introduced just before the imaging lens to compensate for a small and generally negligible path difference delta sub l between the axial and marginal rays.

Finkel, Mitchell W.

Stable-relaxation-oscillation Nd lasers for long-pulse generation

A simple method to produce high-energy neodymium (Nd) laser pulses with durations on the order of 1 microsec is described. Solid-state lasers can be pumped by relatively short flashlamp pulses to produce well-behaved relaxation oscillations in a diffraction-limited beam. Under the right conditons, each output laser pulse consists of a series of discrete subpulses that are ideally suited to efficient, high-energy amplification. Experimental results for an Nd:LiYF4 oscillator/amplifier system are presented along with numerical simulations. These demonstrate that the system operation is predictable and well behaved. Data are also included for a Nd:glass tunable oscillator based on this concept.

Harrison, James

Assessment of structural analysis technology for elastic shell collapse problems

The prediction of the ultimate load carrying capability for compressively loaded shell structures is a challenging nonlinear analysis problem. Selected areas of finite element technology research and nonlinear solution technology are assessed. Herein, a finite element analysis procedure is applied to four shell collapse problems which have been used by computational structural mechanics researchers in the past. This assessment will focus on a number of different shell element formulations and on different approaches used to account for geometric nonlinearities. The results presented confirm that these aspects of nonlinear shell analysis can have a significant effect on the predicted nonlinear structural response. All analyses were performed using the CSM Testbed software system which allowed a convenient assessment of different element formulations with a consistent approach to solving the discretized nonlinear equations.

Knight, N. F., Jr.

Infrared-Proximity-Sensor Modules For Robot

Collision-avoidance system for articulated robot manipulators uses infrared proximity sensors grouped together in array of sensor modules. Sensor modules, called "sensorCells," distributed processing board-level products for acquiring data from proximity-sensors strategically mounted on robot manipulators. Each sensorCell self-contained and consists of multiple sensing elements, discrete electronics, microcontroller and communications components. Modules connected to central control computer by redundant serial digital communication subsystem including both serial and a multi-drop bus. Detects objects made of various materials at distance of up to 50 cm. For some materials, such as thermal protection system tiles, detection range reduced to approximately 20 cm.

Parton, William

Electron and Proton Auroral Dynamics

Data from the Wide-band Imaging Camera (WIC) sensitive to far ultraviolet auroras and from the Spectrographic Imager (SI) channel SI12, sensitive to proton precipitation induced Lyman alpha were analyzed during a high altitude orbit segment of the IMAGE spacecraft. This segment began during the expansive phase of a substorm. The aurora changed into a double oval configuration, consisting of a set of discrete pole-ward forms and a separate diffuse auroral oval equatorwards, Although IMF Bz was strongly southward considerable activity could be seen poleward of the discrete auroras in the region that was considered to be the polar cap. The SI12 Doppler shifted Lyman alpha signature of precipitating protons show that the proton aurora is on the equatorward side of the diffuse aurora. In the following several hours the IMF Bz field changed signed. Although the general character of the proton and electron aurora did not change, the dayside aurora moved equatorward when the Bz was negative and more bright aurora was seen in the central polar cap during periods of positive Bz.

Mende, S. B.

Assessing Organic Preservation and the Implications for Potential Biosignatures in the Bastide Member of the Séítah Formation, Jezero Crater

Introduction: Olivine has the highest CO2 trapping potential of ultramafic minerals, due to its rapid rate of dissolution and high percentage of divalent cations/unpolymerized silicate tetrahedra [1]. It generates divalent carbonates from CO2 and sequesters CO2 into the mineral matrix of the target lithology. The co-occurrence of olivine and carbonate within abraded targets from the Bastide Member of the Séítah Formation (Fm) of the Jezero crater floor [2] suggests the carbonation of olivine occurred within the mineral matrix hosted in a subsurface system. Hydrothermal origins for the subsurface system are hypothesized from orbital data [3]. This is supported by the detection of hydration features within the rock as well as the carbonate features are solely detected within the abrasion patch but not the rock’s surface. Organic preservation potential of abrasion patches: We incorporated the SHERLOC/WATSON results acquired from the Dourbes, Garde, and Quartier abrasion patches in the Séítah Fm to investigate the organomineral associations, and determine the biosignature preservation potential of these rocks. Dourbes is dominated by olivine and has minor amounts of carbonate, hydrated Ca-sulfate, and amorphous or microcrystalline silicate. Fluorescence features (330-340 nm) are detected in discrete locales and could be consistent with double ring aromatic organic molecules; yet, these features do not appear to be associated with an identified mineral phase. Dark subhedral to euhedral olivine grains within the Garde abrasion patch often co-occur with carbonate-consistent spectral signatures in all analyzed scans. The availability of Fe2+ is a known influence on olivine dissolution rates [1] and SuperCam estimates of the olivine composition (Forsterite-60 average for Sols 202-234) may thereby provide a constraint on the carbonation extent. Carbonated olivine phenocrysts within the matrix may be due to aqueous alteration, as the carbonation of nodules is consistent with observations in other hydrothermal systems and within Martian meteorites (ALH84001) [4]. In comparison, the Quartier abrasion patch represents an extensively altered endmember within the Bastide Member of the Séítah Fm. It contains a fluorescence doublet at 305/325 nm that coexists with multiple species of Na-sulfate and Mg-sulfate, Mg-carbonates, olivine [5]. Aqueous alteration and implications for habitability: The identification of primary and secondary mineral phases observed in the Bastide Member suggests the release of cations from primary ultramafic minerals through aqueous alteration. Within Garde, the carbonate detected appears to be Mg-rich and likely formed from the in situ alteration of Mg-rich olivine within the region as carbonate has only been detected within the rock via in situ analysis. Carbonates derived from abiotic and biotic reactions preserve biosignatures (i.e. indices of habitability) on Earth. The detected carbonate phase found in association with fluorescent features within the rock matrix may also indicate potential organic compounds preserved in a putative hydrothermal system. Fluorescent features (~330 nm) are unique to Garde and Dourbes, though they are co-located to carbonate signatures solely within Garde, and between the light toned minerals. Identification of these fluorescent features may be consistent with 1-2 ring aromatic compounds. The limit of detection for Raman is multiple orders of magnitude greater than the limit of detection required for fluorescence [6]. Implications for provenance: Hydrothermal systems represent disequilibrium chemical conditions that are hypothesized to have supported the emergence of life and also preserve ancient carbon within precipitated carbonates [7]. Hydrothermal system associated carbonates are capable of preserving biosignatures up to an estimated ~3.77 – 4.28 Gya [8]. Thus, carbonated olivine found within the Séítah Fm may represent a high-potential biosignature preserving environment on Mars. The formation of carbonates by an aqueous alteration process, such as carbonation of olivine is also consistent with hypotheses for carbonate within the greater regional-olivine bearing unit [2,3], which contains Garde, Dourbes, and Quartier. Acknowledgments: This work was carried out at the Jet Propulsion Laboratory, The California Institute of Technology under a contract from NASA. References: [1] Wood et al., (2019) ES&T, 6, 10. [2] Stack, K. et al., (2020) Space Sci Rev, 216, 127. [3] Tarnas, J. et al., (2021) JGR: Planets, 126, 11. [4] Steele et al., (2007) Meteorit. Planet. Sci., 42, 9. [5] Murphy, A.E. et al., (2022) LPSC [6] Bhartia et al., (2021) Space Sci Rev, 217, 58. [7] Luther (2021) GRL, 48, e2021GL094869. [8] Dodd et al., (2017) Nature, 543, 60-64

E L Cardarelli

Volatile Abundances, Extended Coma Sources, and Nucleus Ice Associations in Comet C/2014 Q2 (Lovejoy)

High-resolution infrared spectra of comet C/2014 Q2 Lovejoy were acquired with NIRSPEC at the W. M. Keck Observatory on two post-perihelion dates(UT 2015 February 2 and 3).H2O was measured simultaneously with CO,CH3OH, H2CO, CH4, C2H6,C2H4, C2H2, HCN, and NH3 on both dates, and rotational temperatures, production rates, relative abundances, H2O ortho-to-para ratios, and spatial distributions in the coma were determined. The first detection of C2H4in a comet from ground-based observations is reported. Abundances relative to H2O for all species were found to be in the typical range compared with values for other comets in the overall population to date. There is evidence of variability in rotational temperatures and production rates on timescales that are small compared with the rotational period of the comet. Spatial distributions of volatiles in the coma suggest complex outgassing behavior.CH3OH, HCN, C2H6, and CH4 spatial distributions in the coma are consistent with direct release from associated ices in the nucleus and are peaked in a more sunward direction compared with co-measured dust. H2O spatial profiles are clearly distinct from these other four species, likely due to a sizable coma contribution from icy grain sublimation.Spatial distributions for C2H2,H2CO, and NH3suggest substantial contributions from extended coma sources, providing further evidence for distinct origins and associations for these species in comets. CO shows a different spatial distribution compared with other volatiles, consistent with jet activity from discrete nucleus ice sources.

Neil Dello Russo

Goal-Oriented Discretization Error Control in Coupled Nearfield-Farfield Low-Boom Simulations

The method of adjoint weighted residuals is used to determine the level of discretization error in loudness predictions of sonic booms on the ground. We analyze the standard nearfied-farfield domain decomposition approach. In the nearfield domain, the three-dimensional Euler equations are solved to obtain a pressure signature generated by the aircraft. In the farfield, this waveform is propagated through the atmosphere to the ground by solving the augmented Burgers’ equation. Loudness is characterized using weighted sound-exposure-level metrics. We formulate discretization error estimates for the ground signature and loudness metrics for this one-way coupled system. Although the nearfield solution is independent of the farfield, the adjoint formulation for the coupled system provides feedback from the farfield to identify high-error regions in the nearfield. The results demonstrate that the discrete adjoint implementation is asymptotically consistent and provides reliable error estimates. Furthermore, we show how the error can be controlled through adaptive refinement of the nearfield mesh. The approach is evaluated on two- and three-dimensional problems, including the X-59 flight demonstration aircraft.

CST

A Multi-Disciplinary Analysis Framework for the Design of Small Launch Vehicles

Between the years of 1995 and 2014 the number of small satellites (1-500 kg) went up from 20 to 180. [1] Out of the 180 launched in 2014 66% were Nano satellites (1-10 kg). [1] With this trend of smaller satellites, one would expect a rise in number of small launch vehicles (SLVs are defined by capability to carry 1-100 kg to orbit), but this has not happened: [3] only 8% of all small satellites are launched on SLVs and the others become secondary payloads on regular launch vehicles. [1] This results in small satellites being placed into either suboptimal orbits or waiting for launch schedules to align with bigger launches, resulting in long waiting times. Dedicated SLVs could improve the responsiveness of small satellite launches, but SLV design is complicated by the large architecture space needed to be explored in order to find efficient and affordable designs. The SLV architecture trade space contains many discrete options, e.g. solid vs. liquid fuels, air launch vs. ground launch, number of stages, etc. [2][3] Performing detailed analysis for all the options at once would be prohibitive. Thus, a sizing environment/framework that is capable of providing necessary information for conceptual-level trade studies and can rapidly explore the vast SLV architecture space is necessary. The framework, illustrated in Figure 1, consists of four disciplines central to the sizing and synthesis of launch vehicles: propulsion, aerodynamics, structures, and trajectory. For the aerodynamics and trajectory disciplines, Missile DATCOM and POST2 are used, respectively. The propulsion and structures disciplines are represented in the framework with tools developed at ASDL Georgia Tech. For propulsion, the Solid Motor Analysis Code (SMAC) is a physics-based conceptual design tool for solid rocket motors. SMAC is capable of geometric burn simulation, ballistic analysis, and prediction of thrust performance. [4] For structures, Launch Vehicle Structural Analysis (LVSA) tool is a physics-based tool that focuses on structural dynamic analysis with sizing capability. These tools are integrated into the framework illustrated in Figure 1 with the corresponding connections described in Table 1. The process flow is as follows: first, SMAC sizes insulation and calculates maximum operating pressures for each vehicle stage. The MEOP and insulation thickness values from SMAC are fed to LVSA which then utilizes this information to size the motor casing. This creates a feedback loop between LVSA and SMAC that converges on the radius available for fuel, casing, and insulation thickness. Once the stage sizing is converged on SMAC creates an engine deck that is passed to POST2. Next, the data from SMAC and LVSA goes into Missile DATCOM which generates an aerodynamics database for POST2. Finally, POST2 performs a targeting optimization while maximizing the payload mass to orbit. Within the framework, POST2 is automated in order to be robust to a wide variety of possible designs by performing a Monte Carlo simulation over the initialization vector of the POST2 optimization variables. To demonstrate the capability of this framework, a sample problem of exploring the design space of an SLV capable of placing satellites into a low Earth orbit (inclination=47 deg, 196.5 by 211.3 nm) is used. This sample problem is a ground-launched SLV consisting of four in-line SRM stages. The multidisciplinary design analysis (MDA) environment is utilized to explore a design space consisting of 22 continuous and 12 discrete variables, shown by Table 1 by blocks 1,2,3. Running a full factorial design of experiments (DOE) would have been prohibitively expensive even with this reduced design space, thus a space filling design with 3,502 and then additional expansion of 2,602 cases was used. The first DOE consisted of 3,502 cases, and all of the variables were varied. These input variables are listed in blocks 1, 2, and 3 in Table 1. Most of the variables are propulsion related with stage lengths determining the delta-V split of the SLV. The expansion consisted of 2,602 cases, and the continuous variables were set to be equal to the most promising designs from the sizing of first DOE. For each set of continuous variables, the discrete variables (grain type, star points, and propellant) were to varied. The results of the DOE can be seen in Figure 2; each of the points in this plot represents a closed launch vehicle that reaches the targeted orbit. For each of the cases, there is information on flown trajectory, structural, and propulsion properties of the SLV. For example, Figure 3 shows changes in altitude and velocity with time for a particular case. The right side of Figure 3 clearly shows the coasting (slow decrease) and burning phases (sharp increase) of the SLV mission. LV mass is positively correlated with the optimized payload mass to orbit because heavier LVs carry more fuel and thus have more stored chemical energy. Furthermore, for any given payload mass to orbit, the most efficient design would result in the smallest LV. The results as visualized in Figure 2 shows this tradeoff, and the Pareto frontier of the efficient designs can be seen along the dotted line. Figure 2 can be divided into regions with the lowest mass vehicles corresponding to the minimum bound on radius, and the highest mass vehicles corresponding to the maximum bound on radius. Within a mass region, the discrete variables, such as propellant type and propellant grain arrangement, have the most effect on payload mass. This paper presents an MDA framework that can perform an automated physics-based sizing of SLV designs and a corresponding methodology to utilize the MDA to explore the design space of SLVs. The proposed methodology was applied to a perform a design space exploration for a four stage SLV. The outputs show the expected pareto frontier forming and provide detailed information about the SLV performance and staging. Using the produced data, it will be possible to select a set of pareto optimal designs that can then be further explored in subsequent design cycles. This demonstration shows that automated design space exploration should be used in the early phase design of future SLV concepts.

Nikita S Birbasov