Results for “Parametric processes” — Page 38

NASA NTRS1994

Attitude stability of a spinning spacecraft during appendage deployment/retraction

The work presented is motivated by the need for a national satellite rescue policy, not the ad hoc policy now in place. In studying different approaches for a national policy, the issue of capture and stabilization of a tumbling spacecraft must be addressed. For a rescue mission involving a tumbling spacecraft, it may be advantageous to have a rescue vehicle which is compact and 'rigid' during the rendezvous/capture phase. After capture, passive stabilization techniques could be utilized as an efficient means of detumbling the resulting system (i.e., both the rescue vehicle and captures spacecraft). Since the rescue vehicle is initially compact and 'rigid,' significant passive stabilization through energy dissipation can only be achieved through the deployment of flexible appendages. Once stabilization is accomplished, retraction of the appendages before maneuvering the system to its final destination may also prove advantageous. It is therefore of paramount interest that we study the effect of appendage deployment/retraction on the attitude stability of a spacecraft. Particular interest should be paid to appendage retraction, since if this process is destabilizing, passive stabilization as proposed may not be useful. Over the past three decades, it has been an 'on-again-off-again affair' with the problem of spacecraft appendage deployment. In most instances, these studies have been numerical simulations of specific spacecraft configurations for which there were specific concerns. The primary focus of these studies was the behavior of the appendage during deployment; the effects of appendage retraction was considered only in one of these studies. What is missing in the literature is a thorough study of the effects of appendage deployment/retraction on the attitude stability of a spacecraft. This paper presents a rigorous analysis of the stability of a spinning spacecraft during the deployment or the retraction of an appendage. The analysis is simplified such that meaningful insights into the problem can be inferred; it is not overly simplified such that critical dynamical behavior is neglected. The system is analyzed assuming that the spacecraft hub is rigid. The appendage deployment mechanism is modeled as a point mass on a massless rod whose length undergoes prescribed changes. Simplified flexibility effects of the appendage are included. The system is examined for stability by linearizing the equations in terms of small deviations from steady, noninterfering coning motion. Routh's procedure for analyzing small deviations from steady motion in dynamical systems is utilized in the analysis. The system of equations are nondimensionalized to facilitate parametric studies. The results are presented in terms of a reduced number of nondimensional parameters so that some general conclusions may be drawn. Verification of the linear analysis is presented through numerical simulations of the complete nonlinear, nonautonomous, coupled equations.

Fitz-Coy, Norman↗

NASA NTRSDate not supplied

Regional and Model-Specific Response Types in A Global Gridded Crop Model Ensemble

