Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Regularization”

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 667 records · Page 37

Charge Transfer Screening and Energy Level Alignment at Complex Organic–Inorganic Interfaces: A Tractable Ab Initio GW Approach

Complex organic–inorganic interfaces are important for device and sensing applications. Charge transfer doping is prevalent in such applications and can affect the interfacial energy level alignments (ELA), which are determined by many-body interactions. We develop an approximate ab initio many-body GW approach that can capture many-body interactions due to interfacial charge transfer. The approach uses significantly less resources than a regular GW calculation but gives excellent agreement with benchmark GW calculations on an F4TCNQ/graphene interface. We find that many-body interactions due to charge transfer screening result in gate-tunable F4TCNQ HOMO–LUMO gaps. We further predict the ELA of a large system of experimental interest—4,4'-bis(dimethylamino)bipyridine (DMAP-OED) on monolayer MoS 2 , where charge transfer screening results in an ~1 eV reduction of the molecular HOMO–LUMO gap. Comparison with a two-dimensional electron gas model reveals the importance of explicitly considering the intraband transitions in determining the charge transfer screening in organic–inorganic interface systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser calibration of the ATLAS Tile Calorimeter during LHC Run 2

This article reports the laser calibration of the hadronic Tile Calorimeter of the ATLAS experiment in the LHC Run 2 data campaign. The upgraded Laser II calibration system is described. The system was commissioned during the first LHC Long Shutdown, exhibiting a stability better than 0.8% for the laser light monitoring. The methods employed to derive the detector calibration factors with data from the laser calibration runs are also detailed. These allowed to correct for the response fluctuations of the 9852 photomultiplier tubes of the Tile Calorimeter with a total uncertainty of 0.5% plus a luminosity-dependent sub-dominant term. Finally, we report the regular monitoring and performance studies using laser events in both standalone runs and during proton collisions. These studies include channel timing and quality inspection, and photomultiplier linearity and response dependence on anode current.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Microbial Anomalies Encountered on the International Space Station

Microorganisms in our living environments are unavoidable. A community of microbes arrived in space with the delivery of the first element of the International Space Station (ISS), attached to hardware and on the bodies of the humans tasked with the initial assembly missions. The risk that microorganisms could cause adverse effects in the health of both the human occupants of the ISS as well as the physical integrity of the station environment and life support systems has been both a driver and a function of engineering and operational controls. Scientists and engineers at NASA have gone to extensive measures to control microbial growth at levels safe for the crewmembers and the spacecraft environment. Many of these measures were initiated with the design of the spacecraft and its systems. Materials used in the ISS were tested for resistance to fungi, such as mold and a paint with a fungus-killing chemical was also used. Controlling the humidity of the air in the Station is also an effective way of discouraging microbe growth. The breathing air is reconditioned by the Environmental Control Life Support System (ECLSS) prior to distribution, utilizing High Efficiency Particulate Air (HEPA) filtration. Requirements restricting the accumulation of water condensate in the air handlers and habitable volume of the ISS were other safeguards added. Water for drinking and food rehydration is disinfected or filtered. A robust in-flight housekeeping regimen for the ISS significantly reduces inappropriate growth of microorganisms and includes a regular cleaning of accessible surfaces with disinfectant wipes. Most of these requirements were suggested by microbiologists to mitigate and possibly prevent many microbiological risks. In addition to these controls, before flight monitoring and analyses of the cabin air, exposed surfaces, water and food, consumables, and crew members are conducted to mitigate microbial risk to the crew and spacecraft. Many microbial risks are much easier to identify and resolve before launch than during space flight. Although the focus has been on prevention of microbiologically related, not all problems can be anticipated. A number of microbial anomalies have occurred on ISS. This paper will discuss the occurrences, root-cause investigations, and mitigation steps taken to remediate the contamination.

Bruce, Rebekah J.↗

Structure of Diisobutylene Maleic Acid Copolymer (DIBMA) and Its Lipid Particle as a “Stealth” Membrane-Mimetic for Membrane Protein Research

