SEARCH · Search NASA

Results for “Distributed model predictive control”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Computations of Flow over a Hump Model Using Higher Order Method with Turbulence Modeling

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k - omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially. nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as those observed in the experiments. The computed separation regions are much longer than those observed experimentally. However, the percentage reduction in the separation region in the steady suction case is closer to what was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than the measured values pointing towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Balakumar, P.↗

NASA NTRS2018

San Francisco Bar Pilot Fatigue Study

The purpose of this study was to evaluate the extent of fatigue among San Francisco Bar Pilots (Maritime Pilots) and its potential impact on safety, and to make recommendations concerning how the risk of fatigue could be managed. Information was gathered via a literature review, observations of Bar Pilots at work, surveys, a task analysis, and an analysis of dispatch records.The work of San Francisco Bar Pilots involves an unusual mix of activities and job demands. Their work calls for situational awareness, reasoning, communication, and perceptual abilities comparable to those required by airline pilots and air traffic controllers. Errors can have severe consequences for public safety and the environment, as well as significant financial costs. Fatigue is increasingly recognized as a hazard that must be managed by the transportation industry. The reduced sleep quality and quantity experienced by personnel who work at night, in conjunction with human circadian rhythms can lead to an operationally significant level of cognitive impairment. The cognitive impact of fatigue includes functions that are critical to safe maritime piloting, such as vigilance, judgment, reaction time and communication. The surveys distributed to Bar Pilots did not uncover evidence of widespread fatigue. The survey results also suggested that fatigue was not a major concern of Bar Pilots. Bar Pilots had overall low scores on the subjective fatigue measures used in the survey, and generally assessed the safety risk due to fatigue as low. Compared to air traffic controllers, Bar Pilots gave significantly lower ratings on questions concerning the prevalence and impact of fatigue. The application of fatigue modeling software to Bar Pilot dispatch records identified that in most cases, pilot's cognitive effectiveness was predicted to be acceptable during their duty periods. However, these results could not be verified with objective data.The study identified a number of fatigue issues that deserve attention. These include Bar Pilot work periods that frequently infringe on the circadian low, consecutive work periods without a significant break, consecutive periods of night work, unpredictable work schedules, start time variability, the potential for sleep inertia, and the number of pilots on the board at any given time.

fatigue↗

NASA NTRSDate not supplied

Summary of the First AIAA Stability and Control Prediction Workshop

Results from the First AIAA Stability and Control Prediction Workshop are summarized in this paper. The workshop series was developed in support of three primary objectives: (1) to establish best practices for the prediction of stability and control derivatives using industry-standard Computational Fluid Dynamics (CFD) solvers, (2) to provide an impartial forum for evaluating the effectiveness of Reynolds-averaged-Navier-Stokes- and Detached-Eddy-Simulation-based modeling techniques, and (3) to identify areas in need of additional research and development. To address these objectives, the inaugural workshop focused on generating computational aerodynamic predictions for the ONERA version of the NASA/Boeing Common Research Model. The configuration includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. While longitudinal wind tunnel test data for the model had been previously documented, unpublished lateral test data at small sideslip angles provided a unique opportunity for participants to generate blind computational predictions for wind tunnel data comparisons. Additional test cases included assessments of the Mach number effect on static lateral/directional stability and the impact of the wind tunnel sting on longitudinal stability. Participants were invited to generate solutions using two workshop-provided series of structured overset and unstructured grids, in addition to any participant custom grids of interest. Total and component-level breakdowns for the force and moment coefficients for each test case are presented, as well as sectional pressure distribution data on the wing and tail components, to assess the agreement between several different Reynolds-averaged Navier-Stokes CFD solvers.

CFD↗

DOE OSTI2025

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1994

Neonatal rat heart cells cultured in simulated microgravity

In vitro characteristics of cardiac cells cultured in simulated microgravity are reported. Tissue culture methods performed at unit gravity constrain cells to propagate, differentiate, and interact in a two dimensional (2D) plane. Neonatal rat cardiac cells in 2D culture organize predominantly as bundles of cardiomyocytes with the intervening areas filled by non-myocyte cell types. Such cardiac cell cultures respond predictably to the addition of exogenous compounds, and in many ways they represent an excellent in vitro model system. The gravity-induced 2D organization of the cells, however, does not accurately reflect the distribution of cells in the intact tissue. We have begun characterizations of a three-dimensional (3D) culturing system designed to mimic microgravity. The NASA designed High-Aspect-Ratio-Vessel (HARV) bioreactors provide a low shear environment which allows cells to be cultured in static suspension. HARV-3D cultures were prepared on microcarrier beads and compared to control-2D cultures using a combination of microscopic and biochemical techniques. Both systems were uniformly inoculated and medium exchanged at standard intervals. Cells in control cultures adhered to the polystyrene surface of the tissue culture dishes and exhibited typical 2D organization. Cells in cultured in HARV's adhered to microcarrier beads, the beads aggregated into defined clusters containing 8 to 15 beads per cluster, and the clusters exhibited distinct 3D layers: myocytes and fibroblasts appeared attached to the surfaces of beads and were overlaid by an outer cell type. In addition, cultures prepared in HARV's using alternative support matrices also displayed morphological formations not seen in control cultures. Generally, the cells prepared in HARV and control cultures were similar, however, the dramatic alterations in 3D organization recommend the HARV as an ideal vessel for the generation of tissue-like organizations of cardiac cells in simulated microgravity.

