Search NASA⌕ Search

SEARCH · Search NASA

Results for “compressibility factor”

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 307 records · Page 17

Engineered Nanoporous Frameworks for Adsorption Cooling Applications

The energy demand for traditional vapor-compressed technology for space cooling continues to soar year after year due to global warming and the increasing human population’s need to improve living and working conditions. Thus, there is a growing demand for eco-friendly technologies that use sustainable or waste energy resources. Here this review discusses the properties of various refrigerants used for adsorption cooling applications followed by a brief discussion on the thermodynamic cycle. Next, sorbents traditionally used for cooling are reviewed to emphasize the need for advanced capture materials with superior properties to improve refrigerant sorption. The remainder of the review focus on studies using engineered nanoporous frameworks (ENFs) with various refrigerants for adsorption cooling applications. The effects of the various factors that play a role in ENF–refrigerant pair selection, including pore structure/dimension/shape, morphology, open-metal sites, pore chemistry and possible presence of defects, are reviewed. Next, in-depth insights into the sorbent–refrigerant interaction, and pore filling mechanism gained through a combination of characterization techniques and computational modeling are discussed. Finally, we outline the challenges and opportunities related to using ENFs for adsorption cooling applications and provide our views on the future of this technology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Process Development of the Vaporizing Foil Actuator Welding Technique

The industrial focus on continuous improvement for products is leading design engineers to consider multi-material design concepts more frequently. The design concepts use the advantages of each material while simultaneously minimizing the drawbacks. It also allows for the overall weight of the products to be reduced. Some traditional joining methods, though, struggle to join dissimilar material combinations. A technique known as vaporizing foil actuator welding, which was developed around ten years ago, is a solid-state joining method that has the potential to overcome the barriers to dissimilar material joining. The biggest drawbacks of vaporizing foil actuator welding are that it is still exclusively used in laboratory settings and early in the technology development process. As such, there is not enough confidence in the process for it to be transitioned to manufacturing environments yet. This work begins by analyzing the current state of the technology and identifying a roadmap to achieve a transition to manufacturing. The analysis found that increased confidence in the technology requires sample to sample repeatability to be improved. Beyond that, the process also needs to be fully automated before it can be considered as a viable technique for mass production. In order to further develop the technique, a fully automated work cell was constructed. By removing the human element from the process, this cell helped identify iii the aspects which contribute to process variability. The positioning of the vaporizing foil actuator was identified as the most important factor in process stability. A small set of samples were made in the automated cell with the addition of one component epoxy adhesive around the perimeter of the weld. These samples showed improved process repeatability compared with the samples made without adhesive. This indicates that small amounts of adhesive may be critical in using the process in industrial settings. Additionally, work was performed with manual positioning and adhesive dispensing. This work found that a hybrid joint behaves like traditional weld bonding, where the addition of adhesive makes the joint stronger than if either of the two joining methods were used on their own. To realize this performance with vaporizing foil actuator welding, though, the adhesive needs to be protected from burning during the rapid compression of air during the impact process.

36 MATERIALS SCIENCE↗

Effect of a micro-scale dislocation pileup on the atomic-scale multi-variant phase transformation and twinning

In this paper, we perform concurrent atomistic–continuum (CAC) simulations to assess the contribution of the internal stress induced by the microscale dislocation pileup at an atomically structured interface to the atomic-scale phase transformations (PTs), reverse PTs, and twinning. The main novelty of this work is to unify the atomistic description of the interface and the coarse-grained (CG) description of the lagging dislocations away from the interface within one single framework. Our major findings are: (a) the interface dynamically responds to a pileup by forming steps/ledges, the height of which is proportional to the number of dislocations arriving at the interface; (b) the pileup-induced internal stress concentration profile follows neither the classical Eshelby model nor the super-dislocation model alone, but a combination of them; (c) when the pre-sheared sample is compressed, a direct square-to-hexagonal PT occurs ahead of the pileup tip and eventually grows into a wedge shape. The two variants of the hexagonal phases form a twin with respect to each other; (d) upon a further increase of the loading, part of the newly formed hexagonal phase transforms back to the square phase. The square product phase resulting from this reverse PT forms a twin with respect to the initial square phase. All phase boundaries (PBs) and twin boundaries (TBs) are stationary and correspond to zero thermodynamic Eshelby driving forces; and (e) the microscale dislocation pileup-induced internal shear stress and the structural change at the atomic-scale interface reduces the stress required for initiating a PT by a factor of 5.5, comparing with that in the sample containing no dislocations. This work is the first characterization of the behavior of PTs/twinning resulting from the reaction between a microscale dislocation slip and an atomically structured interface. In conclusion, the gained knowledge will advance our understanding of how the multi-phase material behaves in many complex physical processes, such as the synthesis of multi-phase high-entropy alloys or superhard ceramics under high-pressure torsion, deep mantle earthquakes in geophysics, and so on, which all involve dislocation slip, PTs, twinning, and their interactions across from the atomistic to the microscale and beyond.

