Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “residual minimization”

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 55 records · Page 3

Convergence acceleration of iterative solutions of Euler equations for transonic flow computations

In this paper, two acceleration techniques for Euler calculations are investigated. The first technique is an extrapolation procedure based on the Power Method; it is applicable when the iterative matrix has dominant eigenvalues. Both real and complex conjugate roots are allowed. The second technique is a generalization of the Minimal Residual Method, where the extrapolation step consists of a weighted combination of the corrections at different iteration levels and the weights are chosen to minimize the Lz norm of the residual. Numerical results, using Jameson's Runge-Kutta Multigrid Code, are presented. The extra computational work to apply either technique is negligible and the extra storage is not a problem on current supercomputers.

Hafez, M.↗

An extension of command shaping methods for controlling residual vibration using frequency sampling

The authors present an extension to the impulse shaping technique for commanding machines to move with reduced residual vibration. The extension, called frequency sampling, is a method for generating constraints that are used to obtain shaping sequences which minimize residual vibration in systems such as robots whose resonant frequencies change during motion. The authors present a review of impulse shaping methods, a development of the proposed extension, and a comparison of results of tests conducted on a simple model of the space shuttle robot arm. Frequency shaping provides a method for minimizing the impulse sequence duration required to give the desired insensitivity.

Singer, Neil C.↗

MTA Removable Beam Absorber

Once commissioned, the secondary production beamline in the MeV Test Area (MTA) will require a removable beam absorber for continuous beam operation. A 4"x 4" x 4" copper core encased within steel will absorb the full 15 Hz 400 MeV Fermilab Linac primary proton beam. However the residual activation will constrain work activities in the hall, so a 6” marble casing is proposed entombing the core. Marble has minimal residual activation. This presentation describes the initial concept and planned ANSYS analysis.

Blas, Efren↗

Damage Tolerance Assessment of Laser Clad Repairs of Coarse Grain Ti-6Al-4V

Background: Laser cladding offers a promising approach to repairing aviation components, with the ability to restore the geometric shape and the static and fatigue strength. Objective: The damage tolerance of a coarse grain Ti-6Al-4 V material is assessed via three sets of samples: (i) baseline, (ii) laser clad repaired, and (iii) laser cladding followed by heat treatment. Method: The specimens were subjected to fatigue loading using a marker band technique, in which cracks initiated from pre-defined laser slots and the small crack growth rates were measured by quantitative fractography. Additional characterization was performed on these samples, including hardness measurements, microstructure analysis, and multimodal residual stress analyses. Results: The region of the specimen undergoing laser cladding exhibited higher hardness values and a compressive residual stress, based on the slitting method. The samples subjected to heat treatment after cladding displayed negligible differences in grain size and minimal residual stress relaxation. Conclusions: The results demonstrated significant variability, with data demonstrating that cracks in the laser clad region (exhibiting the largest compressive residual stress) grew faster than cracks in the substrate region of the sample or in the baseline cases. Finally, these counterintuitive results are rationalized with additional analyses from synchrotron-based energy dispersive X-ray diffraction and high energy X-ray diffraction, demonstrating significant spatial variability in the residual elastic strains in the millimeter-sized prior β grains within the material.

36 MATERIALS SCIENCE↗

Stress In A Fiber During Curing Of Surrounding Matrix Resin

Experiments run to determine variation in tensile stress in single fiber during curing of matrix resin surrounding fiber. Study part of effort to understand physical mechanisms affecting residual stresses in matrix/fiber composites, with view toward optimizing curing cycles (in particular, optimizing temperature-vs.-time schedules of final cooldowns to ambient temperature) to minimize residual stresses. Results signify primary mechanisms affecting residual stress in fibers are thermal expansion and contraction and cure shrinkage of matrix material.

Bowles, Kenneth J.↗

The alignment-distribution graph

Implementing a data-parallel language such as Fortran 90 on a distributed-memory parallel computer requires distributing aggregate data objects (such as arrays) among the memory modules attached to the processors. The mapping of objects to the machine determines the amount of residual communication needed to bring operands of parallel operations into alignment with each other. We present a program representation called the alignment distribution graph that makes these communication requirements explicit. We describe the details of the representation, show how to model communication cost in this framework, and outline several algorithms for determining object mappings that approximately minimize residual communication.

Chatterjee, Siddhartha↗

The alignment-distribution graph

Implementing a data-parallel language such as Fortran 90 on a distributed-memory parallel computer requires distributing aggregate data objects (such as arrays) among the memory modules attached to the processors. The mapping of objects to the machine determines the amount of residual communication needed to bring operands of parallel operations into alignment with each other. We present a program representation called the alignment-distribution graph that makes these communication requirements explicit. We describe the details of the representation, show how to model communication cost in this framework, and outline several algorithms for determining object mappings that approximately minimize residual communication.

