Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple solutions”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Understanding and improving model representation of aerosol optical properties for a Chinese haze event measured during KORUS-AQ

KORUS-AQ was an international cooperative air quality field study in South Korea that measured local and remote sources of air pollution affecting the Korean Peninsula during May–June 2016. Some of the largest aerosol mass concentrations were measured during a Chinese haze transport event (24 May). Air quality forecasts using the WRF-Chem model with aerosol optical depth (AOD) data assimilation captured AOD during this pollution episode but overpredicted surface particulate matter concentrations in South Korea, especially PM 2.5 , often by a factor of 2 or larger. Analysis revealed multiple sources of model deficiency related to the calculation of optical properties from aerosol mass that explain these discrepancies. Using in situ observations of aerosol size and composition as inputs to the optical properties calculations showed that using a low-resolution size bin representation (four bins) underestimates the efficiency with which aerosols scatter and absorb light (mass extinction efficiency). Besides using finer-resolution size bins (8–16 bins), it was also necessary to increase the refractive indices and hygroscopicity of select aerosol species within the range of values reported in the literature to achieve better consistency with measured values of the mass extinction efficiency (6.7 m 2 g –1 observed average) and light-scattering enhancement factor (ƒ(RH)) due to aerosol hygroscopic growth (2.2 observed average). Furthermore, an evaluation of the optical properties obtained using modeled aerosol properties revealed the inability of sectional and modal aerosol representations in WRF-Chem to properly reproduce the observed size distribution, with the models displaying a much wider accumulation mode. Other model deficiencies included an underestimate of organic aerosol density (1.0 g cm –3 in the model vs. observed average of 1.5 g cm –3 ) and an overprediction of the fractional contribution of submicron inorganic aerosols other than sulfate, ammonium, nitrate, chloride, and sodium corresponding to mostly dust (17 %–28 % modeled vs. 12 % estimated from observations). These results illustrate the complexity of achieving an accurate model representation of optical properties and provide potential solutions that are relevant to multiple disciplines and applications such as air quality forecasts, health impact assessments, climate projections, solar power forecasts, and aerosol data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗
DOE OSTI2021

Explainable Artificial Intelligence in Endocrinological Medical Research

Artificial Intelligence (AI) has been a part of the medical community for decades in the form of Clinical Decision Support Systems (CDSS) to aide physicians in diagnosis and categorization of patients (1). Recent years have seen a shift from expert-derived models to the integration of machine learning (ML) algorithms to drive the output of these AI systems due to the ability of ML models to more accurately make predictions by exploiting higher dimensional and often complex data. In many cases ML models gain their advantage in accuracy by capture complex and often non-linear relationships between features being used to make the prediction. However, the hype and excitement around these methods are tempered by the clinical utility of often black-box solutions driven by skepticism of results that are difficult for practitioners to not only interpret but explain to their patients (1-3). This skepticism is not unfounded as multiple examples of black-box solutions identifying incidental correlates as the key predictors have highlighted the potential bias in a training set, or reward system, that a ML model may exploit; for example a model discerning wolves from huskies based on snow in the background rather than features of the dogs (1, 4, 5).

Endocrinological, artificial intelligence, diabete↗
DOE OSTI2026

A Unified Design Theory for Multi-Port Polyphase Transformers Enabling Scalable Power-Multiplexed EV Fleet Charging Systems

This paper presents a unified analytical design theory for multi-port polyphase transformers, targeting scalable and isolated high-power Electric Vehicle (EV) fleet charging systems with power multiplexing capability. As fleet electrification accelerates, conventional one-to-one charger architectures face significant challenges in infrastructure cost, peak power demand, and low utilization of installed power electronics. Power-multiplexed charging architectures, which dynamically distribute power from a shared pool of converter modules across multiple vehicles, have emerged as a promising solution. However, such architectures require scalable, isolated multi-port power interfaces capable of routing energy among multiple inputs and outputs, whose design remains complex and dependent on iterative modeling. To address this gap, the proposed theory provides closed-form expressions for self-inductance, leakage inductance, and mutual coupling terms for arbitrary multi-phase, multi-port transformer structures. The formulation enables direct synthesis of isolated multi-input and multi-output resonant converter systems without reliance on geometry-specific finite-element analysis or extensive parameter extraction. This capability is particularly critical for power-multiplexed systems, where modular converter structures must interface with multiple vehicles while maintaining galvanic isolation and flexible power allocation. The effectiveness of the proposed framework is demonstrated through the design of a 360 kW multi-phase system operating over a 700–900 VDC input and 400–1250 VDC output range. PLECS simulation results confirm accurate prediction of system behavior and validate the applicability of the approach to multi-port, power-multiplexed charging scenarios. The proposed method significantly reduces design complexity while enabling scalable, cost-effective, and fully utilized EV fleet charging infrastructure.

