Search NASA⌕ Search

SEARCH · Search NASA

Results for “positive triangularity”

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 55 records · Page 3

Changes in oxidation state of chromium during LDEF exposure

The solar collector used for the McDonnell-Douglas Cascade Variable Heat Pipe, Experiment A0076 (Michael Grote - Principal Investigator) was finished with black chromium plating as a thermal control coating. The coating is metallic for low emittance, and is finely microcrystalline to a dimension which yields its high absorptivity. An underplate of nickel was applied to the aluminum absorber plate in order to achieve optimal absorptance characteristics from the black chromium plate surface. Experiment A0076 was located at tray position F9, receiving a projected 8.7 x 10 exp 21 atomic oxygen atoms/sq cm and 11,200 ESH solar radiation. During retrieval, it was observed that the aluminized kapton thermal blankets covering most of the tray were severely eroded by atomic oxygen, and that a 'flap' of aluminum foil was overlaying a roughly triangular shaped portion of the absorber panel. The aluminum foil 'flap' was lost sometime between the Long Duration Exposure Facility (LDEF) retrieval and deintegration. At deintegration, the black chromium was observed to have discolored where it had been covered by the foil 'flap'. A summary of the investigation into the cause of the discoloration is presented.

Golden, Johnny L.↗

The Static Longitudinal Stability and Control Characteristics in the Presence of the Ground of a Model Having a Triangular Wing and Canard

A wind-tunnel investigation was made of the low-speed characteristics of a canard configuration having triangular wing and canard surfaces with an aspect ratio of 2. The exposed area of the canard was 6.9 percent of the total wing area. The canard hinge line was located at 0.35 of its mean aerodynamic chord and was 0.5 wing mean aerodynamic chord lengths forward of the wing apex. The ground effects, which made the lift more positive and the -Pitching moment more negative at a given angle of attack, were unaffected by the canard. The stability of the model at a constant canard hinge-moment coefficient decreased to 0 near a lift coefficient of 1.0. In addition, the maximum lift coefficient at which the canard could provide balance was decreased by ground effects to less than 1.0 if the moment center was as far forward as 0.21 of the wing mean aerodynamic chord. The relative magnitude of interference effects between the canard and the wing and body is presented.

Buell, Donald A.↗

Exact Integrations of Polynomials and Symmetric Quadrature Formulas over Arbitrary Polyhedral Grids

This paper is concerned with two important elements in the high-order accurate spatial discretization of finite volume equations over arbitrary grids. One element is the integration of basis functions over arbitrary domains, which is used in expressing various spatial integrals in terms of discrete unknowns. The other consists of quadrature approximations to those integrals. Only polynomial basis functions applied to polyhedral and polygonal grids are treated here. Non-triangular polygonal faces are subdivided into a union of planar triangular facets, and the resulting triangulated polyhedron is subdivided into a union of tetrahedra. The straight line segment, triangle, and tetrahedron are thus the fundamental shapes that are the building blocks for all integrations and quadrature approximations. Integrals of products up to the fifth order are derived in a unified manner for the three fundamental shapes in terms of the position vectors of vertices. Results are given both in terms of tensor products and products of Cartesian coordinates. The exact polynomial integrals are used to obtain symmetric quadrature approximations of any degree of precision up to five for arbitrary integrals over the three fundamental domains. Using a coordinate-free formulation, simple and rational procedures are developed to derive virtually all quadrature formulas, including some previously unpublished. Four symmetry groups of quadrature points are introduced to derive Gauss formulas, while their limiting forms are used to derive Lobatto formulas. Representative Gauss and Lobatto formulas are tabulated. The relative efficiency of their application to polyhedral and polygonal grids is detailed. The extension to higher degrees of precision is discussed.

Liu, Yen↗

An Experimental Investigation of a Triangular Wing of Aspect Ratio 2 and a Body Warped to be Trimmed at M = 2.24