Akins, Robert E.↗

NASA NTRS2004

Computations of Flow Over a Hump Model Using Higher Order Method With Turbulence Modeling

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k-omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as they are observed in the experiments. The computed separation regions are much longer than that are observed. However, the percentage reduction in the separation region in the steady suction case is closer to that was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than that are measured and it points towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Balakumar, Ponnampalam↗

NASA NTRS2012

A Bayesian Framework for Reliability Analysis of Spacecraft Deployments

Deployable subsystems are essential to mission success of most spacecraft. These subsystems enable critical functions including power, communications and thermal control. The loss of any of these functions will generally result in loss of the mission. These subsystems and their components often consist of unique designs and applications for which various standardized data sources are not applicable for estimating reliability and for assessing risks. In this study, a two stage sequential Bayesian framework for reliability estimation of spacecraft deployment was developed for this purpose. This process was then applied to the James Webb Space Telescope (JWST) Sunshield subsystem, a unique design intended for thermal control of the Optical Telescope Element. Initially, detailed studies of NASA deployment history, "heritage information", were conducted, extending over 45 years of spacecraft launches. This information was then coupled to a non-informative prior and a binomial likelihood function to create a posterior distribution for deployments of various subsystems uSing Monte Carlo Markov Chain sampling. Select distributions were then coupled to a subsequent analysis, using test data and anomaly occurrences on successive ground test deployments of scale model test articles of JWST hardware, to update the NASA heritage data. This allowed for a realistic prediction for the reliability of the complex Sunshield deployment, with credibility limits, within this two stage Bayesian framework.

Evans, John W.↗

DOE OSTI2021

A Triple-Moment Representation of Ice in the Predicted Particle Properties (P3) Microphysics Scheme

In the original Predicted Particle Properties (P3) bulk microphysics scheme, all ice-phase hydrometeors are represented by one or more “free” ice categories, where the physical properties evolve smoothly through changes to the four prognostic variables (per category), and with a two-moment representation of the particle size distribution. As such, the spectral dispersion cannot evolve independently, which thus results in limitations in representation of ice—in particular, hail—due to necessary constraints in the scheme to prevent excessive gravitational size sorting. To overcome this, P3 has been modified to include a three-moment representation of the size distribution of each ice category through the addition of a fifth prognostic variable, the sixth moment of the size distribution. The details of the three-moment ice parameterization in P3 are provided. The behavior of the modified scheme, with the single-ice-category configuration, is illustrated through simulations in a simple 1D kinematic model framework as well as with near large-eddy-resolving (250-m grid spacing) 3D simulations of a hail-producing supercell. It is shown that the three-moment ice configuration controls size sorting in a physically based way and leads to an improved capacity to simulate large, heavily rimed ice (hail), including mean and maximum sizes and reflectivity, and thus an overall improvement in the representation of ice-phase particles in the P3 scheme.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Thermal Modeling and Simulation of the Packed-bed Thermal Energy Storage Combined with INL Thermal Energy Distribution System