Crop models are often employed to project crop yields under changing conditions such as global warming and associated management change for adaptation. Multi-model ensembles are promoted to enhance the robustness of projections, but questions remain on what causes often large differences between projections of individual models. Global Gridded Crop Models (GGCMs) are especially exposed to this question when applied for assessing climate change impacts, adaptation, environmental impacts of agricultural production, because their results are used in downstream analyses, such as in integrated assessment or economic modeling for projecting future land-use change. Even though global gridded crop models are often based on detailed field-scale models or have implemented similar modeling principles in other ecosystem models, global-scale models are subject to substantial uncertainties from both model structure and parametrization as well as from calibration and input data quality. AgMIP’s Global Gridded Crop Model Intercomparison (GGCMI) has thus set out to intercompare GGCMs in order to evaluate model performance, describe model uncertainties, identify inconsistencies within the ensemble and underlying reasons, and to ultimately improve models and modeling capacities. In phase 2 of the GGCMI activities, 12 modeling groups followed a modeling protocol that asked for up to 1404 31-year global simulations at 0.5 arc-degree spatial resolution to assess models’ sensitivities to changes in carbon dioxide (C; 4 different levels) temperature (T; 7 different offset levels), water supply (W; 9 levels), and nitrogen (N; 3 levels), the so-called CTWN experiment (Franke et al. 2020; http://dx.doi.org/10.5194/gmd-13-2315-2020). We here present analyses of model response types using impact response surfaces along the C, T, W, and N dimensions, respectively and collectively. Doing so, we can understand differences in simulated responses per driver rather than aggregated changes in yields. We find that models’ sensitivities to the individual driver dimensions are substantially different and often more different across models than across regions. A cluster analysis finds regional and model-specific patterns. There is some agreement across models with respect to the spatial patterns of response types but strong differences in the distribution of response type clusters across models suggests that models need to undergo further scrutiny. We suggest establishing standards in model process evaluation not only against historical dynamics but also against dedicated experiments across the CTWN dimensions.

crop models↗

NASA NTRS2000

Orbiter Capability for Providing Water to the International Space Station according to the Most Probable Flight Attitudes

Water to be generated by, delivered to, and processed by the International Space Station (ISS) is a critical Environmental Control and Life Support (ECLS) element, especially for the early ISS missions. A significant portion of the water required by the ISS shall be provided by the Shuttle Transportation System (STS) Orbiter. The balance of water generated by the Orbiter Fuel Cells (FC), minus that water consumed by the Orbiter Flash Evaporator System (FES) and crew, is available for transfer to the ISS. During later missions, crew respired and perspired water, as well as effluent water from the Orbiter LiOH canisters, will be collected as condensate and available for transfer to the ISS. Orbiter radiator performance provides the most variance in determining the amount of net Orbiter water available for transfer to the ISS. As radiator performance decreases, the dependence upon the FES (and FC water) increases for rejecting Orbiter waste heat. Generally, radiator performance decreases as the ISS assembly size increases (especially as solar arrays are added), and also as beta angle increases. ISS solar array deployment necessitates the use of models with articulating solar arrays (for Earth local-vertical attitudes), as array position dramatically affects Orbiter radiator performance. Recent developments in the relaxation of beta angle limitations have also increased the complexity and difficulty of providing water to the ISS. Other factors that may hinder the ability to transfer water are the number of empty Contingency Water Containers (CWCs) available, duration of open-hatch time, crew activity timeline, and full CWC storage capability. A parametric study has been accomplished that provides a quick-reference table for determining expected water generation rates for ISS missions 2A.2 through 7A.1. An hourly Orbiter water generation rate is reported according to a matrix that consists of: (1) (six) significant changes in ISS assembly configuration; (2) (four) beta angles (0 deg. , +37 deg., +53 deg. , and +75 deg.); (3) the (three) most representative ISS attitudes (XPOP-O, XPOP-180 and +XVV); (4) (four) Orbiter radiator configurations (both stowed, starboard deployed, port deployed, and both deployed) and (5) the (two) conditions (radiator inlet temperatures and fuel cell power) most consistent with sleep and wake periods. Those permutations of higher probability of occurrence than others have been identified. Another parametric study has been accomplished that provides a quick-reference table for determining expected water generation rates for ISS assembly complete missions. An hourly Orbiter water generation rate is reported according to a matrix that consists of: (1) (seven) beta angles (-75 deg., -60 deg., -30 deg., 0 deg., +30 deg., +60 deg., and +75 deg.); (2) the (nine) PYR angles that define the corners of the envelope; (3) (four) Orbiter radiator configurations (both stowed, starboard deployed, port deployed, and both deployed) and (4) the (two) conditions (radiator inlet temperatures and fuel cell power) most consistent with sleep and wake periods.

Dunaway, Brian↗

NASA NTRS2000

SAO and Kelvin Waves in the EuroGRIPS GCMS and the UK Meteorological Offices Analyses

This work is an intercomparison of four tropospheric-stratospheric climate models, the Unified Model (UM) of the U.K. Meteorological Office (UKMO), the model of the Free University in Berlin (FUB). the ARPEGE-climat model of the National Center for Meteorological Research (CNRM), and the Extended UGAMP GCM (EUGCM) of the Center for Global Atmospheric Modelling (CGAM), against the UKMO analyses. This comparison has been made in the framework of the "GSM-Reality Intercomparison Project for SPARC" (GRIPS). SPARC (Stratospheric Processes and their Role in Climate) aims are to investigate the effects of the middle atmosphere on climate and the GRIPS purpose is to organized a comprehensive assessment of current Middle Atmosphere-Climate Models (MACMs). The models integrations were made without identical contraints e.g. boundary conditions, incoming solar radiation). All models are able to represent the dominant features of the extratropical circulation. In this paper, the structure of the tropical winds and the strengths of the Kelvin waves are examined. Explanations for the differences exhibited. between the models. as well as between models and analyses, are also proposed. In the analyses a rich spectrum of waves (eastward and westward) is present and contributes to drive the SAO (SemiAnnual Oscillation) and the QBO (Quasi-Biennal Oscillation). The amplitude of the Kelvin waves is close to the one observed in UARS (Upper Atmosphere Research Satellite) data. In agreement with observations, the Kelvin waves generated in the models propagate into the middle atmosphere as wave packets which underlines convective forcing origin. In most models, slow Kelvin waves propagate too high and are hence overestimated in the upper stratosphere and in the mesosphere, except for the UM which is more diffusive. These waves are not sufficient to force realistic westerlies of the QBO or SAO westerly phases. If the SAO is represented by all models only two of them are able to generate westerlies between 10 hPa and 50 hPa. The importance of the role played by subgrided gravity waves is more and more recognized. Actually, the EUGCM which includes a parametrization of gravity waves with a non-zero phase speed is able to simulate. with however some unrealistic features, clear easterly to westerly transitions as well as westerlies downward propagations. Thermal damping is also important in the westerlies forcing in the stratosphere. The model ARPEGE-climat shows more westerlies in the stratosphere than tile other three models probably due to the use of a simplified scheme to predict the ozone distribution in the middle atmosphere.