Chatterjee, Siddhartha↗

What is the minimal vestibular function required for compensation?

Living with an uncompensated, abnormal vestibular system requires oppressive modification of life style and often prevents return to work and activities of daily living. Patients with vestibular abnormalities were studied to determine the minimal residual vestibular function required to achieve compensation. Three groups of patients with (a) complete unilateral loss of vestibular function with normal horizontal canal-vestibulo-ocular (HCVOR) function in the opposite ear, (b) complete unilateral loss with abnormal HCVOR function in the opposite ear, and (c) bilateral reduction of vestibular function from aminoglycoside toxicity underwent vestibuloocular (VOR), optokinetic (OKN), visual-VOR (VVOR), and computerized dynamic posturography (CDP) tests before and after therapeutic procedures. Results suggest that a minimal VOR response amplitude must be present for compensation of VVOR function to occur. The roles of VOR and OKN phase shifts in vestibular compensation are more complicated and require further study. Compensation of vestibulospinal function does not necessarily accompany VOR or VVOR compensation. Ascending and descending vestibular compensatory mechanisms may involve different spatial sensory inputs. Results of these studies have important implications for the diagnosis and treatment of patients with vestibular disorders, including selection and monitoring of patients for therapeutic regimens such as vestibular nerve section and streptomycin therapy.

Non-NASA Center↗

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗

Hybrid eigensolvers for nuclear configuration interaction calculations

We examine and compare several iterative methods for solving large-scale eigenvalue problems arising from nuclear structure calculations. In particular, we discuss the possibility of using block Lanczos method, a Chebyshev filtering based subspace iterations and the residual minimization method accelerated by direct inversion of iterative subspace (RMM-DIIS) and describe how these algorithms compare with the standard Lanczos algorithm and the locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm. Although the RMM-DIIS method does not exhibit rapid convergence when the initial approximations to the desired eigenvectors are not sufficiently accurate, it can be effectively combined with either the block Lanczos or the LOBPCG method to yield a hybrid eigensolver that has several desirable properties. We will describe a few practical issues that need to be addressed to make the hybrid solver efficient and robust.

97 MATHEMATICS AND COMPUTING↗

Influence of silicon infiltration conditions on microstructure and mechanical properties in binder jet 3D printed SiC

Atmosphere, temperature, and hold time were varied for silicon infiltrated binder jet 3D printed SiC preforms. A two-step infiltration process optimized infiltration. The first step at 1450 °C under vacuum facilitated improved wetting. The second step at elevated temperature in argon advanced infiltration while suppressing vaporization. Variation of the second-step temperature (1550–1750 °C) and hold time revealed that 1550 °C for 1 h in argon following initial vacuum treatment yielded the highest silicon uptake ratio (∼96 %) with minimal residual porosity and strengths of 150–180 MPa. Flexural strength was either stable or improved with elevated temperatures up to 1000°C. In contrast, higher temperatures or extended hold times led to increased silicon vaporization and pores. The optimized process was applied to formulations with extra carbon to enhance SiC content and interconnectedness through reaction bonding, which corresponded to increased flexural strength. These findings provide guidance for silicon infiltration in additively manufactured SiC.

Yoon, Bola [ORNL] (ORCID:0000000260875373)↗

Additive manufacturing of AISI M2 tool steel by binder jetting (BJ): Investigation of microstructural and mechanical properties

The presented research demonstrates for the first time the successful processing of AISI M2 tool steel by binder jetting, a promising additive manufacturing technique capable of producing complex shapes with minimal residual stresses and isotropic properties. The optimal printing parameters were explored by varying processing parameters such as the binder saturation (45 %–105 %), binder set time (0 to 10 s), targeted bed temperature (50–60 °C), oscillator (2600–2750 rpm), recoater (20–28 mm/s), and roller speeds (200–300 rpm). Microstructural characterization and evaluation of mechanical properties of binder jetted parts were performed using x-ray diffraction (XRD), scanning electron microscopy (SEM), and energy dispersive spectroscopy (EDS) to study their chemical composition, powder morphology, microstructure, carbide morphologies, relative density, hardness, compressive strength, and ductility. Two powder sizes (5 and 10 μm) were used, and sintering was performed at varying temperatures (1270, 1280, and 1300 °C) and durations (60 and 120 min), followed by a furnace, air, and water cooling. An optimum hardness of ~970 HV was obtained when parts were sintered at 1270 °C for 60 min, followed by water quenching. Impressive compressive strength of ~ 3580 MPa was observed in the sample sintered at 1280 °C for 60 min duration, followed by air cooling. Furnace-cooled parts showed the highest density of ~95 %, whereas the relative density of air- and water-cooled parts varied between ~91 to 93.50 %, respectively. The microstructure of sintered samples revealed the formation of M 6 C stable carbide, M 2 C metastable carbide, MC as a secondary carbide, and a-Fe matrix, which contributed to the observed increase in mechanical properties.

