SEARCH · Engineering Papers

Results for “optimal start”

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 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS2009

jFuzz: A Concolic Whitebox Fuzzer for Java

We present jFuzz, a automatic testing tool for Java programs. jFuzz is a concolic whitebox fuzzer, built on the NASA Java PathFinder, an explicit-state Java model checker, and a framework for developing reliability and analysis tools for Java. Starting from a seed input, jFuzz automatically and systematically generates inputs that exercise new program paths. jFuzz uses a combination of concrete and symbolic execution, and constraint solving. Time spent on solving constraints can be significant. We implemented several well-known optimizations and name-independent caching, which aggressively normalizes the constraints to reduce the number of calls to the constraint solver. We present preliminary results due to the optimizations, and demonstrate the effectiveness of jFuzz in creating good test inputs. The source code of jFuzz is available as part of the NASA Java PathFinder. jFuzz is intended to be a research testbed for investigating new testing and analysis techniques based on concrete and symbolic execution. The source code of jFuzz is available as part of the NASA Java PathFinder.

Jayaraman, Karthick↗

NASA NTRS2000

Development of a Battery-Free Solar Refrigerator

Recent technology developments and a systems engineering design approach have led to the development of a practical battery-free solar refrigerator as a spin-off of NASA's aerospace refrigeration research. Off-grid refrigeration is a good application of solar photovoltaic (PV) power if thermal storage is incorporated and a direct connection is made between the cooling system and the PV panel. This was accomplished by integrating water as a phase-change material into a well insulated refrigerator cabinet and by developing a microprocessor based control system that allows direct connection of a PV panel to a variable speed compressor. This second innovation also allowed peak power-point tracking from the PV panel and elimination of batteries from the system. First a laboratory unit was developed to prove the concept and then a commercial unit was produced and deployed in a field test. The laboratory unit was used to test many different configurations including thermoelectric, Stirling and vapor compression cooling systems. The final configuration used a vapor compression cooling cycle, vacuum insulation, a passive condenser, an integral evaporator/ thermal storage tank, two 77 watt PV panels and the novel controller mentioned above. The system's only moving part was the variable speed BD35 compressor made by Danfoss. The 365 liter cabinet stayed cold with as little as 274 watt-hours per day average PV power. Battery-free testing was conducted for several months with very good results. The amount of thermal storage, size of compressor and power of PV panels connected can all be adjusted to optimize the design for a given application and climate. In the commercial unit, the high cost of the vacuum insulated refrigerator cabinet and the stainless steel thermal storage tank were addressed in an effort to make the technology commercially viable. This unit started with a 142 liter, mass-produced chest freezer cabinet that had the evaporator integrated into its inner walls. Its compressor was replaced with a Danfoss DC compressor slightly larger than the one used in the laboratory unit. The control system was integrated onto a single electronics card and packaged with its starting capacitors. The water for thermal storage was placed behind a liner that was made to fit inside the original factory liner. The original condenser was also augmented with additional surface area to improve performance. PV panels with a total rated power of 180 watts were used. The unit was tested with very successful results in an outside ambient environment, demonstrating its potential for widespread use in many off-grid applications for solar refrigeration.

Ewert, Michael K.↗

NASA NTRS2011

Predicting Spacecraft Trajectories by the WeavEncke Method

A combination of methods is proposed of predicting spacecraft trajectories that possibly include multiple maneuvers and/or perturbing accelerations, with greater speed, accuracy, and repeatability than were heretofore achievable. The combination is denoted the WeavEncke method because it is based on unpublished studies by Jonathan Weaver of the orbit-prediction formulation of the noted astronomer Johann Franz Encke. Weaver evaluated a number of alternatives that arise within that formulation, arriving at an orbit-predicting algorithm optimized for complex trajectory operations. In the WeavEncke method, Encke's method of prediction of perturbed orbits is enhanced by application of modern numerical methods. Among these methods are efficient Kepler s-equation time-of-flight solutions and self-starting numerical integration with time as the independent variable. Self-starting numerical integration satisfies the requirements for accuracy, reproducibility, and efficiency (and, hence, speed). Self-starting numerical integration also supports fully analytic regulation of integration step sizes, thereby further increasing speed while maintaining accuracy.

Weaver, Jonathan K.↗

DOE OSTI2024

