Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational 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 523 records · Page 29

A High Energy X-ray Diffraction Study of Amorphous Indomethacin

Amorphous pharmaceuticals often possess a wide range of molecular conformations and bonding arrangements. The x-ray pair distribution function (PDF) method is a powerful technique for the characterization of variations in both intra-molecular and inter-molecular packing arrangements. Here, the x-ray PDF of amorphous Indomethacin is shown to be particularly sensitive to the preferred orientations of the chlorobenzyl ring found in isomers in the crystalline state. In some cases, the chlorobenzyl ring has no preferred torsional angle in the amorphous form, while in others evidence of distinct isomer orientations are observed. Amorphous samples with no preferred torsion angles of the chlorobenzyl ring are found to favor enhanced inter-molecular hydrogen bonding, and this is reflected in the intensity of the first sharp diffraction peak. These significant variations in structure rule out amorphous Indomethacin as a possible standard for x-ray PDF measurements. At high humidity, time resolved PDF's for >40 h reveal water molecules forming hydrogen bonds with Indomethacin molecules. Finally, a simple linear hydrogen bond model indicates that water molecules in the wet amorphous form have similar hydrogen bond strengths to those found between Indomethacin dimers or chains in the dry amorphous form.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detection and imaging of chemicals and hidden explosives using terahertz time-domain spectroscopy and deep learning

Detecting concealed chemicals and explosives remains a critical challenge in global security. Terahertz time-domain spectroscopy (THz-TDS) offers a promising non-invasive and stand-off detection technique owing to its ability to penetrate optically opaque materials without causing ionization damage. While many chemicals exhibit distinct spectral features in the terahertz range, conventional terahertz-based detection methods often struggle in real-world environments, where variations in sample geometry, thickness, and packaging can lead to inconsistent spectral responses. In this study, we present a chemical imaging system that integrates THz-TDS with deep learning to enable accurate pixel-level identification and classification of different explosives. Operating in reflection mode and enhanced with plasmonic nanoantenna arrays, our THz-TDS system achieves a peak dynamic range of 96 dB and a detection bandwidth of 4.5 THz, supporting practical, stand-off operation. By analyzing individual time-domain pulses with deep neural networks, the system exhibits strong resilience to environmental variations and sample inconsistencies. Blind testing across eight chemicals—including pharmaceutical excipients and explosive compounds—resulted in an average classification accuracy of 99.42% at the pixel level. Notably, the system maintained an average accuracy of 88.83% when detecting explosives concealed under opaque paper coverings, demonstrating its robust generalization capability. These results highlight the potential of combining advanced terahertz spectroscopy with neural networks for highly sensitive and specific chemical and explosive detection in diverse and operationally relevant scenarios.

Imaging and sensing↗

Skillful multiyear predictions of ocean acidification in the California Current System

The California Current System (CCS) sustains economically valuable fisheries and is particularly vulnerable to ocean acidification, due to its natural upwelling of carbon-enriched waters that generate corrosive conditions for local ecosystems. Here we use a novel suite of retrospective, initialized ensemble forecasts with an Earth system model (ESM) to predict the evolution of surface pH anomalies in the CCS. We show that the forecast system skillfully predicts observed surface pH variations a year in advance over a naive forecasting method, with the potential for skillful prediction up to five years in advance. Skillful predictions of surface pH are mainly derived from the initialization of dissolved inorganic carbon anomalies that are subsequently transported into the CCS. Our results demonstrate the potential for ESMs to provide skillful predictions of ocean acidification on large scales in the CCS. Initialized ESMs could also provide boundary conditions to improve high-resolution regional forecasting systems.

54 ENVIRONMENTAL SCIENCES↗

Forward sensitivity approach for estimating eddy viscosity closures in nonlinear model reduction

In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.

42 ENGINEERING↗

Pairing Directional Solar Inputs From Ray Tracing to Solar Receiver/Reactor Heat Transfer Models on Unstructured Meshes: Development and Case Studies

