Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization methods”

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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fusion of 12 C+ 24 Mg far below the barrier: Evidence for the hindrance effect

Background: The phenomenon of fusion hindrance may have important consequences on the nuclear processes occurring in astrophysical scenarios, if it is a general behavior of heavy-ion fusion at extreme subbarrier energies, including reactions involving lighter systems, e.g., reactions in the carbon and oxygen burning stages of heavy stars. The hindrance is generally identified by the observation of a maximum of the S-factor vs energy. Whether there is an S-factor maximum at very low energies for systems with a positive fusion Q value is an experimentally challenging question. Purpose: Our aim has been to search for evidence of fusion hindrance in 12 C+ 24 Mg which is a medium-light m system with positive Q value for fusion, besides the heavier cases where hindrance is recognized to be a general phenomenon. 12 C+ 24 Mg is very close to the 16 O + 16 O and 12 C + 12 C systems that are important for the late evolution of heavy stars. Methods: The experiment has been performed in inverse kinematics using the 24 Mg beam from the XTU Tandem accelerator of LNL in the energy range 26-52 MeV with an intensity of 4-8 pnA. The targets were 12 C evaporations 50 mu g/cm 2 thick, isotopically enriched to 99.9%. The fusion-evaporation residues were detected at small angles by a E-Delta E-ToF detector telescope following an electrostatic beam deflector. Results: Previous measurements of fusion cross section for 12 C+ 24 Mg were limited to above-barrier energies. In the present experiment the excitation function has been extended down to similar or equal to 15 mu b and it appears that the S factor develops a clear maximum vs energy, indicating the presence of hindrance. This is the first convincing evidence of an S factor maximum in a medium-light system with a positive fusion Q value. These results have been fitted following a recently suggested method and a detailed analysis within the coupled-channels model that has been performed using a Woods-Saxon potential and including the ground state rotational band of 24 Mg. The coupled-channels calculations give a good account of the data near and above the barrier but overpredict the cross sections at very low energies. Conclusions: The hindrance phenomenon is clearly observed in 12 C+ 24 Mg, and its energy threshold is in reasonable agreement with the systematics observed for several medium-light systems. The fusion cross sections at the hindrance threshold show that the highest value (sigma(s) = 1.6 mb) is indeed found for this system. Therefore it may even be possible to extend the measurements further down in energy to better establish the position of the S-factor maximum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Teaching Freight Mode Choice Models New Tricks Using Interpretable Machine Learning Methods

Understanding and forecasting the intricate freight mode choice behavior under various industry, policy, and technology contexts is essential in freight planning and policymaking. Numerous models have been developed in prior studies to provide insights into freight mode selection, the majority of which use discrete choice models such as multinomial logit (MNL) models. However, logit models often rely on linear specifications of independent variables, despite potential nonlinear relationships in the data. Moreover, there often lacks a heuristic and efficient approach to identify such complex relationships to define the logit model specifications. To fill this gap, we developed an MNL model for freight mode choice using the insights from state-of-the- art machine learning (ML) models. ML models can capture the nonlinear nature of the complex decision-making process, and recent advances in 'explainable AI' have greatly improved their interpretability. The interpretable ML methods help enhance the performance of MNL models and advance knowledge of freight mode choice. Specifically, the influential factors and their relationship with individual modes are identified using SHapley Additive exPlanations (SHAP) to improve the MNL's performance. The workflow is demonstrated in a case study of Austin, Texas, and the SHAP results reveal multiple nonlinear relationships predicted by ML models. Incorporating those relationships into MNL model specifications improves the interpretability and accuracy of the MNL model compared to a conventional MNL model. Findings from this study can be used to guide freight planning and inform policymakers and practitioners on how key factors affect freight decision-making.

ADVANCED PROPULSION SYSTEMS,MATHEMATICS AND COMPUT↗
DOE OSTI2024

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Extracting the femtometer structure of strange baryons using the vacuum polarization effect

