Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “target variability”

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 649 records · Page 36

Control-oriented core-SOL-divertor model to address integrated burn and divertor control challenges in ITER

The real-time regulation of a burning plasma’s temperature and density, or burn control, will be necessary to produce high fusion power in future tokamaks like ITER. This is made more challenging due to the plasma’s nonlinear characteristics and the interdependence between the core-plasma and edge-plasma regions. For example, a raising plasma temperature leads to increasing reactivity and therefore to more alpha-particle heating, which further increases temperature. Furthermore, a raise of the fusion power increases the heat flow through the scrape-off-layer (SOL), which can compromise the integrity of the divertor without proper safeguards. For control design, a model-based approach is attractive because it can directly incorporate the nonlinear, coupled, burning-plasma dynamics into the design. To facilitate this design approach, a control-oriented core-SOL-divertor (CSD) model is presented in this work. In this CSD model, a core-plasma model captures the nonlinear dynamics of the core’s density and temperature, and a SOL-divertor model defines the plasma conditions at the separatrix and divertor including the heat load on the target plates. The core-plasma and SOL-divertor models are coupled through the exchange of various variables. In particular, the SOL-divertor model yields the separatrix temperature and the influx of recycled particles into the core-plasma. Further, these variables influence the power and particle balances captured by the core-plasma model. In return, the core-plasma model determines the intensity of the heat and particles fluxes across the separatrix, and this outflow strongly impacts the SOL-divertor model. Therefore, the power and density of the core-plasma, which can be readily modulated through external heating systems and pellet injection, can be viewed as control knobs for the SOL-divertor region in addition to the gas puffing. In simulations of the CSD model, it is demonstrated how external actuation can be utilized to meet burn control and divertor control objectives simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Distribution and diversity of enzymes for polysaccharide degradation in fungi

Fungi are important polysaccharide degraders in the environment and for biotechnology. Here, the increasing number of sequenced fungal genomes allowed for systematic identification of genes and proteins involved in polysaccharide degradation in 218 fungi. Globally, 9,003 sequences for glycoside hydrolases and lytic polysaccharide mono-oxygenases targeting cellulose, xylan, and chitin, were identified. Although abundant in most lineages, the distribution of these enzymes is variable even between organisms from the same genus. However, most fungi are generalists possessing several enzymes for polysaccharide deconstruction. Most identified enzymes were small proteins with simple domain organization or eventually consisted of one catalytic domain associated with a non-catalytic accessory domain. Thus unlike bacteria, fungi's ability to degrade polysaccharides relies on apparent redundancy in functional traits and the high frequency of lytic polysaccharide mono-oxygenases, as well as other physiological adaptation such as hyphal growth. Globally, this study provides a comprehensive framework to further identify enzymes for polysaccharide deconstruction in fungal genomes and will help identify new strains and enzymes with potential for biotechnological application.

59 BASIC BIOLOGICAL SCIENCES↗

Marginal Cost Pricing in a World without Perfect Competition: Implications for Electricity Markets with High Shares of Low Marginal Cost Resources

A common approach to regulating electricity is through auction-based competitive wholesale markets. The goal of this approach is to provide a reliable supply of power at the lowest reasonable cost to the consumer. This necessitates market structures and operating rules that ensure revenue sufficiency for all generators needed for resource adequacy purposes. Wholesale electricity markets employ marginal-cost pricing to provide cost-effective dispatch such that resources are compensated for their operational costs. However, marginal-cost pricing alone cannot guarantee cost recovery outside of perfect competition, and electricity markets have at least six attributes that preclude them from functioning as perfectly competitive markets. These attributes include market power, externalities, public good attributes, lack of storage, wholesale price caps, and ineffective demand curve. Until (and unless) these failures are ameliorated, some form of corrective action(s) will be necessary to improve market efficiency so that prices can correctly reflect the needed level of system reliability. Many of these options necessarily involve some form of administrative or out-of-market actions, such as scarcity pricing, capacity payments, bilateral or other out-of-market contracts, or some hybrid combination. A key focus with these options is to create a connection between the electricity market and long-term reliability/loss-of-load expectation targets, which are inherently disconnected in the native markets because of the aforementioned market failures. The addition of variable generation resources can exacerbate revenue sufficiency and resource adequacy concerns caused by these underlying market failures. Because variable generation resources have near-zero marginal costs, they effectively suppress energy prices and reduce the capacity factors of conventional generators through the merit-order effect in the simplest case of a convex market; non-convexities can also suppress prices.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Final Report: Optimal Model Complexity in Geological Carbon Sequestration: A Response Surface Uncertainty Analysis

