Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical solution”

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,297 records · Page 72

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗

Metallic Technetium Sequestration in Nickel Core/Shell Microstructure during Fe(OH)2 Transformation with Ni doping

This study investigates the impacts of Ni doping on technetium-99 (Tc) sequestration in aqueous solutions through transformation of Fe(OH)2(s) to iron spinel (magnetite) under alkaline conditions. Extensive solid characterization was performed for the mineral phases produced, as well as the Tc/Ni speciation and distribution within these phases. X-ray diffraction results show that iron spinel was the dominant mineral product without detectable Ni incorporation. The doped Ni ions mainly precipitated as fine Fe/Ni oxide/hydroxide particles, including strongly reduced nanometer?sized spheroidal Ni-rich and metallic Ni phases. High-resolution analytical scanning transmission electron microscopy using energy dispersive X-ray spectroscopy and electron energy loss spectroscopy on the produced solid samples (focused ion beam-prepared specimens) revealed three Tc distribution domains dominated by nanocrystals and, especially, a Tc-rich metallic phase. Instances of metallic Tc were specifically found in spheroidal, Ni-rich and metallic nanoparticles exhibiting a core/shell microstructure that suggests strong reduction and sequential precipitation of Ni-Tc-Ni. Mass balance analysis showed nearly 100% Tc removal from the 4.8 × 10-4 M Tc solutions. The finding of the metallic Tc encapsulation indicates that Tc sequestration through Ni-doped Fe(OH)2(s)?to?iron spinel transformation process likely provides an alternative treatment pathway for Tc removal and could be combined into further waste treatment approaches.

Wang, Guohui↗

Scattering in black hole backgrounds and higher-spin amplitudes. Part II

We continue to investigate correspondences between, on the one hand, scattering amplitudes for massive higher-spin particles and gravitons in appropriate quantum-to-classical limits, and on the other hand, classical gravitational interactions of spinning black holes according to general relativity. We first construct an ansatz for a gravitational Compton amplitude, at tree level, constrained only by locality, crossing symmetry, unitarity and consistency with the linearized-Kerr 3-point amplitude, to all orders in the black hole’s spin. We then explore the extent to which a unique classical Compton amplitude can be identified by comparing with the results of the classical process of scattering long-wavelength gravitational waves off an exact Kerr black hole, determined by appropriate solutions of the Teukolsky equation. Up to fourth order in spin, we find complete agreement with a previously conjectured exponential form of the tree-level Compton amplitude. At higher orders, we extract tree-level contributions from the Teukolsky amplitude by an analytic continuation from a physical (a/GM < 1) to a particle-like (a/GM > 1) regime. Up to the sixth order in spin, we identify a unique conservative part of the amplitude which is insensitive both to the choice of boundary conditions at the black hole horizon and to branch choices in the analytic continuation. The remainder of the amplitude is determined modulo an overall sign from a branch choice, with the sign flipping under exchanging purely ingoing and purely outgoing boundary conditions at the horizon. Along the way, we make contact with novel applications of massive spinor-helicity variables pertaining to their relation to EFT operators and (spinning) partial amplitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

COMPOFF: A Compiler Cost model using Machine Learning to predict the Cost of OpenMP Offloading

The HPC industry is inexorably moving towards an era of extremely heterogeneous architectures, with more devices configured on any given HPC platform and potentially more kinds of devices, some of them highly specialized. Writing a separate code suitable for each target system for a given HPC application is not practical. The better solution is to use directive-based parallel programming models such as OpenMP. OpenMP provides a number of options for offloading a piece of code to devices like GPUs. To select the best option from such options during compilation, most modern compilers use analytical models to estimate the cost of executing the original code and the different offloading code variants. Building such an analytical model for compilers is a difficult task that necessitates a lot of effort on the part of a compiler engineer. Recently, machine learning techniques have been successfully applied to build cost models for a variety of compiler optimization problems. In this paper, we present COMPOFF, a cost model which uses the multi-layer perceptrons to statically estimates the Cost of OpenMP OFFloading. We used six different transformations on a parallel code of Wilson Dslash Operator to support GPU offloading, and we predicted their cost of execution on different GPUs using COMPOFF during compile time. Our results show that this model can predict offloading costs with a root mean squared error in prediction of less than 0.5 seconds. Our preliminary findings indicate that this work will make it much easier and faster for scientists and compiler developers to port legacy HPC applications that use OpenMP to new heterogeneous computing environment.

97 MATHEMATICS AND COMPUTING↗

Improving and Automating Building Model Data Exchange