Amodei, M.↗

NASA NTRSDate not supplied

A New Concurrent Engineering Tool for the New Mission Design Center at NASA Ames Research Center

The NASA Ames Mission Design Center (MDC) has undergone a significant transformation process in recent years, culminating in becoming a true Concurrent Engineering Center. The underlying goal of this transformation is to better serve the Ames PI community in a cost-effective and rapidly-responsive way, providing quality products pivotal to the decision-making proposal strategy of our center. The four pillars of the change focus on four areas: personnel, physical places, tools, and training. Personnel: the MDC switched from using a dedicated pool of engineers—mostly entry-level career—for a matrixed approach with more experienced Subject Matter Experts, usually involved in spaceflight projects. These SMEs get detailed to the MDC for the duration of the study on a part-time or full-time basis. Physical places: The physical aspect pertains mostly to the layout of the new Engineering building and the dedicated concurrent engineering layout to foster collaboration. The current plan includes a dedicated space for concurrent engineering sessions as well as dedicated project rooms for concept studies. Tool: the MDC team is developing a new in-house concurrent engineering tool to facilitate concurrent engineering sessions. The team leveraged the many lessons learned and insight gained from using the current concurrent engineering tool, in use since 2007. SMEs have been involved throughout its design and development, ensuring accuracy and providing validity to the implementation approach. This new tool is database-driven and allows a team to work concurrently on the same model of a mission concept. “Workspaces” are provided for each of the traditional spacecraft mission design disciplines, including cost accounting and systems engineering. There are multiple tiers of fidelity available within each workspace, which can be completed parametrically or independently from the other disciplines as needed or synchronously with the other disciplines. A “commit” step synchronizes a workspace with the mission database and identifies any potential conflicts, along with the user who entered the conflicting data. A Mission Summary workspace enables a Study Lead to run an effective CE session, displaying key graphics, workspace status information, and technical and monetary budget information. This new CE tool will be undergoing Beta Testing with its primary user base and a limited Parts database in the summer of 2020. The team plans to release the first fully operational version in December 2020. Training: working in early concept development and maturation for a space mission requires a unique mindset and being comfortable with uncertainties in an environment where requirements are not yet fully developed, and where changes and trade-offs happen at a high tempo. The MDC started a weekly training program to ensure a common knowledge base on how do develop and mature early concept studies, as well as how to operate in a concurrent engineering environment. In conclusion, the MDC transformation will provide the overall PI community focused on cost-effective small satellite scientific missions with a powerful tool to address and mature early concepts more effectively and efficiently.

David Mauro↗

NASA NTRS2015

Advanced Concept Studies for Supersonic Commercial Transports Entering Service in the 2018-2020 Period Phase 2