The critical component of a risk assessment study in evaluating GCS is an analysis of uncertainty in CO2 modeling. In such analyses, direct numerical simulation of CO2 flow and leakage requires many time-consuming model runs. Alternatively, analytical methods have been developed which allow fast and efficient estimation of CO2 storage and leakage, although restrictive assumptions on formation rock and fluid properties are employed. In this study, an intermediate approach is proposed based on the Design of Experiment and Response Surface methodology, which consists of using a limited number of numerical simulations to estimate a prediction outcome as a combination of the most influential uncertain site properties. The methodology can be implemented within a Monte Carlo framework to efficiently assess parameter and prediction uncertainty while honoring the accuracy of numerical simulations. The choice of the uncertain properties is flexible and can include geologic parameters that influence reservoir heterogeneity, engineering parameters that influence gas trapping and migration, and reactive parameters that influence the extent of fluid/rock reactions. The method was tested and verified on modeling long-term CO2 flow, non-isothermal heat transport, and CO2 dissolution storage by coupling two-phase flow with explicit miscibility calculation using an accurate equation of state that gives rise to convective mixing of formation brine variably saturated with CO2. All simulations were performed using three-dimensional high-resolution models including a target deep saline aquifer, overlying caprock, and a shallow aquifer. To evaluate the uncertainty in representing reservoir permeability, sediment hierarchy of a heterogeneous digital stratigraphy was mapped to create multiple irregularly shape stratigraphic models of decreasing geologic resolutions: heterogeneous (reference), lithofacies, depositional environment, and a (homogeneous) geologic formation. To ensure model equivalency, all the stratigraphic models were successfully upscaled from the reference heterogeneous model for bulk flow and transport predictions (Zhang & Zhang, 2015). GCS simulation was then simulated with all models, yielding insights into the level of parameterization complexity that is needed for the accurate simulation of reservoir pore pressure, CO2 storage, leakage, footprint, and dissolution over both short (i.e., injection) and longer (monitoring) time scales. Important uncertainty parameters that impact these key performance metrics were identified for the stratigraphic models as well as for the heterogeneous model, leading to the development of reduced/simplified models at lower characterization cost that can be used for the reservoir uncertainty analysis. All the CO2 modeling was conducted using PFLOTRAN – a massively parallel, multiphase, multi-component, and reactive transport simulator developed by a multi-laboratory DOE/SciDAC (Scientific Discovery through Advanced Computing) project (Zhang et al., 2017, in review). Within the uncertainty analysis framework, increasing reservoir depth were investigated to explore its effect on the uncertainty outcomes and the potential for developing gravity-stable injection with increased storage security (Dai et al., 20126; Dai et al., 2017, in review). Finally, to accurately model CO2 fluid-rock reactions and resulting long-term storage as secondary carbonate minerals, a modified kinetic rate law for general mineral dissolution and precipitation was proposed and verified that is invariant to a scale transformation of the mineral formula weight. This new formulation will lead to more accurate assessment of mineral storage over geologic time scales (Lichtner, 2016).

58 GEOSCIENCES↗

Specification of EBR-II Gamma Scan Data from the Precision Gamma Scanner

The gamma scan data stored in the FIPD was originally measured using the Precision Gamma Scanner (PGS) at the Hot Fuel Examination Facility (HFEF) at Idaho National Laboratory (INL). This system was originally manufactured in 1971 (when the facilities, now at INL, were operated by Argonne National Laboratory) and was upgraded on several occasions. It consists of a germanium gamma ray detector (the exact nature of which changed over time as technology improved), a variable-slit collimator that blocks the view of gamma emissions from all but the selected area of the target, and a computer-controlled system for translating a sample past the slit and logging the data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

