Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Input Sensitivity”

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 1,603 records · Page 89

A Technique for Increasing the Sensitivity of a Solid-State Fission Probe

A small silicon p-n junction wafer, when coated with uranium 235, can be used as a compact fission probe for low power flux and power mapping. Because of the inverse relation between the magnitude of a neutron-induced fission pulse and the inherent capacitance of the detecting element (capacitance is proportional to area), the size, and hence the sensitivity, of the semiconductor detector has been limited. New developments in the field of semiconductor detectors have made it possible to fabricate large area detectors which are essentially free from the capacitance effect. However, preliminary results indicate that they are much more susceptible to radiation damage than the detectors described in this report and as such may not be suitable for flux mapping. increasing the sensitivity cannot be accomplished by simply fabricating a larger detector. It has been observed that by combining the silicon p-n junction wafers in a series configuration the capacitance effect can be bypassed, and a fission probe can be made with a resultant increase in sensitivity by a factor of ten while sustaining only a minor decrease in pulse height. Analysis further indicates that for n silicon wafers in series, if n(C(sub i)) + C(sub c)/C(sub b) less than 0.1 where C(sub i) and C(sub c) are the preamplifier input and cable capacitances, respectively, and C(sub b) is the junction capacitance of a single silicon wafer, there should be no substantial reduction in pulse height due to series circuitry.

Steinberg, Robert↗

Phase-Shifting Zernike Interferometer Wavefront Sensor

The canonical Zernike phase-contrast technique1,2,3,4 transforms a phase object in one plane into an intensity object in the conjugate plane. This is done by applying a static pi/2 phase shift to the central core (approx. lambda/D) of the PSF which is intermediate between the input and output planes. Here we present a new architecture for this sensor. First, the optical system is simple and all reflective. Second, the phase shift in the central core of the PSF is dynamic and or arbitrary size. This common-path, all-reflective design makes it minimally sensitive to vibration, polarization and wavelength. We review the theory of operation, describe the optical system, summarize numerical simulations and sensitivities and review results from a laboratory demonstration of this novel instrument

wavefront sensing↗

Periodic Material-Based Seismic Base Isolators for Small Modular Reactors

Important structures such as nuclear power plants should experience very little vibrations in order to maintain the safety of the nuclear facilities during seismic events. Protection of these structures by seismic isolators requires the isolators to provide seismic isolation irrespective of the direction of excitations and to develop a stable response without rocking. Conventional seismic isolation systems, however, can only fulfill partial requirements for the seismic protection. This study adopts the principle of periodic materials for seismic isolation. These materials exhibit unique properties of frequency band gaps, where incoming waves with frequencies falling inside the frequency band gaps are forbidden. The new isolation systems, known as periodic foundations, function both as a structural foundation to support the gravitational weight of the superstructure and also as a seismic isolator. The isolation mechanism, hypothetically, can easily fulfill all the requirements for seismic protection of the critical facilities. This study focuses on the design of one-dimensional (1D) and three-dimensional (3D) periodic foundations for seismic isolation of small modular reactor (SMR) buildings. The theoretical study was first conducted to reveal the fundamental behavior of periodic foundations. Then global sensitivity analysis was utilized to study the effect of design parameters on the frequency band gaps. In addition, simple straight-forward design equations based on the sensitivity analysis are proposed for the design of periodic foundations. Utilizing the proposed equations, prototypes of 1D and 3D periodic foundations were designed to isolate an SMR building. Scaled models were fabricated and experimentally tested to validate the hypothesis and design. The periodic foundation structural systems were tested under various input waves, including seismic waves, in the horizontal and vertical directions and the torsional mode. The shake table test results show that the periodic foundations can reduce the acceleration response of the SMR building up to 90% in the horizontal direction and the torsional mode. As much as a 40% response reduction in the vertical direction is also observed. Moreover, the periodic foundation-isolated structure exhibited stable response with negligible rocking on the structure. This study proved the capability of periodic foundations to enhance the seismic safety of critical structures.

42 ENGINEERING↗

Performance Assessment Program for the Savannah River Site Liquid Waste Facilities - 13610