A cambered and twisted triangular wing of aspect ratio 2 in combination with a cambered body was investigated experimentally to determine the effectiveness of the camber in reducing the drag due to lift at trim at supersonic speeds. Four arrangements were tested comprising all combinations of a symmetrical and a cambered wing with a symmetrical and a cambered body. The camber shape investigated was derived by linearized lifting surface theory for triangular wings with sonic leading edges and satisfied the requirement that the wing be trimmed at the design Mach number and lift coefficient. The experimental results for the cambered wing and cambered body showed that the drag coefficient at trim was always greater, at the same lift coefficient, than that for the untrimmed symmetrical wing and body. The trim lift coefficient was positive and decreased with increasing Mach number. At the design Mach number of 2.24, the trim lift coefficient was somewhat lower and the drag coefficient was higher than values predicted by linearized lifting surface theory for the wing alone. A comparison of the trim lift-drag ratio of the cambered wing and cambered body with values obtained by trimming the symmetrical wing and symmetrical body either with a canard or a trailing-edge flap showed that, at approximately the design Mach number the cambered configuration developed a somewhat higher value than the trailing-edge flap configuration but a lower value than the canard configuration.

Adams, Gaynor J.↗

Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a tubine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of independence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

Choo, Yung K.↗

Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a turbine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of dependence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

Choo, Yung K.↗

Space station architecture, module, berthing hub, shell assembly, berthing mechanism and utility connection channel

The geometric form of a space station is presented that includes a description of a plurality of modules and berthing hubs, joined by interconnections which are sideways connectable. The modules and hubs are fastened together in a triangular configuration in three dimensions. The interconnections include a pair of opposed, axially aligned, flanged ports and a clamp latch formed from a plurality of sections hinged along their length and extending circumferentially around the flanged ports. A hermetic seal is formed between the ports. A utilities connection channel extends between the ports. The channel has a shell with utilities connectors movable between an extended position to mating connectors in the modules and a withdrawn position. Assembly sequence and common module shell structure is detailed.

Cohen, M. M.↗

Resolution of the buoyancy in the 8-foot high temperature tunnel combustor

Currently, the 8-Foot High Temperature Tunnel (8-Ft. HTT) combustor produces a good profile at only one point (2000 psia and 3650 R with oxygen enrichment). Air is enriched with oxygen (liquid) so that the combustor product gas will contain the volumetric amount of oxygen normally found in air. The oxygen enriched air has a large fraction that is not reacted and flows through the outer periphery of the fuel injector. This ring of cold air in addition to the relatively cold walls of the combustor set up buoyancy forces that produce a segregation of relatively cool gases at the bottom of the combustor exit. The basic problem is to produce a test gas that has uniform properties at all combustor conditions. The combustor temperature may be as high as 3700 R or as low as 2000 R. Combustor pressures can be as high as 3500 psia (no oxygen enrichment) and as low as 600 psia. The segregation is most severe with oxygen enriched air, since its temperature is lower and its density is high. The combustor is lined with nickel 201 and can be operated at about 1600 R maximum. A global mixing process is desired that produces an acceptable profile of temperature, species, and velocity at the exit of the combustor. The ultimate goal is a temperature profile with about 100 R variance and about 2 percent variance in oxygen. The exit total temperature must not be lowered significantly by the mixing apparatus or mechanisms employed. If immersed bodies are used, they must also be kept very hot. All combustor wall modifications must be able to survive the heat and structural conditions of the varied operating conditions. Our approach to resolving this issue is being conducted in three stages: (1) Consider mixing exclusively, (2) Resolve the heat transfer concerns resulting from the chosen mixing strategy, and (3) Solve the material and structural problems resulting from stages (1) and (2). Since the 8-Ft. HTT is unavailable for experimentation, the study is conducted exclusively with computational fluid dynamic (CPD) codes (Fluent/Uns and Rampant 3.1) using unstructured grid through body fitted coordinates. Both CFD codes are general purpose Navier-Stokes solution packages that can solve integral conservation equations for conservation of mass, momentum and energy. The governing equations are discretized using a control-volume finite-element method on unstructured triangular 2-D grids. In the interest of time, a 3-D tetrahedral grid was used to check the 2-D results on one mixing strategy and the 2-D results were confirmed. Preliminary results indicate that excellent mixing can be achieved with a body placed in the center of the flowing hot fluid with a minor modification to the combustor wall similar to a model positioned in the test section of a wind tunnel. The concept here, is to create longitudinal vortices strong enough to bring both fluids into intimate contact with each other near the wall where the cold fluid resides. However, there may be a trade-off in the choice of mixing strategy, heat transfer and structural requirements. Currently, we are examining the more promising geometries for heat transfer concerns and developing strategies for the material of construction of the center body.

Loney, Norman W.↗

Space station architecture, module, berthing hub, shell assembly, berthing mechanism and utility connection channel