Abstract One of the fundamental goals of particle physics is to gain a microscopic understanding of the strong interaction. Electromagnetic form factors quantify the structure of hadrons in terms of charge and magnetization distributions. While the nucleon structure has been investigated extensively, data on hyperons are still scarce. It has recently been demonstrated that electron-positron annihilations into hyperon-antihyperon pairs provide a powerful tool to investigate their inner structure. We present a method useful for hyperon-antihyperon pairs of different types which exploits the cross section enhancement due to the effect of vacuum polarization at theJ/ψresonance. Using the 10 billionJ/ψevents collected with the BESIII detector, this allows a precise determination of the hyperon structure function. The result is essentially a precise snapshot of the$$\bar{\Lambda }{\Sigma }^{0}\,(\Lambda {\bar{\Sigma }}^{0})$$ Λ ¯ Σ 0 ( Λ Σ ¯ 0 ) transition process, encoded in the transition form factor ratio and phase. Their values are measured to beR = 0.860 ± 0.029(stat.) ± 0.015(syst.),$$\Delta {\Phi }_{\bar{\Lambda }{\Sigma }^{0}}=(1.011\pm 0.094({{\rm{stat.}}})\pm 0.010({{\rm{syst.}}}))\,{{\rm{r}}}ad$$ Δ Φ Λ ¯ Σ 0 = ( 1.011 ± 0.094 ( stat. ) ± 0.010 ( syst. ) ) r a d and$$\Delta {\Phi }_{\Lambda {\bar{\Sigma }}^{0}}=(2.128\pm 0.094({{\rm{stat.}}})\pm 0.010({{\rm{syst.}}}))\,{{\rm{r}}}ad$$ Δ Φ Λ Σ ¯ 0 = ( 2.128 ± 0.094 ( stat. ) ± 0.010 ( syst. ) ) r a d . Furthermore, charge-parity (CP) breaking is investigated in this reaction and found to be consistent with CP symmetry.

Science & Technology - Other Topics↗
DOE OSTI2024

Parton distributions from boosted fields in the Coulomb gauge

We propose a new method to calculate parton distribution functions (PDFs) from lattice correlations of boosted quarks and gluons in the Coulomb gauge. Compared to the widely used gauge-invariant Wilson-line operators, these correlations greatly simplify the renormalization thanks to the absence of linear power divergence. Additionally, they enable access to larger off-axis momenta under preserved 3D rotational symmetry, as well as enhanced long-range precision that facilitates the Fourier transform. We verify the factorization formula that relates this new observable to the quark PDF at one-loop order in perturbation theory. Moreover, through a lattice calculation of the pion valence quark PDF, we demonstrate the aforementioned advantage and features of the Coulomb gauge correlation and show that it yields consistent results with the gauge-invariant method. This opens the door to a more efficient way to calculate parton physics on the lattice.

Baryons↗
DOE OSTI2023

Efficient Multigrid Reduction-in-Time for Method-of-Lines Discretizations of Linear Advection

Parallel-in-time methods for partial differential equations (PDEs) have been the subject of intense development over recent decades, particularly for diffusion-dominated problems. It has been widely reported in the literature, however, that many of these methods perform quite poorly for advection-dominated problems. In this report we analyze the particular iterative parallel-in-time algorithm of multigrid reduction-in-time (MGRIT) for discretizations of constant-wave-speed linear advection problems. We focus on common method-of-lines discretizations that employ upwind finite differences in space and Runge-Kutta methods in time. Using a convergence framework we developed in previous work, we prove for a subclass of these discretizations that, if using the standard approach of rediscretizing the fine-grid problem on the coarse grid, robust MGRIT convergence with respect to CFL number and coarsening factor is not possible. This poor convergence and non-robustness is caused, at least in part, by an inadequate coarse-grid correction for smooth Fourier modes in space-time known as characteristic components. We propose an alternative coarse-grid operator that provides a better correction of these modes. This coarse-grid operator is related to previous work and uses a semi-Lagrangian discretization combined with an implicitly treated truncation error correction. Theory and numerical experiments show the proposed coarse-grid operator yields fast MGRIT convergence for many of the method-of-lines discretizations considered, including for both implicit and explicit discretizations of high order. Parallel results demonstrate speed-up over sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Invariant discovery of features across multiple length scales: Applications in microscopy and autonomous materials characterization

Physical imaging is a foundational characterization method in areas from condensed matter physics and chemistry to astronomy and spans length scales from atomic to universe. Images encapsulate crucial data regarding atomic bonding, materials microstructures, and dynamic phenomena such as microstructural evolution and turbulence, among other phenomena. The challenge lies in effectively extracting and interpreting this information. Variational Autoencoders (VAEs) have emerged as powerful tools for identifying the underlying factors of variation in image data, providing a systematic approach to distilling meaningful patterns from complex data sets. However, a significant hurdle in their application is the definition and selection of appropriate descriptors reflecting local structures. Here, we introduce the scale-invariant VAE approach (SI-VAE) based on the progressive training of the VAE with the descriptors sampled at different length scales. The SI-VAE allows the discovery of the length scale-dependent factors of variation in the system. Here, we illustrate this approach using the ferroelectric domain images and generalize it to the movies of the electron-beam induced phenomena in graphene and topography evolution across combinatorial libraries. This approach can further be used to initialize the decision making in automated experiments including structure–property discovery and can be applied across a broad range of imaging methods. This approach is universal and can be applied to any spatially resolved data including both experimental imaging studies and simulations, and can be particularly useful for exploration of phenomena such as turbulence and scale-invariant transformation fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of thermal conditioning on the initiation threshold of secondary high‐explosives