The Liquid Waste facilities at the U.S. Department of Energy's (DOE) Savannah River Site (SRS) are operated by Liquid Waste Operations contractor Savannah River Remediation LLC (SRR). A separate Performance Assessment (PA) is prepared to support disposal operations at the Saltstone Disposal Facility and closure evaluations for the two liquid waste tank farm facilities at SRS, F-Tank Farm and H-Tank Farm. A PA provides the technical basis and results to be used in subsequent documents to demonstrate compliance with the pertinent requirements identified in operations and closure regulatory guidance. The Saltstone Disposal Facility is subject to a State of South Carolina industrial solid waste landfill permit and the tank farms are subject to a state industrial waste water permit. The three Liquid Waste facilities are also subject to a Federal Facility Agreement approved by the State, DOE and the Environmental Protection Agency (EPA). Due to the regulatory structure, a PA is a key technical document reviewed by the DOE, the State of South Carolina and the EPA. As the waste material disposed of in the Saltstone Disposal Facility and the residual material in the closed tank farms is also subject to reclassification prior to closure via a waste determination pursuant to Section 3116 of the Ronald W. Reagan National Defense Authorization Act of Fiscal Year 2005, the U.S. Nuclear Regulatory Commission (NRC) is also a reviewing agency for the PAs. Pursuant to the Act, the NRC also has a continuing role to monitor disposal actions to assess compliance with stated performance objectives. The Liquid Waste PA program at SRS represents a continual process over the life of the disposal and closure operations. When the need for a PA or PA revision is identified, the first step is to develop a conceptual model to best represent the facility conditions. The conceptual model will include physical dimensions of the closed system, both the engineered and natural system, and modeling input parameters associated with the modeled features, both initial values (at the time of facility closure) and degradation rates/values. During the development of the PA, evaluations are conducted to reflect not only the results associated with the best available information at the time but also to evaluate potential uncertainties and sensitivities associated with the modeled system. While the PA will reflect the modeled system results from the best available information, it will also identify areas for future work to reduce overall PA uncertainties moving forward. DOE requires a PA Maintenance Program such that work continues to reduce model uncertainties, thus bolstering confidence in PA results that support regulatory decisions. This maintenance work may include new Research and Development activities or modeling as informed by previous PA results and other new information that becomes available. As new information becomes available, it is evaluated against previous PAs and appropriate actions are taken to ensure continued confidence in the regulatory decisions. Therefore, the PA program is a continual process that is not just the development of a PA but seeks to incorporate new information to reduce overall model uncertainty and provide continuing confidence in regulatory decisions. (author)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Extensible Perturbed Parameter Ensemble for the Community Atmosphere Model Version 6

This paper documents the methodology and preliminary results from a Perturbed Parameter Ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2-CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations are analyzed and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks. The use of various emulators is explored in the multi- dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that using PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2-CAM6 PPE is publicly available, and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g. coupling to the land surface).

Machine learning↗

Exploring Sensitivity of ICF Outputs to Design Parameters in Experiments Using Machine Learning

We report building a sustainable burn platform in inertial confinement fusion (ICF) requires an understanding of the complex coupling of physical processes and the effects that key experimental design changes have on implosion performance. While simulation codes are used to model ICF implosions, incomplete physics and the need for approximations deteriorate their predictive capability. Identification of relationships between controllable design inputs and measurable outcomes can help guide the future design of experiments and development of simulation codes, which can potentially improve the accuracy of the computational models used to simulate ICF implosions. In this article, we leverage developments in machine learning (ML) and methods for ML feature importance/sensitivity analysis to identify complex relationships in ways that are difficult to process using expert judgment alone. We present work using random forest (RF) regression for prediction of yield, velocity, and other experimental outcomes given a suite of design parameters, along with an assessment of important relationships and uncertainties in the prediction model. We show that RF models are capable of learning and predicting on ICF experimental data with high accuracy, and we extract feature importance metrics that provide insight into the physical significance of different controllable design inputs for various ICF design configurations. These results can be used to augment expert intuition and simulation results for optimal design of future ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗

On approaches to analyze the sensitivity of simulated hydrologic fluxes to model parameters in the community land model

Here, effective sensitivity analysis approaches are needed to identify important parameters or factors and their uncertainties in complex Earth system models composed of multi-phase multi-component phenomena and multiple biogeophysical-biogeochemical processes. In this study, the impacts of 10 hydrologic parameters in the Community Land Model on simulations of runoff and latent heat flux are evaluated using data from a watershed. Different metrics, including residual statistics, the Nash-Sutcliffe coefficient, and log mean square error, are used as alternative measures of the deviations between the simulated and field observed values. Four sensitivity analysis (SA) approaches, including analysis of variance based on the generalized linear model, generalized cross validation based on the multivariate adaptive regression splines model, standardized regression coefficients based on a linear regression model, and analysis of variance based on support vector machine, are investigated. Results suggest that these approaches show consistent measurement of the impacts of major hydrologic parameters on response variables, but with differences in the relative contributions, particularly for the secondary parameters. The convergence behaviors of the SA with respect to the number of sampling points are also examined with different combinations of input parameter sets and output response variables and their alternative metrics. This study helps identify the optimal SA approach, provides guidance for the calibration of the Community Land Model parameters to improve the model simulations of land surface fluxes, and approximates the magnitudes to be adjusted in the parameter values during parametric model optimization.

58 GEOSCIENCES↗

Nuclear Data Uncertainty Quantification for Adjoint-Weighted Quantities via XGPT in Monte Carlo Codes

The Total Monte Carlo method (TMC) is an innovative approach to nuclear data uncertainty propagation. In TMC, a large number of independent ENDF files are randomly generated starting from nuclear models parameters and their uncertainties. These files are then processed with NJOY or similar codes to produce a set of formatted continuous-energy cross sections (i.e., ACE files), and independent Monte Carlo neutron transport simulations are run with these different ACE files as input. Finally the distributions of the considered response functions are obtained directly from the results of the Monte Carlo runs. Despite several improvements of the TMC method, this approach is often considered an inefficient way of propagating nuclear data uncertainties, as it requires several computationally expensive Monte Carlo runs. Covariance-based methods present important advantages for several applications, especially when different sources of uncertainties are to be investigated. For this reason, sensitivity analysis and uncertainty quantification capabilities are being implemented in several Monte Carlo neutron transport codes. On the other hand, the fact that no artificial multi-group covariance matrices are required and that higher moments of the responses distributions can be easily obtained, are definite advantages of the TMC method. In an effort to overcome the limitations of both approaches, a new method, called eXtended Generalized Perturbation Theory (XGPT) is under development. When performing nuclear data uncertainty propagation, such XGPT approach does not require generating discretized (i.e., multigroup) covariance matrices, and allows to estimate higher moments of the response probability distributions, via linear perturbation theory. One of the promising applications of the XGPT method involves its adoption in conjunction with the Proper Orthogonal Decomposition (POD) of nuclear data, to build a Reduced Order Model (ROM) for the approximation of thousands of independent TMC results in a single criticality source simulation. In the present work, this approach is employed for the first time to evaluate the impact of uncertainty in nuclear data on adjoint-weighted quantities. The method has been implemented in an extended version of the Serpent-2.1.24 Monte Carlo code and is demonstrated on two fast spectrum benchmarks, focusing on {sup 208}Pb nuclear data. A new method for propagating nuclear data uncertainties in criticality source Monte Carlo simulations has been developed. This eXtended Generalized Perturbation Theory (XGPT) approach allows to accurately reproduce reference Total Monte Carlo (TMC) uncertainty estimates in a single criticality source simulation. The present work aims to demonstrate that the proposed method can quantify nuclear data uncertainty based on adjoint-weighted quantities. The XGPT approach was applied to evaluate the impact of {sup 208}Pb nuclear data uncertainties on the effective prompt neutron lifetime (l{sub eff}) in two fast spectrum systems. The obtained XGPT results show a good agreement with the reference TMC results, obtained from independent Monte Carlo runs with {sup 208}Pb nuclear data from random TENDL-2013 evaluations. The XGPT method, implemented in a extended version of the Serpent Monte Carlo code and adopted along with the Proper Orthogonal Decomposition (POD) of the nuclear data, is able to provide an accurate reduced order model to approximate thousands TMC results within a single run, ensuring a large saving in computational requirements. The effect of {sup 208}Pb nuclear data uncertainties appears to be rather large. This is consistent with previous studies that adopted the same input random evaluations data from TENDL. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radionuclides in bats using a contaminated pond on the Nevada National Security Site, USA