Toward Accelerating Discovery via Physics-Driven and Interactive Multifidelity Bayesian Optimization

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and often nondifferentiable parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, processing spaces, and molecular embedding spaces. Often these systems are expensive or time consuming to evaluate a single instance, and hence classical approaches based on exhaustive grid or random search are too data intensive. This resulted in strong interest toward active learning methods such as Bayesian optimization (BO) where the adaptive exploration occurs based on human learning (discovery) objective. However, classical BO is based on a predefined optimization target, and policies balancing exploration and exploitation are purely data driven. In practical settings, the domain expert can pose prior knowledge of the system in the form of partially known physics laws and exploration policies often vary during the experiment. Here, we propose an interactive workflow building on multifidelity BO (MFBO), starting with classical (data-driven) MFBO, then expand to a proposed structured (physics-driven) structured MFBO (sMFBO), and finally extend it to allow human-in-the-loop interactive interactive MFBO (iMFBO) workflows for adaptive and domain expert aligned exploration. These approaches are demonstrated over highly nonsmooth multifidelity simulation data generated from an Ising model, considering spin–spin interaction as parameter space, lattice sizes as fidelity spaces, and the objective as maximizing heat capacity. Detailed analysis and comparison show the impact of physics knowledge injection and real-time human decisions for improved exploration with increased alignment to ground truth. Here, the associated notebooks allow to reproduce the reported analyses and apply them to other systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Adaptive sampling for accelerating neutron diffraction-based strain mapping *

Abstract Neutron diffraction is a useful technique for mapping residual strains in dense metal objects. The technique works by placing an object in the path of a neutron beam, measuring the diffracted signals and inferring the local lattice strain values from the measurement. In order to map the strains across the entire object, the object is stepped one position at a time in the path of the neutron beam, typically in raster order, and at each position a strain value is estimated. Typical dwell times at neutron diffraction instruments result in an overall measurement that can take several hours to map an object that is several tens of centimeters in each dimension at a resolution of a few millimeters, during which the end users do not have an estimate of the global strain features and are at risk of incomplete information in case of instruments outages. In this paper, we propose an object adaptive sampling strategy to measure the significant points first. We start with a small initial uniform set of measurement points across the object to be mapped, compute the strain in those positions and use a machine learning technique to predict the next position to measure in the object. Specifically, we use a Bayesian optimization based on a Gaussian process regression method to infer the underlying strain field from a sparse set of measurements and predict the next most informative positions to measure based on estimates of the mean and variance in the strain fields estimated from the previously measured points. We demonstrate our real-time measure-infer-predict workflow on additively manufactured steel parts—demonstrating that we can get an accurate strain estimate even with 30%–40% of the typical number of measurements—leading the path to faster strain mapping with useful real-time feedback. We emphasize that the proposed method is general and can be used for fast mapping of other material properties such as phase fractions from time-consuming point-wise neutron measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Building Initial Dynamic System Models for Digital Twins of the Cryogenic Moderator System at the ORNL Spallation Neutron Source

This work describes the initial development of dynamic system models of the cryogenic moderator system (CMS) of the Spallation Neutron Source (SNS) at ORNL as a part of the ORNL LDRD funded project Building TRANSFORM to Accelerate Digital Twin Applications for Nuclear Systems, LOIS 10563. The goal of the work is to start the dynamic system modeling effort with the end goal of using them for real-time applications as digital twins. The CMS is a cryogenic liquid hydrogen flow loop that provides moderation of the neutrons that are generated by the SNS. For optimal neutron production, the CMS needs to maintain a steady and controlled density of cryogenic hydrogen in the moderator section thus requiring precise temperature and pressure control. Due to the varied time scales and system characteristics, control of the system is complex, and diagnostics are also difficult. Difficulty in accessing the flow loop during operations, limited instrumentation and unknown design details of the equipment combine to make the case for having sophisticated digital twin models of the system. Operationally the CMS also provides a strong use case for digital twins due to the constant need of optimization and for troubleshooting/diagnostics. The large amount of data collected which are freely available for using in building the model and verifying and validating the model also makes it a great candidate for a proof-of-concept for digital twins. The project extends ORNL's capacity of development and implementation of the open-source dynamic system modeling tool TRANSFORM for engineering design and digital twin/real-time applications. Specific system configuration data for the CMS have been gathered and an initial dynamic model was created in the TRANSFORM library using Dymola as the solution platform. Models of increasing complexity are created to demonstrate the need for a multi-layered approach in digital twin modeling depending on the scale and phenomena being focused on. The dynamic modeling is shown to bring the dynamic operational aspects to the design process for systems as well as serve as a digital twin to the hardware and allow for models to be tuned and compared against real time operational data. These aims should help to push forward strategic goals of application of digital twins and increase the impact of ORNL systems modeling capabilities with TRANSFORM/Modelica for various advanced energy systems.