36 MATERIALS SCIENCE↗

Low Peripheral Blood Counts and Elevated Proinflammatory Cytokines Signal a Poor CD19 Chimeric Antigen Receptor T-cell Response in Acute Lymphoblastic Leukemia

CD19 chimeric antigen receptor T-cell (CAR-T) therapy has significantly improved outcomes for patients with relapsed/refractory B-cell acute lymphoblastic leukemia (R/R B-ALL). However, approximately 20% of patients fail to achieve a complete remission (CR), and some develop severe, life-threatening toxicities. Understanding the biological mechanisms underlying both dysfunctional responses and severe toxicity is essential for optimizing patient management and improving therapeutic efficacy. This study aimed to (1) characterize cytokine profiles associated with dysfunctional responses and severe toxicity following CAR-T infusion, (2) examine the timing and trajectory of cytokine changes in relation to treatment outcomes, and evaluate potential strategies for mitigating toxicity and treatment failure. We conducted a comprehensive analysis of serum cytokine profiles in 86 adult and pediatric patients undergoing autologous CD19 CAR-T therapy for B-ALL. Patients were categorized into three groups: (1) Dysfunctional response—Patients who failed to achieve a minimal residual disease-negative CR (MRD-CR) by Day 63 or who experienced recurrence of CD19+ disease in the setting ongoing CAR-T cell detection before Day 63. (2) Functional response with severe cytokine release syndrome (CRS) and/or neurotoxicity (NTX)—Patients with best response of MRD-CR by Day 63 who experienced grade 3 or higher CRS or NTX. (3) Functional response without severe CRS or NTX—Patients with best response of MRD-CR by Day 63 who did not experience grade =3 CRS or NTX. Cytokine levels were measured during the first-week postinfusion and correlated with treatment efficacy, toxicity outcomes, complete blood counts, and CAR-T expansion dynamics. This analysis aimed to better understand how cytokine profiles relate to patient outcomes and immune responses in CAR-T therapy. Patients with dysfunctional response exhibited decreased neutrophils, platelets, and levels of granulocytic cytokines (suggestive of low bone marrow reserve) alongside elevated pro-inflammatory cytokines by Day 1. Functional response with severe toxicity patients showed a progressive rise in proinflammatory cytokines, reaching similar levels to dysfunctional response patients by Day 7. We observed that high cytokines at both the Day 1 and Day 7 time points were associated with poor survival. These findings remained significant when adjusting for high disease burden, a known predictor of severe inflammatory toxicity and lack of response. Early post-CAR-T infusion inflammation is associated with both dysfunctional response and severe toxicity—even after adjusting for disease burden. This suggests that inflammation, in addition to disease burden, plays a role in determining patient outcome. Therefore, strategies aimed at reducing the pro-inflammatory state prior to or early after CAR-T cell infusion may improve outcomes for R/R B-ALL patients.

Serum cytokines↗

Experimental and simulated heating in nanocrystalline diamond charge exchange injection films using an electron beam to mimic SNS conditions

Nanocrystalline diamond (NCD) foils have been in use at the Spallation Neutron Source as primary stripper foils for charge exchange injection since the machine’s inception. NCD’s low atomic number leads to reduced beam scattering while still being easy to handle, as opposed to even lighter elements. NCD foils also have the benefit of being free-standing, and rigid when grown under the right conditions to minimize residual stresses. This study overviews a method to mimic the Spallation Neutron Source (SNS) beam and characterize thermally driven failure of NCD foils. To do this, a foil test stand was developed with in situ diagnostics that tracks signs of foil sublimation and thinning. The foil test stand’s electron beam is equated to the SNS beamline by comparing peak deposited energy densities for both beams. Post-mortem testing is also conducted to help elucidate changes the NCD foils undergo during exposure to the beam. Further, a COMSOL simulation was also developed as a method to assess potential future changes to the NCD film and predict how changes to beam conditions will effect foil temperatures. This method allows for examination of the NCD foil’s thermal limits and can be used to assess future changes to SNS stripper foils.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Iterative methods in GPU-resident linear solvers for nonlinear constrained optimization