36 MATERIALS SCIENCE↗

Effect of laser scanning speed on the microstructure, phase transformation and mechanical property of NiTi alloys fabricated by $\mathrm{LPBF}$

This study investigated the effects of laser scanning speed on the microstructure, phase transformation and properties of NiTi alloys fabricated by laser powder bed fusion (LPBF). In this study, the contributions of metallurgical factors under different scanning speeds, such as Ni evaporation, Ni4Ti3 precipitation, dislocations and internal stress, to the transformation temperature and transformation latent heat were clarified through specially designed experiments. Ni evaporation is found to have the most profound effect, followed by precipitation. Increasing scanning speed is found to reduce Ni loss, thus cause less increase in the transformation temperature and transformation heat of the LPBF-NiTi alloys. Increasing scanning speed also increases the microstructure non-uniformity and thus widens the transformation temperature interval. The orientations of residual stress exhibit strong crystallographic stiffness dependence. The LPBF-NiTi alloys with different scanning speeds all exhibited high strains (>13.4%) and excellent shape memory effect. A LPBF-NiTi honeycomb structure exhibited 96% shape recovery rate after a 60% pre-compressive deformation. Besides, there is an optimum scanning speed for minimum porosity and smallest average pore size. However, the pore structure is found to have weak influence on the tensile behaviour of the LPBF-NiTi alloys, possibly due to the high defect tolerance of the martensitic transformation.

36 MATERIALS SCIENCE↗

Multiscale Richtmyer-Meshkov instability experiments to isolate the strain rate dependence of strength

Theoretical analysis of Richtmyer-Meshkov instability (RMI) experiments for solid strength shows that the strain rate for a given shock should be inversely proportional to the length scale of the sine wave perturbations when η 0 k , the nondimensional amplitude to wavelength ratio, is held fixed. To isolate the effect of strain rate on strength, free-surface RMI specimens of annealed copper were prepared with three perturbation regions with the same η 0 k but different length scales, characterized by the wavelength λ varying by a factor of 4.9 from 65 to 130 to 320 µ m . Three such targets with different fixed η 0 k ′ s were impacted to a shock pressure of 25 GPa, and the instability evolution was measured with photon Doppler velocimetry. Strengths estimated by comparing hydrocode simulation to the data increased from 700 to 1200 MPa as λ decreased. The different η 0 k targets exercised increasing amounts of plastic strain yet showed no evidence of strain hardening. Physical regime sensitivity analysis determined that for 320 − 65 µ m wavelength perturbations, the effective strain rates increased from 8.7 × 10 6 to 3.3 × 10 7 s − 1 , a factor of 3.8. Thus, the predicted strain rate scaling was mostly achieved but slightly suppressed by increased strength at higher rates. The RMI strength estimates were plotted against constitutive testing data on copper from the literature to show striking evidence of the strength upturn at higher strain rates. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Collective Summary of sCO2 Materials Development (Supercritical Transformational Electric Power Generation (STEP) Level 2 Milestone Report) (Parts I - II)