42 ENGINEERING↗

DOE OSTI2023

Grid Resilient, Self-Powered, Fuel Flexible, High Efficiency Heating System (Final Report)

The Grid Independent High Efficiency System (GIHES) delivers a heating system that operates free from the power grid and provides ultra-low emissions and high efficiency heating for residential and commercial buildings. The project demonstrated two weeks of continuous grid independent operation of a storage water heater equipped with a powered damper. The Thermoelectric Generator (TEG) integrated GIHES technology revolutionizes heating systems to provide grid resilient hot water production, while delivering significant value to the end user by eliminating impacts of power interruption and enhancing comfort, convenience and productivity. The system maximizes the thermal-to-electric conversion efficiency by optimizing location and orientation of the TEGs and durability while lowering system costs. The TEGs generate enough power to charge a battery and run the storage water heater unit with a powered damper uninterrupted by providing the required parasitic electrical power for start-up, shutdown and during idling. Several TEGs were analytically and experimentally evaluated based on the size, surface area and the power being generated before designing and fabricating a complete multi-TEG assembly for integration with the water heater. The GIHES technology can be extended to boilers, furnaces and tankless water heaters and integrated with heat pump to further increase the overall equipment efficiency. The technology has the potential to reduce the amount of dispatch power and peak power plant operation. An advanced burner was designed, developed and 3D printed for fuel flexible operation. The burner was designed to enable tighter integration with TEGs for potentially higher TEG output, should this be necessary. A bench-scale test rig to perform testing of a 200,000 Btu/h burner and capable of handling up to 10% H2 (by volume) was setup at GTI’s laboratory with appropriate safety and controls to ensure smooth and safe operation. Testing of the 3D burner showed < 5 ppm NOx emissions can be achieved for the entire firing range and with different levels of hydrogen blends with natural gas. The burner will enable a more integrated TEG-burner module design for improved system performance, should this be required in the future. Techno-economic analysis (TEA) of the TEG integrated with the water heaters was performed. The testing data and the power output were used to develop the TEA. The analysis for the storage unit with damper showed that the TEG integrated system costs $95. This is much lower than a one-time charge needed to wire and install dedicated power. Discussions with TEG manufacturers and OEM’s will provide more detailed information and methods to further lower these costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1975

A random search algorithm for laboratory computers

The small laboratory computer is ideal for experimental control and data acquisition. Postexperimental data processing is often performed on large computers because of the availability of sophisticated programs, but costs and data compatibility are negative factors. Parameter optimization can be accomplished on the small computer, offering ease of programming, data compatibility, and low cost. A previously proposed random-search algorithm ('random creep') was found to be very slow in convergence. A method is proposed (the 'random leap' algorithm) which starts in a global search mode and automatically adjusts step size to speed convergence. A FORTRAN executive program for the random-leap algorithm is presented which calls a user-supplied function subroutine. An example of a function subroutine is given which calculates maximum-likelihood estimates of receiver operating-characteristic parameters from binary response data. Other applications in parameter estimation, generalized least squares, and matrix inversion are discussed.

Curry, R. E.↗

NASA NTRS1989

The potential application of the blackboard model of problem solving to multidisciplinary design

The potential application of the blackboard model of problem solving to multidisciplinary design is discussed. Multidisciplinary design problems are complex, poorly structured, and lack a predetermined decision path from the initial starting point to the final solution. The final solution is achieved using data from different engineering disciplines. Ideally, for the final solution to be the optimum solution, there must be a significant amount of communication among the different disciplines plus intradisciplinary and interdisciplinary optimization. In reality, this is not what happens in today's sequential approach to multidisciplinary design. Therefore it is highly unlikely that the final solution is the true optimum solution from an interdisciplinary optimization standpoint. A multilevel decomposition approach is suggested as a technique to overcome the problems associated with the sequential approach, but no tool currently exists with which to fully implement this technique. A system based on the blackboard model of problem solving appears to be an ideal tool for implementing this technique because it offers an incremental problem solving approach that requires no a priori determined reasoning path. Thus it has the potential of finding a more optimum solution for the multidisciplinary design problems found in today's aerospace industries.