Lockheed Martin Aeronautics Company (LM), working in conjunction with General Electric Global Research (GE GR) and Stanford University, executed a 19 month program responsive to the NASA sponsored "N+2 Supersonic Validation: Advanced Concept Studies for Supersonic Commercial Transports Entering Service in the 2018-2020 Period" contract. The key technical objective of this effort was to validate integrated airframe and propulsion technologies and design methodologies necessary to realize a supersonic vehicle capable of meeting the N+2 environmental and performance goals. The N+2 program is aligned with NASA's Supersonic Project and is focused on providing system level solutions capable of overcoming the efficiency, environmental, and performance barriers to practical supersonic flight. The N+2 environmental and performance goals are outlined in the technical paper, AIAA-2014-2138 (Ref. 1) along with the validated N+2 Phase 2 results. Our Phase 2 efforts built upon our Phase 1 studies (Ref. 2) and successfully demonstrated the ability to design and test realistic configurations capable of shaped sonic booms over the width of the sonic boom carpet. Developing a shaped boom configuration capable of meeting the N+2 shaped boom targets is a key goal for the N+2 program. During the LM Phase 1 effort, LM successfully designed and tested a shaped boom trijet configuration (1021) capable of achieving 85 PLdB under track (forward and aft shock) and up to 28 deg off-track at Mach 1.6. In Phase 2 we developed a refined configuration (1044-2) that extended the under 85 PLdB sonic boom level over the entire carpet of 52 deg off-track at a cruise Mach number of 1.7. Further, the loudness level of the configuration throughout operational conditions calculates to an average of 79 PLdB. These calculations rely on propagation employing Burger's (sBOOM) rounding methodology, and there are indications that the configuration average loudness would actually be 75 PLdB. We also added significant fidelity to the design of the configuration in this phase by performing a low speed wind tunnel test at our LTWT facility in Palmdale, by more complete modelling of propulsion effects in our sonic boom analysis, and by refining our configuration packaging and performance assessments. Working with General Electric, LM performed an assessment of the impact of inlet and nozzle effects on the sonic boom signature of the LM N+2 configurations. Our results indicate that inlet/exhaust streamtube boundary conditions are adequate for conceptual design studies, but realistic propulsion modeling at similar stream-tube conditions does have a small but measurable impact on the sonic boom signature. Previous supersonic transport studies have identified aeroelastic effects as one of the major challenges associated with the long, slender vehicles particularly common with shaped boom aircraft (Ref. 3). Under the Phase 2 effort, we have developed a detailed structural analysis model to evaluate the impact of flexibility and structural considerations on the feasibility of future quiet supersonic transports. We looked in particular at dynamic structural modes and flutter as a failure that must be avoided. We found that for our N+2 design in particular, adequate flutter margin existed. Our flutter margin is large enough to cover uncertainties like large increases in engine weight and the margin is relatively easy to increase with additional stiffening mass. The lack of major aeroelastic problems probably derives somewhat from an early design bias. While shaped boom aircraft require long length, they are not required to be thin. We intentionally developed our structural depths to avoid major flexibility problems. So at the end of Phase 2, we have validated that aeroelastic problems are not necessarily endemic to shaped boom designs. Experimental validation of sonic boom design and analysis techniques was the primary objective of the N+2 Supersonic Validations contract; and in this Phase, LM participated in four high speed wind tunnel tests. The first so-called Parametric Test in the Ames 9x7 tunnel did an exhaustive look at variation effects of the parameters: humidity, total pressure, sample time, spatial averaging distance and number of measurement locations, and more. From the results we learned to obtain data faster and more accurately, and made test condition tolerances easy to meet (eliminating earlier 60 percent wasted time when condition tolerances could not be held). The next two tests used different tunnels. The Ames 11 ft tunnel was used to test lower Mach numbers of 1.2 and 1.4. There were several difficulties using this tunnel for the first time for sonic boom including having to shift the measurement Mach numbers to 1.15 and 1.3 to avoid flow problems. It is believed that the 11 ft could be used successfully to measure sonic boom but there are likely to be a number of test condition restrictions. The Glenn 8x6 ft tunnel was used next and the tunnel has a number of desirable features for sonic boom measurement. While the Ames 9x7 can only test Mach 1.55 to 2.55 and the 11 ft can only test Mach 1.3 and lower, the Glenn 8x6 can test continuously from Mach 0.3 to 2.0. Unfortunately test measurement accuracy was compromised by a reference pressure drift. Post-test analysis revealed that the drift occurred when Mach number drifted slightly. Test measurements indicated that if Mach number drift is eliminated, results from the 8x6 would be more accurate, especially at longer distances, than results from the 9x7. The fourth test in the 9x7, called LM4, used everything we learned to comprehensively and accurately measure our new 1044-02 configuration with a full-carpet shaped signature design. Productivity was 8 times greater than our Phase 1 LM3 test. Measurement accuracy and repeatability was excellent out to 42 in. However, measurements at greater distances require the rail in the aft position and become substantially less accurate. Further signature processing or measurement improvements are needed for beyond near-field signature validation.