There are many instances throughout a project’s lifecycle where there arises a need for quick and accurate risk assessment of building designs. For example, an unexpected design change during construction may necessitate structural engineers to perform a seismic risk assessment on analytical models of the updated building design using high fidelity structural analysis software, such as ANSYS or Abaqus. However, the efficiency of such workflows often depends upon the interoperability of architectural design software and structural analysis software. When the quality of this interoperability is lacking or even non-existent, the efficiency of virtual engineering workflows is hampered, which increases project costs. A McGraw Hill industry survey of professional users of Building Information Modeling (BIM) technologies found that there is high demand for BIM interoperability for structural analysis, but that the value/difficulty ratio is currently too low for practical use. There have been efforts by the academic community to facilitate model data exchange between the architectural design and structural analysis domains, but such solutions have not been widely adopted by industry, face technical challenges, and oftentimes are limited in applicability for users of various BIM software. Therefore, INL is developing capabilities to improve, automate, and generalize model data exchange between architectural BIM software (e.g., Revit) and structural analysis software (e.g., SAP2000, ANSYS). The goal is to help expedite and automate as much of the pre-processing step for creating analytical models in finite element analysis software as reasonably as possible. Such a "BIM-to-FEA" conversion tool should provide direct benefit to end-users through accuracy, automation, quick turn-around, and wide applicability. To generalize the application of this BIM-to-FEA conversion tool and increase its useability among the many different commercial BIM software currently used by industry, the program is being developed with the concept of openBIM. OpenBIM is the application of non-proprietary, open data standards that allow for BIM model data exchange in a format that is accessible, retainable, and useable for all users. The most widely used open, non-proprietary data exchange format for BIM is the Industry Foundation Classes (IFC) schema. IFC is developed by buildingSMART international and is ISO certified (ISO 16739-1:2018). The BIM-to-FEA conversion tool is being developed for compatibility with typical commercial building designs of steel framed structures. The tool is currently capable of importing architectural BIM data of framed building structures, recognizing and extracting the aspects of the model that are required for structural analysis, adjusting the connectivity of frame members, and finally exporting to an analytical model stored in the IFC format. The exported IFC analytical model can then be imported into various openBIM compliant software, such as SAP2000. Such capabilities have already been tested on commercial software, as shown above, and continue to be improved. Work is underway to test the conversion on various commercial BIM software, develop a user-friendly interface, incorporate the program into the broader DeepLynx data warehouse project being developed by INL, and to eventually open-source the tool for the benefit of the community. Future development of the tool envisions the ability for efficient iterative risk assessment of generative building designs, all within a workflow utilizing open-source tools. One such open-source tool will be MOOSE, an advanced finite element analysis tool developed at INL. The conversion tool will also branch out from typical commercial building designs and will aim to incorporate nuclear construction. The aim will be to convert both structural and non-structural components of nuclear facilities, such as curved concrete containment structures and piping systems, respectively.

97 MATHEMATICS AND COMPUTING↗

In‐Situ Product Removal for the Enzymatic Depolymerization of Poly(ethylene terephthalate) via a Membrane Reactor

Poly(ethylene terephthalate) (PET) is a common single-use plastic and a major contributor to plastic waste. PET upcycling through enzymatic depolymerization has drawn significant interests, but lack of robust enzymes in acidic environments remains a challenge. This study investigates in-situ product removal (ISPR) of protons and monomers from enzymatic PET depolymerization via a membrane reactor, focusing on the ICCG variant of leaf branch compost cutinase. More than two-fold improvements in overall PET depolymerization and terephthalic acid yields were achieved employing ISPR for an initial PET loading of 10 mgPET ml buffer −1 . The benefit of ISPR was reduced for a lower initial loading of 1 mgPET ml buffer −1 due to decreased need for pH stabilization of the enzyme-containing solutions. A back-of-envelop analysis suggests that at a modest dilution ratio, ISPR could help achieve savings on caustic base solutions used for pH control in a bioreactor. Our study provides valuable insights for future ISPR developments for enzymatic PET depolymerization, addressing the pressing need for more sustainable solutions towards plastic recycling and environmental conservation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanoscale serration characteristics of additively manufactured superalloys