Rogers, James L.↗

NASA NTRS2009

Precision Cleaning of Oxygen Systems and Components

Currently, NASA uses Dichloropentafluoropropane (HCFC-225), a Class II ozone depleting substance (ODs), to clean contaminated oxygen systems. Starting in 20 15, the Montreal Protocols and Clean Air Act prohibit the production and importation of all hydrochlorofluorocarbons (HCFC), except for limited use in refrigeration applications. Thus, a new non-ozone depleting solvent needs to be developed for use in cleaning. Optimally, such a solvent should also be environmentally benign or green to avoid needing to replace the new solvent with yet another solvent in the future due to other environmental concerns. Work for the first year consisted of two parts. The first part was developing a method of testing the cleaning efficiency of potential solvents. Stainless steel coupons were contaminated with a known weight of various contaminants and contaminant combinations and then immersed in solvent for ten minutes. The coupons were then removed and dried in an oven until all solvent had evaporated. Once dry, the coupons were weighed and the mass of the non-volatile residue (NVR) left on the coupon was determined. The cleaning efficiency of the solvents is reported as percent cleaning, with 100% cleaning being zero NVR left on the test coupon. The second half of the first goal was to use the develop method to perform baseline testing on current solvents. The second part of the work was to begin exploring alternative cleaning solvents. A variety of hydrofluorocarbons (HFCs) were tested. Preliminary testing was also performed with ionic liquids and aqueous surfactant solutions. Once potential solvents were identified, an analysis of the performance and environmental characteristics of each was to be conducted. Four contaminants were specified for use in testing. These are Mil-Spec-H-5606 (5606), a hydraulic fluid, Mil-H-83282B (83282), another hydraulic fluid, diethylhexyl sebacate (Sebacate), and WD-40. The structures of these contaminants are all similar, with long aliphatic hydrocarbon chains of lengths between fifteen and fifty. Contaminants were tested both individually and in combination. All combinations were done with equal weights of the constituent contaminants.

McLaughlin, Russell↗

NASA NTRSDate not supplied

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

DOE OSTI2022

Dynamic Learning of Correlation Potentials for a Time-Dependent Kohn-Sham System

We develop methods to learn the correlation potential for a time-dependent Kohn-Sham (TDKS) system in one spatial dimension. We start from a low-dimensional two-electron system for which we can numerically solve the time-dependent Schr¨odinger equation; this yields electron densities suitable for training models of the correlation potential. We frame the learning problem as one of optimizing a least-squares objective subject to the constraint that the dynamics obey the TDKS equation. Applying adjoints, we develop efficient methods to compute gradients and thereby learn models of the correlation potential. Our results show that it is possible to learn values of the correlation potential such that the resulting electron densities match ground truth densities. We also show how to learn correlation potential functionals with memory, demonstrating one such model that yields reasonable results for trajectories outside the training set.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

An Optimization Approach to Support Science Decision Making for Lunar Surface Exploration