Polymers such as PTFE (polytetrafluorethylene or Teflon), EPDM (ethylene propylene diene monomer) rubber, FKM fluoroelastomer (Viton), Nylon 11, Nitrile butadiene (NBR) rubber, hydrogenated nitrile rubber (HNBR) and perfluoroelastomers (FF_202) are commonly employed in super critical CO 2 (sCO2) energy conversion systems. O-rings and gaskets made from these polymers face stringent performance conditions such as elevated temperatures, high pressures, pollutants, and corrosive humid environments. In FY 2019, we conducted experiments at high temperatures (100°C and 120°C) under isobaric conditions (20 MPa). Findings showed that elevated temperatures accelerated degradation of polymers in sCO2, and that certain polymer microstructures are more susceptible to degradation over others. In FY 2020, the focus was to understand the effect of sCO2 on polymers at low (10 MPa) and high pressures (40 MPa) under isothermal conditions (100°C). It was clear that the same selectivity was observed in these experiments wherein certain polymeric functionalities showed more propensity to failure over others. Fast diffusion, supported by higher pressures and long exposure times (1000 hours) at the test temperature, caused increased damage in sCO2 environments to even the most robust polymers. We also looked at polymers under compression in sCO2 at 100°C and 20 MPa pressure to imitate actual sealing performance required of these materials in sCO2 systems. Compression worsened the physical damage that resulted from chemical attack of the polymers under these test conditions. In FY 2021, the effect of cycling temperature (from 50°C to 150°C to 50°C) for polymers under a steady sCO2 pressure of 20 MPa was studied. The aim was to understand the influence of cycling temperatures of sCO2 for typical polymers under isobaric (20 MPa) conditions. Thermoplastic polymers (Nylon, and PTFE) and elastomers (EPDM, Viton, Buna N, Neoprene, FF202, and HNBR) were subjected to 20 MPa sCO2 pressure for 50 cycles and 100 cycles in separate experiments. Samples were extracted for ex-situ characterization at 50 cycles and upon the completion of 100 cycles. Each cycle constituted of 175 minutes of cycling from 50°C to 150°C. The polymer samples were examined for physical and chemical changes by Dynamic Mechanical and Thermal Analysis (DMTA), Fourier Transform Infrared (FTIR) spectroscopy, and compression set. Density and mass changes immediately after removal from test were measured for degree of swell comparisons. Optical microscopy techniques and micro computer tomography (micro CT) images were collected on select specimens. Evaluations conducted showed that exposures to super-critical CO2 environments resulted in combinations of physical and/or chemical changes. For each polymer, the dominance of cycling temperatures under sCO2 pressures, were evaluated. Attempts were made to qualitatively link the permanent sCO2 effects to polymer micro- structure, free volume, backbone substitutions, presence of polar groups, and degree of crystallinity differences. This study has established that soft polymeric materials are conducive to failure in sCO2 through mechanisms of failure that are dependent on polymer microstructure and chemistry. Polar pendant groups, large atom substitutions on the backbone are some of the factors that are influential structural factors.

36 MATERIALS SCIENCE↗

DeepAdversaries: examining the robustness of deep learning models for galaxy morphology classification

With increased adoption of supervised deep learning methods for work with cosmological survey data, the assessment of data perturbation effects (that can naturally occur in the data processing and analysis pipelines) and the development of methods that increase model robustness are increasingly important. In the context of morphological classification of galaxies, we study the effects of perturbations in imaging data. In particular, we examine the consequences of using neural networks when training on baseline data and testing on perturbed data. We consider perturbations associated with two primary sources: (a) increased observational noise as represented by higher levels of Poisson noise and (b) data processing noise incurred by steps such as image compression or telescope errors as represented by one-pixel adversarial attacks. We also test the efficacy of domain adaptation techniques in mitigating the perturbation-driven errors. We use classification accuracy, latent space visualizations, and latent space distance to assess model robustness in the face of these perturbations. For deep learning models without domain adaptation, we find that processing pixel-level errors easily flip the classification into an incorrect class and that higher observational noise makes the model trained on low-noise data unable to classify galaxy morphologies. On the other hand, we show that training with domain adaptation improves model robustness and mitigates the effects of these perturbations, improving the classification accuracy up to 23% on data with higher observational noise. Domain adaptation also increases up to a factor of ${\approx}2.3$ the latent space distance between the baseline and the incorrectly classified one-pixel perturbed image, making the model more robust to inadvertent perturbations. Successful development and implementation of methods that increase model robustness in astronomical survey pipelines will help pave the way for many more uses of deep learning for astronomy.

