Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical 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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Predicting large-scale pool fire dynamics using an unsteady flamelet- and large-eddy simulation-based model suite

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. Here, the coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier–Stokes/eddy dissipation concept-based methodology.

42 ENGINEERING↗
DOE OSTI2024

On the Statistical Mechanics of Mass Accommodation at Liquid–Vapor Interfaces

Here we propose a framework for describing the dynamics associated with the adsorption of small molecules to liquid-vapor interfaces using an intermediate resolution between traditional continuum theories that are bereft of molecular detail and molecular dynamics simulations that are replete with them. In particular, we develop an effective single particle equation of motion capable of describing the physical processes that determine thermal and mass accommodation probabilities. The effective equation is parametrized with quantities that vary through space away from the liquid-vapor interface. Of particular importance in describing the early time dynamics is the spatially dependent friction, for which we propose a numerical scheme to evaluate from molecular simulation. Taken together with potentials of mean force computable with importance sampling methods, we illustrate how to compute the mass accommodation coefficient and residence time distribution. Throughout, we highlight the case of ozone adsorption in aqueous solutions and its dependence on electrolyte composition.

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

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Solving Coupled Cluster Equations by the Newton Krylov Method

We describe using the Newton Krylov method to solve the coupled cluster equation. The method uses a Krylov iterative method to compute the Newton correction to the approximate coupled cluster amplitude. The multiplication of the Jacobian with a vector, which is required in each step of a Krylov iterative method such as the Generalized Minimum Residual (GMRES) method, is carried out through a finite difference approximation, and requires an additional residual evaluation. The overall cost of the method is determined by the sum of the inner Krylov and outer Newton iterations. We discuss the termination criterion used for the inner iteration and show how to apply pre-conditioners to accelerate convergence. We will also examine the use of regularization technique to improve the stability of convergence and compare the method with the widely used direct inversion of iterative subspace (DIIS) methods through numerical examples.

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

Analysis of persistent contaminants and personal care products by dispersive liquid-liquid microextraction using hydrophobic magnetic deep eutectic solvents

Here, in this work, hydrophobic magnetic deep eutectic solvents (HMDESs) were used in the development of a simple and rapid dispersive liquid-liquid microextraction (DLLME) approach coupled to high performance liquid chromatography with UV detection (HPLC-UV) for the determination of ten organic contaminants including five polycyclic aromatic hydrocarbons, four UV filters, and a pesticide from water at trace levels. The HMDESs were prepared by mixing a hydrogen bond acceptor, metal halide salt, and hydrogen bond donor in suitable molar ratios. Two HMDESs, 2 tetraoctylammonium bromide ([N 8888 + ][Br - ]): cobalt chloride (CoCl 2 ): 4 octanoic acid (OA) and 3 trioctylphosphine oxide (TOPO): neodymium chloride (NdCl 3 ): 3 OA, offered the highest analyte extraction efficiency overall and were chosen as suitable solvents for validation of the microextraction method. Under optimized extraction conditions, the method required 30 µL of HMDES as extraction solvent, acetone (87.5 µL) as disperser solvent, a NaCl concentration of 30% (w/v), and an extraction time of 120 s at 20°C. Enrichment factors of the analytes ranged from 44.6 for 3-(4-methylbenzylindene) camphor to 66.0 for 2-ethylhexyl-4-(dimethyl)aminobenzoate. The method provided low limits of detection (LODs) ranging from 0.5 to 4.5 µg L -1 , and acceptable precision, with RSD values lower than 9.6%. Furthermore, the validated method was successfully applied for tap and lake water analysis, resulting in relative recoveries of spiked samples ranging between 94.7 and 119.2%.

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

Modeling and implementation of vertical excursion FFA in the Zgoubi ray-tracing code

Vertical Fixed Field Accelerators (vFFAs) feature complex and highly non-linear magnetic fields, which require simulation codes allowing step-wise particle tracking. Methods to model the 3D magnetic field of scaling vFFAs have been developed in the ray-tracing code Zgoubi. The field modeling and particle tracking methods include the field non-linearities, the fringe fields, and the field superposition of neighboring magnets. The procedure implements the vFFA analytical field expressions, allowing design studies and parameter optimizations using the Zgoubi built-in fit method. The vFFA procedure has been applied to a ten-fold symmetry ring with a triplet focusing structure designed to accelerate protons from 3 MeV to 12 MeV and studied under the ISIS-II proton driver prototype project. Results from particle tracking in externally generated 3D semi-analytical field maps and the developed vFFA analytical model are shown to be in excellent agreement.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

