Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗
NASA NTRS1978

Lunar occultation of Saturn. II - The normal reflectance of Rhea, Titan, and Iapetus

An inversion procedure to obtain the reflectance of the central region of a satellite's disk from lunar occultation data is presented. The scheme assumes that the limb darkening of the satellite depends only on the radial distance from the center of the disk. Given this assumption, normal reflectances can be derived that are essentially independent of the limb darkening and the diameter of the satellite. The procedure has been applied to our observations of the March 1974 lunar occultation of Tethys, Dione, Rhea, Titan, and Iapetus. In the V passband we derive the following normal reflectances: Rhea (0.97 plus or minus 0.20), Titan (0.24 plus or minus 0.03), Iapetus, bright face (0.60 plus or minus 0.14). For Tethys and Dione the values derived have large uncertainties, but are consistent with our result for Rhea.

Elliot, J. L.↗
NASA NTRS1987

A natural-language interface to a mobile robot

The present work on robot instructability is based on an ongoing effort to apply modern manipulation technology to serve the needs of the handicapped. The Stanford/VA Robotic Aid is a mobile manipulation system that is being developed to assist severely disabled persons (quadriplegics) in performing simple activities of everyday living in a homelike, unstructured environment. It consists of two major components: a nine degree-of-freedom manipulator and a stationary control console. In the work presented here, only the motions of the Robotic Aid's omnidirectional motion base have been considered, i.e., the six degrees of freedom of the arm and gripper have been ignored. The goal has been to develop some basic software tools for commanding the robot's motions in an enclosed room containing a few objects such as tables, chairs, and rugs. In the present work, the environmental model takes the form of a two-dimensional map with objects represented by polygons. Admittedly, such a highly simplified scheme bears little resemblance to the elaborate cognitive models of reality that are used in normal human discourse. In particular, the polygonal model is given a priori and does not contain any perceptual elements: there is no polygon sensor on board the mobile robot.

Michalowski, S.↗
NASA NTRS2004

Relative Sensor with 4(pi) Coverage for Formation Flying Missions

The Terrestrial Planet Finder (TPF) pre-project, an element of NASA's Origins program, is currently developing two architectures for a mission to search for earth-like planets around nearby stars. One of the architectures being developed is the Formation Flying Interferometer (FFI). The FFI is envisioned to consist of up to seven spacecraft (as many as six 'collectors' with IR telescopes, and a 'combiner') flying in precise formation within f 1 cm of pre-determined trajectories for synchronized observations. The spacecraft-to-spacecraft separations are variable between 20 m and 100 m or more during observations to support various configurations of the interferometer in the planet-finding mode. The challenges involved with TPF autonomous operations, ranging from formation acquisition and formation maneuvering to high precision formation control during science observations, are unprecedented. In this paper we discuss the development of the formation acquisition sensor, which uses novel modulation and duplexing schemes to enable fast signal acquisition, multiple-spacecraft operation, and mitigation of inherent jamming conditions, while providing precise formation sensing and integrated radar capability. This approach performs delay synthesis and carrier cycle ambiguity resolution to improve range measurement, and uses differential carrier cycle ambiguity resolution to make precise bearing angle measurements without calibration maneuvers.

formation flying sensors↗
NASA NTRS1989

Launch packaging options for the photovoltaic power module cargo element

The National Aeronautics and Space Administration recently embarked on the Space Station Freedom program, which will utilize the Shuttle Orbiter for transportation to orbit. This task will be accomplished with a number of flights over several years. Each flight is unique in terms of the hardware that is manifested and the method by which it is integrated to form viable cargo elements. Work Package 4 is responsible for the electric power system for Space Station Freedom, and was delegated the authority to develop a photovoltaic (PV) power module cargo element. The PV power module consists of several unique assemblies. The first of these is the combined solar array/beta gimbal assembly. The remaining assemblies form the single combined integrated equipment assembly for each PV power module. These three combined assemblies are packaged into a launch cradle to form the PV power module cargo element, which is placed in the cargo bay of the Shuttle Orbiter for transportation to orbit. Various constraints determine the packaging options for the three PV power module combined assemblies. The size and shape of the combined assemblies in relation to the Shuttle Orbiter cargo bay dimensions and other manifested hardware are ultimately a factor in determining the acceptable packaging schemes for the PV power module cargo element. Several packaging options for the PV power module cargo element are presented. These options are discussed in terms of their impact on the overall flight hardware manifest as determined by the various constraints.