Dynamic Energy Transport and Integration Lab (DETAIL) at Idaho National Laboratory is to support experimental demonstration and validation research on Nuclear-Renewable Hybrid Energy System [1]. The Thermal Energy Distribution System (TEDS) is a thermal-hydraulic flow loop in DETAIL with its own dedicated control system to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit. Among various TES options, the packed-bed TES is adopted in TEDS because it offers a low-cost single-tank thermal storage option compared to the traditional two-tank TES. However, since the TES tank is filled with granular fillers having different thermophysical properties from those of TES tank wall, there is a potential thermo-mechanical issue to be carefully addressed like thermal ratcheting which may pose a significant design concern for the packed-bed TES tank. Thermal ratcheting is a thermomechanical process caused by the repeated rearrangement of granular filler inside a TES tank during continuous thermal cycling operation of the packed-bed TES system. If the thermally induced stress exceeds the yield strength of a TES tank wall, catastrophic consequences may happen like rupture of the TES tank. Thus, it is important to understand the phenomenon to ensure the robust operation of the packed-bed TES tank. Given that thermal ratcheting is caused by complex interaction of thermal transport in the porous bed and solid mechanics, the accurate prediction of transient thermal behavior of the packed-bed TES tank, which is the focus of this paper, is critical to the reliable thermal ratcheting analysis. This paper discusses the numerical modeling, simulation, and validation studies that are ongoing at INL to investigate the transient thermal behavior of the packed-bed TES. Of particular concern is the transient thermal process occurring in the packed-bed TES unit that is operated in conjunction with the INL TEDS. The main goal of this research is three-fold: (i) provide preliminary insights into the transient thermal behavior of the TEDS TES tank, (ii) support the thermal measurement and validation plan for TEDS experiment, and (iii) provide transient thermal boundary conditions to support the reliable thermal ratcheting analysis of the TEDS TES tank. For the transient thermal modeling and analysis, a CFD model was developed, and the validity of the modeling approach was examined via comparing the numerical simulation results with the experimental data obtained from various design characteristics of packed-bed TES tanks. Then, the present modeling method was applied for the transient thermal analysis of the TEDS TES tank, and the results are discussed along with the potential improvement of data acquisition strategy for the future TEDS experiments for more precise validation study.

25 ENERGY STORAGE↗

DOE OSTI2024

A Modeling framework for flocculated cohesive sediment transport in the current bottom boundary layer

Cohesive sediment transport, where its settling velocity is controlled by the flocculation process, is a crucial component in determining biochemical cycles, fate of pollutants, and morphodynamics in many aquatic ecosystems. In this study, a modeling framework is presented to investigate how flocculation influences cohesive sediment transport in the current bottom boundary layer in dilute conditions, consistent with the calibration range of the flocculation model. From a local analysis of floc dynamics in homogenous turbulence, we identify that the floc size distribution is mainly controlled by floc cohesion and yield strength. The uncertainty in fractal dimension plays a minor role for the floc size but it influences the resulting floc density and settling velocity. The transport analysis in the current boundary layer shows that the flocculation process alters the vertical distribution of the settling velocity and hence the sediment concentration with a strong dependence on cohesion, floc yield strength, and floc structure. When the flocs are more susceptible to breaking, a well-mixed concentration profile is obtained. In contrast, for flocs with higher cohesion or yield strength, higher concentration with a sharp gradient is observed close to the bed. Overall, the settling velocity exhibits a low vertical variability within 20 % of the depth-averaged value except near the bed. Further, this suggests that using a depth-averaged settling velocity yields acceptable predictions of the sediment concentration profiles, especially for flocs with lower cohesion.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Event-Driven Predictive Approach for Real-Time Volt/VAR control with CVR in solar PV rich Active Distribution Network

The focus of this paper is on analyzing the impact of conservation voltage reduction in the presence of active devices such as solar photovoltaic (PV) and developing controls that leverage these distributed energy resources. An event-driven predictive approach for real-time volt/volt-ampere reactive (VAR) optimization, along with local two-level adaptive volt/VAR droop-based control algorithm for advanced distribution management systems, is introduced. The methodology covers aggregated and autonomous controls under different timescale operations, including the impact and effect of unpredicted events such as cloud transients on PV power production. In addition, the control schemes include the uncertainties in PV power generation and load power demand. The proposed methodology is validated in a real-time framework using the real-time digital simulator platform through co-simulation with models based on Python and OpenDSS (Open Distribution System Simulator). The developed methodology is tested on the modified IEEE 123-feeder test system. The results reveal that the proposed methodology works well in the presence of high penetrations of PV power, produces significant energy savings, and mitigates over-/undervoltage problems.

14 SOLAR ENERGY↗

DOE OSTI2022

Modeling the hydrodynamic impact on the tool influence function during hemispherical subaperture optical polishing

To fabricate high-precision and accurate optics relative to the optical design surface, a high level of deterministic control of material removal (i.e., the tool influence function, TIF) during subaperture tool polishing is required. In this study, a detailed analysis of the pressure distribution, which is a key component of the TIF, has been performed using finite element analysis to couple together solid mechanics and fluid dynamics. Modeling experimental parameters of recently published work reveals that, when considering tool deformation, which in turn influences the fluid film thickness between the tool and workpiece, the effective pressure profile has a flat-top distribution. This flat-top pressure profile differs from the parabolic pressure distributions predicted by Hertzian mechanics. Furthermore, the shear contribution is shown here to be a key contributor to material removal, inducing the removal at the periphery of the contact edge and even outside the generally accepted contact area. Finally, the simulated fluid velocities provide evidence of mixed-mode contact polishing, supporting recent experimental findings that also suggest that onset of hydroplaning contributions lead to material removal drop-off.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Load Flexibility Analysis of Residential HVAC and Water Heating and Commercial Refrigeration