A space station (20) includes a plurality of modules (24) and berthing hubs (22), joined by interconnections (26) which are sideways connectable. The modules (24) and hubs (22) are fastened together in a triangular configuration in three dimensions. The interconnections (26) include a pair of opposed, axially aligned, flanged ports (50) and a clamp latch (52) formed from a plurality of sections (54, 56 and 58) hinged along their length and extending circumferentially around the flanged ports (50). A hermetic seal (63) is formed between the ports (50). A utilities connection channel (68) extends between the ports (50). The channel (68) has a shell (70) with utilities connectors (74) movable between an extended position to mating connectors in the modules (24) and a withdrawn position. Assembly sequence and common module shell structure is detailed.

Cohen, Marc M.↗

X-ray investigation of molten crystal hydrates H2SO4(nH2O) and HNO3(nH2O)

Integral analysis of the intensity of the electron density distribution curve in molten crystal hydrates provided by X-ray analysis, permits the following conclusions on the structure of the complex SO and NO ions, and the short-range order in the structure of the solution. The SO4 ion in the solution has a tetrahedral structure with an S to O distance equal to 1.5 A. For the NO3 in the solution, a planar triangular shape is probable, with an N to O distance equal to 1.2 A. Preferential distances between each of the oxygens of the SO ion and the nearest molecules of water proved near to the corresponding distances in solid crystal hydrates. For an (H2SO4)(H2O) solution, the average number of water molecules surrounding each oxygen atom of the SO4 (--) ion was on the order of 1.3 molecules. Hence the preferential distances between the water molecules and the oxygen atoms of the SO ion, and the preference of their mutual position, correspond to the fixed position of these same elements of the structure in the solid crystal hydrate.

Romanova, A. V.↗

Estimation of Rotary Stability Derivatives at Subsonic and Transonic Speeds

The first part of this paper pertains to the estimation of subsonic rotary stability derivatives of wings. The unsteady potential flow problem is solved by a superposition of steady flow solutions. Numerical results for the damping coefficients of triangular wings are presented as functions of aspect ratio and Mach number, and are compared with experimental results over the Mach number range 0 to 1. In the second part, experimental results are used. to point out a close correlation between the nonlinear variations with angle of attack of the static pitching-moment curve slope and the damping-in-pitch coefficient. The underlying basis for the correlation is found as a result of an analysis in which the indicial function concept and. the principle of super-position are adapted to apply to the nonlinear problem. The form of the result suggests a method of estimating nonlinear damping coefficients from results of static wind-tunnel measurements.

Tobak, Murray↗

Lift, Drag, and Pitching Moment of an Aspect-Ratio-2 Triangular Wing with Leading-Edge Flaps Designed to Simulate Conical Camber

An investigation was conducted to determine the effectiveness of leading-edge flaps in reducing the drag at lifting conditions of a triangular wing of aspect ratio 2.0. The flaps, deflected to simulate conically cambered wings having a wide range of design lift coefficients, were tested over a Mach number range of 0.70 to 2.22 through an angle-of-attack variation from -6 deg to +18 deg at a constant Reynolds number of 3.68 million based on the wing mean aerodynamic chord. A symmetrical wing of the same plan form and aspect ratio was also tested to provide a basis for comparison. The experimental results showed that with the flaps in the undeflected position, a small amount of fixed leading-edge droop incorporated over the outboard 5 percent of the wing semispan was as effective at high subsonic speeds as conical camber in improving the maximum lift-drag ratio above that of the symmetrical wing. At supersonic speeds, the penalty in minimum drag above that of the symmetrical wing was less than that incurred by conical camber. Deflecting the leading-edge flaps about the hinge line through 80 percent of the wing semispan resulted in further improvements of the drag characteristics at lift coefficients above 0.20 throughout the Mach number range investigated. The lift and pitching-moment characteristics were not significantly affected by the leading-edge flaps.

Menees, Gene P.↗

An Experimental Investigation of the Hinge-moment Characteristics of a Constant-chord Control Surface Oscillating at High Frequency