79 ASTRONOMY AND ASTROPHYSICS↗

Electrolyte Assisted Hydrogen Storage Reactions (Final Technical Report)

The goal of this project is to address critical deficiencies of hydrogen storage systems design based on hydride materials, as originally identified by the DOE Hydrogen Storage Engineering Center of Excellence (HSECoE). Baseline hydrogen storage technology presently relies on compressed gas operating at ~700 bar pressure, which imposes huge demands on fuel delivery, fuel storage and system cost. For onboard storage applications, Type IV composite overwrapped pressure vessels and associated balance-of-plant components are necessary to ensure safe and effective fuel delivery. However, such compressed gas technology falls well-short of volumetric targets even at 700 bar, given the density of gaseous molecular hydrogen is only 40 g·H2/l at ambient temperatures. One alternative is to utilize hydride materials which accommodate hydrogen in atomic form. Certain hydrides can attain volumetric densities that exceed the density of liquid H2 (71 g·H2/l) while also offering advantageous thermodynamic properties. However, such material systems presently rely on solid-state diffusion for hydrogen release, which has a very high activation barrier for atom mobility and, thus, requires impractically-high temperatures for operation. The initial focus of our research effort is to employ and demonstrate an electrolyte system to mediate the diffusion of species at lower temperatures relevant to transportation applications, with the goal of establishing the critical factors necessary to obviate the need for high-temperature release of hydrogen. A parallel goal of this exploratory effort is to determine the effectiveness of modest electrochemical potentials in overcoming any endothermic requirements for hydrogen release in a similar electrolyte-promoted scenario.

08 HYDROGEN↗

Valley-Tunable Even-Denominator Fractional Quantum Hall State in the Lowest Landau Level of an Anisotropic System

Fractional quantum Hall states (FQHSs) at even-denominator Landau level filling factors (v) are of prime interest as they are predicted to host exotic, topological states of matter. Here, we report here the observation of a FQHS at v = 1/2 in a two-dimensional electron system of exceptionally high quality, confined to a wide AlAs quantum well, where the electrons can occupy multiple conduction-band valleys with an anisotropic effective mass. The anisotropy and multivalley degree of freedom offer an unprecedented tunability of the v = 1/2 FQHS as we can control both the valley occupancy via the application of in-plane strain, and the ratio between the strengths of the short- and long-range Coulomb interaction by tilting the sample in the magnetic field to change the electron charge distribution. Thanks to this tunability, we observe phase transitions from a compressible Fermi liquid to an incompressible FQHS and then to an insulating phase as a function of tilt angle. We find that this evolution and the energy gap of the v = 1/2 FQHS depend strongly on valley occupancy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Carbon dioxide pipeline network transportation cost model: evaluating economic and geographic factors for efficient carbon capture, storage, and utilization

This study presents a comprehensive pipeline network modeling framework to estimate the CO 2 delivery cost for CO 2 utilization and geologic CO 2 storage across the United States. We developed a Python-based CO 2 pipeline transportation cost model leveraging Argonne National Laboratory’s pipeline engineering expertise and detailed natural gas transmission pipeline cost data across U.S. regions. Using existing road corridors as practical routing guides, the model designs pipeline networks that aggregate CO 2 from one or multiple sources and deliver it to selected destinations. It then minimizes the total transportation cost by optimizing pipeline diameters and incorporating booster pumps. A key contribution is the incorporation of up-to-date, region-specific cost factors with itemized components for materials, labor, miscellaneous construction expenses, and right-of-way acquisition. Results emphasize that regional variation and economies of scale associated with CO 2 pipeline costs are significant and should be explicitly accounted for in screening and planning studies. By combining realistic routing constraints with regionalized cost inputs, the model provides transparent design methodology and location-specific insights into source–destination delivery costs, including the effects of routing complexity along existing road networks. We demonstrate the model with two illustrative case studies – one for CO 2 storage and one for CO 2 utilization – in which the model designs pipeline networks spanning hundreds of miles across the states, collecting CO 2 from multiple sources and delivering it to designated endpoints while minimizing levelized cost of delivery via diameter and compression optimization. The model offers a practical, scalable approach for alternative design option screening and early-stage CO 2 transportation planning.