While most performance metrics of high–explosive (HE) based devices like detonation velocity, detonation pressure, and energy output are expected to degrade over time, the evolution of initiation threshold appears less clear, with claims of both increasing and decreasing trends in threshold having been made in the literature. This work analyzes D–optimally designed sequential binary test data for a few thermally conditioned porous–powder and polymer–bonded HE initiator systems using a Bayesian likelihood method employing the probit regression model. Here we find that in most cases the initiation threshold decreases (i.e., sensitivity increases) upon accelerated thermal conditioning. However, such results are nuanced and influenced by factors like the contact area of initiating stimulus, HE characteristics like density and specific surface area, as well as possible thermally induced changes to other materials and interfaces involved.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2023

Unraveling capacity fading in lithium-ion batteries using advanced cyclic tests: A real-world approach

Battery lifespan estimation is essential for effective battery management systems, aiding users and manufacturers in strategic planning. However, accurately estimating battery capacity is complex, owing to diverse capacity fading phenomena tied to factors such as temperature, charge-discharge rate, and rest period duration. In this work, we present an innovative approach that integrates real-world driving behaviors into cyclic testing. Unlike conventional methods that lack rest periods and involve fixed charge-discharge rates, our approach involves 1000 unique test cycles tailored to specific objectives and applications, capturing the nuanced effects of temperature, charge-discharge rate, and rest duration on capacity fading. This yields comprehensive insights into cell-level battery degradation, unveiling growth patterns of the solid electrolyte interface (SEI) layer and lithium plating, influenced by cyclic test parameters. Here, the results yield critical empirical relations for evaluating capacity fading under specific testing conditions.

25 ENERGY STORAGE↗
DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗
DOE OSTI2023

Deep reinforcement learning assisted co-optimization of Volt-VAR grid service in distribution networks

With the increasing penetration of distributed energy resources in distribution networks, Volt-VAR control and optimization (VVC/VVO) have become very important to ensure an acceptable quality of service to all customers. System operators can rely on slow-responding utility devices, including capacitor banks and on-load tap changing transformers, along with fast-responding battery and photovoltaic (PV) inverters for the VVC/VVO implementation. Because of variations in response time of these two classes of devices, and different control actions (discrete versus continuous), coordinated and optimal scheduling and operation have become of utmost importance. Here, this paper develops a look-ahead deep reinforcement learning (DRL)-based multi-objective VVO technique to improve the voltage profile of active distribution networks, decrease network and inverter power loss, and save the operational cost of the grid. It proposes a deep deterministic policy gradient (DDPG)-based approach to schedule the optimal reactive and/or active power set-points of fast-responding inverters, and a deep Q-network (DQN)-based DRL agent to schedule the discrete decisions variables of slow-responding assets. The reactive power output of PV and battery smart inverters are scheduled at 30-minute intervals and the capacitors’ commitment status is scheduled with several hour intervals. The proposed framework is validated on the modified IEEE 34-bus and 123-bus test cases with embedded PV and PV-plus-storage. To validate the efficacy of the proposed VVO, it is compared with several scenarios, including the base case without VVO, localized droop control of DERs, DDPG-only, and twin delayed DDPG (TD3) agent-based DRL techniques. The results justify the superior performance of the proposed method to improve the voltage profile, reduce network power loss, and minimize the look-ahead grid operational cost while minimizing the undesirable power losses in inverters as a result of power factor adjustments.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Fractal analysis on Ag 2 O thin film using a data-driven approach

The synthesis of fractal Ag oxide (Ag 2 O) on the surface of Ag thin film has been achieved at room temperature by using Synchrotron X-ray irradiation. We have performed an automated quantitative analysis of a batch of 1879 fractal Ag 2 O patterns in a scanning electron microscopy (SEM) image within a radius of 2000 mm outward from the center of the X-ray beam. The morphology is similar to that of the diffusion-limited cluster aggregation (DLCA) model. The fractal dimension (D) of Ag 2 O is between 1.7 and 1.5 from the center to the edge. The area distribution density of fractal Ag 2 O follows a quadratic function with radius R. It is found that the branches’ number of fractal Ag 2 O is a key factor affecting the fractal dimension. The more branches the fractal has, the greater the D is. This is the first time that Ag fractal has been investigated by combining automated data analysis methods with batch experimental data. This data-driven approach provides a new research perspective for rationally regulating materials’ fractal morphology and performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sub-Sampled Imaging for STEM: Maximising Image Speed, Resolution and Precision Through Reconstruction Parameter Refinement