In this study, perched groundwater percolating through radionuclide contamination in the E Tunnel Complex on the Nevada National Security Site, formerly the Nevada Test Site, emerges and is stored in a series of ponds making it available to wildlife, including bats. Since many bat species using the ponds are considered sensitive or protected/regulated and little information is available on dose to bats from radioactive water sources, bats were sampled to determine if the dose they were receiving exceeded the United States Department of Energy dose limit of 1.0E-3 Gy/day. Radionuclide concentrations in water, sediment, and flying insects were also measured as input parameters to the dose rate model and to examine trophic level relationships. The RESRAD-Biota model was used to calculate dose rates to bats using different screening levels. Efficacy of RESRAD-Biota and suggested improvements are discussed. Finally, dose to bats foraging and drinking at these ponds is well below the dose limit set to protect terrestrial biota populations.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Newly Released Capabilities in the Distributed-Memory SuperLU Sparse Direct Solver

We present the new features available in the recent release of SuperLU_DIST, Version 8.1.1. SuperLU_DIST is a distributed-memory parallel sparse direct solver. The new features include (1) a 3D communication-avoiding algorithm framework that trades off inter-process communication for selective memory duplication, (2) multi-GPU support for both NVIDIA GPUs and AMD GPUs, and (3) mixed-precision routines that perform single-precision LU factorization and double-precision iterative refinement. Apart from the algorithm improvements, we also modernized the software build system to use CMake and Spack package installation tools to simplify the installation procedure. Throughout the article, we describe in detail the pertinent performance-sensitive parameters associated with each new algorithmic feature, show how they are exposed to the users, and give general guidance of how to set these parameters. We illustrate that the solver’s performance both in time and memory can be greatly improved after systematic tuning of the parameters, depending on the input sparse matrix and underlying hardware.

97 MATHEMATICS AND COMPUTING↗

Synthesis techniques for insensitive aircraft control systems

Comparative evaluation of seven insensitive synthesis techniques applied to the control of the longitudinal axis and symmetric flexure modes of the C-5A aircraft is described. These techniques have the common goal of achieving insensitivity to parameter uncertainties and unmodeled dynamics. For this comparison, the 'complete' model was assumed to consist of rigid body dynamics, fifteen flexure modes, unsteady aerodynamics, gust dynamics, and first order actuators for two control inputs. Thirteen flexure modes and the unsteady aerodynamics were treated as neglected dynamics. Dynamic pressure, structural frequency and damping, and the stability derivative were treated as uncertain parameters. Design specifications consisted of a transient response criterion, statistical response criteria, eigenvalue criteria and stability margins. A nominal controller designed without regard to sensitivity served as a benchmark for performance comparisons of the insensitive controllers.

Harvey, C. A.↗

Absolute calibration technique for broadband ultrasonic transducers

Calibrating an ultrasonic transducer can be performed with a reduced number of calculations and testing. A wide-band pulser is connected to an ultrasonic transducer under test to generate ultrasonic waves in a liquid. A single frequency is transmitted to the electrostatic acoustic transducer (ESAT) and the voltage change produced is monitored. Then a broadband ultrasonic pulse is generated by the ultrasonic transducer and received by the ESAT. The output of the ESAT is amplified and input to a digitized oscilloscope for fast Fourier transform. The resulting plot is normalized with the monitored signal from the single frequency pulse. The plot is then corrected for characteristics of the membrane and diffraction effects. The transfer function of the final plot is determined. The transfer function gives the final sensitivity of the ultrasonic transducer as a function of frequency. The advantage of the system is the speed of calibrating the transducer by a reduced number of measurements and removal of the membrane and diffraction effects.