Hoberecht, Mark A.↗
DOE OSTI2021

Towards an Agent-Based Blackboard System for Reactor Design Optimization

The field of reactor design is rich with opportunities for applications of computational optimization algorithms; these applications can range from preliminary core design to reactor shuffling patterns. Many of these schemes rely on sets of previously generated solutions (sometimes referred to as “generations”) to inform future decisions. While it is important to build upon prior knowledge, this process requires a full generation of solutions to be formed before future solutions can be examined. Rather than relying on a generational scheme to perform an optimization, we propose using an agent-based approach in conjunction with a blackboard framework for performing reactor design optimizations. Utilizing an agent-based approach allows agents to perform tasks independently, while retaining the ability to build off of previous solutions. We develop an agent-based blackboard system (ABBS) for determining the Pareto front (PF) in sodium fast reactor design optimization problems and compared this with the Non-Dominated Sorting Genetic Algorithm II (NSGA-II). Our goal is to evaluate the viability of the ABBS in producing a PF that is comparable with the NSGA-II algorithm. The design space consists of the fuel height, fuel smear, and plutonium fraction in the core, and we seek to minimize the reactivity swing and plutonium mass, while maximizing the burnup. The diversity, coverage, and spread of the PFs generated by the two methods are examined, and the ABBS is able to converge to the same PF as the NSGA-II algorithm. These results show that the ABBS is able to find optimal designs that are similar to those found by the NSGA-II algorithm. We conclude our study by applying the ABBS to the design of a sodium-cooled fast reactor to dispose of weapons-grade plutonium. The ABBS finds a core design that can burn upwards of 17.5 kg of weapons-grade plutonium per year and degrade an additional 195 kg of weapons-grade plutonium per year into non-weapons-grade material.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1983

Numerical solution of the Navier-Stokes equations for compressible turbulent two/three dimensional flows in the terminal shock region of an inlet/diffuser

The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique have been used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations have been carried out: (1) 2-D, steady, subsonic; (2) 2-D, steady, transonic with normal shock, (3) 2-D, steady, supersonic with terminal shock, (4) 2-D, transient process of shock development and (5) 3-D, steady, transonic with normal shock. The numerical results obtained for the 2-D and 3-D transonic shocked flows have been compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data. The predicted transient responses of the flow to externally applied impulsive disturbance are consistent with results obtained via asymptotic analysis.

Liu, N.-S.↗
NASA NTRS2002

Investigation of Aerodynamic and Radiometric Land Surface Temperatures

The surface temperature, T(sub s), of a land surface measured by a radiometer, T(sub s,r), and the temperature "felt" by the air, T(sub aero), often differ significantly and are difficult if not impossible to define rigorously. However, recent studies conducted by the principal investigators with several land surface models suggest that this problem can be largely resolved. The main goal of this project was to use model-based and empirical studies to improve understanding and reconcile the difference between T(sub s,r) and T(sub aero), while maintaining consistency within the models and with theory and data. The results from this effort have contributed progress towards the effective use of remotely sensed surface temperature measurements taken from an arbitrary view angle over a partial canopy cover for producing high quality sensible and latent heat flux estimates. In addition, we have developed parameterizations that are designed to improve the representation of the roughness length for heat in climate and mesoscale models. This not only provides improved representation of surface energy balance in such models, but should also facilitate the use of surface temperature measurements for validating or updating the surface temperature produced by SVATs (soil-vegetation-atmosphere schemes) in climate or mesoscale models.