Introduction: Scientific exploration is one of the three pillars of NASA’s Moon2Mars architecture, with crew surface extra vehicular activities (EVA) serving a critical enabling function. Development of surface EVA operational planning and execution, specifically integrating science and flight control teams (FCT), is currently being explored through analog scenarios. This integration, exercised, for example, through the Joint EVA and Hu-man Surface Mobility Test Team (JETT), allows for science input on EVA activities in near real-time through a Science Evaluation Room (SER), or Arte-mis science backroom, which integrates with the broader FCT through the Science Officer. The SER works within the FCT to support dynamic EVA planning in response to changes in operational constraints as well as science opportunities and re-prioritization, increasing the mission science return and accelerating the accomplishment of the Moon2Mars science objectives. The SER works within the FCT to provide recommendations to traverse execution in near real-time. One challenge is the requirement to deliver SER inputs to the FCT on operationally relevant timelines. Failure to do so may result in suboptimal execution of science exploration EVAs or even loss of key science objectives. To close this gap, we present a network optimization tool to allow the SER to provide rapid input to the FCT in response to changes in operational constraints or science opportunities. Inputs are predicated on approved science objectives, and clear rationale must be provided to the FCT for any requested change. Accordingly, this tool incorporates the Science Traceability Matrix (STM), SER prioritization scheme, and station characterization and action planning with operational constraints such as duration, traverse speed, and distance to maximize science objectives based on SER priorities, consistent with FCT operational requirements. Method: As a proof of concept, we used an existing linear programing software package used to simulate optimal routes through cellular metabolism. We built a Demonstrative Model with three STM objectives and four stations on a region of the Moon. The objectives were given an arbitrary prioritization and mapped to the stations through four possible crew actions. (Figs. 1 and 2). This station to STM mapping is consistent with the method used by the JETT5 Science Team to develop analog surface EVA science planning. We used a grid system with the landing site at the origin and the four stations placed across the positive x,y quadrant. Actions were assigned to each station and the accomplishment of those actions resulted in a numerical “reward” based on the ability of that action to achieve science objectives. The aggregate reward from each individual STM objective contributes to a global score (Science Yield), weighted by its priority. Operational constraints included a requirement to start and end at the landing site, 5 minutes each for initial station characterization and “clean up,” and variable total EVA time, traverse rate (fixed to 0.5 meters per second in our example), and time to perform each action (10, 5, 7, and 15 min for actions 1, 2, 3, and 4, respectively). Additional constraints and variables will be added in the future (e.g., sample mass, number of stations, traverse route constraints, illumination). Optimization. We converted the connections (arcs) between these stations (nodes) into a mixed integer linear programming optimization problem (arcs = constraints, nodes = variables) with the objective to maximize Science Yield. For any action, the Science Yield is equal to the relevance of that action to an STM objective [3, 2, and 1 point(s) for High, Med., and Low relevance, respectively], multiplied by the STM Objective Priority [3, 2, and 1 point(s) for High, Med., and Low priority, respectively]. This resulted in a model that computes the optimal station and action combination to maximize the Science Yield. These weightings can be adjusted by the SER as desired. Results: We explored three test cases for the Demonstrative Model. First, we set the maximum EVA duration to 120 minutes and computed the optimal route (Fig. 3A). The model suggested per-forming Actions 1 and 2 at Station P01, followed by Actions 1 and 2 at Station P02, and finally Actions 1 and 3 at Station P04 before returning to the Landing Site. Second, we adjusted the STM Objective Priori-ty order and computed the new optimal route (Fig. 3B). Under this situation, the model suggested per-forming all Actions at Station P02 followed by all Actions at Station P03. The previous test cases were relevant to SER planning activities. Next, we explored providing mid-EVA replanning input to the FCT. Scenario: While executing the Route in Fig. 3A the crew finishes at Station P01 and FCT decides that the EVA needs to finish in 45 minutes back at the Landing Site. FCT asks SER to recommend changes to the plan to accommodate this operation-al change. Using the model and incorporating these new constraints (start at Station P01, max. time of 45 min), the model suggested performing Actions 2 and 4 at Station P03 (Fig. 4), requiring 41 minutes to complete and return to the Landing Site. Interestingly, Station 3 was not part of the original route. Using the model, we determined the EVA would need 66 minutes, instead of 45, in order for the original Station P04 to yield a larger Science Yield than Station P03. The parametrization and simulation was per-formed in less than a minute, demonstrating the operational relevance of the approach. Future Efforts: The results from the Demonstrative Model suggest this tool can accelerate SER decision making on operationally relevant timelines. Use in analog activities, such as JETT5 or follow-ons, which have over a dozen stations for a crew to explore and over a dozen actions per station, will provide needed validation of the utility of this tool for planning EVAs, replanning mid-EVA, or planning follow-on EVAs based on previous results. Further integration with FCT execution monitoring tools may provide additional efficiency gains, al-lowing rapid and iterative exploration of operation-al and science decision space by the FCT and SER.

Science Operations↗

NASA NTRSDate not supplied

Design and Optimization of a 3-Stage Axial Supercritical CO2 Compressor