The results of an experimental investigation of the hinge-moment characteristics of a constant-chord control surface oscillating at high frequency is presented. The control surface was mounted on an aspect-ratio-2 triangular wing. The aerodynamic restoring-moment coefficient and damping-moment coefficient were determined at a frequency of 260 cycles per second for a Mach number range of 0.6 to 0.8 and 1.3 to 1.9 at angles of attack of 5 degrees and 10 degrees. The test results showed linear theory to be a reliable guide to the prediction of the trend of the restoring-moment coefficient with Mach number for the supersonic speed range of the investigation but overestimated the magnitude of the coefficient. The experimental values of the damping-moment coefficient were, for the most part, more positive than those indicated by the theory and, for some conditions, could lead to instability of the control surface. Comparison of the results of this investigation with those of previous investigations at 0 and 50 cycles per second showed that frequency had little effect on the restoring-moment coefficient. The damping-moment coefficient was similarly insensitive to frequency at an oscillation amplitude of plus-or-minus 1.0 degrees but at an amplitude ofplus-or-minus 2.5 degrees the results showed a destabilizing shift with increasing frequency.

Reese, David E JR↗

Longitudinal Stability and Drag Characteristics at Mach Numbers from 0.70 to 1.37 of Rocket-propelled Models Having a Modified Triangular Wing

A modified triangular wing of aspect ratio 2.53 having an airfoil section 3.7 percent thick at the root and 5.98 percent thick at the tip was designed in an attempt to improve the lift and drag characteristics of triangular wings. Free-flight drag and stability tests were made using rocket-propelled models equipped with the modified wing. The Mach number range of the test was from 0.70 to 1.37. Test results indicated the following: The lift-curve slope of wing plus fuselage approaches the theoretical value of wing alone at supersonic Mach numbers. The drag coefficient, based on total wing area, for wing plus interference was approximately 0.0035 at subsonic Mach numbers and 0.0080 at supersonic Mach numbers. The maximum shift in aerodynamic center for the complete configuration was 14 percent in the rearward direction from the forward position of 51.5 percent of mean aerodynamic chord at subsonic Mach numbers. The variation of lift and moment with angle of attack was linear at supersonic Mach numbers for the range of coefficients covered in the test. The high value of lift-curve slope was considered to be a significant result attributable to the wing modifications.

Chapman, Rowe, Jr↗

Distortion Representation of Forecast Errors for Model Skill Assessment and Objective Analysis

We study a novel characterization of errors for numerical weather predictions. In its simplest form we decompose the error into a part attributable to phase errors and a remainder. The phase error is represented in the same fashion as a velocity field and will be required to vary slowly and smoothly with position. A general distortion representation allows for the displacement and a bias correction of forecast anomalies. In brief, the distortion is determined by minimizing the objective function by varying the displacement and bias correction fields. In the present project we use a global or hemispheric domain, and spherical harmonics to represent these fields. In this project we are initially focusing on the assessment application, restricted to a realistic but univariate 2-dimensional situation. Specifically we study the forecast errors of the 500 hPa geopotential height field for forecasts of the short and medium range. The forecasts are those of the Goddard Earth Observing System data assimilation system. Results presented show that the methodology works, that a large part of the total error may be explained by a distortion limited to triangular truncation at wavenumber 10, and that the remaining residual error contains mostly small spatial scales.

Hoffman, Ross N.↗

Dividing the Concentrator Target From the Genesis Mission