Linear solvers are major computational bottlenecks in a wide range of decision support and optimization computations. The challenges become even more pronounced on heterogeneous hardware, where traditional sparse numerical linear algebra methods are often inefficient. For example, methods for solving ill-conditioned linear systems have relied on conditional branching, which degrades performance on hardware accelerators such as graphical processing units (GPUs). To improve the efficiency of solving ill-conditioned systems, our computational strategy separates computations that are efficient on GPUs from those that need to run on traditional central processing units (CPUs). Our strategy maximizes the reuse of expensive CPU computations. Iterative methods, which thus far have not been broadly used for ill-conditioned linear systems, play an important role in our approach. In particular, we extend ideas from Arioli et al., (2007) to implement iterative refinement using inexact LU factors and flexible generalized minimal residual (FGMRES), with the aim of efficient performance on GPUs. In conclusion, we focus on solutions that are effective within broader application contexts, and discuss how early performance tests could be improved to be more predictive of the performance in a realistic environment.

97 MATHEMATICS AND COMPUTING↗

Atomic-scale Imaging of PGM-free Catalyst Active Sites by 30 keV 4D-STEM

Platinum group metal (PGM)-free catalysts have attracted a large amount of attention due to their potential for enabling low-cost, commercially viable hydrogen fuel cells and electrolyzers [1]. Degradation of cell performance remains a significant challenge, however, currently limiting implementation of devices that utilize PGM-free materials [2]. Controlling degradation in these materials requires a better understanding of the atomic structure of active sites, thought to be FeN4 structures within a graphitic carbon lattice, which would enable more accurate prediction of the associated degradation pathways. The exact structural arrangement of active sites in these materials is still debated since a range of structures have been computationally predicted and methods to directly validate these models are still needed [1,3-5].While scanning transmission electron microscopy (STEM) has provided initial glimpses into the nature of the proposed active sites, detailed evaluation of the structure of these sites is challenging since they involve defects, edges, and nitrogen dopants in the graphitic lattice, which are susceptible to beam damage [6,7]. Conventional STEM imaging and spectroscopy methods exacerbate this problem with dose-inefficiency or nonideal contrast characteristics for imaging light and heavy elements simultaneously. A method that minimizes damage while generating dose-efficient, easily-interpretable contrast for both light and heavy elements is therefore needed to facilitate imaging of sensitive active site atomic structure.Here, we demonstrate direct atomic-scale imaging of PGM-free catalyst active sites by low-voltage four-dimensional (4D)-STEM. To accomplish this, we pair a 30 keV aberration-corrected probe, which minimizes knock-on damage, with a fast pixelated detector that has optimal performance at low beam energies [8]. We show how this setup enables relatively simple center-of-mass (CoM) techniques [9,10] to produce images with an increased signal-to-noise ratio (SNR) and light-element contrast over conventional imaging modes, allowing the entire active site structure to be imaged at the atomic scale. Moreover, we show how electron ptychography [11,12] enables images with further improved characteristics to be obtained, for example by minimizing residual aberrations, which provides a more accurate structural representation. The increased understanding that these low-voltage 4D-STEM techniques will provide about the atomic structure of PGM-free active sites and their associated degradation pathways will facilitate rational design of next-generation materials, promoting development of low-cost hydrogen fuel cells, electrolyzers, and other energy conversion devices [13]. References:[1] U Martinez et al., Adv. Mater. 31 (2019), p. 1806545.[2] Y Shao et al., Adv. Mater. 31 (2019), p. 1807615.[3] J Kneebone et al., J. Phys. Chem. C 121, (2017), p.16283.[4] T Mineva et al., ACS Catal. 9 (2019), p. 9359.[5] J Li et al., Nat. Catal. 4 (2021), p. 10.[6] T Susi et al., ACS Nano 6 (2012), p. 8837.[7] P Zelenay and DJ Myers, US Department of Energy Hydrogen and Fuel Cells Program 2017 Annual Merit Review and Peer Evaluation Meeting, Washington, DC (2017).[8] H Ryll et al., J. Inst. 11 (2016), p. P04006.[9] K Müller et al., Nat. Commun. 5 (2014), p. 5653.[10] I Lazic et al., Ultramicroscopy 160 (2016), p. 265.[11] H Yang et al., Nat. Commun. 7 (2016), p. 12532.[12] Y Jiang et al., Nature 559 (2018), p. 343. [13] Research sponsored by the Hydrogen and Fuel Cell Technologies Office, Office of Energy Efficiency and Renewable Energy, US Department of Energy (DOE). Research was conducted at the Center for Nanophase Materials Sciences, which is a DOE Office of Science User Facility.

Zachman, Michael↗