Asa, Erdem [ORNL] (ORCID:0000000190884812)↗
DOE OSTI2022

RDPM: An Extensible Tool for Resilience Design Patterns Modelling

Resilience to faults, errors, and failures in extreme-scale high-performance computing (HPC) systems is a critical challenge. Resilience design patterns offer a new, structured hardware and software design approach for improving resilience. While prior work focused on developing performance, reliability, and availability models for resilience design patterns, this paper extends it by providing a Resilience Design Patterns Modeling (RDPM) tool which allows (1) exploring performance, reliability, and availability of each resilience design pattern, (2) offering customization of parameters to optimize performance, reliability, and availability, and (3) allowing investigation of trade-off models for combining multiple patterns for practical resilience solutions.

Kumar, Mohit↗
DOE OSTI2024

L-PBF High-Throughput Data Pipeline Approach for Multi-modal Integration

Abstract Metal-based additive manufacturing requires active monitoring solutions for assessing part quality. Multiple sensors and data streams, however, generate large heterogeneous data sets that are impractical for manual assessment and characterization. In this work, an automated pipeline is developed that enables feature extraction from high-speed camera video and multi-modal data analysis. The framework removes the need for manual assessment through the utilization of deep learning techniques and training models in a weakly supervised paradigm. We demonstrate this pipeline’s capability over 700,000 high-speed camera frames. The pipeline successfully extracts melt pool and spatter geometries and links them to corresponding pyrometry, radiography, and processparameter information. 715 individual prints are examined to reveal melt pool areas that exceeds 0.07 mm 2 and pyrometry signal over a threshold (375 pyrometry units) were more likely to have defects. These automated processes enable massive throughput of characterization techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nickel hydroxide–nickel carbonate competitive growth on carbonate surfaces

The thermodynamic and kinetic factors controlling the competitive heterogeneous nucleation and growth of ubiquitous metal carbonate and hydroxide phases are poorly understood. In this work, calcite (CaCO 3 ) and magnesite (MgCO 3 ) powders were reacted with NiCl 2 (0–600 μM) for 7 days at 22 °C and 5 °C. The reacted powders were analyzed with X-ray photoelectron spectroscopy (XPS), scanning electron microscopy, and energy-dispersive X-ray spectroscopy to characterize the Ni surface precipitates formed. Evidence from these techniques pointed to the formation of mixed Ni carbonate-Ni hydroxide amorphous surface precipitates. On calcite, XPS detected primarily Ni(OH) 2 despite the initial solutions being more supersaturated with respect to gaspéite (NiCO 3 ) than to theophrastite (Ni(OH) 2 ) by a factor of 17–18. In contrast, NiCO 3 was the dominant component detected by XPS on magnesite in the same conditions. Decreasing the temperature had the effect of increasing the proportion of NiCO 3 to the detriment of Ni(OH) 2 . The experimental observations were consistent with low lattice/cation size mismatch favoring NiCO 3 nucleation and temperature most influencing Ni(OH) 2 nucleation. Comparison to previous work on Co-reacted powders indicated the differences in lattice/cation size mismatch and/or water exchange rate impacted the composition of the surface precipitates more than the relative thermodynamic stabilities of the competing minerals. This work explored the heterogeneous growth regime of Ni carbonate and hydroxide phases on carbonate surfaces and shed light on the factors that control the competition between surface precipitates when mineral surfaces are in contact with aqueous solutions supersaturated with respect to multiple mineral phases. In conclusion, these results contribute to geochemists’ efforts toward interpreting data from geochemical systems with elevated Ni concentrations, improving Ni environmental remediation and recovery strategies, and predicting the fate and transport of Ni in geochemical systems.