A novel method for pairing surface irradiation and volumetric absorption from Monte Carlo ray tracing to computational heat transfer models is presented. The method is well-suited to directionally and spatially complex concentrated radiative inputs (e.g., solar receivers and reactors). The method employs a generalized algorithm for directly mapping absorbed rays from a Monte Carlo ray tracing model to boundary or volumetric source terms in the computational mesh. The algorithm is compatible with unstructured, two and three-dimensional meshes with varying element shapes. Four case studies were performed on a directly irradiated, windowed solar thermochemical reactor model to validate the method. The method was shown to conserve energy and preserve spatial variation when mapping rays from a Monte Carlo ray tracing model to a computational heat transfer model in ansys fluent.

14 SOLAR ENERGY↗

Development of systematic uncertainty-aware neural network trainings for binned-likelihood analyses at the LHC

We propose a neural network training method capable of accounting for the effects of systematic variations of the data model in the training process and describe its extension towards neural network multiclass classification. The procedure is evaluated on the realistic case of the measurement of Higgs boson production via gluon fusion and vector boson fusion in the τ τ decay channel at the CMS experiment. The neural network output functions are used to infer the signal strengths for inclusive production of Higgs bosons as well as for their production via gluon fusion and vector boson fusion. We observe improvements of 12 and 16% in the uncertainty in the signal strengths for gluon and vector-boson fusion, respectively, compared with a conventional neural network training based on cross-entropy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Variance reduction via simultaneous importance sampling and control variates techniques using vegas

Monte Carlo (MC) integration is an important calculational technique in the physical sciences. Practical considerations require that the calculations are performed as accurately as possible for a given set of computational resources. To improve the accuracy of MC integration, a number of useful variance reduction algorithms have been developed, including importance sampling and control variates. In this work, we demonstrate how these two methods can be applied simultaneously, thus combining their benefits. We provide a python wrapper, named CoVVVR, which implements our approach in the VEGAS program. The improvements are quantified with several benchmark examples from the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Codebase release r1.4 for CoVVVR

Monte Carlo (MC) integration is an important calculational technique in the physical sciences. Practical considerations require that the calculations are performed as accurately as possible for a given set of computational resources. To improve the accuracy of MC integration, a number of useful variance reduction algorithms have been developed, including importance sampling and control variates. In this work, we demonstrate how these two methods can be applied simultaneously, thus combining their benefits. We provide a python wrapper, named COVVVR, which implements our approach in the VEGAS program. The improvements are quantified with several benchmark examples from the literature.

Shyamsundar, Prasanth↗

Ultrasonic Characterization of Ethylene Vinyl Acetate (EVA) Crosslinking for Quality Assurance and Lamination Process Control (US-Xlink)

Module makers strive to cut lamination processing times to enhance production throughput and reduce costs. If overdone, this might lead to poor EVA quality due to incomplete EVA crosslinking and a large concentration of potentially harmful reactants. Those deficiencies are frequently missed during production quality testing as initially their impact on power output is small. Additionally, the crosslinking agent is frequently inhomogeneously dispersed across the EVA foils such that local destructive characterization procedures such as differential scanning calorimetry (DSC), Soxhlet extraction, swelling methods, or mechanical tests provide only a limited amount of information. However, these defects and inhomogeneities can become considerably more important during field operation, dramatically reducing long-term power yield, and increasing LCOE. Examples of typical long-term module degradation modes relating to poor lamination process conditions include cell breakage, corrosion of the metallization, delamination, and local quality deficiencies. To provide reliable material characterization in a manufacturing setting, we devised a non-destructive technology that uses ultrasound to evaluate the quality of interface adhesion and the degree of crosslinking. As a calibration reference, DSC measurements were employed. Although the potential of ultrasonic approaches for this purpose has already been noted, those previous methods were generally limited to local qualitative measurements. Furthermore, variations in EVA thickness and temperature had a substantial impact on them. The solution we propose solves these limitations by employing ultrasonic absorption rather than relying solely on sound velocity data.

14 SOLAR ENERGY↗

Impact of Humidity and Temperature on the Stability of the Optical Properties and Structure of MAPbI3, MA0.7FA0.3PbI3 and (FAPbI3)0.95(MAPbBr3)0.05 Perovskite Thin Films