The emergence of connected devices and the Internet of Things has opened the door to the concept of widespread use of building equipment for load-side flexibility. This stems from the increased use of distributed and clean energy resources for electric generation and the growing demand for a more resilient and reliable grid. A significant challenge associated with the use of building equipment as a distributed energy resource is quantifying the magnitude of power or energy flexibility. Unlike batteries that have well-defined and relatively constant capacity and maximum charge/discharge rates, building equipment often has characteristics that vary based on external factors (e.g., weather or occupant behavior). This makes quantifying the availability and flexibility of these resources difficult. Controlling a fleet of equipment while minimizing any adverse impacts to the occupants (e.g., discomfort or food spoilage) is a nontrivial matter. Although there are many different control methodologies available, this study focuses on two: priority-based control (PBC) and model-predictive control (MPC). These were selected because they represent two ends of the control spectrum. PBC is a relatively simple control strategy that requires only the current state of the devices to be controlled and an estimate of how much power each device consumes when turned on. MPC is much more complex and requires a model of the system and forecasts of any disturbances that influence the system. Because MPC has a predictive element, it can prepare for a desired response before the event and therefore is expected to provide a more optimal control relative to PBC. Both MPC and PBC are compared with a baseline simulation that uses feedback control typical of the thermostatically controlled devices being evaluated. This report investigates the load flexibility of three building systems: residential air conditioners, residential water heaters, and commercial refrigeration. Air heating and cooling and water heating account for 45% of residential electricity use, making them good candidates for providing load flexibility. Refrigeration in commercial buildings accounts for 16% of commercial building electricity use, making it another good candidate. Simplified models of each system were generated and validated based on measured data from actual equipment. These models were used to simulate the thermal and power response of the systems using PBC and MPC for a fleet of devices trying to achieve a variety of different responses, including hourly load shed, hourly load up, and daily peak demand reduction. Significant findings from each piece of equipment follow.

42 ENGINEERING↗

NASA NTRS1997

Neonatal rat heart cells cultured in simulated microgravity

In vitro characteristics of cardiac cells cultured in simulated microgravity are reported. Tissue culture methods performed at unit gravity constrain cells to propagate, differentiate, and interact in a two-dimensional (2D) plane. Neonatal rat cardiac cells in 2D culture organize predominantly as bundles of cardiomyocytes with the intervening areas filled by nonmyocyte cell types. Such cardiac cell cultures respond predictably to the addition of exogenous compounds, and in many ways they represent an excellent in vitro model system. The gravity-induced 2D organization of the cells, however, does not accurately reflect the distribution of cells in the intact tissue. We have begun characterizations of a three-dimensional (3D) culturing system designed to mimic microgravity. The NASA-designed High-Aspect Ratio Vessel (HARV) bioreactors provide a low shear environment that allows cells to be cultured in static suspension. HARV-3D cultures were prepared on microcarrier beads and compared to control-2D cultures using a combination of microscopic and biochemical techniques. Both systems were uniformly inoculated and medium exchanged at standard intervals. Cells in control cultures adhered to the polystyrene surface of the tissue culture dishes and exhibited typical 2D organization. Cells cultured in HARVs adhered to microcarrier beads, the beads aggregated into defined clusters containing 8 to 15 beads per cluster, and the clusters exhibited distinct 3D layers: myocytes and fibroblasts appeared attached to the surfaces of beads and were overlaid by an outer cell type. In addition, cultures prepared in HARVs using alternative support matrices also displayed morphological formations not seen in control cultures. Generally, the cells prepared in HARV and control cultures were similar; however, the dramatic alterations in 3D organization recommend the HARV as an ideal vessel for the generation of tissuelike organization of cardiac cells in vitro.

Non-NASA Center↗

DOE OSTI2025

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗

DOE OSTI2020

Better together: Elements of successful scientific software development in a distributed collaborative community