This paper presents the detailed design and optimization of a three-stage axial supercritical carbon dioxide (sCO2) compressor to support the advancement of sCO 2 power cycles, recognized for their efficiency and compactness in energy conversion. The three-stage design is a scaled model of the original nine-stage 100 MW design. The aerodynamic design and structural analysis were performed at the University of Cincinnati, and the results are presented. The first stage of the design has been meticulously manufactured and experimentally tested at University of Notre Dame Turbomachinery Laboratory. The design process starts with the preliminary design process, considering the required operating boundary conditions and geometrical constraints for full stage design. The design is then scaled down to 3-stage based on the power limitation. The definition for design parameters has been inspired from the EEE HPC design and tweaked to fit the sCO 2 operation. Axisymmetric analysis is performed next and the parametric geometry modeler -Tblade3 has been used to create 3D blade geometry from the results which can then be exported for the 3D CFD analysis and optimization. An efficiency of 89.85% is predicted for the final design at the design point.

Compressor↗

DOE OSTI2021

Climate Drives Modeled Forest Carbon Cycling Resistance and Resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond (resistance) and recover from (resilience) to rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the ED-2.2 model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900, and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes, and interacted with one another. We then calculated metrics of functional ecosystem stability. Modeled gross primary production exhibited a mean resistance of -0.13, -0.23, and -0.42 (unitless) in response to increasing disturbance severity of 45%, 65%, and 85% mortality, respectively; and resilience of 0.002, 0.004, 0.008 (yr-1) respectively. This pattern held true for other carbon fluxes, meaning that in general simulated GPP declined more initially but recovered more quickly from higher mortality. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Data compression and covariance matrix inspection: Cosmic shear

Covariance matrices are among the most difficult pieces of end-to-end cosmological analyses. In principle, for two-point functions, each component involves a four-point function, and the resulting covariance often has hundreds of thousands of elements. We investigate various compression mechanisms capable of vastly reducing the size of the covariance matrix in the context of cosmic shear statistics. This helps identify which of its parts are most crucial to parameter estimation. We start with simple compression methods, by isolating and “removing” 200 modes associated with the lowest eigenvalues, then those with the lowest signal-to-noise ratio, before moving on to more sophisticated schemes like compression at the tomographic level and, finally, with the massively optimized parameter estimation and data compression (MOPED). We find that, while most of these approaches prove useful for a few parameters of interest, like Ω m , the simplest yield a loss of constraining power on the intrinsic alignment (IA) parameters as well as S 8 . For the case considered—cosmic shear from the first year of data from the Dark Energy Survey—only MOPED was able to replicate the original constraints in the 16-parameter space. Finally, we apply a tolerance test to the elements of the compressed covariance matrix obtained with MOPED and confirm that the IA parameter A IA is the most susceptible to inaccuracies in the covariance matrix.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Neutrino Theory in the Precision Era

This document summarises discussions on future directions in theoretical neutrino physics, which are the outcome of a neutrino theory workshop held at CERN in February 2025. The starting point is the realisation that neutrino physics offers unique opportunities to address some of the most fundamental questions in physics. This motivates a vigorous experimental programme which the theory community fully supports. \textbf{A strong effort in theoretical neutrino physics is paramount to optimally take advantage of upcoming neutrino experiments and to explore the synergies with other areas of particle, astroparticle, and nuclear physics, as well as cosmology.} Progress on the theory side has the potential to significantly boost the physics reach of experiments, as well as go well beyond their original scope. Strong collaboration between theory and experiment is essential in the precision era. To foster such collaboration, \textbf{we propose to establish a CERN Neutrino Physics Centre.} Taking inspiration from the highly successful LHC Physics Center at Fermilab, the CERN Neutrino Physics Centre would be the European hub of the neutrino community, covering experimental and theoretical activities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1978

The GISS sounding temperature impact test

The impact of DST 5 and DST 6 satellite sounding data on mid-range forecasting was studied. The GISS temperature sounding technique, the GISS time-continuous four-dimensional assimilation procedure based on optimal statistical analysis, the GISS forecast model, and the verification techniques developed, including impact on local precipitation forecasts are described. It is found that the impact of sounding data was substantial and beneficial for the winter test period, Jan. 29 - Feb. 21. 1976. Forecasts started from initial state obtained with the aid of satellite data showed a mean improvement of about 4 points in the 48 and 772 hours Sub 1 scores as verified over North America and Europe. This corresponds to an 8 to 12 hour forecast improvement in the forecast range at 48 hours. An automated local precipitation forecast model applied to 128 cities in the United States showed on an average 15% improvement when satellite data was used for numerical forecasts. The improvement was 75% in the midwest.

Halem, M.↗