CCS↗

Interstitial Atoms and the Frustrated and Allowed Structural Transitions Principle: Tunability in the Electronic Structure of AuCu 3 ‐type Frameworks in Dy 4 T 1− x Ga 12 (T = Ag, Ir)

In this Article, we explore how the chemical pressure (CP) features of an intermetallic phase may provide opportunities to couple perturbations in electron count with the stabilization of the underlying geometrical structure. AuCu 3 ‐type LnGa 3 (Ln = lanthanide or group 3 metal) phases contain octahedral cavities of negative CP held open by overly compressed Ln–Ga interactions, leading to a series of transition metal‐stuffed derivatives. We present new additions to this family with the synthesis and crystal structures of Dy 4 T 1−x Ga 12 with (T, x) = (Ag, 0.29) and (Ir, 0.15), adopting Y 4 PdGa 12 ‐type superstructures of the AuCu 3 ‐type. density functional Ttheory (DFT)‐CP calculations, when adjusted to avoid dipolar CP features, affirm that T atom incorporation provides a mechanism for the relief of packing tensions, while electronic density of states distributions illustrate that the T atoms serve largely as electron or hole donors to the band structure, as needed for them to attain d 10 configurations. The maximum obtainable value for x may be limited by a mismatch between the Fermi energy and pseudogap, in line with the balance of factors envisioned by the frustrated and allowed structural transitions principle. Furthermore, trends in resistivity measurements on T = Ir, Pd, and Ag compounds are interpretable in terms of the varying degrees of disorder arising from x < 1.0.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon Fiber Prepreg Recycling- Automated Perform Manufacturing Equipment (Final Technical Report)

Discontinuous fiber composites have a compelling business case as related to light metals such as aluminum and magnesium machined or die cast components. DFC’s have been in commercial production using carbon fiber/epoxy prepreg for both commercial aircraft and for military applications for a number of years and their usage is expanding into automotive and consumer goods components. Reasons for this include significantly lower cost as compared to continuous fiber composites easier design and understanding of the material by metal traditionalists, higher degree of part complexity with relatively simple and low cost molds when compared to high-pressure die casting, and potential for near zero waste. IACMI’s 10-year goals of reduction in carbon fiber reinforced plastics manufacturing cost by 25%, reduction in their embodied energy by 50%, and recyclability of 95% can all be advanced through robust recycling. With recycling We have shown that a 15% contribution in manufacturing cost reduction is feasible, that a significant reduction in embodied energy (over 90% for the recycled fraction) is achievable, and that when combined with primary pyrolysis recycling, the 95% target can be met. Large-scale application of advanced composites began with the aerospace industry, and CRTC’s location in Washington State provided access to a waste stream of nearly 900 tonnes/year that was going to landfill. As the aerospace industry has automated parts production, their material forms have become significantly easier to recycle, and some evidence of this exists for the future automotive and wind turbine components as well. CRTC decided to focus on re-purposing the carbon fiber/epoxy aerospace scrap streams through materials reformatting and reuse, and not on primary recycling (pyrolysis mainly) as this technology is commercially available. The main reasons for this were that it offers a route to effective use of the most widely available and consistent feedstock called edge-trim, as well as allows creation of complex and low-cost parts using rapid and cost-effective compression molding technologies. The reasons also include leading the way for potential automotive pre-preg recycling as the combination of roughly 85% virgin pre-preg and 15% recycled pre-preg directly supports IACMI’s goals and enables significant part cost reduction and a path to zero-waste composite manufacturing at the OEMs or Tier 1’s. The purpose of the project was to develop equipment and technology to transform continuous fiber pre-preg in roll forms into smaller chip-format DFC feedstock, and to determine impacts on moldability and performance based on lessons learned along the path. The availability of scrap from pre-preg manufacturers was a significant value to the program even though the aerospace pre-preg is designed to be autoclave molded thus has handling, tackiness, and processing factors that are not ideal for DFCs. The scope of the project was to develop a Phase 1 system that could identify and resolve some of the major hurdles, determine where to focus effort on building a commercial-scale Phase 2 system, and develop technology for handling a wide variety of incoming scrap materials to produce the most usable product. A chopping system was developed and reduced to practice via an experimental approach that enabled all forms of scrap coming from the pre-preg manufacturers and some forms of the scrap coming from OEM production lines to be very effectively re-purposed. The systems that were developed were not considered to be production robust but rather identified and resolved some of the major technical hurdles, and demonstrated the path towards effective full-scale commercialization.