Structural elements made of nickel-based superalloys usually operate at high temperatures. Many surface failure mechanisms such as creep, fatigue, fretting fatigue, corrosion, and stress corrosion cracking start from the surface. Thus, studying the surface mechanical properties and deformation behavior of materials is vital in order to draw a correlation between the surface properties and failures. Herein, a nanoindentation technique was implemented to characterize surface properties of a widely used additively manufactured (AM) superalloy, i.e., Inconel 625. The specimens were tested in a range of temperatures in a vacuum chamber. Serrated plastic flow characterized by pop-in events during nanoindentation, known as Portevin-Le Chatelier (PLC) effect, was observed in all three tested elevated temperatures. This phenomenon depicts itself as bursts of plasticity in the loading section of the load-displacement curves. These bursts were studied comprehensively to explore incipient plasticity. Hertzian contact mechanics were implemented to extract the maximum shear stress beneath the indenter from the pop-in loads. The results show that the initial pop-in load increases as the temperature increases. The average initial pop-in load increases by almost four times from 300 °C to 650 °C. This was attributed to the formation of strengthening precipitates. The average size of the serrations and the magnitude of dislocation nucleation increase from 300 °C to 500 °C and decrease from 500 °C to 650 °C. Finally, the serrations are attributed to the dislocation generation and movement as well as their interaction with the solute atoms and precipitates.

36 MATERIALS SCIENCE↗

Overview and technology opportunities for thermochemically-produced bio-blendstocks

Global demand for transportation fuels is projected to increase 40% by 2040, and biomass-derived fuels (biofuels) play a crucial role in substituting fossil fuels and mitigating greenhouse gas emissions. Currently, biofuels are mainly consumed as blendstocks combined with petroleum-based fuels, and effective conversion technologies can address the quality challenges for offering standalone biofuels. Thermochemical conversion process is one of the most promising pathways among existing technologies for biofuel production. However, the major barriers are unwanted characteristics (e.g., thermal instability) of intermediate products, such as bio-oil, and required upgrading treatments for producing compatible fuels. Here, this study highlights the merits and critical challenges of thermochemical conversion and physicochemical upgrading technologies for bio-blendstock production from lignocellulosic biomass. The novelty of this study lies in potential directions for future research through both critical and systematic literature reviews, and the proposed intensified process for lignocellulosic-based fuel blendstocks production. It is concluded that recovery and fractionation strategies (e.g., quenching and stripping) can maximize process yields and add values in the efficient conversion pathways. Effective quenching can stop secondary free radical reactions and improve liquid yields over gas and solid yields. Stripping process can improve process yield, catalyst lifespan, and thermal stability. It is further concluded that physicochemical treatments are not as effective as thermochemical treatments, but have advantages of mild operating conditions and potential for integrated solutions in conjunction with other treatments.

09 BIOMASS FUELS↗

Sorption kinetics and stability of conventional adsorbents for mercury remediation

In-situ remediation of mercury at numerous contaminated sites worldwide is a challenging and costly endeavor due to the persistency of this contaminant. In this study, we evaluated eight commercially available sorbent media ranging from carbon-, clays- and silica-based materials (PBC– Biochar, eSorb – Sorbster, nsPAC – Powdered Activated Carbon, fsPAC – Powdered Activated Carbon with Mackinawite, F300 – Filtrasorb 300, Si-SH – Silica Thiol, eBind – RemBind, Q-Clay – Organoclay PM-199), for their effectiveness in sorbing mercury (Hg 2+ ) and mercury complexed with dissolved organic matter (Hg-DOM). Under the chosen experimental conditions of this study, results showed that in the absence of DOM, the kinetic rates of Hg 2+ sorption onto the evaluated sorbents were in the order of 0.31 min-1 (Si-SH) to 2.98 min -1 (nsPAC), whereas in the presence of DOM, the rates varied from 0.16 min-1 (F300) to 0.95 min -1 (nsPAC). The measured sorption capacity for Hg 2+ in the absence of DOM varied from 3.02 mg/g (Q-Clay) to 35.15 mg/g (Si-SH), whereas in the presence of DOM, calculated partition coefficient (KD) ranged from 69.7 mL/g (Q-Clay) to 41,510 mL/g (Si-SH). Furthermore, kinetic data suggest liquid film diffusion was the rate-limiting steps governing mercury sorption onto the studied media. Overall, the obtained study parameters (kinetics/isotherm) are particularly important in informing robust engineering designs for deployment of the vast majority of evaluated sorbents. Thus, sorbent-based strategies offer viable solutions for cost-effective cleanup of mercury at industrially contaminated sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrochemical, thermodynamic, and physical properties of tetradecyltrihexylphosphonium ([P 6,6, 6,14 ] + ) and methyl-propyl piperidinium containing ionic liquids and their propylene carbonate solutions