Many scientific disciplines rely on computational methods for data analysis, model generation, and prediction. Implementing these methods is often accomplished by researchers with domain expertise but without formal training in software engineering or computer science. This arrangement has led to underappreciation of sustainability and maintainability of scientific software tools developed in academic environments. Some software tools have avoided this fate, including the scientific library Rosetta. We use this software and its community as a case study to show how modern software development can be accomplished successfully, irrespective of subject area. Rosetta is one of the largest software suites for macromolecular modeling, with 3.1 million lines of code and many state-of-the-art applications. Since the mid 1990s, the software has been developed collaboratively by the RosettaCommons, a community of academics from over 60 institutions worldwide with diverse backgrounds including chemistry, biology, physiology, physics, engineering, mathematics, and computer science. Developing this software suite has provided us with more than two decades of experience in how to effectively develop advanced scientific software in a global community with hundreds of contributors. Here we illustrate the functioning of this development community by addressing technical aspects (like version control, testing, and maintenance), community-building strategies, diversity efforts, software dissemination, and user support. We demonstrate how modern computational research can thrive in a distributed collaborative community. The practices described here are independent of subject area and can be readily adopted by other software development communities

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2024

Computational Analysis of the X-57 Maxwell Airplane, the Landing Configuration with High-Lift Blowing and Aileron Deflections (Preliminary Fuselage)

The X-57 Maxwell is an all-electric airplane with a distributed electric propulsion system used for a high-lift system at takeoff and landing conditions. The Kestrel and USM3D flow solvers were used at NASA Langley to investigate the performance of the X-57 Maxwell in the development of an aerodynamic database. The configuration investigated in this paper had a 30° flap deflection, and in addition the pilot’s right aileron were deflected and undeflected in different cases. The solutions were computed at an airspeed of 58 KEAS, for an altitude of 6000 feet, and a flight Reynolds number of 0.588e+06 per foot. To evaluate the high-lift distributed electric propulsion system without aileron deflections, the solutions were computed for an angle-of-attack sweep from −2° to 20°, with the high-lift propellers blowing and the cruise propellers excluded from the simulation. To investigate the aileron effectiveness, the solutions were computed at angles of attack of −2° and 14°, for aileron deflections from −25° to 18° with the high-lift propellers blowing and the cruise propellers windmilling at idle-power. The high-lift propellers and the idle-power cruise propellers were modeled with an actuator disk. Results show negligible differences in lift, drag and pitching moment whether the idle-powered cruise propellers were included or excluded from the simulation. In general, the Kestrel and USM3D codes compared well for lift, drag and pitching moment for the landing configuration with no aileron control. The Kestrel code, using the Spalart-Allmaras turbulence model with the rotation corrections terms, predicted an increasing lift coefficient up to maximum lift coefficient of 4.65 at a 15° angle of attack. The USM3D code, using the Spalart-Allmaras turbulence model turbulence model with the Quadratic Constitutive Relation predicted an increasing lift coefficient up to maximum lift coefficient of 4.5 at a 12° angle of attack, with the lift remaining constant through a 14◦ angle of attack. A possible difference in lift coefficient between the codes for high angles of attack may result from the different available options used with the standard Spalart-Allmaras turbulence model. The Kestrel code predicted more drag across the range of angle of attack than the USM3D code. The codes compared well for pitching moment coefficient across the range of angle of attack. The Kestrel and USM3D codes compared well for aileron effectiveness at a 2° angle of attack. The Kestrel code predicts better aileron effectiveness than USM3D at a 14° angle of attack. There is more flow separation in the region outboard of the last high-lift nacelle for the USM3D solutions at a 14° angle of attack, which diminishes the ability of the aileron to be effective.

X-57 Maxwell Airplane↗

DOE OSTI2022

Predicted Impacts of Pt and Ionomer Distributions on Low-Pt-Loaded PEMFC Performance

Low-cost, high performance proton exchange membrane fuel cells (PEMFCs) have been difficult to develop due to limited understanding of coupled processes in the cathode catalyst layer (CCL). Low-Pt-loaded PEMFCs suffer losses beyond those predicted solely due to reduced catalyst area. Although consensus links these losses to thin ionomer films in the CCL, a precise mechanistic explanation remains elusive. In this publication, we present a physically based PEMFC model with novel structure-property relationships for thin-film Nafion, validated against PEMFC data with low Pt loading. Results suggest that flooding exacerbates kinetic limitations in low-loaded PEMFCs, shifting the Faradaic current distribution. As current density increases, protons travel further into the CCL, resulting in higher Ohmic overpotentials. We also present a parametric study of CCL design parameters. We find that graded Pt and ionomer loadings reduce Ohmic losses and flooding, but individually do not provide significant improvements. However, a dual-graded CCL (i.e., graded Pt and ionomer) is predicted to significantly improve the maximum power density and limiting current compared to uniformly loaded CCLs. This work highlights the importance of accurate transport parameters for thin-film Nafion and provides a pathway to low-cost PEMFCs via precise control of CCL microstructures.

08 HYDROGEN↗