Friedl, Mark↗
DOE OSTI2025

Angular Correlations of Gamma Rays from Fission Product Solution

Fission products from photon-induced fission are routinely used at Idaho National Laboratory as part of nuclear forensics activities. Typically, the focus is on fission products with relatively large fission yields (Ba-140, Mo-99, etc.). The nuclear data associated with the well-known isotopes have small uncertainties, and many have well-known angular correlations between emitted gamma-rays, if applicable. However, for isotopes with smaller fission yields or less intense gamma-ray emission intensities, the nuclear decay schemes and gamma-ray emission intensities may be less well-known. The Germanium Rotational Measurements for Angular Correlation (GeRMAC) system was designed to use two high-purity germanium (HPGe) detectors in a coincidence counting set-up but allows one of the HPGe detectors to rotate at a constant distance from a radioactive source, changing the relative angle between the two HPGe detectors. Angular correlations can then be determined from the relative photopeak areas between two gamma-rays emitted as part of the same nuclear decay. By recording the HPGe spectra in list mode, the angular correlations of gamma-rays emitted from complex sources, consisting of multiple isotopes, can be measured simultaneously and the individual isotope emissions can be isolated. Using the GeRMAC system, a multi-isotope fission product source, produced with photon-induced fission of U-238, was measured using a single rotational arc (180 degrees) to obtain the angular correlation between gamma-rays emitted by the large number of fission products contained in the solution with an angular accuracy of 1 degree. In particular, La-140 (38 gamma rays) and Nd-147 (26 gamma rays) are of interest to the nuclear forensics community.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
NASA NTRS2024

Launch Complex 34, SWMU CC054 2023 DNAPL Source Zone Operations, Maintenance, and Monitoring, Site-Wide Long-Term Monitoring, and Hot Spot 6 Air Sparge System Annual Performance Monitoring and Phase Two Expansion Construction Completion Report Cape Canaveral Space Force Station, Florida