Here, the viscosity, conductivity, electrochemical window and related thermodynamic properties such as excess volume and dynamic viscosity deviation of two phosphonium ionic liquids were measured and calculated for “as-supplied” and dried neat liquids and also v/v mixtures of 99/1, 95/5, 90/10, 75/25, and 50/50 of the ionic liquids with propylene carbonate (PC). Tetradecyltrihexylphosphonium dicyanamide ([P 6,6,6,14 ] + dicyanamide), Cyphos 105, and tetradecyltrihexylphosphonium bis(trifluororomethanesulfonyl)imide (called bistriflimide) ([P 6,6,6,14 ] + bistriflimide), Cyphos 109, ionic liquids were studied. Generally speaking, there were slight differences in the measured properties of the wet and dried IL solutions which were reflected in differences if the calculated properties. The measured and calculated properties were compared with those collected and derived from a piperidinium-based ionic liquid, methyl-propyl piperidinium bistriflimide. Arrhenius plots for both the viscosity and conductivity were linear, with the linearity of the Litovitz plots being slightly higher. The viscosity and conductivity of the phosphonium solutions yielded Walden Plots that differed from Walden plots of piperidinium solutions and other ionic liquid solutions. For the phosphonium based ionic liquids, the ionicity was found to increase with increasing dilution in propylene carbonate. The electrochemical windows of the ionic liquids were determined as a function of the concentration of the ionic liquid present in the solution. The size of the window generally initially decreased and then increased with the addition of PC for the [P 6,6,6,14 ] + dicyanamide ionic liquid increasing from 3.4 V in the 95/5 v/v (volume Cyphos 105) to volume of PC solvent) to 3.7 in the 50/50 dried IL case. The electrochemical windows of the dried liquids were slightly larger 3.5 V for the 95/5 v/v and 4.2 V for the 50/50 v/v Cyphos 105/PC solutions. The electrochemical window was larger and the conductivity was higher for the piperidinium IL solutions, but since the ionicity increases with dilution, the phosphonium based IL solutions may prove favorable for processes such as electrochemical reduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Does Water Enhance Mg Intercalation in Oxides? The Case of a Tunnel Framework

The presence of H 2 O has been linked to enhancements in the reactivity of cathodes for Mg 2+ electrochemistry. If the enhancements were mimicked by nonaqueous solvents, they could enable Mg batteries with transformational energy density. However, the extent to which H 2 O may boost actual intercalation of Mg 2+ , as opposed to competing reactions, has not been elucidated. Here, in this paper, we evaluate its role as additive in the electrochemistry of a tunnel polymorph of V 2 O 5 in a nonaqueous Mg 2+ electrolyte. The electrochemical response and V reduction in the cathodes positively correlated with H 2 O concentration, but it was not concurrent with commensurate changes in cell volume and Mg content. These observations indicate that H 2 O does not enhance Mg 2+ intercalation, but rather, it promotes competing pathways. This work shows the importance of accurately probing reactions in multivalent electrolytes. Importantly, it indicates that H 2 O is not a universal solution to the challenge of Mg 2+ intercalation in oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bridging molecular-scale interfacial science with continuum-scale models

Solid–water interfaces are crucial for clean water, conventional and renewable energy, and effective nuclear waste management. However, reflecting the complexity of reactive interfaces in continuum-scale models is a challenge, leading to oversimplified representations that often fail to predict real-world behavior. This is because these models use fixed parameters derived by averaging across a wide physicochemical range observed at the molecular scale. Recent studies have revealed the stochastic nature of molecular-level surface sites that define a variety of reaction mechanisms, rates, and products even across a single surface. To bridge the molecular knowledge and predictive continuum-scale models, we propose to represent surface properties with probability distributions rather than with discrete constant values derived by averaging across a heterogeneous surface. This conceptual shift in continuum-scale modeling requires exponentially rising computational power. By incorporating our molecular-scale understanding of solid–water interfaces into continuum-scale models we can pave the way for next generation critical technologies and novel environmental solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Study of perovskite CsPbBr 3 detector polarization and its mitigation with ultrahigh x-ray flux