NanoSIP: NanoSIMS Applications for Microbial Biology

High-resolution imaging with secondary ion mass spectrometry (nanoSIMS) has become a standard method in systems biology and environmental biogeochemistry and is broadly used to decipher ecophysiological traits of environmental microorganisms, metabolic processes in plant and animal tissues, and cross-kingdom symbioses. When combined with stable isotope-labeling—an approach we refer to as nanoSIP—nanoSIMS imaging offers a distinctive means to quantify net assimilation rates and stoichiometry of individual cell-sized particles in both low- and high-complexity environments. While the majority of nanoSIP studies in environmental and microbial biology have focused on nitrogen and carbon metabolism (using 15 N and 13 C tracers), multiple advances have pushed the capabilities of this approach in the past decade. The development of a high-brightness oxygen ion source has enabled high-resolution metal analyses that are easier to perform, allowing quantification of metal distribution in cells and environmental particles. New preparation methods, tools for automated data extraction from large data sets, and analytical approaches that push the limits of sensitivity and spatial resolution have allowed for more robust characterization of populations ranging from marine archaea to fungi and viruses. Further, NanoSIMS studies continue to be enhanced by correlation with orthogonal imaging and ‘omics approaches; when linked to molecular visualization methods, such as in situ hybridization and antibody labeling, these techniques enable in situ function to be linked to microbial identity and gene expression. Here we present an updated description of the primary materials, methods, and calculations used for nanoSIP, with an emphasis on recent advances in nanoSIMS applications, key methodological steps, and potential pitfalls.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Bridging Equipment Reliability Data and Risk Informed Decisions in a Plant Operation Context

Industry equipment reliability and asset management programs are essential elements that help ensure the safe and economical operation of nuclear power plants. The effectiveness of these programs is addressed in several industry-developed and regulatory programs. The Risk-Informed Asset Management (RIAM) project is tasked to develop tools in support of the equipment reliability and asset management programs at nuclear power plants. These tools are designed to create a direct bridge between component health/lifecycle data and decision making (e.g., maintenance scheduling and project prioritization). The goal of this article is to provide a guide for specific use cases that the RIAM project is targeting. We have grouped uses cases into three main areas. The first area focuses on the analysis of equipment reliability data with a particular emphasis on condition-based data, such as test/surveillance reports and component monitoring data. The second area focuses on the integration of equipment reliability into system/plant reliability models to determine system/plant health and identify the components that are critical to maintain an operational system. Lastly, the third area manages plant resources, such as maintenance activities and replacement scheduling using optimization methods. Here the primary focus is on supporting typical system engineer decisions regarding maintenance activity scheduling and component aging management. This is performed in a risk-informed context where the term “risk” is broadly constructed to include both plant reliability and economics. This framework combines data analytics tools to analyze equipment reliability data with risk-informed methods designed to support system engineer decisions (e.g., maintenance and replacement schedules, optimal maintenance posture) in a customizable workflow.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Analysis of the SiMPL Method for Density-Based Topology Optimization

We present a rigorous convergence analysis of a new method for density-based topology optimization that provides pointwise bound-preserving design updates and faster convergence than other popular first-order topology optimization methods. Due to its strong bound preservation, the method is exceptionally robust, as demonstrated in numerous examples here and in the companion article [D. Kim et al., Struct. Multidiscip. Optim., 68 (2025), 74]. Furthermore, it is easy to implement with clear structure and analytical expressions for the updates. Our analysis covers two versions of the method, characterized by the employed line search strategies. We consider a modified Armijo backtracking line search and a Bregman backtracking line search. For both line search algorithms, our algorithm delivers a strict monotone decrease in the objective function and further intuitive convergence properties, e.g., strong and pointwise convergence of the density variables on the active sets, norm convergence to zero of the increments, convergence of the Lagrange multipliers, and more. In addition, the numerical experiments demonstrate apparent mesh-independent convergence of the algorithm. Here, we refer to the new algorithm as the SiMPL method (pronounced “simple”), which stands for Sigmoidal Mirror descent with a Projected Latent variable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

From anti-Arrhenius to Arrhenius behavior in a dislocation-obstacle bypass: Atomistic simulations and theoretical investigation

Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Integrating in situ environmental covariates in an American lobster catch model to improve impact assessment