This Annual Performance Monitoring Report (PMR) for the Dense Non-Aqueous Phase Liquid (DNAPL) Source Zone (DSZ), Site-Wide Long-Term Monitoring (LTM), and Hot Spot 6 (HS 6) Air Sparge (AS) System presents the results of Year 14 operations and performance monitoring of the hydraulic containment (HC) Interim Measure (IM), details associated with construction and implementation of the HS 6 AS system expansion (Phase Two), and the results of operations and performance sampling of the HS 6 AS IM at Launch Complex 34 (LC34), located at Cape Canaveral Space Force Station (CCSFS), Florida. The timeframe for activities documented in this PMR extends from April 1, 2023 to March 31, 2024. LC34 has been designated Solid Waste Management Unit CC054 under the Kennedy Space Center (KSC) Resource Conservation and Recovery Act Corrective Action Program. The objective of the HC IM at LC34 is to contain the shallow and deep DSZ and surrounding dissolved-phase trichloroethene (TCE) high concentration plume via operation of a hydraulic containment system (HCS). The pre-IM design 300 micrograms per liter (μg/L) TCE groundwater contour was used to establish the deep zone capture area for deep recovery wells, and the shallow zone capture area was defined by the DSZ. The system began operating in 2010, and in 2015, the system was expanded to provide HC for areas within the 300 μg/L TCE groundwater isocontours of HS 3 and 4. In 2018 and 2019, an investigation was conducted to recharacterize the DSZ, which included investigating TCE mass in Layer 7. This data was subsequently used to optimize the pumping rates of the HCS and install additional recovery wells in Layer 7 to more adequately capture residual contaminant mass. The operational period for Year 14 of the HCS was from April 1, 2023 to March 31, 2024. Operational runtime for the system was 94 percent during Year 14, with downtime events attributed to planned maintenance, system repairs, and power outages. As of March 31, 2024, a total of 344,849,634 cumulative gallons of groundwater containing 94,656 pounds of chlorinated volatile organic compounds (CVOCs) have been removed by the HCS. During the reporting period covered under this report, the HCS recovered 31,176,393 gallons and approximately 6,319 pounds of CVOC mass. Total combined influent concentrations of TCE have decreased since startup from approximately 280,000 µg/L (January 2010) to 25,000 µg/L (March 2024). During the reporting period, all effluent concentrations from the HCS (aqueous and vapor) were below regulatory reporting limits, indicating the system continues to operate as intended. Performance monitoring was conducted in January 2024 within the DSZ to evaluate TCE contamination. Groundwater samples were collected via DPT at nine locations, consistent with previous events between 2017 and 2022. Full vertical profile sampling was completed at each DPT from 8 to 98 feet below land surface (bls), at 5-foot intervals. The DPT performance monitoring results are summarized in this PMR. The results revealed TCE remains at concentrations greater than 11,000 µg/L in the DSZ (1-percent solubility, indicative of DNAPL) at eight of the nine DPT locations and at depths ranging from 28 to 98 feet bls. An overall decreasing trend of TCE concentrations was observed in DPT samples during this reporting period, which is a reduction from the previous event (December 2022) and the peak event in December 2021, where TCE percentages appeared to increase in all depth zones because several recovery wells were turned off during the AS pilot study in the DSZ. The maximum TCE concentration in January 2024 was 1,600,000 µg/L in the 53 feet bls depth interval at DPT594 (previous maximum result in 2022 was 1,800,000 µg/L in the 48 feet bls depth interval at DPT599). This maximum concentration in the 53 feet bls depth interval is in the deep capture zone. During the January 2024 DPT event, the largest portion of TCE mass was observed in the 48 feet bls interval above/within Layer 4. This trend remains consistent with previous years and appears to indicate continued mass discharge from Layer 4 (fine-grained unit). In addition to DPT sampling, annual groundwater samples were collected from 11 deep monitoring wells in the DSZ area (Layers 7 and 8) in December 2023 to verify vertical and horizontal delineation. Three of the wells were also sampled biweekly to evaluate operations of recovery well RW21D (screened 86 to 106 feet bls), which was installed in January 2023. Of the Layer 7/8 monitoring sampled only annually, results were non-detect or less than groundwater cleanup target levels GCTLs in December 2023, with the exception of one well, IW45D2, which had a cis-1,2-dichloroethene (cDCE), detection greater than the GCTL. Of the three wells sampled biweekly during the operational period, the well located closest to Layer 7 recovery well RW21D (IW44D2, screened 105 to 115 feet bls) had concentrations of TCE, cDCE and vinyl chloride (VC) greater than GCTLs throughout the operational period, but displayed a decreasing trend since the peak concentrations in September 2023. The maximum TCE concentration during this operational period was 190,000 µg/L at IW44D2 in September 2023, but reduced to 700 µg/L in March 2024, indicating the HCS is still effectively removing mass from the source area. Expansion of the HCS and addition of new recovery wells is ongoing and will continue to be evaluated as the groundwater recovery scheme is optimized. Details of the expansion and optimization will be provided in a future PMR. The HS 6 AS IM was initiated in 2018 with 160 AS wells and expanded in 2019 with another 140 AS wells. An additional expansion of the HS 6 AS IM was completed during the reporting period covered under this report and details of the construction implementation and startup of the expansion are detailed in Section III of this report. The new expansion, referred to as Phase Two, was implemented between August 17, 2022 and August 28, 2023, and included the installation of 190 air sparge wells to treat an additional 11.2 acres. The original configuration (referred to as Phase One) operated until Phase Two came online, then all but 52 AS wells were turned off so the components could be moved and utilized in the Phase Two area. The 52 AS wells that remain on are in a barrier configuration preventing contaminated groundwater from impacting the treated area. The HS 6 AS system (both Phase One and Two) operated normally during the reporting period covered under this report. Semi-annual performance monitoring of the Phase One configuration was conducted in April and November 2023, consistent with previous years. For the Phase Two configuration, 21 new monitoring wells were installed and sampled quarterly, with a baseline event in July 2023, and quarterly events in November 2023 and February 2024 summarized in this report. Semi-annual monitoring results collected in April and October 2023 show concentrations of contaminants of concern (cDCE, trans-1,2-dichloroethene, and VC) have decreased to less than GCTLs in nearly all wells and not impacting the surface water drainage canal, indicating the HS 6 IM continues to meet objectives. The baseline and quarterly sampling for the Phase Two configuration indicate generally decreasing concentrations in wells within and around the perimeter of the treatment area. At least two more quarters of monitoring will be conducted and once those results are evaluated a reduced the sampling frequency may be considered. Overall, the tasks associated with Year 14 operation of the HC IM and operation of the HS 6 AS IM were performed in accordance with recommendations included in the previous 2022 LC34 (Year 13) PMR. Evaluation of results from the HC IM and HS 6 IM show that these systems are operating as designed and meeting performance objectives.