In situ real-time spectroscopic ellipsometry (RTSE) measurements have been conducted on MAPbI3, MA0.7FA0.3PbI3, and (FAPbI3)0.95(MAPbBr3)0.05 perovskite thin films when exposed to different levels of relative humidity at given temperatures over time. Analysis of RTSE measurements track changes in the complex dielectric function spectra and structure, which indicate variations in stability influenced by the underlying material, preparation method, and perovskite composition. MAPbI3 and MA0.7FA0.3PbI3 films deposited on commercial fluorine-doped tin oxide coated glass are more stable than corresponding films deposited on soda lime glass directly. (FAPbI3)0.95(MAPbBr3)0.05 films on soda lime glass showed improved stability over the other compositions regardless of the substrate, and this is attributed to the preparation method as well as the final composition.

36 MATERIALS SCIENCE↗

Radial Profiles of the Binary Fraction in Elliptical Galaxies

The radial profile of the binary fraction may vary with environment and is of significant importance for studying the formation mechanisms of binary stars and their dynamical evolution within globular clusters and galaxies. However, existing studies remain limited to the Milky Way and its neighboring galaxies. Leveraging the method proposed by Zhang et al. for estimating the variation in the binary fraction from integrated spectral features, we analyze a sample of 513 elliptical galaxies drawn from the Mapping Nearby Galaxies at Apache Point Observatory survey to measure their radial binary fraction profiles. Our results show that after accounting for the effect induced by radial variations in the stellar population (SP), the median SP-subtracted binary fraction, $r^{\textrm{med}}_{\textrm{b,sub}}$, becomes approximately flat. For nearly all elliptical galaxies in our sample, the variation in the binary fraction relative to the galaxy center at 1R e is less than 5%. No clear correlation is found between the ​​​​​gradient of the binary fraction and the gradients of the SP properties. Moreover, we also compare differences between ultraviolet (UV) upturn and non-UV upturn galaxies. The overall binary fraction profiles and the SP properties of the non-UV upturn galaxies in our sample are comparable to those of the UV upturn galaxies. This similarity may arise from the presence of residual star formation in the non-UV upturn systems.