Yost, William T.↗

Fitting a Two-Component Scattering Model to Polarimetric SAR Data from Forests

Two simple scattering mechanisms are fitted to polarimetric synthetic aperture radar (SAR) observations of forests. The mechanisms are canopy scatter from a reciprocal medium with azimuthal symmetry and a ground scatter term that can represent double-bounce scatter from a pair of orthogonal surfaces with different dielectric constants or Bragg scatter from a moderately rough surface, which is seen through a layer of vertically oriented scatterers. The model is shown to represent the behavior of polarimetric backscatter from a tropical forest and two temperate forest sites by applying it to data from the National Aeronautic and Space Agency/Jet Propulsion Laboratory's Airborne SAR (AIRSAR) system. Scattering contributions from the two basic scattering mechanisms are estimated for clusters of pixels in polarimetric SAR images. The solution involves the estimation of four parameters from four separate equations. This model fit approach is justified as a simplification of more complicated scattering models, which require many inputs to solve the forward scattering problem. The model is used to develop an understanding of the ground-trunk double-bounce scattering that is present in the data, which is seen to vary considerably as a function of incidence angle. Two parameters in the model fit appear to exhibit sensitivity to vegetation canopy structure, which is worth further exploration. Results from the model fit for the ground scattering term are compared with estimates from a forward model and shown to be in good agreement. The behavior of the scattering from the ground-trunk interaction is consistent with the presence of a pseudo-Brewster angle effect for the air-trunk scattering interface. If the Brewster angle is known, it is possible to directly estimate the real part of the dielectric constant of the trunks, a key variable in forward modeling of backscatter from forests. It is also shown how, with a priori knowledge of the forest height, an estimate for the attenuation coefficient of the canopy can be obtained directly from the multi-incidence-angle polarimetric observations. This attenuation coefficient is another key variable in forward models and is generally related to the canopy density.

temperate forests↗

Observing System Variations Effect on Reanalyses

Reanalyses integrate multitudes of satellite and conventional observations data assimilation and numerical weather prediction. The result is that many disparate observation platforms, discontinuous in space and time, lead to complete and consistent representations the state of the weather. The component also provides physical fields rarely or never observed. However, the numerical model bias is continuously being corrected by the observational analysis, and this bias changes as variations in the observations occur. NASA's Modern Era Retrospective-analysis for Research and Applications (MERRA) sensitivity to variations in the observing systems are explored. Specifically, we will evaluate the water budget and transport processes as they relate to the advent of SSM/I and AMSU-A radiance assimilation, and an additional case of radiosonde station that exhibits a dramatic shift in mean water states. The MERRA input observation data, now available online, is used to explore these variations.

Bosilovich, Michael G.↗

Long‐Term Support of an Active Subglacial Hydrologic System in Southeast Greenland by Firn Aquifers

The state of the subglacial hydrologic system, which can modify ice motion, is sensitive to the volume and rate of meltwater reaching it. Bare-ice regions rapidly transport meltwater to the bed via moulins, while in certain accumulation-zone regions, meltwater first flows through firn aquifers, which can introduce a substantial delay. We use a subglacial hydrological model forced with idealized meltwater input scenarios to test the effect of this delay on subglacial hydrology. We find that addition of firn-aquifer water to the subglacial system elevates the inland subglacial water pressure while reducing water pressure and enhancing subglacial channelization near the terminus. This effect dampens seasonal variations in subglacial water pressure and may explain regionally anomalous ice-velocity patterns observed in Southeast Greenland. As surface melt rates increase and firn aquifers expand inland, it is crucial to understand how inland drainage of meltwater affects the evolution of the subglacial hydrologic system.

Kristin Poinar↗

Estimating Future Changes of Energy Demand for Heating and Cooling Buildings at NASA Centers GC23J-1198