groundwater remediation↗
NASA NTRS1990

Numerical simulations of interacting disk galaxies

Galaxy-galaxy interactions have long attracted many extragalactic astronomers in various aspects. A number of computer simulations performed in the 1970s have successfully reproduced the peculiar morphologies observed in interacting disk galaxies and clarified that tidal deformation explains most of the observed global peculiarities. However, most of these simulations have used test particles in modelling the disk component. Tidal response of a self-gravitating disk remains to be further clarified. Another topic which is intensely discussed at present is the relation between galaxy-galaxy interactions and activity. Many observations suggest that interactions trigger strong starbursts and possibly active galactic nuclei (AGN). However, the detailed mechanism of triggering is not yet clear. It is vital here to understand the dynamics of interstellar gas. In order to understand various phenomena related to galaxy-galaxy interactions (mainly for disk galaxies), the author performed a series of numerical simulations on close galaxy encounters which includes both interstellar gas and self-gravitating disk components. In these simulations, the galaxy model to be perturbed (target galaxy) consists of a halo and a disk. The halo was treated as a rigid spherical gravitational field which is assumed to remain fixed during the interaction. The disk is composed of stars and gas. The stellar disk was constructed by 20000 collisionless particles of the same mass. Those particles move in the halo gravitational field, interacting with each other and with the perturber. Therefore, the self-graviy of the disk is properly taken into account. Stellar particles were initially given circular velocities with small random motions required to stabilize the disk against local axisymmetric disturbances. The gravitational field of the stellar disk was calculated by the particle-mesh scheme (e.g. Hockney and Eastwood 1981). The gaseous component was modelled by the cloud-particle scheme (e.g. Roberts and Hausman 1984). Here, the authors represent the gas as an ensemble of small spheres (i.e. clouds) and include the creation of an OB star in a cloud-cloud collision and subsequent velocity push on nearby clouds due to a supernova explosion.

Noguchi, Masafumi↗
DOE OSTI2020

A non-perturbative pairwise-additive analysis of charge transfer contributions to intermolecular interaction energies

Energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs) decomposes the interaction energy between molecules into physically interpretable components like geometry distortion, frozen interactions, polarization, and charge transfer (CT, also sometimes called charge delocalization) interactions. In this work, a numerically exact scheme to decompose the CT interaction energy into pairwise additive terms is introduced for the ALMO-EDA using density functional theory. Unlike perturbative pairwise charge-decomposition analysis, the new approach does not break down for strongly interacting systems, or show significant exchange–correlation functional dependence in the decomposed energy components. Both the energy lowering and the charge flow associated with CT can be decomposed. Complementary occupied–virtual orbital pairs (COVPs) that capture the dominant donor and acceptor CT orbitals are obtained for the new decomposition. It is applied to systems with different types of interactions including DNA base-pairs, borane-ammonia adducts, and transition metal hexacarbonyls. While consistent with most existing understanding of the nature of CT in these systems, we report the results also reveal some new insights into the origin of trends in donor–acceptor interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2015