Morgenstern, John↗

NASA NTRS2016

Satellite Constellation Cost Modeling: An Aggregate Model

Satellite constellations and Distributed Spacecraft Mission (DSM) architectures offer unique benefits to Earth observation scientists and unique challenges to cost estimators. The Cost and Risk (CR) module of the Tradespace Analysis Tool for Constellations (TAT-C) being developed by NASA Goddard seeks to address some of these challenges by providing a new approach to cost modeling, which aggregates existing Cost Estimating Relationships (CER) from respected sources, cost estimating best practices, and data from existing and proposed satellite designs. Cost estimation through this tool is approached from two perspectives: parametric cost estimating relationships and analogous cost estimation techniques. The dual approach utilized within the TAT-C CR module is intended to address prevailing concerns regarding early design stage cost estimates, and offer increased transparency and fidelity by offering two preliminary perspectives on mission cost. This work outlines the existing cost model, details assumptions built into the model, and explains what measures have been taken to address the particular challenges of constellation cost estimating. The risk estimation portion of the TAT-C CR module is still in development and will be presented in future work. The cost estimate produced by the CR module is not intended to be an exact mission valuation, but rather a comparative tool to assist in the exploration of the constellation design tradespace. Previous work has noted that estimating the cost of satellite constellations is difficult given that no comprehensive model for constellation cost estimation has yet been developed, and as such, quantitative assessment of multiple spacecraft missions has many remaining areas of uncertainty. By incorporating well-established CERs with preliminary approaches to approaching these uncertainties, the CR module offers more complete approach to constellation costing than has previously been available to mission architects or Earth scientists seeking to leverage the capabilities of multiple spacecraft working in support of a common goal.

Science Data Processing↗

NASA NTRS2016

Streamlining the Design Tradespace for Earth Imaging Constellations

Satellite constellations and Distributed Spacecraft Mission (DSM) architectures offer unique benefits to Earth observation scientists and unique challenges to cost estimators. The Cost and Risk (CR) module of the Tradespace Analysis Tool for Constellations (TAT-C) being developed by NASA Goddard seeks to address some of these challenges by providing a new approach to cost modeling, which aggregates existing Cost Estimating Relationships (CER) from respected sources, cost estimating best practices, and data from existing and proposed satellite designs. Cost estimation through this tool is approached from two perspectives: parametric cost estimating relationships and analogous cost estimation techniques. The dual approach utilized within the TAT-C CR module is intended to address prevailing concerns regarding early design stage cost estimates, and offer increased transparency and fidelity by offering two preliminary perspectives on mission cost. This work outlines the existing cost model, details assumptions built into the model, and explains what measures have been taken to address the particular challenges of constellation cost estimating. The risk estimation portion of the TAT-C CR module is still in development and will be presented in future work. The cost estimate produced by the CR module is not intended to be an exact mission valuation, but rather a comparative tool to assist in the exploration of the constellation design tradespace. Previous work has noted that estimating the cost of satellite constellations is difficult given that no comprehensive model for constellation cost estimation has yet been developed, and as such, quantitative assessment of multiple spacecraft missions has many remaining areas of uncertainty. By incorporating well-established CERs with preliminary approaches to approaching these uncertainties, the CR module offers more complete approach to constellation costing than has previously been available to mission architects or Earth scientists seeking to leverage the capabilities of multiple spacecraft working in support of a common goal.

Science Data Processing↗

NASA NTRS2016

A Survey of Cost Estimating Methodologies for Distributed Spacecraft Missions