With its unique and trusted earth observations, NASA is a critical source in informing decisions that will help achieve the U.S. goal of Net-Zero Greenhouse Gas (GHG) Emissions by 2050. NASA’s Prediction of Worldwide Energy Resource (POWER) project facilitates the use of NASA Earth Science data holdings within the energy, agricultural, and building heating/cooling design industries. POWER packages solar and meteorological data from several NASA projects in a user friendly GIS-enabled web services system (https://power.larc.nasa.gov). As part of the development of new data products to support the energy and building heating/cooling design communities, we estimate the changes in energy required to heat and cool buildings in the future climate at 14 different NASA site locations spread throughout the continental United States, as projected by CMIP6 climate models under different emissions scenarios. Bias-corrected downscaled time series of meteorological variables are taken from NASA Earth Exchange (NEX) Global Daily Downscaled Projections (GDDP-CMIP6) downscaled climate model data. The spread between the different model projections is accounted for by analyzing both the ensemble average of 22 CMIP6 models and 6 representative models with different climate sensitivities and different interannual variability. Changes in energy use are estimated in two ways. First, changes in the total annual heating and cooling degree days (HDD and CDD, respectively) are calculated relative to the current climate. This is done at all 14 sites. Second, the downscaled time series are used as inputs into RETScreen(R), a clean energy management decision tool, to give an estimate of heating/cooling energy use for a typical office building. We use this estimation method with model data at Langley Research Center. In the next 50 years, the annual total of HDD (CDD) is projected to decrease by 8-38% (increase by 5-28%) at all sites, with the increase in CDD typically a larger magnitude the decrease in HDD. For a typical small office building at Langley Research Center, the amount of energy needed to cool increases by 33-54% and the amount of energy to heat decreases by 29-40%. POWER is working to develop long term climate data services based on these results to include in future data products to provide to users.

Bradley M. Hegyi↗

Mathematical Model of a Regenerative Fuel Cell for System Optimization

This thesis developed a system-level optimization model of a regenerative fuel cell (RFC) system for long-duration, off-world energy storage applications. Prior RFC design studies have typically been limited to reduced parameter sets and simplified constraints due to computational limitations relative to the number of relevant degrees of freedom. As a result, important nonlinear interactions between subsystems have not been fully captured. This work began to address that gap by developing a higher-fidelity, nonlinear optimization framework that incorporates a broader set of design variables and coupled constraints, enabling a multidimensional model that captures the coupled behavior of RFC subsystems and demonstrates the feasibility of applying optimization to such systems. An expanded system-level optimization approach was established that captures interactions between electrochemical performance, structural requirements, and storage design. This enabled a more comprehensive evaluation of trade-offs than conventional formulations. The model integrates four coupled subsystems: a fuel cell, an electrolyzer, reactant gas, and high-pressure storage tanks, and was formulated to accommodate a wide range of mission parameters, including operational time and required output power. It incorporates constraints on available solar array power, reactant mass balance between production and consumption, and pressure-dependent storage requirements. To enable reliable convergence, the optimization problem was reformulated to reduce dimensionality and improve numerical stability, with subsystem models organized for efficient evaluation. Problem dimensionality was reduced by consolidating lower-level design variables into higher-level representative quantities, and subsystem behavior was evaluated within the optimization loop. A multi-start initialization strategy was employed to mitigate sensitivity to local minima and improve solution quality, while nonlinear relationships were solved using robust numerical methods. The results showed that convergence was achieved across a range of required output power values. Specific energy reached a maximum at a critical mission power level, where the electrolyzer power matched the available solar input and operated near its voltage and current density limits. Beyond this point, further increases in required power resulted in less mass-efficient operation, increasing total system mass and reducing overall performance. The developed model represents an advancement in RFC system-level optimization by enabling analysis of a broader and more tightly coupled design space than previous considerations. While convergence behavior and computational cost remain challenges, the methods introduced improve solvability and allow inclusion of additional design variables with minimal loss of physical fidelity. However, the numerical results should not be interpreted as definitive design recommendations, as the model includes simplifying assumptions and omits several higher-order effects. Future work should extend this framework by incorporating additional subsystems and loss mechanisms, such as thermal management, parasitic power consumption, and reactant losses, to improve fidelity and ensure more representative design conclusions.

Electrochemistry↗