Advanced Stirling Convertor Dual Convertor Controller Testing at NASA Glenn Research Center in the Radioisotope Power Systems System Integration Laboratory

NASA Glenn Research Center developed a nonnuclear representation of a Radioisotope Power System (RPS) consisting of a pair of Advanced Stirling Convertors (ASCs), Dual Convertor Controller (DCC) EMs (engineering models) 2 and 3, and associated support equipment, which were tested in the Radioisotope Power Systems System Integration Laboratory (RSIL). The DCC was designed by the Johns Hopkins University Applied Physics Laboratory (JHU/APL) to actively control a pair of ASCs. The first phase of testing included a Dual Advanced Stirling Convertor Simulator (DASCS), which was developed by JHU/APL and simulates the operation and electrical behavior of a pair of ASCs in real time via a combination of hardware and software. RSIL provides insight into the electrical interactions between a representative radioisotope power generator, its associated control schemes, and realistic electric system loads. The first phase of integration testing included the following spacecraft bus configurations: capacitive, battery, and super-capacitor. A load profile, created based on data from several missions, tested the RPS's and RSIL's ability to maintain operation during load demands above and below the power provided by the RPS. The integration testing also confirmed the DCC's ability to disconnect from the spacecraft when the bus voltage dipped below 22 volts or exceeded 36 volts. Once operation was verified with the DASCS, the tests were repeated with actual operating ASCs. The goal of this integration testing was to verify operation of the DCC when connected to a spacecraft and to verify the functionality of the newly designed RSIL. The results of these tests are presented in this paper.

Controllers↗
DOE OSTI2025

Blinding scheme for the scale-dependence bias signature of local primordial non-Gaussianity for DESI 2024

The next generation of spectroscopic surveys is expected to achieve an unprecedented level of accuracy in the measurement of cosmological parameters. To avoid confirmation bias and thereby improve the reliability of these results, blinding procedures become a standard practice in the cosmological analyses of such surveys. Blinding is especially crucial when the impact of observational systematics is important relative to the cosmological signal, and a detection of that signal would have significant implications. This is the case for local primordial non-gaussianity, as probed by the scale-dependent bias of the galaxy power spectrum at large scales that are heavily sensitive to the dependence of the target selection on the imaging quality, known as imaging systematics. We propose a blinding method for the scale-dependent bias signature of local primordial non-gaussianity at the density field level which consists in generating a set of weights for the data that replicate the scale-dependent bias. The applied blinding is predictable, and can be straightforwardly combined with other catalog-level blinding procedures that have been designed for the baryon acoustic oscillation and redshift space distortion signals. The procedure is validated through simulations that replicate data from the first year of observation of the Dark Energy Spectroscopic Instrument, but may find applications to other upcoming spectroscopic surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture

Meshfree discretizations of state-based peridynamic models are attractive due to their ability to naturally describe fracture of general materials. However, two factors conspire to prevent meshfree discretizations of state-based peridynamics from converging to corresponding local solutions as resolution is increased: quadrature error prevents an accurate prediction of bulk mechanics, and the lack of an explicit boundary representation presents challenges when applying traction loads. Herein, we develop a reformulation of the linear peridynamic solid (LPS) model to address these shortcomings, using improved meshfree quadrature, a reformulation of the nonlocal dilatation, and a consistent handling of the nonlocal traction condition to construct a model with rigorous accuracy guarantees. In particular, these improvements are designed to enforce discrete consistency in the presence of evolving fractures, whose a priori unknown location render consistent treatment difficult. In the absence of fracture, when a corresponding classical continuum mechanics model exists, our improvements provide asymptotically compatible convergence to corresponding local solutions, eliminating surface effects and issues with traction loading which have historically plagued peridynamic discretizations. When fracture occurs, our formulation automatically provides a sharp representation of the fracture surface by breaking bonds, avoiding the loss of mass. We provide rigorous error analysis and demonstrate convergence for a number of benchmarks, including manufactured solutions, free-surface, nonhomogeneous traction loading, and composite material problems. Finally, we validate simulations of brittle fracture against a recent experiment of dynamic crack branching in soda-lime glass, providing evidence that the scheme yields accurate predictions for practical engineering problems.