The installation and operation of floating offshore wind power is an integral component of societal transition to renewable energy generation where fixed bottom offshore wind is not possible. However, it will cause unique ecosystem changes. To disentangle the effects of offshore wind installations from the concurrent effects of climate change and the fishing practices on commercially significant resources, we must develop detailed characterizations of the resources before development occurs. In the Gulf of Maine, American lobster is the most commercially and culturally important fishery. At the time of writing, this is the largest fishery by value in North America. Our understanding of baseline localized parameters (such as catch per trap at the spatial scale of individual turbines) should be informed by relationships to environmental, biological, and survey-specific functional drivers of catch. A more mechanistic understanding of catch will allow for strategic adjustments to Post- Deployment fishery responses and ultimately, the development of research- and commercial-scale floating offshore wind development. Here, we used survey data from the New England Aqua Ventus Pre-Construction Commercial Trapping Survey to develop Generalized Additive Models describing seasonal catch per trap for legal and sublegal lobsters. We found fall catch to be nearly twice that of spring. Bottom temperature dynamics could be used to predict catch, and the Fall survey was associated with a warmer temperature regime. By using analytical tools that incorporate environmental heterogeneity, we developed monitoring methods from preconstruction baseline data that will be applicable over the post-construction operating period of an offshore wind farm.

BACI↗
DOE OSTI2025

Dehydrogenation vs Apparent Hydrogenation: Unraveling the Mechanisms of He and O 2 Plasma Etching on Colloidal Nanocrystal Films

Removing organic ligands from colloidal nanoparticles is critical for fabricating solid-state devices, yet accurately quantifying this removal remains a significant analytical challenge. Here, we establish a robust and accessible method for this quantification by calibrating Raman spectroscopy against precise ion beam analysis (IBA) for nanoparticle assemblies (CNAs) processed by helium (He) and oxygen (O 2 ) plasmas. We demonstrate that the calibration curves are remarkably independent of plasma power and pressure, depending critically only on the choice of feed gas. He plasma induces rapid dehydrogenation and cross-linking, evidenced by a much faster decrease in the C–H Raman signal relative to the actual carbon loss. Conversely, O 2 plasma leads to a surprising “apparent hydrogenation”, where the carbon backbone is removed significantly faster than the C–H signal diminishes. This counterintuitive effect is explained by a serial mechanism of oxidative fragmentation; β-scission cleaves the alkyl chains, and subsequent stabilization steps enrich the remaining film with hydrogen-rich methyl-terminated fragments, while carbon is efficiently removed as volatile CO. This work provides calibrated functions that enable the rapid determination of absolute carbon content in processed CNAs using simple Raman spectroscopy with uncertainties of ∼8% for O 2 and ∼12% for He plasma, offering a vital tool for both process diagnostics and fundamental studies of plasma–matter interactions in colloidal nanocrystal films.

Animal feed↗
DOE OSTI2024

Verification and Demonstration of One-Dimensional Freezing Model in SAM for Salt-Cooled Reactor Analysis Applications

This work presented the development and implementation of the one-dimensional freezing model in system analysis code, SAM, as well as code verification, and code demonstration during a postulated overcooling transient, for fluoride salt-cooled high-temperature reactor (FHR) system and safety analysis applications. The paper at first summarized the freezing model, finite element numerical method, and special numerical treatment for handling phase appearance/disappearance. Analytical solutions were derived for two cases (with and without solid walls) for code verifications purpose. As expected, numerical results predicted by the SAM code agreed very well with the analytical solution. A code demonstration was then performed on a postulated protected overcooling event transient of a generic reference PB-FHR design. The code was found to successfully predict salt freezing during such a postulated event. However, due to lack of salt freezing testing data, code validation has not been performed in this work, which will be pursued in later studies when such data becomes available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Structure of native four-repeat satellite III sequence with non-canonical base interactions