High-flux capable semiconductor x-ray detectors are essential in various applications, but the detrimental effects of detector polarization limit their use in many cases. Here, for this work, we studied the polarization of perovskite CsPbBr 3 semiconductor detectors using ultrahigh flux synchrotron x rays (10 6 –10 12 photons s –1 mm –2 at 58.61 keV). The CsPbBr 3 detectors did not show immediate polarization prominently until a flux higher than 10 10 photons s –1 mm –2 . Using the pump-and-probe technique, we visualized the spatial and temporal effects of polarization. The polarized region, represented by reduced photocurrent, extended beyond the area under direct irradiation, and the reduced photocurrent persisted after potential de-polarization treatments. We found that stronger applied electric fields and fewer carrier traps can mitigate polarization, represented by less photocurrent deficit. By examining the detectors’ current response under controlled ambient light, low, and high-flux x rays, we studied the trap filling and release behavior of CsPbBr 3 . We discovered that the polarization is caused by partial detector damage due to deep defects generated by ultrahigh flux x-ray irradiation. Our work provides insight into CsPbBr 3 polarization under extremely intense x-ray radiation and shows that reducing crystal defects and increasing detector bias voltage are effective solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-speed two-dimensional terahertz spectroscopy with echelon-based shot-to-shot balanced detection

By using a reflective-echelon-based electro-optic sampling technique and a fast detector, we develop a two-dimensional terahertz (THz) spectrometer capable of shot-to-shot balanced readout of THz waveforms at a full 1-kHz repetition rate. To demonstrate the capabilities of this new detection scheme for high-throughput applications, we use gas-phase acetonitrile as a model system to acquire two-dimensional THz rotational spectra. The results show a two-order-of-magnitude speedup in the acquisition of multidimensional THz spectra when compared to conventional delay-scan methods while maintaining accurate retrieval of the nonlinear THz signal. Our report presents a feasible solution for bringing the technique of multidimensional THz spectroscopy into widespread practice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of distribution coefficients of 41 elements in nitric and hydrochloric acids for actinide resin

The distribution coefficients of 41 elements were determined on Eichrom’s Actinide resin in HCl and HNO 3 solutions of molarities between 0.01 M and 9 M. Herein a separation scheme for these 41 elements was designed utilizing these distribution coefficients that combines the Actinide resin with Eichrom’s TRU resin. Overall, these values improve the ability to create separations targeting specific analytes utilizing Actinide resin.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Bubble departure and sliding in high-pressure flow boiling of water

Bubble growth, departure and sliding in low-pressure flow boiling has received considerable attention in the past. However, most applications of boiling heat transfer rely on high-pressure flow boiling, for which very little is known, as experimental data are scarce and very difficult to obtain. In this work, we conduct an experiment using high-resolution optical techniques. By combining backlit shadowgraphy and phase-detection imaging, we track bubble shape and physical footprint with high spatial ($6\,\mathrm {\mu }{\rm m}$) and temporal ($33\,\mathrm {\mu }{\rm s}$) resolutions, as well as bubble size and position as bubbles nucleate and slide on top of the heated surface. We show that at pressures above 1 MPa bubbles retain a spherical shape throughout the growth and sliding process. We analytically derive non-dimensional numbers to correlate bubble velocity and liquid velocity throughout the turbulent boundary layer and predict the sliding of bubbles on the surface, solely from physical properties and the bubble growth rate. We also show that these non-dimensional solutions can be leveraged to formulate elementary criteria that predict the effect of pressure and flow rate on bubble departure diameter and growth time.

Mechanics↗

Analytical Tool to evaluate Heterogeneous Neuromorphic Architectures (ATHENA)

SAND2022-7062 O The Analytical Tool to Evaluate Heterogeneous Neuromorphic Architectures (ATHENA) quickly evaluates performance metrics like energy, area, and latency for an AI/ML network. It currently supports the evaluation of analog neural networks. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Hahn, Andrew↗

SPARTAN (Scalable Probabilistic Application Reconfigurable Tensor Autonomous Network)

The technical founder of Ludwig Computing Inc has been competitively selected for support by Cyclotron Road, a U.S. Department of Energy (DOE) Advanced Manufacturing Office (AMO) Lab-Embedded Entrepreneurship Program (LEEP) through an approved merit review process. Ludwig Computing Inc, supported by the U.S. Department of Energy's Advanced Manufacturing Office through the Cyclotron Road program, has investigated the advantages of probabilistic computing for real-world compute-intensive applications. This research adds to the understanding of alternative computing paradigms by exploring a unique hardware-software co-design that integrates quantum computing methods with nature-inspired problem-solving techniques. The project's focus on areas such as combinatorial optimization, graph analytics, and machine learning demonstrates the potential for significant advancements in computational efficiency and performance. By harnessing natural randomness to streamline large circuits into fewer devices, Ludwig's approach enables massive parallelism, potentially offering higher throughput, speed, and energy efficiency compared to conventional hardware solutions. This work benefits the public by paving the way for more efficient computing solutions that could address complex real-world problems while potentially reducing energy consumption in data-intensive industries.

97 MATHEMATICS AND COMPUTING↗