Satellite constellations and Distributed Spacecraft Mission (DSM) architectures offer unique benefits to Earth observation scientists and unique challenges to cost estimators. The Cost and Risk (CR) module of the Tradespace Analysis Tool for Constellations (TAT-C) being developed by NASA Goddard seeks to address some of these challenges by providing a new approach to cost modeling, which aggregates existing Cost Estimating Relationships (CER) from respected sources, cost estimating best practices, and data from existing and proposed satellite designs. Cost estimation through this tool is approached from two perspectives: parametric cost estimating relationships and analogous cost estimation techniques. The dual approach utilized within the TAT-C CR module is intended to address prevailing concerns regarding early design stage cost estimates, and offer increased transparency and fidelity by offering two preliminary perspectives on mission cost. This work outlines the existing cost model, details assumptions built into the model, and explains what measures have been taken to address the particular challenges of constellation cost estimating. The risk estimation portion of the TAT-C CR module is still in development and will be presented in future work. The cost estimate produced by the CR module is not intended to be an exact mission valuation, but rather a comparative tool to assist in the exploration of the constellation design tradespace. Previous work has noted that estimating the cost of satellite constellations is difficult given that no comprehensive model for constellation cost estimation has yet been developed, and as such, quantitative assessment of multiple spacecraft missions has many remaining areas of uncertainty. By incorporating well-established CERs with preliminary approaches to approaching these uncertainties, the CR module offers more complete approach to constellation costing than has previously been available to mission architects or Earth scientists seeking to leverage the capabilities of multiple spacecraft working in support of a common goal.

0↗

NASA NTRSDate not supplied

Development of a Cold-Walled Molten Regolith Electrolysis Reactor for Lunar Oxygen Production

On the lunar surface, production of commodi-ties to support human presence, such as water, food and oxygen, and sustain the growth of a per-manent outpostwill likely require the use oflocal resources. The moon is covered almost entirely withfragmented oxide minerals known as regolithhundreds of meter thick.As a resource, it is rich in oxygen (> 42 wt.%) bound in a solid state with a variety of metals. The molten regolith electrolysis (MRE)reactor is a promising technology for the production of gaseous oxygen from the lunar reg-olithin a simple, single-stepreaction that requires minimalconsumable materials, produces oxygen and metals with high electrical efficiency and high yields from any regolithcomposition.This process involvesmelting regolith to~1600°C then electro-lyzing the molten pool to separate metal and oxy-gen ions that are then collected as liquid metal and gaseous oxygen at the respective electrodes. Lab-scale demonstrations of the MRE technology have previously reliedon external heating sources to bring the entirety of the reactor up to the operating temperaturewhich creates corrosive interfaces be-tween the molten regolith and the containment ma-terial in the reactor, limiting the overall lifespan of a reactor[1]. The Gaseous Lunar Oxygen from Regolith Electrolysis (GaLORE) project is focused on the development of a “cold-walled” or “Joule-heated” reactor design in which an internal heating source is used to selectively melt a pool of regolith between the electrodes of the reactor, leaving a shell of solidified regolith between the molten pool and the containment vessel of the reactor. This next generation reactor concept has been under development as molten oxide electrolysis (MOE) by MIT and Boston Metal for the production of iron from pure ores for terrestrial application [2]. The GaLORE project in engaged in early development of the technology for use with varying lunar regolith compositionsin the lunar environment. Thermal modelling of a proposed cold-walledreac-tor design were used as a scaffold to develop pa-rameters for a feasible reactor shape and size as well as target energy consumption[3]. The current development effort for the cold-walled reactor de-sign will be presented as a trade study of the most promising techniques for melting regolithwithin the constraints imposed by the lunar environment.Heater devices are designed to accommodate lim-ited electrical power availabilityon the moon, a wide range of regolith compositions that may be seen on the moon, limited metals available for re-placing consumed parts,and the low thermal con-ductivity of granular regolith in vacuum.Heater de-vices will be down-selectedbased on performance measurements within the above operationalcon-straints,and selected devices will be integrated into a reactor with electrodes to begin producing oxygen. [1]Sibille,L.,Sadoway, D.R.,Sirk, A., Tripathy,P.,Melendez, O., Standish, E., Dominguez, J. A., Stefanescu, D.M.,Curreri, P.A., Poizeau,S.,2009. “Recent Advancesin Scale-up Development of Molten Regolith Electrolysis for Oxygen Production in support of a Lunar Base.”AIAA 2009-659, 47th AIAA Aerospace Sciences Meeting, 5 -8 January 2009, Orlando, FL. [2] Boston Metal, https://www.boston-metal.com/moe-technology/#moe-process[3]Schreiner, S.S.,Sibille, L., Dominguez, J.A., Hoffman, J.A., 2016. "A parametric sizing model for Molten Regolith Electrolysis reactors to pro-duce oxygen on the Moon." Advances in Space Research 57.,7,1585-1603.

K D Grossman↗