Abstract Tandem-repetitive DNA (where two or more DNA bases are repeated numerous times) can adopt non-canonical secondary structures. Many of these structures are implicated in important biological processes. Human Satellite III (HSat3) is enriched for tandem repeats of the sequence ATGGA and is located in pericentromeric heterochromatin in many human chromosomes. Here, we investigate the secondary structure of the four-repeat HSat3 sequence 5′-ATGGA ATGGA ATGGA ATGGA-3′ using X-ray crystallography, NMR, and biophysical methods. Circular dichroism spectroscopy, thermal stability, native PAGE, and analytical ultracentrifugation indicate that this sequence folds into a monomolecular hairpin with non-canonical base pairing and B-DNA characteristics at concentrations below 0.9 mM. NMR studies at 0.05–0.5 mM indicate that the hairpin is likely folded-over into a compact structure with high dynamics. Crystallographic studies at 2.5 mM reveal an antiparallel self-complementary duplex with the same base pairing as in the hairpin, extended into an infinite polymer. The non-canonical base pairing includes a G–G intercalation sandwiched by sheared A–G base pairs, leading to a cross-strand four guanine stack, so called guanine zipper. The guanine zippers are spaced throughout the structure by A–T/T–A base pairs. Our findings lend further insight into recurring structural motifs associated with the HSat3 and their potential biological functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A glimpse into the fungal metabolomic abyss: Novel network analysis reveals relationships between exogenous compounds and their outputs

Fungal specialized metabolites are a major source of beneficial compounds that are routinely isolated, characterized, and manufactured as pharmaceuticals, agrochemical agents, and industrial chemicals. The production of these metabolites is encoded by biosynthetic gene clusters that are often silent under standard growth conditions. There are limited resources for characterizing the direct link between abiotic stimuli and metabolite production. Herein, we introduce a network analysis-based, data-driven algorithm comprising two routes to characterize the production of specialized fungal metabolites triggered by different exogenous compounds: the direct route and the auxiliary route. Both routes elucidate the influence of treatments on the production of specialized metabolites from experimental data. The direct route determines known and putative metabolites induced by treatments and provides additional insight over traditional comparison methods. The auxiliary route is specific for discovering unknown analytes, and further identification can be curated through online bioinformatic resources. We validated our algorithm by applying chitooligosaccharides and lipids at two different temperatures to the fungal pathogen Aspergillus fumigatus. After liquid chromatography–mass spectrometry quantification of significantly produced analytes, we used network centrality measures to rank the treatments’ ability to elucidate these analytes and confirmed their identity through fragmentation patterns or in silico spiking with commercially available standards. Later, we examined the transcriptional regulation of these metabolites through real-time quantitative polymerase chain reaction. Our data-driven techniques can complement existing metabolomic network analysis by providing an approach to track the influence of any exogenous stimuli on metabolite production. Our experimental-based algorithm can overcome the bottlenecks in elucidating novel fungal compounds used in drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Code associated with Publication “Analytic Solutions and Field-Scale Application for Verification of Coupled Thermo-Hydro-Mechanical Processes in Subsurface Fractured Media”

As part of a submitted paper, which is collection of previously published analytical solutions to coupled thermo-hydro-mechanical problems in subsurface flow and transport, we have prepared a collection of python scripts to compute and plot those analytic solutions. All code to be released implements existing methods; there are no novel algorithms nor any major innovations to existing software.

Hyman, Jeffrey↗
DOE OSTI2025

The Measurement of Grain Size by Electron Beam Backscattered Diffraction

The method of electron back-scatter diffraction provides an analytic means to identify crystalline phase and the grain size of structural materials. The need to utilize imaging parameters consistent with revealing the contributing microstructural scale is often overlooked through packaged imaging scripts. It is shown that imaging parameters consistent with a post processing assessment of low misorientation angles between boundaries are achievable as using pixel sizes on the order of 10 1 nm.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Characterizing the Reproducibility of Noisy Quantum Circuits

The ability of a quantum computer to reproduce or replicate the results of a quantum circuit is a key concern for verifying and validating applications of quantum computing. Statistical variations in circuit outcomes that arise from ill-characterized fluctuations in device noise may lead to computational errors and irreproducible results. While device characterization offers a direct assessment of noise, an outstanding concern is how such metrics bound the reproducibility of a given quantum circuit. Here, we first directly assess the reproducibility of a noisy quantum circuit, in terms of the Hellinger distance between the computational results, and then we show that device characterization offers an analytic bound on the observed variability. We validate the method using an ensemble of single qubit test circuits, executed on a superconducting transmon processor with well-characterized readout and gate error rates. The resulting description for circuit reproducibility, in terms of a composite device parameter, is confirmed to define an upper bound on the observed Hellinger distance, across the variable test circuits. This predictive correlation between circuit outcomes and device characterization offers an efficient method for assessing the reproducibility of noisy quantum circuits.

97 MATHEMATICS AND COMPUTING↗