JVLA Observations of Young Brown Dwarfs

We present sensitive 3.0 cm JVLA radio continuum observations of six regions of low-mass star formation that include twelve young brown dwarfs (BDs) and four young BD candidates. We detect a total of 49 compact radio sources in the fields observed, of which 24 have no reported counterparts and are considered new detections. Twelve of the radio sources show variability in timescales of weeks to months, suggesting gyrosynchrotron emission produced in active magnetospheres. Only one of the target BDs, FU Tau A, was detected. However, we detected radio emission associated with two of the BD candidates, WL 20S and CHLT 2. The radio flux densities of the sources associated with these BD candidates are more than an order of magnitude larger than expected for a BD and suggest a revision of their classification. In contrast, FU Tau A falls on the well-known correlation between radio luminosity and bolometric luminosity, suggesting that the emission comes from a thermal jet and that this BD seems to be forming as a scaled-down version of low-mass stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Heterodyne efficiency for a coherent laser radar with diffuse or aerosol targets

The performance of a Coherent Laser Radar is determined by the statistics of the coherent Doppler signal. The heterodyne efficiency is an excellent indication of performance because it is an absolute measure of beam alignment and is independent of the transmitter power, the target backscatter coefficient, the atmospheric attenuation, and the detector quantum efficiency and gain. The theoretical calculation of heterodyne efficiency for an optimal monostatic lidar with a circular aperture and Gaussian transmit laser is presented including beam misalignment in the far-field and near-field regimes. The statistical behavior of estimates of the heterodyne efficiency using a calibration hard target are considered. For space based applications, a biased estimate of heterodyne efficiency is proposed that removes the variability due to the random surface return but retains the sensitivity to misalignment. Physical insight is provided by simulation of the fields on the detector surface. The required detector calibration is also discussed.

Frehlich, R. G.↗

Epitaxial niobium dioxide thin films by reactive-biased target ion beam deposition

Epitaxial NbO{sub 2} thin films were synthesized on Al{sub 2}O{sub 3} (0001) substrates via reactive bias target ion beam deposition. X-ray diffraction and Raman spectra were used to confirm the tetragonal phase of pure NbO{sub 2}. Through XPS, it was found that there was a ∼1.3 nm thick Nb{sub 2}O{sub 5} layer on the surface and the bulk of the thin film was NbO{sub 2}. The epitaxial relationship between the NbO{sub 2} film and the substrate was determined. Electrical transport measurement was measured up to 400 K, and the conduction mechanism was discussed.

36 MATERIALS SCIENCE↗

First Measurement of Differential Charged Current Quasielastic-like $\nu_\mu$–Argon Scattering Cross Sections in Kinematic Imbalance Variables with the MicroBooNE Detector

We report the first measurement of flux-integrated multi-differential cross sections for charged-current events with muon neutrinos scattering on argon with solely a muon and a single proton in the final state as a function of kinematic imbalance variables. The measurement was carried out using the Booster Neutrino Beam at Fermi National Accelerator Laboratory within the MicroBooNE Liquid Argon Time Projection Chamber detector with an exposure of 6.79 × 10 20 protons on target. Events were selected to enhance the contribution of charged-current mesonless interactions with one proton detected in the final state. The data discussed here are reported in terms of multi-differential cross sections in kinematic imbalance variables, which are generally sensitive to nuclear effects. The double-differential results in these variables can provide an excellent handle to disentagle specific nuclear aspects not easily isolated via single differential cross sections. Our results pave a path towards identifying regions of the phase-space where future interaction modeling development and Monte Carlo neutrino generator tuning efforts should concentrate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The effect of irradiation and reflectance variability on vegetation condition assessment

Surface reflectance and global irradiance variations across an imaged area determine the minimum signal difference from groups of pixels in different targets needed to discriminate them. The factors affecting such variations and their magnitudes are considered. The sensor output also depends upon the interaction of the spectral response of the sensor and the spectral upwelling radiance from the target. From a consideration of the spectral instrument responses of detectors in the Landsat multispectral scanners (MS), together with measured reflectance factors of a vegetative canopy stressed to different levels of severity, the feasibility of mapping and quantifying disease stress witlh Landsat vegetative indices (VINs) is demonstrated.