The study of membrane proteins remains challenging, especially in a native membrane environment. Recently, major progress has been made using maleic acid copolymers, such as styrene maleic acid, to purify membrane proteins and study them directly with native lipids associated with the membrane. Additional maleic acid copolymers, such as diisobutylene maleic acid (DIBMA) membrane-mimetic systems, are being developed and found to have improved spectroscopic properties and pH stability. We studied DIBMA and its lipid particles in solution to better understand its assembly, without and with the lipids, to provide an insight regarding how to use it in solution for better membrane extraction. Using small-angle neutron and X-ray scattering (SANS/SAXS), we show that DIBMA organizes into structures of different size scales at various concentrations and ionic strengths. The polymer performed reasonably well under most solvent conditions except in very low concentrations and high-salt conditions that could result in limited interaction with lipids. To explore DIBMA lipid particles as a suitable membrane-mimetic system for neutron scattering studies of membrane proteins, we measured and determined the contrast-matching point of DIBMA to be ~12% (v/v) D 2 O — similar to that of most protiated lipid molecules but distinct from that of regular protiated proteins — providing a natural contrast for separating their neutron scattering signals. Using SANS contrast variation, we demonstrated that the scattering from the whole lipid particle can be annihilated. Further, we determined that a well-defined lipid nanodisc structure with DIBMA was contrast-matched. These results demonstrate that the DIBMA lipid particle is an outstanding “stealth” membrane-mimetic for membrane proteins. The results provide a structural framework for understanding the organization and assembly process of the polymer itself and the lipid molecules. Such an understanding is imperative for structural techniques such as cryo-electron microscopy, nuclear magnetic resonance, small-angle scattering, and other biophysical techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetic moments of A = 3 nuclei obtained from chiral effective field theory operators

Chiral effective field theory (χEFT) provides a framework for obtaining internucleon interactions in a systematically improvable fashion from first principles, while also providing for the derivation of consistent electroweak current operators. In this study, we apply consistently derived interactions and currents towards calculating the magnetic dipole moments of the A = 3 systems 3 H and 3 He. We focus here on LENPIC interactions obtained using semilocal coordinate-space (SCS) regularization. Starting from the momentum-space representation of the LENPIC χEFT vector current, we derive the SCS-regularized magnetic dipole operator up through next-to-next-to-leading order (N 2 LO). We then carry out no-core shell-model calculations for 3 H and 3 He systems using the SCS LENPIC interaction at N 2 LO in χEFT and evaluate the magnetic dipole moments obtained using the consistently derived one-nucleon and two-nucleon electromagnetic currents. As anticipated by prior results with χEFT currents, the current corrections through N 2 LO provide improved, but not yet complete, agreement with experiment for the 3 H and 3 He magnetic dipole moments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NGC 6845 Revisited: Stellar Kinematics on the Early Type Member Galaxies

We observed the compact group of galaxies NGC 6845 with GMOS-S imaging and long-slit spectroscopy. The group itself is comprised of two spiral galaxies and two S0 galaxies, and shows signs of gravitational interaction among them. While the spiral galaxies have been already well studied, the two S0 galaxies have received little attention. We investigate the morphology and stellar kinematics of the two S0 galaxies in detail in order to derive the value for their mass and subsequently incorporate them into the mass model for the group. It was found that the internal kinematics of the largest S0 (NGC 6845 C) is mostly regular, while the smallest S0 (NGC 6845 D) shows some hint of interaction, most likely with NGC 6845 C and NGC 6845 A. The new determination of the dynamical mass of the group implies that the group mass has to be 10 times the combined mass of the individual members according to their internal kinematics. The results point to an early stage of the group interaction and also to the possibility that it has still not reached virial equilibrium. From the galaxy evolution viewpoint, this group represents an example of the case of the transformation of Spirals into S0 galaxies in dense environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Study the Vertical Structure and Transportation of the Extreme African Dust Storms using MERRA-2 Data

The Modern-Era Retrospective Analysis for Research and Application, Version 2 (MERRA-2) provides the first long-term global reanalysis to assimilate space-based observations of aerosols and represents their interactions with other physical processes in the climate system. In this study, we have examined the variations of atmospheric aerosols for the last 20 years since 2002 using the sub-daily MERRA-2 data and found seven extreme African dust storms that were transported westward across the Atlantic Ocean from the Sahara, crossing 70o-80oW. In particular, the well-known ‘Godzilla’ dust storm occurred in June 2020, and its dust cloud, with the highest-on-record aerosol optical depths, was transported toward the Americas. This storm greatly degraded air quality over large areas of the Caribbean Basin and the United States. The air quality index reached unhealthy levels for sensitive groups in more than ten U.S. states. In our study, the vertical structure and transport characteristics of the dust layers during this extreme dust event are investigated. The geopotential height and temperature were found anomalously low (around 600 hPa) over the Atlantic Ocean off northwest Africa before the June 2020 dust storm. This anomalous circulation pattern was persistent for more than 12 days starting from around May 30, breaking the regular easterly waves that transport dust from the Sahara Desert to the west. To verify the data quality, daily MERRA-2 PM2.5 data were calculated and compared with PM2.5 observations from the U.S. Environmental Protection Agency (EPA) at several selected ground stations in Florida. We provide this case study to illustrate how to effectively use various MERRA-2 data services at Goddard Earth Sciences Data and Information Services Center (GES DISC) where MERRA-2 data are archived, hoping to help data users in exploring their own topics of interest using data services at GES DISC.

