SEARCH · Engineering Papers
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.
Acute esophagitis (AE) is a common dose-limiting toxicity in radiation therapy of locally advanced non-small cell lung cancer (LA-NSCLC). We developed an early AE prediction model from weekly accumulated esophagus dose and its associated local volumetric change.
Definitive, nonsurgical management of gynecologic malignancies involves external beam radiation therapy (EBRT) and/or brachytherapy (BT). Summation of the cumulative dose is critical to assess the total biologic effective dose to targets and organs at risk. Cumulative dose calculation from EBRT and BT can be performed with or without image registration (IR) and biologic dose summation. Among these dose summation strategies, linear addition of dose-volume histogram (DVH) parameters without IR is the global standard for composite dose reporting. This approach stems from an era without image guidance and simple external beam and brachytherapy treatment approaches. With technological advances, EBRT and high-dose-rate BT have evolved to allow for volume-based treatment planning and delivery. Modern conformal therapeutic radiation involves volumetric or intensity modulated EBRT, capable of simultaneously treating multiple targets at different specified dose levels. Therefore, given the complexity of modern radiation treatment, the linear addition of DVH parameters from EBRT and high-dose-rate BT is challenging to represent the combined dose distribution. Deformable image registration (DIR) between EBRT and image guided brachytherapy (IGBT) data sets may provide a more nuanced calculation of multimodal dose accumulation. However, DIR is still nascent in this regard, and needs further development for accuracy and efficiency for clinical use. Biologic dose summation can combine physical dose maps from EBRT and each IGBT fraction, thereby generating a composite DVH from the biologic effective dose. However, accurate radiobiologic parameters are tissue-dependent and not well characterized. A combination of voxel-based DIR and biologic weighted dose maps may be the best approximation of dose accumulation but remains invalidated. The purpose of this report is to review dose summation strategies for EBRT and BT, including conventional equivalent dose in 2-Gy fractions dose summation without image registration, physical dose summation using 3-dimensional rigid IR and DIR, and biologic dose summation. We also provide general clinical workflows for IGBT with a focus on cervical cancer.
The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.
Conventional multi-cylinder drying of paper and board typically relies on both conductive drying from steam-heated cylinders and convective drying, where heated air flows over the paper web surface. Conduction primarily contributes to heat transfer, while convection is the main driver of mass transfer. However, conventional drying systems are heavily dependent on steam, usually powered by fossil fuels, and are often energy-inefficient with high levels of waste. Additionally, these surface-driven processes result in a lower percentage of energy absorption compared to the energy supplied, leading to significant energy losses. To improve this long-standing process, an experimental system was developed to investigate a process intensification approach involving the integration of Radio Frequency (RF) heating, a volumetric electromagnetic technology, alongside traditional conduction and convection drying methods. This study also emphasizes the use of sensors to continuously monitor key parameters such as moisture content, supply system temperatures, sample temperatures, air flows, and drying rates in real-time. The effect of RF as an auxiliary energy source in localized environments at varying moisture levels was explored to optimize industrial drying systems, quantify potential improvements, and provide insights for future studies. Experimental results from trials combining RF with convection and with the base case alternating conduction-convection drying processes are presented. It is shown that RF is a viable process intensification approach for paper drying improving the drying rate and energy intensity at higher moisture contents. These findings offer valuable insights for process intensification and contribute to the process development, modeling, and simulation of advanced paper drying techniques.
The thermal evolution of lattice parameters coupled with heat capacity data provide insight into tailorable magnetism-structure attributes in the orthorhombic compound AlFe 2 B 2 that was synthesized with and without small additions of gallium. Temperature-dependent X-ray powder diffraction experiments conducted through the magnetic phase transition reveal that the a- and b-parameters of both samples increase with increasing temperature while the c-parameter decreases. While a weak volumetric thermal expansion is noted over a range of temperatures well below and above the magnetic phase transition, anomalous behavior was observed within the phase transition region itself to reveal a magnetostructural phase transition with borderline first-order character in the Ga-modified sample but of more second-order character in the Ga-free sample. It is established that the nearest-neighbor Fe-Fe interatomic distance within the (ab)-plane plays a dominant role in influencing the magneto-functional response of these compounds. The magnetocaloric properties are discussed in the context of temperature-induced changes of the interatomic bonding that are influenced by the hypothesized presence of iron antisite defects in the AlFe 2 B 2 lattice.
Intermetallic metal hydrides are critical materials for hydrogen storage applications, however, metal hydrides with greater storage capacities are still needed. Within metal hydrides, the volumetric storage capacities are limited by the number of hydrogen-accommodating interstitial sites which can be simultaneously occupied given a minimum hydride nearest-neighbor distance of ~2.1 A, according to the Switendick-Westlake criterion. To date, violations of this criterion are rare. Perhaps the most well studied compounds violating this criterion are the RNiInHx compounds (R = Ce, La, Nd). Previous neutron diffraction studies on the deuterated species revealed the presence of Ni-D∙∙∙D-Ni-D∙∙∙D-Ni chains with anomalously close D∙∙∙D contacts of ~1.6 A. Yet there are no neutron vibrational spectroscopic investigations reported for these atypical hydrides. Here we use neutron vibrational spectroscopy (NVS) measurements to probe the hydrogen dynamics in LaNiInHx (x = 0.67, 1.6) and CeNiInH 1.4 . For x > 0.67, the presence of close H...H contacts yields two related features in the vibrational spectrum centered near ~90 meV corresponding to the oscillations of paired H atoms simultaneously occupying neighboring R 3 Ni tetrahedra. Notably, these features are energetically distinct from comparable vibrational motions for "unpaired" H atoms when x = 0.67. To compare, we also present powder neutron diffraction and NVS measurements for the newly characterized, chemically similar Sn compounds CeNiSnH, CeNiSnH 2 , and CeNiSnD 2 . These compounds also contain R 3 Ni tetrahedra, however, the H-occupied tetrahedra are well separated from each other with the closest H∙∙∙H distances exceeding 2.1 A, and the Switendick-Westlake criterion is not violated. Consequently, the spectral signature of the close H∙∙∙H contacts is absent in these hydrides.
Here we present a study on geophysical inverse modelling for subsurface structural properties of an unconventional hydrocarbon site that was monitored previously by Interferometric Synthetic Aperture Radar (InSAR) technology for surface deformation. A static three-dimensional geomodel along with extracted property maps replicates the depth of each underlying stratigraphic unit and structural feature with the density of each geological layer. We examine the hypothesis that integration of elastic properties of each formation layer with InSAR observations in a stratified elastic medium will lead to a viscoelastic geophysical inverse problem that can be solved to estimate fractional volume change at the reservoir level. Moreover, we examine synthetic scenarios in which the elastic properties of the formations are perturbed before determining the resulting impact on the rate of surface deformation. The results show that although the slope of underlying formations, their density and depth can define the extent and pattern of a deformation signal, their properties have a marginal impact on volumetric change compared to the dense network of shallow depth Coal Seam Gas (CSG) mining wells. Besides, it is also demonstrated that the inversion of InSAR deformation maps can resolve the uncertainties associated with low-resolution seismic interpretation as well as filling the data gaps within seismic acquisitions. A significant contribution of this investigation to the geological basin modelling involves a) introducing a remote and non-invasive technology such as InSAR to improve geophysical mapping of subsurface structures such as faults in areas with sparse or no reflective seismic information, and b) applying a multi-layer viscoelastic geophysical source model for an unconventional hydrocarbon reservoir such as CSG.
This paper presents a novel sol-gel method to synthesize large and thick silica aerogel monoliths at near ambient conditions using a commercial aqueous solution of colloidal silica nanoparticles as building blocks. To achieve slabs with high visible transmittance and low thermal conductivity, the method combines the strategies of (i) synthesizing gels on an omniphobic perfluorocarbon liquid substrate, (ii) aging at temperatures above room temperature, and (iii) performing solvent exchange with a low-surface-tension organic solvent prior to ambient drying. The omniphobic liquid substrates were used to prevent cracking and ensure an optically-smooth surface, while nanoparticle building blocks were small (<10 nm) to limit volumetric light scattering. Gels were aged at temperatures between 25 and 80 °C for up to 21 days to make them stronger and stiffer and to reduce shrinkage and cracking during ambient drying. Ambient drying was achieved by first exchanging water in the gel pores for octane, followed by drying in an octane-rich atmosphere to decrease capillary forces. The synthesized nanoparticle-based silica aerogel monoliths had thicknesses up to 5 mm, diameters up to 10 cm, porosities exceeding 80%, and thermal conductivities as low as 0.08 W m –1 K –1 . Notably, the slabs featured visible transmittance exceeding 75% even for slabs as thick as 5 mm. The as-synthesized aerogel monoliths were exposed to TMCS vapor to induce hydrophobic properties resulting in a water contact angle of 140° that prevented water infiltration into the pores and protected the aerogels from water damage. This simple synthesis route conducted at near ambient conditions produces hydrophobic aerogel monoliths with promising optically transparent and thermally insulating properties that can be adhered to glass panes for window insulation and solar-thermal energy conversion applications.
Graphene oxide (GO)-based energy storage faces dual bottlenecks: unsustainable reduction methods and sluggish proton transfer kinetics. Here, we introduce a groundbreaking green strategy using lignin-derived vanillyl alcohol (VA) as a dual-function monomer to simultaneously address these challenges. By thermally annealing GO/VA films at mild temperatures (<100 °C), VA triggers an interface-confined reduction of GO while self-polymerizing into redox-active oligomers (P-VA) that intercalate between graphene layers. This dual role-reducing agent and proton highway enables a 3D conductive network with minimized graphene restacking, abundant redox sites, and rapid H + transport pathways. Density Functional Theory (DFT) reveals how P-VA optimizes proton dynamics, while the resulting rGO-P-VA4-T90 electrode achieves a record volumetric capacitance of 311.1F/cm 3 (777.8F/cm 2 ) and retains 87.8 % capacity after 10,000 cycles. Flexible solid-state supercapacitors deliver 94.2 μWh/cm 2 energy density at 63.8 μW/cm 2 , rivaling state-of-the-art devices. Furthermore, this work redefines sustainable graphene engineering, merging biomass valorization with high-performance energy storage in a scalable, eco-friendly paradigm.
Electro-fuels can be produced from concentrated sources of carbon dioxide and hydrogen using electricity generated from renewable sources; this process enables energy storage at high volumetric energy density. Among the electro-fuels options, FT (Fischer-Tropsch) fuel is attractive for heavy-duty trucks and non-road transportation applications. This study conducts a techno-economic analysis of FT liquid fuel production from H 2 and CO 2 using a detailed performance analysis. Minimum fuel selling price is estimated for a broad range of H 2 and CO 2 prices and for a range of potential CO 2 credits. The analysis indicates that H 2 price has the largest impact on the minimum selling price of FT fuel. FT fuel production with a CO 2 price of $17.3/metric ton requires an H 2 price of $0.8/kg to be cost-competitive with the pre-tax petroleum diesel price of $3.1/gal in 2050 (before the application of any CO 2 credits). When the H 2 price is $2.0/kg from central water electrolysis (2020 target), the minimum selling price of the FT fuel is $5.4–5.9/gal. A sensitivity analysis shows that future system optimization of FT fuel production could focus on improving the H 2 and CO 2 recycle contributions and FT fuel conversion ratio. The analysis results can be combined with various upstream systems for H 2 and CO 2 production.
Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.
This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.
The versatility of hydrogen will provide opportunities in large scale, long duration energy storage, and the decarbonization of industry. Further, hydrogen energy carriers are molecules and materials that store hydrogen at higher volumetric density than possible using gas phase hydrogen.
Agricultural production is a major source of carbon dioxide (CO 2 ) and nitrous oxide (N 2 O) globally. The effects of conservation practices on soil CO 2 and N 2 O emissions remain a high degree of uncertainty. In this study, soil CO 2 and N 2 O emissions under different residue and tillage practices in an irrigated, continuous corn system, were investigated using the Root Zone Water Quality Model (RZWQM2). Combinations of no/high stover removal (NR and HR, respectively) and no-till/conventional tillage (NT and CT, respectively) field experiments were tested over the four crop-years (Apr. 2011–Apr. 2015). The model was calibrated using the NRCT, and validated with other treatments. The simulation results showed that soil volumetric water content (VWC) in the NR treatments (i.e., NRCT and NRNT) was 1.3%–1.9% higher than that in the HR treatments (i.e., HRCT and HRNT) averaged across the four years. A higher amount of CO 2 and N 2 O emissions were simulated in the NRCT across the four years (annual average: 7034 kg C/ha/yr for CO 2 and 3.8 kg N/ha/yr for N 2 O), and lower emissions were in the HRNT (annual average: 6329 kg C/ha/yr and 3.7 kg N/ha/yr for N 2 O). A long-term simulation (2001–2015) suggested that the CO 2 and N 2 O emissions were closely correlated with the stover removal degree (SRD), tillage, VWC, soil temperature (ST), years in management (Y), and fertilizer application. Stover and tillage practices had cumulative effects on CO 2 emissions. The simulated annual CO 2 emissions in 1st year from NRCT, NRNT, and HRCT were 7.8%, 0.0%, and 7.7% higher than that from HRNT, respectively; then the emissions in 15th year were 63.6%, 47.7%, and 29.1% higher, respectively. Meanwhile, there were no cumulative effects on N 2 O emissions. Furthermore, the results also demonstrated that the RZWQM2 is a promising tool for evaluating the long-term effects of CO 2 and N 2 O emissions on different conservation practices.
We report the chemical and isotopic compositions of fallout melt glasses from nuclear tests contain a range of information constraining the physical conditions within the fireball and the mechanisms of fallout formation but historic studies tended to exclude the behavior of stable major and trace elements. Here, we present a large study specifically focused on major and trace element relationships within a population of macroscale fallout samples from a single event. We interpret these data to better constrain how fallout melt glass formation in near surface environments is influenced by that environment and demonstrate how major and trace element abundances can provide useful insights into chemical processes within the fireball. Data confirm that the uranium in the fallout glass population derives from two isotopically distinct endmembers: isotopically enriched uranium (presumably from the weapon), and natural composition uranium that may be a combination of anthropogenic and environmental materials from within the blast zone. The similarity between major and trace element concentrations in fallout and corresponding local soils from the event site confirm the local soils as the most probable source of entrained material into the fireball and the source of carrier material into which the bomb vapor was incorporated. The lack of correlation between major and trace element abundances with size indicates that volatility driven processes, such as condensation from the fireball, do not control the composition of macroscale fallout melt glass. Although the fallout has major and trace element chemical characteristics broadly similar to those of the local, associated soils, some systematic differences are observed between the two populations. Fallout melt glass is depleted in volatile elements such as K, Na, Tl and Pb, consistent with heating to temperatures above ~1000 °C for 3-10 s. This is supported by the results of laser heating experiments performed on rhyolitic soil at temperatures (1600-2200 °C) and timescales (1-120 s) that are broadly relevant to fallout formation conditions. Relative enrichments of metals such as Cu and Co do not correlate with the abundance of uranium, suggesting that fallout also records input of near field anthropogenic materials. Our observations suggest that major chemical features can be related to processing in the fireball and used to inform the thermal-chemical evolution of the system. Ultimately, these data are consistent with a fallout formation mechanism that involves rapid melting of surface materials to form carrier material melts with minor incorporation of bomb vapor and a degree of volumetric volatile loss due to heating.
Directed energy deposition (DED) based remanufacturing leverages the flexibility of additive manufacturing to add value to broken or worn components. DED offers the ability to repair cast iron, a material difficult to repair with traditional welding techniques. Despite this, development of appropriate DED process conditions for bimetallic cast iron structures lags low- and medium-carbon steel repair. Thermal stresses and porosity generated by high-temperature deposition on cast iron often lower mechanical properties and hinder the qualification process. In this report scanning speed, powder mass flow rate, and stepover width are studied in multiple-track structures deposited on gray cast iron. Dilution and residual stresses are found to be highly dependent on the selected process parameters. Samples with a higher volumetric energy input, e.g., slower scanning speeds and higher powder feed rates, showed improved density and lower residual stresses but suffered lower dilution into the substrate. The presented conditions further the development of additive manufacturing technologies for automotive repair.
Although bimetallic structures can offer unique solutions to engineering problems with varying functionality, joining two dissimilar alloys using additive manufacturing involves challenges such as cracking and segregation. In this work, bimetallic structures of Inconel 740H superalloy and P91 steel are processed using wire-arc additive manufacturing (WAAM) for the first time. P91 steel is directly deposited over 740H with a cooling time of 20 minutes in between the two depositions. During microstructure characterization, a large gradient zone with coarse grains spanning for ~2200 μm is observed. Moreover, the microhardness of the gradient zone is the least in comparison with the P91 and 740H deposits. No deleterious phase formation is found; however, intergranular cracks are observed along the gradient zone. In order to deduce the reason for cracking, non-equilibrium solidification simulations are performed for the mixed compositions of P91 and 740H in the gradient zone. Here, the reasons for crack formation in the gradient zone are the sudden change in the volumetric coefficient of thermal expansion due to the formation of MC carbides as well as the development of local strains from the thermal residual stresses. Therefore, altering the deposition sequence and introducing gradient layers with mixed composition between the P91 and 740H deposits will be beneficial in achieving defect-free builds. Through the integration of experiments and simulation, this work gains valuable insights into the processability of dissimilar alloy manufacturing using the WAAM technique.
Friction stir welding (FSW) is a solid-state joining process that minimizes the heat-affected zone (HAZ) compared with fusion-based arc welding, making it well suited for joining martensitic armor steels where hardness and ballistic resistance are critical. This study investigates residual stress formation in three defect-free FSW butt-joint configurations relevant to armored-vehicle fabrication: similar rolled homogeneous armor (RHA–RHA, Case 1), similar high-hardness armor (HHA–HHA, Case 2), and dissimilar HHA–RHA (Case 3) joints produced under temperature-controlled conditions (770 °C). Neutron diffraction was employed to quantify the magnitude and spatial distribution of residual stresses in the longitudinal, transverse, and normal directions and to correlate them with weld microstructure and hardness. Tensile residual stresses were concentrated in the softened HAZ, reaching approximately 300 MPa for Case 2 and 400 MPa for Case 1 (≈50–70 % of the base-metal yield strength; ∼581 MPa for RHA and ∼566 MPa for HHA), while compressive residual stresses dominated the stir zone. The spatial extent of tensile stresses scaled with the width of the softened HAZ, which was largest in the dissimilar HHA–RHA joint and smallest in the HHA–HHA joint. Full-width-at-half-maximum (FWHM) analysis revealed low microstrain in overtempered HAZ regions and high microstrain in the stir zone associated with severe plastic deformation and fresh martensite formation. This work demonstrates that residual stress evolution in FSW of martensitic armor steels is governed not primarily by peak temperature or thermal contraction, as inferred from fusion-welding analogies, but by the competition between transformation-induced volumetric expansion and tempering-induced stress relaxation. The relative dominance of these mechanisms is shown to depend on alloy hardenability and local thermal history, leading to more extensive HAZ softening and broader tensile stress regions in the lower-hardenability RHA steel. These findings establish a transferable mechanistic framework for optimizing solid-state joining strategies in high-strength steels and other transformation-hardening alloys beyond armor applications.