Duggin, M. J.↗

Leveraging knowledge engineering and machine learning for microbial bio-manufacturing

Genome scale modeling (GSM) predicts the performance of microbial workhorses and helps identify beneficial gene targets. GSM integrated with intracellular flux dynamics, omics, and thermodynamics have shown remarkable progress in both elucidating complex cellular phenomena and computational strain design (CSD). Nonetheless, these models still show high uncertainty due to a poor understanding of innate pathway regulations, metabolic burdens, and other factors (such as stress tolerance and metabolite channeling). Besides, the engineered hosts may have genetic mutations or non-genetic variations in bioreactor conditions and thus CSD rarely foresees fermentation rate and titer. Metabolic models play important role in design-build-test-learn cycles for strain improvement, and machine learning (ML) may provide a viable complementary approach for driving strain design and deciphering cellular processes. In order to develop quality ML models, knowledge engineering leverages and standardizes the wealth of information in literature (e.g., genomic/phenomic data, synthetic biology strategies, and bioprocess variables). Data driven frameworks can offer new constraints for mechanistic models to describe cellular regulations, to design pathways, to search gene targets, and to estimate fermentation titer/rate/yield under specified growth conditions (e.g., mixing, nutrients, and O 2 ). Furthermore this review highlights the scope of information collections, database constructions, and machine learning techniques (such as deep learning and transfer learning), which may facilitate "Learn and Design" for strain development.

59 BASIC BIOLOGICAL SCIENCES↗

HIV-1 Neutralizing Antibody Signatures and Application to Epitope-Targeted Vaccine Design

Eliciting HIV-1-specific broadly neutralizing antibodies (bNAbs) remains a challenge for vaccine development, and the potential of passively delivered bNAbs for prophylaxis and therapeutics is being explored. We used neutralization data from four large virus panels to comprehensively map viral signatures associated with bNAb sensitivity, including amino acids, hypervariable region characteristics, and clade effects across four different classes of bNAbs. The bNAb signatures defined for the variable loop 2 (V2) epitope region of HIV-1 Env were then employed to inform immunogen design in a proof-of-concept exploration of signature-based epitope targeted (SET) vaccines. V2 bNAb signature-guided mutations were introduced into Env 459C to create a trivalent vaccine, and immunization of guinea pigs with V2- SET vaccines resulted in increased breadth of NAb responses compared with Env 459C alone. These data demonstrate that bNAb signatures can be utilized to engineer HIV-1 Env vaccine immunogens capable of eliciting antibody responses with greater neutralization breadth.

60 APPLIED LIFE SCIENCES↗

Sensor-induced temporal variability of Landsat MSS data

Landsat-1 and Landsat-2 multispectral scanner (MSS) data were studied to determine the consistency of the calculated reflectance values over time. Data from six spectrally stable targets were collected over a 3-year period (1975-1977). Reflectance values calculated from the digital numbers were regressed against time to note any long term changes. Results indicate that, over a 1000-day period beginning 1 January 1975, MSS 2 reflectances were stable. MSS 1 reflectances over that same period decreased approximately 25-32 percent. The observed decreases are most likely due to a deterioration of the transmissive qualities of the Landsat-1 scanner's internal optical path. Although MSS temporal variation may be unique to MSS 1, time (days in orbit) should be considered as an independent variable when intersatellite calibration equations are computed. A case is made for a stable, monitored calibration system which would permit the calculation of true top-of-the-atmosphere reflectance measures.

Nelson, R. F.↗

SU-E-T-230: Creating a Large Number of Focused Beams with Variable Patient Head Tilt to Improve Dose Fall-Off for Brain Radiosurgery