42 ENGINEERING↗
DOE OSTI2021

Development of Enabling Technologies for Chemical Looping Combustion and Chemical Looping with Oxygen Uncoupling (Final Report)

This report summarizes results from the project, “Development of Enabling Technologies for Chemical Looping Combustion and Chemical Looping with Oxygen Uncoupling,” which evaluated several aspects of dual fluidized bed chemical looping combustion and chemical looping with oxygen uncoupling (CLOU). The objective was to provide tools and enabling technologies to help advance fluidized bed chemical looping technology to pilot, demonstration and commercial scale. One focus area is oxygen carriers, which are key to chemical looping combustion. The copper oxygen carrier-coal ash system was systematically evaluated through a combination of thermodynamic modeling and lab-scale experiments, taking into consideration different oxygen carrier support materials and coal types. A method of mapping “safe” and “risk” zones for different combinations was established, and recommendations for target conditions are provided. In addition, a simple solution for limiting negative influence of some coal ashes, namely adding small amounts of calcium to the system, is proposed. In addition, a novel process for recovering and recycling copper from spent oxygen carriers is proposed. A new approach for design and operation of loop seals in a CLOU system was developed, and involves distinct gas injection points and a short horizontal section to help control the fate of the fluidizing gases. For the air-to-fuel reactor loop seal, the upstream side is fluidized with air, ideally input into the side rather than into the bottom, which prevents uncoupling (reduction) of the oxygen carrier before entering the fuel reactor. Also, alternative for separating oxygen carrier particles and gas in circulating fluidized bed-based chemical looping system is proposed. A new reduced reaction scheme for conversion of coal in copper-based CLOU was developed and implemented into CPFD Software’s Barracuda VR package. The resulting model provides higher fidelity than the baseline, especially when it comes to minor reactions that are part of the overall combustion environment. For in-reactor heat extraction, it was determined in this project that the best way to do that in a dual circulating fluidized-bed system is through heat exchange low in the air reactor. The air reactor is the hotter of the two and for a fast-fluidized circulating bed, the lower, more dense section has a higher heat transfer coefficient and offers more consistent particle-wall contact since it doesn’t have the splashing behavior the top of the bed does. It was determined that just the wall surface area in the lower quarter to third of the air reactor is sufficient to control temperatures in both reactors. Finally, consideration was given to a new concept for CLC, which involves a staged fuel reactor with a different type of oxygen carrier in each stage.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

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 and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks are analyzed. 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 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).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Ultra-high-energy cosmic rays in a galactic wind and its termination shock

Results are reported from numerical modeling of the acceleration and transport of ultra-high-energy cosmic rays in a galactic wind and its termination shock. A two-dimensional (azimuthally symmetric) wind and spiral magnetic field, with a spherical termination shock, where the velocity drops suddenly, is assumed. The time-dependent cosmic-ray transport equation, including all major transport effects is solved using an implicit finite-difference scheme. Particles are injected as the shock of low energy, and the subsequent evolution of the distribution function is followed. Iron nuclei are readily accelerated at the shock to energies up to 100 billion GeV, and protons to 10 billion GeV. A major effect aiding the acceleration of these particles is the spiral of the magnetic field carried out by the wind, caused by the rotation of the Galaxy, with the result that the shock is nearly normal over most of its area. Increasing the magnetic field or rotation rate increases the maximum energy attainable. Anisotropies and energy densities of the particles are also discussed. It is concluded that the process is consistent with observations of ultra-high-energy cosmic rays.

Jokipii, J. R.↗