Heterogeneous nucleation↗
DOE OSTI2021

The effect of hydrogen-bonding on flavin’s infrared absorption spectrum

Cluster and continuum solvation computational models are employed to model the effect of hydrogen bonding interactions on the vibrational modes of lumiflavin. Calculated spectra were compared to experimental Fourier-transform infrared (FTIR) spectra in the diagnostic 1450–1800 cm –1 range, where intense V C=C , V C=N , V C 2 =O , and V C 4 =O stretching modes of flavin’s isoalloxazine ring are found. Local mode analysis is used to describe the strength of hydrogen-bonding in cluster models. The computations indicate that V C=C and V C=N mode frequencies are relatively insensitive to intermolecular interactions while the V C 2 =O and V C 4 =O modes are sensitive to direct (and also indirect for V C 2 =O ) hydrogen-bonding interactions. Although flavin is neutral, basis sets without the diffuse functions provide incorrect relative frequencies and intensities. The 6-31+G* basis set is found to be adequate for this system, and there is limited benefit to considering larger basis sets. Calculated vibrational mode frequencies agree with experimentally determined frequencies in solution when cluster models with multiple water molecules are used. Accurate simulation of relative FTIR band intensities, on the other hand, requires a continuum (or possibly quantum mechanical/molecular mechanical) model that accounts for long-range electrostatic effects. Lastly, an experimental peak at ca. 1624 cm –1 that is typically assigned to the V C 2 =O vibrational stretching mode has a complicated shape that suggests multiple underlying contributions. Our calculations show that this band has contributions from both the C 6 –C 7 and C 2 = O stretching vibrations.

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

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Binding of Sulfates and Water to Monovalent Cations

The binding of the sulfate ligand group to monovalent cations in the presence of water is important for many systems. To understand the structure and energetics of sulfate complexes, we use density functional theory to study ethyl sulfate binding to the monovalent cations Li + , Na + , and K + , and to water. The free energies of binding and optimal structures are calculated for a range of the number of ethyl sulfates and waters. Without water, the most optimal structure for all the cations is bidentate binding by two ethyl sulfates, yielding a 4-fold coordination. With water, the lowest free energy structures also have two ethyl sulfates, but the coordination varies with cations. For complexes with water, the four oxygen atoms in the sulfate group enable multiple binding geometries for the cations and for hydrogen bonding with water. Many of these geometries differ in free energy by only a small amount (1–2 kcal/mol), meaning there will be multiple binding configurations in bulk solution. In comparison to the optimal structures for binding to the carboxylate group, there is more variation for binding to the sulfate group as a function of cation type and the number of waters. Further, the polarization of the atoms is significant and varies among the sulfate oxygen atoms. The water oxygen charge is often larger than that of sulfate oxygen, which plays a role in the preference for monodentate ligand binding to cations in the presence of water.

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

Scalable nanomanufacturing of chalcogenide inks: a case study on thermoelectric V–VI nanoplates

Solution-processed semiconducting main-group chalcogenides (MMCs) have attracted increasing research interest for next-generation device technologies owing to their unique nanostructures and superior properties. To achieve the full potential of MMCs, the development of highly universal, scalable, and sustainable synthesis and processing methods of chalcogenide particles is thus becoming progressively more important. Here we studied scalable factors for the synthesis of two-dimensional (2D) V–VI chalcogenide nanoplates (M 2 Q 3 : M = Sb, Bi; Q = Se, Te) and systematically investigated their colloidal behaviour and chemical stability. Based on a solvent engineering technique, we demonstrated scale-up syntheses of MMCs up to a 900% increase of batch size compared with conventional hydrazine-based gram-level syntheses, and such a scalable approach is highly applicable to various binary and ternary MMCs. Furthermore, we studied the stability of printable chalcogenide nanoparticle inks with several formulation factors including solvents, additives, and pH values, resulting in inks with high chemical stability (>4 months). As a proof of concept, we applied our solution-processed chalcogenide particles to multiple additive manufacturing methods, confirming the high printability and processability of MMC inks. Furthermore, the ability to combine the top-down designing freedom of additive manufacturing with bottom-up scalable synthesis of chalcogenide particles promises great opportunities for large-scale design and manufacturing of chalcogenide-based functional devices for broad application.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Shifts in valence states in bimetallic MXenes revealed by electron energy-loss spectroscopy (EELS)