Purpose: To develop a treatment delivery and planning strategy by increasing the number of beams to minimize dose to brain tissue surrounding a target, while maximizing dose coverage to the target. Methods: We analyzed 14 different treatment plans via Leksell PFX and 4C. For standardization, single tumor cases were chosen. Original treatment plans were compared with two optimized plans. The number of beams was increased in treatment plans by varying tilt angles of the patient head, while maintaining original isocenter and the beam positions in the x-, y- and z-axes, collimator size, and beam blocking. PFX optimized plans increased beam numbers with three pre-set tilt angles, 70, 90, 110, and 4C optimized plans increased beam numbers with tilt angles increasing arbitrarily from range of 30 to 150 degrees. Optimized treatment plans were compared dosimetrically with original treatment plans. Results: Comparing total normal tissue isodose volumes between original and optimized plans, the low-level percentage isodose volumes decreased in all plans. Despite the addition of multiple beams up to a factor of 25, beam-on times for 1 tilt angle versus 3 or more tilt angles were comparable (<1 min.). In 64% (9/14) of the studied cases, the volume percentage decrease by >5%, with the highest value reaching 19%. The addition of more tilt angles correlates to a greater decrease in normal brain irradiated volume. Selectivity and coverage for original and optimized plans remained comparable. Conclusion: Adding large number of additional focused beams with variable patient head tilt shows improvement for dose fall-off for brain radiosurgery. The study demonstrates technical feasibility of adding beams to decrease target volume.

60 APPLIED LIFE SCIENCES↗

X-Ray Radar Imaging Technique Using a 2 Mev Linear Electron Accelerator

X-ray radar imaging, patented in 2013 by James R. Wood, combines standard radar techniques with the penetration power of x-rays to image scenes. Our project strives to demonstrate the technique using a 2 MeV linear electron accelerator (linac) to generate the S-band–modulated x-ray signals. X-ray detectors such as photodiodes and scintillators are used to detect the signals in backscatter and transmission detection schemes. The S-band microstructure is imposed on the variable-width electron pulse, and this modulation carries over to the bremsstrahlung x-rays after the electron beam is incident upon a copper-tungsten alloy target. Using phase/distance calculations and a low-jitter timing system, we expect to detect different object distances by comparing the measured phase differences. The experimental setup, which meets strict jitter requirements, and preliminary experimental results are presented.

43 PARTICLE ACCELERATORS↗

PERIODIC AND APERIODIC VARIABILITY IN THE MOLECULAR CLOUD ρ OPHIUCHUS

Presented are the results of a near-IR photometric survey of 1678 stars in the direction of the ρ Ophiuchus (ρ Oph) star forming region using data from the 2MASS Calibration Database. For each target in this sample, up to 1584 individual J-, H-, and K{sub s} -band photometric measurements with a cadence of ∼1 day are obtained over three observing seasons spanning ∼2.5 yr; it is the most intensive survey of stars in this region to date. This survey identifies 101 variable stars with ΔK{sub s} -band amplitudes from 0.044 to 2.31 mag and Δ(J – K{sub s} ) color amplitudes ranging from 0.053 to 1.47 mag. Of the 72 young ρ Oph star cluster members included in this survey, 79% are variable; in addition, 22 variable stars are identified as candidate members. Based on the temporal behavior of the K{sub s} time-series, the variability is distinguished as either periodic, long time-scale or irregular. This temporal behavior coupled with the behavior of stellar colors is used to assign a dominant variability mechanism. A new period-searching algorithm finds periodic signals in 32 variable stars with periods between 0.49 to 92 days. The chief mechanism driving the periodic variability for 18 stars is rotational modulation of cool starspots while 3 periodically vary due to accretion-induced hot spots. The time-series for six variable stars contains discrete periodic ''eclipse-like'' features with periods ranging from 3 to 8 days. These features may be asymmetries in the circumstellar disk, potentially sustained or driven by a proto-planet at or near the co-rotation radius. Aperiodic, long time-scale variations in stellar flux are identified in the time-series for 31 variable stars with time-scales ranging from 64 to 790 days. The chief mechanism driving long time-scale variability is variable extinction or mass accretion rates. The majority of the variable stars (40) exhibit sporadic, aperiodic variability over no discernable time-scale. No chief variability mechanism could be identified for these variable stars.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