Sub-sampling during image acquisition in scanning transmission electron microscopy (STEM) has been shown to provide a means to increase the overall speed of acquisition while at the same time providing an efficient means to control the dose, dose rate and dose overlap delivered to the sample. In this paper, we discuss specifically the parameters used to reconstruct sub-sampled images and highlight their effect on inpainting using the beta-process factor analysis (BPFA) methodology. The selection of the main control parameters can have a significant effect on the resolution, precision and sensitivity of the final inpainted images, and here we demonstrate a method by which these parameters can be optimised for any image in STEM. As part of this paper, we also provide a link to open source code and a tutorial on its use, whereby these parameters can be tested for any datasets. When coupled with the hardware necessary to rapidly sub-sample images in STEM, this approach can have significant implications for imaging beam sensitive materials and dynamic processes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A Solid-State Support for Separating Astatine-211 from Bismuth

Increasing access to the short-lived α-emitting radionuclide astatine-211 ( 211 At) has the potential to advance targeted α-therapeutic treatment of disease and to solve challenges facing the medical community. For example, there are numerous technical needs associated with advancing the use of 211 At in targeted α-therapy, e.g., improving 211 At chelates, developing more effective 211 At targeting, and characterizing in vivo 211 At behavior. There is an insufficient understanding of astatine chemistry to support these efforts. The chemistry of astatine is one of the least developed of all elements on the periodic table, owing to its limited supply and short half-life. Increasing access to 211 At could help address these issues and advance understanding of 211 At chemistry in general. Here, we contribute an extraction chromatographic processing method that simplifies 211 At production in terms of purification. It utilizes the commercially available Pre-Filter resin to rapidly (<1.5 h) isolate 211 At from irradiated bismuth targets (Bi decontamination factors ≥876 000), in reasonable yield (68–55%) and in a form that is compatible for subsequent in vivo study. We are excited about the potential of this procedure to address 211 At supply and processing/purification problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Synthesis of Narrow SnTe Nanowires Using Alloy Nanoparticles

Topological crystalline insulator tin telluride (SnTe) provides a rich playground to examine interactions of correlated electronic states, such as ferroelectricity, topological surface states, and superconductivity. The study of SnTe nanowires may lead to even richer physics owing to the one-dimensional (1D) confinement effect and an increased contribution from the topological surface states. Thus, for transport measurements, SnTe nanowires must be synthesized with reduced diameters and high crystalline quality to ensure 1D confinement and phase coherence of the topological surface electrons. Here, we report a facile growth method to produce narrow SnTe nanowires with a high yield using alloy nanoparticles as growth catalysts. The average diameter of the SnTe nanowires grown using alloy nanoparticles is 85 nm, nearly a factor of three reduction compared to the average diameter of 240 nm when using gold nanoparticles as growth catalysts. Transport measurements reveal the effect of the nanowire diameter on the residual resistance ratio and magnetoresistance. Particularly, the ferroelectric transition temperature for SnTe evolves systematically with the nanowire diameter. In situ cryogenic cooling of narrow SnTe nanowires in a transmission electron microscope directly reveals the cubic to rhombohedral structural transition, which is associated with the ferroelectric transition. Thus, these narrow SnTe nanowires represent a model system to study electronic states arising from 1D confinement, such as 1D topological superconductivity and potential multiband superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Wavelength modulation laser-induced fluorescence for plasma characterization

Laser-Induced Fluorescence (LIF) spectroscopy is an essential tool for probing ion and atom velocity distribution functions (VDFs) in complex plasmas. VDFs carry information about the kinetic properties of species that is critical for plasma characterization. Accurate interpretation of these functions is challenging due to factors such as multicomponent distributions, broadening effects, and background emissions. Our research investigates the use of Wavelength Modulation (WM) LIF to enhance the sensitivity of VDF measurements. Unlike standard Amplitude Modulation (AM) methods, WM–LIF measures the derivative of the LIF signal. This approach makes variations in VDF shape more pronounced. VDF measurements with WM–LIF were investigated with both numerical modeling and experimental measurements. The developed model enables the generation of both WM and AM signals, facilitating comparative analysis of fitting outcomes. Experiments were conducted in a weakly collisional argon plasma with magnetized electrons and non-magnetized ions. Measurements of the argon ion VDFs employed a narrow-band tunable diode laser, which scanned the 4p 4 D 7/2 –3d 4 F 9/2 transition centered at 664.553 nm in vacuum. A lock-in amplifier detected the second harmonic WM signal, which was generated by modulating the laser wavelength with an externally controlled piezo-driven mirror of the diode laser. Finally, our findings indicate that the WM–LIF signal is more sensitive to fitting parameters, allowing for better identification of VDF parameters such as the number of distribution components, their temperatures, and velocities. In addition, WM–LIF can serve as an independent method to verify AM measurements and is particularly beneficial in environments with substantial light noise or background emissions, such as those involving thermionic cathodes and reflective surfaces.

47 OTHER INSTRUMENTATION↗