36 MATERIALS SCIENCE↗

Elucidating the mechanical and thermal response of nanotwinned Ni alloys

Transformative advances in nanoscale materials synthesis, characterization and modeling are enabling the synthesis of materials with unprecedent properties and greater understanding of the nanoscale mechanisms that underpin these properties. Nanotwinned Cu alloys have received considerable attention due to their impressive balance of strength and ductility, but they have limited microstructural stability. Recent instantiation of nanotwins in sputter deposited Ni-Mo-W, and several commercial Ni-based superalloys, point to the potential development of a new class of high temperature materials with a very beneficial suite of mechanical and physical properties. The experimental study described in this report was undertaken to elucidate the nanoscale origins of the thermal and mechanical behavior of nanotwinned Ni alloys. Micropillar compression experiments have shown nanotwinned Ni 85 Mo 15-x W x alloys to possess unusually high strengths above 3.5GPa and enhanced microstructural stability. The strength of these nanotwinned alloys is determined by the abrupt formation of shear bands. This study focused on identifying the nanoscale trigger, or triggers, for shear banding in nanotwinned materials using in situ experiments and atomic-scale postmortem analysis. Contrasting and comparing the mechanical response and postmortem nanostructures of “Mo-rich” and “W-rich” nanotwinned specimens elucidated the importance of the lateral motion of easy glide defects that results in erasure of twins and local coarsening of the nanotwinned microstructure. This twin coarsening was then associated with very localized plasticity, strain softening, and shear band formation. The difference between alloys was further studied by considering the role of various material factors: stacking fault energy, coherent twin boundary spacing and flatness, grain size, and the interactions of solute atoms with twin and grain boundaries. The effect of alloy composition and sputtering power and temperature were also investigated and used to control twin spacing and geometry. Combined with atomic-scale simulations, these experiments insights should allow us to identify strategies for achieving concomitant ultrahigh strength and ductility in nanotwinned Ni alloys. In parallel but synergistic studies, our work was buoyed by collaborative state-of-the-art nanoscale orientation and strain mapping at the DOE National Center for Electron Microscopy (NCEM). For example, novel 4D-STEM techniques were used to acquire nanoscale strain maps. The extremely fine twin spacing limited our measurements of atomic-scale thermal expansion within twins, but recent results from NCEM suggest that it will soon be possible to conduct such measurements and to unravel the origins of the novel thermal expansion characteristics of nanotwinned Ni alloys.

36 MATERIALS SCIENCE↗

Carbon Fiber Reinforced Plastic with Substantially Improved Through-plane Thermal Conductivity

Thermal conductivity of carbon fiber support structures that also serve as cooling substrates is of increasing importance for high energy physics detectors. One of the shortcomings of current carbon fiber laminates is that they exhibit excellent thermal conductivity of order several hundred W/(m*K) along the fiber direction, whereas the through plane thermal conductivity perpendicular to the fibers is orders of magnitude lower, and at about 1 W/(m*K) more in the range of thermal insulators. This proposal aims to significantly improve the through plane thermal conductivity, by at least a factor five, changing CFRP from being a thermal insulator to a thermal conductor in the through plane direction. Between June 2021 and December 2022, a total of 36 carbon fiber laminate test samples and two carbon loaded epoxy test samples were manufactured and analyzed for thermal performance. Unfortunately, the promising results obtained with a similar sample manufactured in 2017 were not reproduced. We believe this to be mainly due the excessive age of the carbon fiber prepreg material used for this SBIR - the same material batch, produced in 2015, was used in 2017 and in 2021/2022, i.e. the material was seven years old when the final samples were produced - , and to a lesser degree to the difference in pressure attainable with the 2017 and 2021/2022 laminate curing setups. While disappointing, the results are not completely discouraging, since for the samples with the highest compression the through-plane thermal conductivities plotted as a function of per-ply thickness extrapolate well to the data point from the 2017 sample.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization

Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. Here, we demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss–Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.