data management, reanalysis↗

Rapid data-driven model reduction of nonlinear dynamical systems including chemical reaction networks using ℓ 1 -regularization

We develop a new data-driven paradigm for efficient model reduction of a broad class of nonlinear dynamical systems. Our model reduction method directly enables the interpretation of key components of the dynamical system, unlike traditional projection-based model reduction methods that focus on reducing computational complexity more than interpretability. Our method is not application specific and is simple to implement on nonlinear dynamical systems arising from a variety of different fields. It requires minimal parameterization using a single parameter to trade-off between model complexity and estimation error. We use a data-driven paradigm to formulate model reduction as an efficient convex optimization problem that scales polynomially in the original size of the complex system, enabling systems with as many as thousands of components to be reduced in a matter of minutes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quasinormal modes for non-minimally coupled scalar fields in regular black hole spacetimes: Grey-body factors, area spectrum and shadow radius

Black hole perturbation theory is a useful approach to study interactions between black holes and fundamental fields. A particular class of black hole solutions arising out of modification of Einstein’s general theory of relativity are regular black holes (RBHs) which can be constructed using a nonlinear electrodynamic Lagrangian. Because of their importance, we are interested in studying the behavior of three kinds of such RBHs under perturbations generated by an external field. Indeed, we investigate the quasinormal modes (QNMs) of a massive scalar field propagating near the RBHs which is non-minimally coupled to the Ricci scalar-tensor of background geometry. We will attempt to find the low-lying quasinormal frequencies of the perturbations by using WKB approximation. We shall also study the relationship between the QNMs and some characteristic properties of black holes such as grey-body factors, area quantization, and shadow radius.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Easy-plane anisotropic-exchange magnets on a honeycomb lattice: Quantum effects and dealing with them

We provide analytical and numerical insights into the phase diagram and other properties of the extended Kitaev-Heisenberg model on the honeycomb lattice in the easy-plane limit, in which interactions are only between spin components that belong to the plane of magnetic ions. This parameter subspace allows for a much-needed systematic quantitative investigation of spin excitations in the ordered phases and of their generic features. Specifically, we demonstrate that in this limit one can consistently take into account magnon interactions in both zero-field zigzag and field-polarized phases. For the nominally polarized phase, we propose a regularization of the unphysical divergences that occur at the critical field and are plaguing the 1/S approximation in this class of models. For the explored parameter subspace, all symmetry-allowed terms of the standard parametrization of the extended Kitaev-Heisenberg model, such as K, J, and Γ, are significant, making the offered consideration relevant to a much wider parameter space. Furthermore, the dynamical structure factor near the paramagnetic critical point illustrates this relevance by showing features that are reminiscent of the ones observed in α–RuCl 3 , underscoring that they are not unique and should be common to a wide range of parameters of the model and, by extension, to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Energetics of hydroxylbastnäsite solid solutions, La 1-$\chi$ Nd $\chi$ CO 3 OH

Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-scale, Multi-disciplinary, and Multi-agent Explainable AI with Koopman-Undergirded Learning, Prediction, and Analysis (M3EA KULPA) (Project Closeout Report)

The goal of this project was to develop and use domain-aware machine learning formulations, based on the Koopman Operator (KO), for modelling multi-scale, multi-disciplinary (e.g., multi-physics), and/or multi-agent systems. The project developed these formulations for the following cases: • Systems with dynamics at two separate time scales, • Systems with a bi-level hierarchical control structure, • Systems with bi-level hierarchical control and dynamics at two separate time scales (the lower level controls operating at the faster time scale), and • Systems with n separate but interacting agents/disciplines (with/without control, respectively); the controls for each agent could include bi-level hierarchical control and dynamics at two separate time scales as described above. The project then defined a set of dynamical systems consisting of different nonlinear oscillators that could be used to test these different formulations and then subsequently learned the KO models for those systems. With the KO models, we were able to do the following: • Quantify system stability, including both long-term and transient behavior, • Quantify the effects of feedbacks between the different time scales and agents/disciplines in terms of those feedbacks’ effects on system stability, • Replace a standard Proportional-Integral (PI) control in the hierarchical control structure with a KO-based Linear-Quadratic Regular (LQR), a form of optimal control, • Calculate optimal supervisory control policies a) with and without time scale separated dynamics at the lower level control levels and b) with both PI and KO-based LQR lower level control policies, and • Calculate dynamic Nash equilibria for multi-agent systems where each agent makes its own control decisions.

97 MATHEMATICS AND COMPUTING↗