The Genesis spacecraft, launched in 2001, traveled to a Lagrangian point between the Earth and Sun to collect particles from the solar wind and return them to Earth. However, during the return of the spacecraft in 2004, the parachute failed to open during descent, and the Genesis spacecraft crashed into the Utah desert. Many of the solar wind collectors were broken into smaller pieces, and the field team rapidly collected the capsule and collector pieces for later assessment. On each of the next few days, the team discovered that various collectors had survived intact, including three of four concentrator targets. Within a month, the team had imaged more than 10,000 fragments and packed them for transport to the Astromaterials Acquisition and Curation Office within the ARES Directorate at JSC. Currently, the Genesis samples are curated along with the other extraterrestrial sample collections within ARES. Although they were broken and dirty, the Genesis solar wind collectors still offered the science community the opportunity to better understand our Sun and the solar system as a whole. One of the more highly prized concentrator collectors survived the crash almost completely intact. The Genesis Concentrator was designed to concentrate the solar wind by a factor of at least 20 so that solar oxygen and nitrogen isotopes could be measured. One of these materials was the Diamond-on-Silicon (DoS) concentrator target. Unfortunately, the DoS concentrator broke on impact. Nevertheless, the scientific value of the DoS concentrator target was high. The Genesis Allocation Committee received a request for approximately 1 cm(sup 2) of the DoS specimen taken near the focal point of the concentrator for the analysis of solar wind nitrogen isotopes. The largest fragment, Genesis sample 60000, was designated for this allocation and needed to be precisely cut. The requirement was to subdivide the designated sample in a manner that prevented contamination of the sample and minimized the risk of losing or breaking the precious requested sample fragment. The Genesis curator determined that the use of laser scribing techniques to "cut" a precise line and subsequently cleave the sample (in a controlled break of the sample along that line) was the best method for accomplishing the sample subdivision. However, there were risks, including excess heating of the sample, that could cause some of the implanted solar wind to be lost via thermal diffusion. Accidentally breaking the sample during the handling and cleaving process was an additional risk. Early in fiscal year 2013, to address this delicate, complicated task, the ARES Directorate assembled its top scientists to develop a cutting plan that would ensure success when applied to the actual concentrator target wafer; i.e., to produce an approximately 1 cm(sup 2) piece from the requested area of the wafer. The team, subsequently referred to as the JSC Genesis Tiger Team, spent months researching and testing parameters and techniques related to scribing, cleaving, transporting, handling, and holding (i.e., mounting) the specimen. The investigation required considerable "thinking outside the box," and many, many trials using nonflight wafer analogs. After all preliminary testing, the following method was adopted as the final cutting plan. It was used in two final end-to-end practice runs before being used on the actual flight target wafer. The wafer was oriented on the laser cutting stage with the 100 and 010 directions of the wafer parallel to the corresponding X and Y directions of the cutting stage. The laser was programed to scribe 31 lines of the appropriate length along the Y stage direction. The programed scribe lines were separated by 5 micron in the X direction. The laser parameters were set as follows: (1) The laser power was 0.5 watts; (2) each line consisted of 50 passes, with the Z position being advanced 5 micron per pass; and (3) 30 s would elapse before the next line was scribed to allow for wafer cool down from any possible heating via the laser. The ablated material that "stuck" in the "scribe-cut" was removed from the "cut" using an ultrasonic micro-tool. After all the ablated silicon was removed from the wafer, the wafer was repositioned in exactly the same orientation on the laser stage. The laser was focused using the bottom of the wafer channel, and the 31-line scribing pattern described above was reprogrammed using the Z position of the groove bottom as the starting Z value instead of the top wafer surface, which was used previously. Upon completion of the second set of scribes, the ultrasonic micro-tool was again used to clean out the cut. The wafer was remounted on the stage in exactly the same orientation as before. The laser was again focused on the bottom of the groove. This time, however, the laser was.programed to scribe only one line down the exact center of the channel. The final scribe line consisted of 100 passes with a Z advance of 5 micron per pass and with the laser power set at 0.5 watts. As mentioned above, the final cutting plan was practiced in two end-to-end trials using non-flight, triangular-shaped silicon wafers similar in size and orientation to the actual DOS 60000 target sample. The actual scribing of the triangular-shaped wafers required scribing two lines and cleaving (i.e. scribe-cleave, then scribe-cleave) to obtain the piece requested for allocation. Early in December 2012, after many months of experiments and practicing and perfecting the techniques and procedures, the team successfully subdivided the Genesis DoS 60000 target sample, one of the most scientifically important samples from the Genesis mission (figure 2). On December 17, 2012, the allocated piece of concentrator target sample was delivered to the requesting principal investigator.The cutting plan developed for the subdivision of this sample will be used as the model for subdividing future requested Genesis flight wafers (appropriately modified for different wafer types).

Lauer, H. V., Jr.↗

Wind-Tunnel Investigation of a Small-Scale Model of an Aerial Vehicle Supported by Tilting Ducted Fans

A wind-tunnel investigation has been made to study the static longitudinal and lateral stability characteristics of a simplified aerial vehicle supported by ducted fans that tilt relative to the airframe. The ducts were in a triangular arrangement with one duct in front and two at the rear in order to minimize the influence of the downwash of the front duct on the rear ducts. The results of the investigation were compared with those of a similar investigation for a tandem two-duct arrangement in which the ducts were fixed (rather than tiltable) relative to the airframe, since the three-duct configuration had been devised in an attempt to avoid some of the deficiencies of the tandem fixed-duct configuration. The results of the investigation indicated that the tilting-duct arrangement had less noseup pitching moment for a given forward speed than the tandem fixed-duct arrangement. The model had less angle-of-attack instability than the tandem fixed-duct arrangement. The model was directionally unstable but had a positive dihedral effect throughout the test speed range.

Smith, Charles C., Jr.↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