97 MATHEMATICS AND COMPUTING↗

An indirect approach to optimize the reaction rates of thermal NO formation for diesel engines

With stringent emission regulations, it has become more important for modern diesel engine manufacturers to accurately predict engine-out nitrogen oxide (NO x ) emissions across a wide range of operating conditions. Thermal NO is the major source of engine-out NO x in modern diesel engines. For thermal NO formation, several earlier studies have recommended the forward and reverse reaction rate coefficients of the rate-limiting reaction (O + N 2 ⇌ NO + N). However, due to deficiencies in sub-models and inadequacies of reduced chemical mechanisms to represent diesel combustion, these recommended values more often than not need to be adjusted in reduced order combustion models to accurately predict engine-out NO x . Hence, in this work a systematic and computationally efficient approach has been proposed to streamline the process of determining the optimum reaction rate coefficients. Here, to develop the optimization approach, four different production diesel engines with different operating conditions in terms of speed, load, and exhaust-gas recirculation have been considered. Numerical simulations have been performed using a detailed zero-dimensional velocity-composition-frequency transported probability density function (0D-VCF-tPDF) model that uses hundreds of notional particles to capture in-cylinder stratification. Four different combinations of hydrocarbon and NO x chemical mechanisms were used to represent chemistry. It was found that for the rate-limiting reaction, the pre-exponent factors (A f1 , A r1 ) and activation energies (E A,f1 , E A,r1 ) of the forward and reverse reaction rates follow a linear band in A f1 - E A,f1 and A r1 - E A,r1 space where predicted engine-out NO x match the measured values closely. By encompassing such bands from different engines and considering constraints on activation energies, a reduced search domain of pre-exponent factors and activation energies was constructed that is expected to be applicable to any diesel engine. Eventually, computationally efficient three-line and one-line search approaches were proposed to determine the optimum values of the pre-exponent factors and activation energies that led to a minimum error between measured and predicted engine-out NO x . Finally, these three-line and one-line NO x optimization approaches were applied to a fifth production diesel engine for which the 0D-VCF-tPDF model showed a very good predictive performance in terms of predicting peak pressure, 50% burn rate, and engine-out NO x when compared to measured and 3D-CFD values.

33 ADVANCED PROPULSION SYSTEMS↗

Dynamic recrystallization of a wrought magnesium alloy: Grain size and texture maps and their application for mechanical behavior predictions

Effect of dynamic recrystallization (DRX) on the grain refinement and texture modification was studied by conducting a series of hot compression on AZ31B Mg alloy. Processing-grain size-texture maps were established as a function of temperature, strain rate, strain, and the Zener-Hollomon parameter (Z). Moreover, influence of simultaneous changes in the grain size and texture on tensile yield strength and ductility of the hot-worked Mg alloy was studied using the grain size and texture maps established. The processing-grain size map showed that the DRX grain size decreased with the increase in Z, also revealing various characteristics ranging from a grain growth, grain refinement, to a bimodal distribution of ultrafine grains and partially recrystallized grains. The effect of twinning on the grain refinement was also evident at high Z conditions. The processing-texture map revealed that the initial fiber texture was altered significantly to a shear, off-normal, or extension-twin texture with the increase in Z. The Schmid-factor maps were calculated to identify a dominant deformation mode during a subsequent tensile deformation of the hot-worked samples. Finally, a corresponding Hall-Petch relationship for each dominant deformation mode was used to establish tensile yield-strength maps, which agree well with the measured data.

36 MATERIALS SCIENCE↗

Beryllium–tungsten graded density inner shells in double shell capsules for improved hydrodynamic stability

The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