Li, Xiejin [Chinese Academy of Sciences; Universit↗

In-line, High-Throughput Quality Monitoring for Fuel Cell and Electrolyzer Components Based on Transmission and Reflection Imaging

During the manufacturing of fuel cell and electrolyzer membranes and membrane electrode assemblies (MEAs), real-time, in-line, high-throughput optical-based quality monitoring methods are essential for detecting defects and monitoring thickness variations, thus improving the performance and increasing the durability of fuel cell and electrolyzer in the hydrogen industry. For the MEAs with very opaque coatings, optical transmission-based imaging has been developed and applied in the Roll-to-Roll system using a flashlight and a high-sensitivity CCD camera. We observed high signal-to-noise ratio images while the Roll-to-Roll system ran at 5 ft/min. The entire sample image could quickly be recovered from the discrete frames using customized Python codes for automatic frame cropping and stitching. We detected significant non-uniformities in our experimental MEAs specimen. For fuel cell and low-temperature electrolysis (LTE) transparent membranes, we used optical reflectance hyperspectral imaging with interference fringe-based thickness mapping. We set up a hyperspectral camera to measure various rolls of commercial membranes. The measurement results are analyzed to find the thickness distribution of each roll and to check for defects. Transmission and reflection imaging-based quality monitoring techniques demonstrated in this project can be widely used in the mass production environment to improve the production yield and performance of hydrogen devices.

DIRECT ENERGY CONVERSION,ENGINEERING↗

Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation

This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: (i) selected configuration interaction (sCI) realized both via an a posteriori (ap-sCI) and variational selection of the most important determinants, (ii) an a posteriori singular value decomposition (SVD) of the FCI tensor (SVD-FCI), and (iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that the DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods (except DMRG) perform best with a delocalized basis, and the efficiency of sCI and SVD-FCI is closer to that of DMRG. For larger analogs of the models, the DMRG consistently requires the fewest parameters but still scales exponentially in 2D and 3D systems, and the performance of SVD-FCI is essentially equivalent to that of ap-sCI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗

Tempering kinetics during multilayer laser additive manufacturing of a ferritic steel

Grade 91 steel forms martensite during additive manufacturing and the extent of tempering of martensite significantly affects the mechanical properties of parts. Currently, there is a lack of quantitative understanding of the tempering kinetics for Grade 91 steel, and as a result, the effects of repeated thermal cycles on properties for different processing conditions cannot be determined. Here we evaluate the tempering kinetics by determining the constant terms in the Johnson Mehl Avrami kinetic equation from the tempering data available in the literature and the thermal cycles computed using a rigorously-tested heat and fluid flow model of multi-layer additive manufacturing. The raw tempering data are cleaned using a neural network to enhance accuracy. The lower layers experience repeating cycles of heating and cooling when the upper layers are added. As a result, the hardness is reduced owing to the tempering of martensite. In contrast, martensite formed in the upper layers is not tempered to the same extent and the hardness remains high. Therefore, the hardness of the part increases with the distance from the substrate. Variations in the heat input at different laser powers and scanning speeds significantly affect the extent of tempering. Finally, since the method used here can provide a quantitative understanding of the tempering of martensite and the spatial variation in hardness, it can be used to tailor the microstructure and hardness of heat treatable printed metallic parts.

36 MATERIALS SCIENCE↗

Smoothed aggregation for difficult stretched mesh and coefficient variation problems

Abstract Four adaptations of the smoothed aggregation algebraic multigrid (SA‐AMG) method are proposed with an eye toward improving the convergence and robustness of the solver in situations when the discretization matrix contains many weak connections. These weak connections can cause higher than expected levels of fill‐in within the coarse discretization matrices and can also give rise to suboptimal smoothing within the prolongator smoothing phase. These smoothing drawbacks are due to the relatively small size of some diagonal entries within the filtered matrix that one obtains after dropping the weak connections. The new algorithms consider modifications to the Jacobi‐like step that defines the prolongator smoother, modifications to the filtered matrix, and also direct modifications to the resulting grid transfer operators. Numerical results are given illustrating the potential benefits of the proposed adaptations.

Hu, Jonathan J.↗

Linearization errors in discrete goal-oriented error estimation

This paper is concerned with goal-oriented a posteriori error estimation for nonlinear functionals in the context of nonlinear variational problems solved with continuous Galerkin finite element discretizations. A two-level, or discrete, adjoint-based approach for error estimation is considered. The traditional method to derive an error estimate in this context requires linearizing both the nonlinear variational form and the nonlinear functional of interest which introduces linearization errors into the error estimate. In this paper, we investigate these linearization errors. In particular, we develop a novel discrete goal-oriented error estimate that accounts for traditionally neglected nonlinear terms at the expense of greater computational cost. We demonstrate how this error estimate can be used to drive mesh adaptivity. Here, we show that accounting for linearization errors in the error estimate can improve its effectivity for several nonlinear model problems and quantities of interest. We also demonstrate that an adaptive strategy based on the newly proposed estimate can lead to more accurate approximations of the nonlinear functional with fewer degrees of freedom when compared to uniform refinement and traditional adjoint-based approaches.

42 ENGINEERING↗

A performant energy-conserving particle reweighting method for Particle-in-Cell simulations

A new particle-based reweighting method is developed and demonstrated in the Aleph Particle-in-Cell with Direct Simulation Monte Carlo (PIC-DSMC) program. Novel splitting and merging algorithms ensure that modified particles maintain physically consistent positions and velocities. This method allows a single reweighting simulation to efficiently model plasma evolution over orders of magnitude variation in density, while accurately preserving energy distribution functions (EDFs). Demonstrations on electrostatic sheath and collisional rate dynamics show that reweighting simulations achieve accuracy comparable to fixed weight simulations with substantial computational time savings. This highly performant reweighting method is recommended for modeling plasma applications that require accurate resolution of EDFs or exhibit significant density variations in time or space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