MXenes are an emergent class of two-dimensional materials with a very wide spectrum of promising applications. The synthesis of multiple MXenes, specifically solid-solution MXenes, allows fine tuning of their properties, expands their range of applications, and leads to enhanced performance. The functionality of solid-solution MXenes is closely related to the valence state of their constituents: transition metals, oxygen, carbon, and nitrogen. However, the impact of changes in the oxidation state of elements in MXenes is not well understood. In this work, three interrelated solid-solution MXene systems (Ti$_{2–y}$Nb$_y$CT$_x$, Nb$_{2–y}$V$_y$CT$_x$, and Ti$_{2–y}$V$_y$CT$_x$) were investigated with scanning transmission electron microscopy and electron energy-loss spectroscopy to determine the localized valence states of metals at the nanoscale. The analysis demonstrates changes in the electronic configuration of V upon modification of the overall composition and within individual MXene flakes. These shifts of oxidation state can explain the nonlinear optical and electronic features of solid-solution MXenes. Vanadium appears to be particularly sensitive to modification of the valence state, while titanium maintains the same oxidation state in Ti–Nb and Ti–V MXenes, regardless of stoichiometry. Here, the study also explains Nb's influential role in the previously observed electronic properties in the Nb–V and Nb–Ti systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CrysFieldExplorer : rapid optimization of the crystal field Hamiltonian

A new approach to the fast optimization of crystal electric field (CEF) parameters to fit experimental data is presented. This approach is implemented in a lightweight Python-based program, CrysFieldExplorer. The main novelty of the method is the development of a unique loss function, referred to as the spectrum characteristic loss (L Spectrum ), which is based on the characteristic polynomial of the Hamiltonian matrix. Particle swarm optimization and a covariance matrix adaptation evolution strategy are used to find the minimum of the total loss function. It is demonstrated that CrysFieldExplorer can perform direct fitting of CEF parameters to any experimental data such as a neutron spectrum, susceptibility or magnetization measurements etc. CrysFieldExplorer can handle a large number of non-zero CEF parameters and reveal multiple local and global minimum solutions. Crystal field theory, the loss function, and the implementation and limitations of the program are discussed within the context of two examples.

42 ENGINEERING↗
DOE OSTI2022

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Predicting Radiation-Induced Plutonium Redox Chemistry using Multiscale Modeling Methods

Over the the last 70 years plutonium (Pu) has been integral in the development of several technologies that have changed the world, yet our fundamental understanding of its chemistry is still far from complete. This is a testament to this element?s unique and complex properties, such as its ability to coexist as multiple oxidation states in aqueous solution. Careful manipulation of plutonium oxidation states is essential in the study and utilization of its rich chemistry. To achieve this level of control, a comprehensive mechanistic understanding of radiation-induced plutonium redox chemistry is critical due to the unavoidable exposure of plutonium to ionizing radiation fields, both inherent and from in-process applications. For this reason, we have developed an experimentally evaluated multi-scale computer model for the prediction of gamma radiation-induced Pu(IV) redox chemistry in concentrated nitric acid solutions (1.0, 3.0, and 6.0 M).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Designing FAIR Workflows at OLCF: Building Scalable and Reusable Ecosystems for HPC Science

High Performance Computing (HPC) centers, such as the Oak Ridge Leadership Computing Facility (OLCF), provide advanced infrastructure that enables scientific research at extreme scale. These centers operate with unique hardware configurations, specialized software environments, and elevated security re quirements that differ substantially from what most users encounter on their local systems. As a result, users often develop customized digital artifacts that are tightly coupled to the specific configuration of a given HPC center. Although necessary, this practice can lead to significant duplication of effort as multiple users independently create similar solutions to common problems.

97 MATHEMATICS AND COMPUTING↗