The National Criticality Experiments Research Center and its role in support of advanced reactor design

The National Criticality Experiments Research Center (NCERC) located at the Nevada National Security Site (NNSS) in the Device Assembly Facility (DAF) and operated by Los Alamos National Laboratory (LANL) is the only general purpose critical experiments facility in the United States. Experiments from subcritical to critical and above prompt critical are carried out at NCERC on a regular basis. In recent years, NCERC has become more involved in experiments related to nuclear energy, including the Kilopower/KRUSTY demonstration and the recent Hypatia experiment. Multiple nuclear energy related projects are currently ongoing at NCERC. This paper discusses NCERC’s role in advanced reactor design and how that role may change in the future.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Temporally-consistent koopman autoencoders for forecasting dynamical systems

Absence of sufficiently high-quality data often poses a key challenge in data-driven modeling of high-dimensional spatio-temporal dynamical systems. Koopman Autoencoders (KAEs) harness the expressivity of deep neural networks (DNNs), the dimension reduction capabilities of autoencoders, and the spectral properties of the Koopman operator to learn a reduced-order feature space with simpler, linear dynamics. However, the effectiveness of KAEs is hindered by limited and noisy training datasets, leading to poor generalizability. To address this, we introduce the Temporally-Consistent Koopman Autoencoder (tcKAE), designed to generate accurate long-term predictions even with limited and noisy training data. This is achieved through a consistency regularization term that enforces prediction coherence across different time steps, thus enhancing the robustness and generalizability of tcKAE over existing models. We provide analytical justification for this approach based on Koopman spectral theory and empirically demonstrate tcKAE’s superior performance over state-of-the-art KAE models across a variety of test cases, including simple pendulum oscillations, kinetic plasma, and fluid flow data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dark photon constraints from CMB temperature anisotropies

The resonant conversion, within the inter-galactic medium, of regular photons into dark photons amplifies the anisotropy observed in the CMB, thereby imposing stringent constraints on the existence of light dark photons. In this study, we investigate the impact of light dark photons, with masses in the range 3 × 10$^{-15}$ eV < mA' < 3 × 10$^{-12}$ eV on the power spectrum of temperature anisotropies within the cosmic microwave background (CMB) radiation utilizing the state-of-the-art large-volume FLAMINGO cosmological simulations. Our results show that using full Planck data, one can expect the existing constraints on the dark photon mixing parameter in this mass range to improve by an order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Calibration sources for the LEGEND-200 experiment

In the search for a monochromatic peak as the signature of neutrinoless double beta decay an excellent energy resolution and an ultra-low background around the Q-value of the decay are essential. The LEGEND-200 experiment performs such a search with high-purity germanium detectors enriched in 76 Ge immersed in liquid argon. To determine and monitor the stability of the energy scale and resolution of the germanium diodes, custom-made, low-neutron emission 228 Th sources are regularly deployed in the vicinity of the crystals. Here we describe the production process of the 17 sources available for installation in the experiment, the measurements of their alpha- and gamma- activities, as well as the determination of the neutron emission rates with a low-background LiI(Eu) detector operated deep underground. With a flux of ( 4.27 ± 0.60 stat ± 0.92 syst ) × 10 -4 n / (kBq·s), approximately one order of magnitude below that of commercial sources, the neutron-induced background rate, mainly from the activation of 76 Ge, is negligible compared to other background sources in LEGEND-200.

47 OTHER INSTRUMENTATION↗

Coexistence of single-particle and collective excitation in Ni 61

Here, the high-spin states in 61 Ni have been studied using the fusion evaporation reaction, 50 Ti( 14 C, 3n) 61 Ni at an incident beam energy of 40 MeV. A Compton suppressed multi-HPGe detector setup, consisting of six clover detectors and three single-crystal HPGe detectors, was used to detect the de-exciting $\gamma$ rays from the excited states. The level scheme has been extended up to an excitation energy of 12.8 MeV and a tentative J π =35/2 + . The low-lying negative parity levels are found to be generated by single-particle excitation within the fp shell and also excitations to the g 9/2 orbitals as explained well with shell model calculations using the GXPF1Br+V MU (modified) interaction. Two rotational structures of regular E2 sequences with small to moderate axial deformation have been established at higher excitation energy. Most interestingly, two sequences of M1 transitions are reported for the first time and described as magnetic rotational bands. The shears mechanism for both the bands can be described satisfactorily by the geometrical model. The shell model calculation involving the cross shell excitation beyond the fp shell well reproduce the M1 and E2 sequences. The shell model predicted B(M1) values for the magnetic rotational band B1 show the decreasing trend with spin as expected with closing of the shears.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