Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local minimum”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Benchmarking image processing techniques for porosity measurement in polymer additive manufacturing: Review and experimental analysis

An image processing workflow is proposed for porosity measurement in polymer additive manufacturing. Various techniques, including global and local thresholding, region growing, and K-means clustering, were applied to microscopic images of carbon fiber reinforced acrylonitrile butadiene styrene (CF-ABS) and benchmarked for their ability to accurately measure porosity. Global methods included Otsu, minimum error, iterative, and entropy-based thresholding, while local methods included Niblack, Bernsen, Sauvola, and Bradley-Roth algorithms. Artificial uneven illumination was introduced to test local adaptive thresholds. Results showed significant differences in porosity values across methods. Otsu, region growing, and K-means clustering excelled under uniform illumination, while Sauvola and Bradley-Roth performed better with uneven illumination. Comparison with X-ray computed tomography (XCT) revealed slightly lower porosity values (2.55 %) than optimized methods (2.73–2.79 %) due to XCT's lower resolution excluding smaller pores. While XCT offers finer pore detection, it limits sample volume and underestimates porosity due to spatial variation. Validation using artificial grayscale images with 5 % porosity confirmed that Otsu, Bradley-Roth, region growing, and Sauvola algorithms produced accurate results. Although tested on a single material system, these methods can be adapted to others with optimization. In conclusion, given XCT's high computational and time costs, this study highlights suitable image processing techniques as cost-effective alternatives for porosity analysis in polymer composites.

Additive manufacturing↗
DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Temperature estimation in a spatially inhomogeneous flame by diode laser absorption spectroscopy

A two-temperature (2T) model is proposed for estimating the temperature of a spatially inhomogeneous hot zone using single-beam diode laser absorption spectrometry. The proposed algorithm is based on fitting experimental absorption spectra by a linear combination of two database-simulated single-temperature spectra with different temperatures. The model efficiency is experimentally demonstrated when determining the temperature in different flame zones of a Wolfgard – Parker slot burner. The proposed 2T algorithm is used to find the maximum and minimum temperatures in different flame sections, which are compared with local temperatures in the same flame zones measured by the method of coherent anti-Stokes Raman scattering (CARS). It is shown that the maximum temperatures determined by the proposed 2T algorithm are in good agreement with CARS data in flame zones where high-temperature regions prevail. A good agreement between the minimum temperatures obtained by these two methods is also observed for predominantly cold zones. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Generating Bright Emissive States by Modulating the Bandgap of Monolayer Tungsten Diselenide

Transition metal dichalcogenides (TMDs) are essential due to their fascinating electronic and optical properties, strong exciton binding energy, and layer-dependent bandgap. They can be tuned to function as a single-photon emitter, but the quantum yield of photoluminescence of single-layer WSe 2 is low. There is some evidence that a bright emissive state can be produced by introducing local defects through functionalization. In this paper we use spin-polarized periodic density functional theory (DFT) to study the effect of functionalization for the specific case of cyclic carbenes. We find that the simplest of these molecules, cyclopentadiene carbene (Cyc), binds to the surface by a covalent bond leading to a substantial change in the bandgap (1.24 eV compared to 1.64 eV for the pristine surface). Moreover, there are semiflat bands below the Fermi level that originate from σ bonding and the interaction between the lone-pair p orbitals on the carbon of Cyc that bonds with Se on the TMD. Cyc is found to form a Type IIa heterojunction before and after contact with the surface, where the “a” refers to ordering of the energy levels where the TMD levels have the larger splitting and “b” to the reverse. The effect of various electron-withdrawing and electron-donating groups on Cyc is investigated, and it is found that the direct bandgap and heterointerfaces can be chemically tuned with covalently bound functional groups. With an electron-withdrawing group (EWG) such as –CHO, –COCl, and –CN attached to Cyc, a Type III//Type IIa interface is formed, whereas all other EWGs used in this study form a Type IIa junction, before and after contact to the surface. The electron-donating groups (EDG) form Type Ib//Type Ia junctions, and in particular, we find that the system WSe 2 + Cyc-Me band structure consists of semiflat bands at valence band maxima localized on Cyc-Me and the conduction band minimum is coupled between the surface and functional group with a direct bandgap of 0.88 eV. Hence, we predict that the Cyc-Me-functionalized monolayer WSe 2 will produce a red-shifted bright emissive state. Furthermore, we also find that carbenes with the Cyc all-carbon ring have (1) a triplet ground state and (2) form covalent bonds to WSe 2 , while heterocyclic carbenes (1) have a singlet ground state and (2) do not form covalent bonds to WSe 2 .

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

Hybrid Density Functional Theory Comparison of Oxygen Release and Solvent Decomposition Kinetics on Li x NiO 2 Surfaces

High-nickel-content layered oxides are among the most promising electric vehicle battery cathode materials. However, their interfacial reactivity with electrolytes and tendency toward oxygen release (possibly yielding reactive 1 O 2 ) remain degradation concerns. Elucidating the most relevant (i.e., fastest) interfacial degradation mechanism will facilitate future mitigation strategies. We apply screened hybrid density functional (HSE06) calculations to compare the reaction kinetics of Li x NiO 2 surfaces with ethylene carbonate (EC) with those of O 2 release. Further, on both the (001) and (104) facets, EC oxidative decomposition exhibits lower activation energies than O 2 release. Our calculations, coupled with previously computed liquid-phase reaction rates of 1 O 2 with EC, strongly question the role of “reactive 1 O 2 ” species in electrolyte oxidative degradation. The possible role of other oxygen species is discussed. To deal with the challenges of modeling Li x NiO 2 surface reactivity, we emphasize a “local structure” approach instead of pursuing the global energy minimum.

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

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

Abstract We develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Geochemistry & Geophysics↗
DOE OSTI2024

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Mixed-Material Scintillators (MMSS Quarterly Report FY20Q2)

This project aims to invent, model, and prototype mixed-material scintillator systems (MMSSs), a new class of radiation detectors that use heterogenous internal structures of different scintillating materials to detect additional properties of radiation. These internal structures can be produced using additive manufacture (3D printing) of scintillator, a currently emerging application of additive manufacture technology. MMSSs combine the low cost and complexity of conventional scintillation detectors with capabilities currently only available in more expensive and complex detectors. By identifying promising MMSS designs, this project will enable a new class of detectors to meet DNN’s mission needs for SNM detection. This quarter, we delivered one the MMSS Particle-ID Modeling report, a major deliverable of this project. This report describes our results on simulations of the particle-ID (PID) class of detectors. We found that these detectors outperform competing approaches for neutron-gamma discrimination, neutron source pointing, and neutron spectroscopy. Further highlights of the report are described below. During this quarter, We alerted HQ to a change in schedule for the remaining simulations planned for this project: we plan to put off further simulation and reporting of the position-resolving (PR) detector class until July. This change in schedule has allowed us to go into additional detail in the PID work that is both very promising and motivates our upcoming proposal. The other active task in this project is prototyping 3D-printed scintillators with the aim of realizing the designs we’ve simulated. On that front, we had three focuses: measuring the light output, increasing the size of printed parts, and combining blue and green scintillators. We succeeded in measuring the light output, showing an output of 30% of commercial standard. This result leaves room for improvement but is within striking distance of where we need to be. We were on track to demonstrate a key requirement for increasing the size of parts and also to show combined green/blue prints before the laboratory moved to a Minimum Safe (MinSafe) Operations posture as a result of local Shelter In Place orders due to Covid-19.

42 ENGINEERING↗
DOE OSTI2021

The Area Localized Coupled Model for Analytical Mean Flow Prediction in Arbitrary Wind Farm Geometries

This work introduces the area localized coupled (ALC) model, which extends the applicability of approaches that couple classical wake superposition models and atmospheric boundary layer models to wind farms with arbitrary layouts. Coupling wake and top–down boundary layer models is particularly challenging since the latter requires averaging over planform areas associated with turbine-specific regions of the flow that need to be specified. The ALC model uses Voronoi tessellation to define this local area around each turbine. A top–down description of a developing internal boundary layer is then applied over Voronoi cells upstream of each turbine to estimate the local mean velocity profile. Coupling between the velocity at hub-height based on this localized top–down model and a wake model is achieved by enforcing a minimum least-square-error in mean velocity in each cell. The wake model in the present implementation takes into account variations in wind farm inflow velocity and represents the wake profile behind each turbine as a super-Gaussian function that smoothly transitions between a top-hat shape in the region immediately following the turbine to a Gaussian profile downstream. Detailed comparisons to large-eddy simulation (LES) data from two different wind farms demonstrate the efficacy of the model in accurately predicting both wind farm power output and local turbine hub-height velocity for different wind farm geometries. These validations using data generated from two different LES codes demonstrate the model's versatility with respect to capturing results from different simulation setups and wind farm configurations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2022

One and Two Dimensional Quantum Lattice Algorithms for Maxwell Equations in Inhomogeneous Scalar Dielectric Media I: Theory

A quantum lattice algorithm (QLA) is developed for Maxwell equations in scalar dielectric media using the Riemann-Silberstein representation on a Cartesian grid. For x-dependent and y-dependent dielectric inhomogeneities, the corresponding QLA requires a minimum of 8 qubits/spatial lattice site. This is because the corresponding Pauli spin matrices have off-diagonal components which permit the local collisional entanglement of these qubits. However, z-dependent inhomogeneities require a QLA with a minimum of 16 qubits/lattice site since the Pauli spin matrix σz is diagonal. For 2 dimensional inhomogeneities, one can readily couple the 8-8 qubit schemes for x-y variations. z-x and y-z variations can be treated by either a 16-8 qubit scheme or a 16-16 qubit representation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Local practically safe extremum seeking with assignable rate of attractivity to the safe set

We present Assignably Safe Extremum Seeking (ASfES), an algorithm designed to minimize a measured, static objective function while maintaining a measured, static metric of safety (a control barrier function or CBF) to be positive in a practical sense. We ensure that for trajectories with safe initial conditions, the violation of safety can be made arbitrarily small through appropriately chosen design constants. We also guarantee an assignable “attractivity” rate: from unsafe initial conditions, the trajectories approach the safe set, in the sense of the measured CBF, at a rate no slower than a user-assigned rate. Similarly, from safe initial conditions, the trajectories approach the unsafe set, in the sense of the CBF, no faster than the assigned attractivity rate. The feature of assignable attractivity is not present in the semiglobal version of safe extremum seeking, where the semiglobality of convergence is achieved by slowing the adaptation. We also demonstrate local convergence of the parameter to a neighborhood of the minimum of a quadratic objective function constrained to the safe set with a linear CBF. The ASfES algorithm and analysis are multivariable, but we also extend the algorithm to a Newton-Based ASfES scheme which we show is only useful in the scalar case. The proven properties of the designs are illustrated through simulation examples.

42 ENGINEERING↗
DOE OSTI2023

Numerical Simulations of High Prandtl Number Liquid Jets Impinging on a Flat Plate

In this work, 3D simulations of oil jets impinging on a flat, heated wall are presented. The numerical setup uses the Volume of Fluid (VoF) method to model the two-phase flow. A careful grid definition across the liquid film, along with the use of the Conjugate Heat Transfer (CHT) approach allowed local heat transfer to be solved with fine resolution at the wall. Variations of liquid flow rate, liquid temperature and surface temperature allow to cover a wide range of local Reynolds and Prandtl numbers (226 < Re < 2850, 77 < Pr < 161). Resulting surface-averaged heat transfer compares very well with experimental measurements conducted in a previous study. In-depth analysis of the flow has identified expected features from the literature. In particular, the impact of jet axial velocity profiles on the heat transfer distribution in the stagnation zone was clearly stated. The increase in heat transfer when warming the liquid film was also reproduced and explained by a decrease in oil viscosity and an increase in film velocity. All those effects were taken into account in correlations for stagnation and local values of Nusselt number. A grid sensitivity study was also conducted, showing that if the grid solving the thermal boundary layer in the stagnation zone can be coarsened without impacting local and surface-averaged predictions of heat transfer, a minimum resolution (2 to 3 cells) within the thermal boundary layer is however required for an accurate prediction of heat transfer.

conjugate heat transfer↗
DOE OSTI2022

Mesoscale Investigation of Dislocation-Grain Boundary Interactions in Metals and Alloys [Slides]

GBs have been subject of intense study for several years, including the effects of microstructure on mechanical response. Until recently these have been primarily focused on a limited number of idealized GB structures (i.e, the minimum energy structure). We are specifically focusing on changes in local GB structure and the effect this has on dislocation-grain boundary interactions. This will help us to better understand both (more) realistic and general GB structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Energy cost of localization of relational quantum information

Entanglement of spatially separated quantum states is usually defined with respect to a reference frame provided by some external observer. Thus, if one wishes to localize the quantum information within a spatially separated entangled state, one must enact an entanglement extraction protocol also defined with respect to that external frame. Entanglement extraction for Gaussian ground states in such an external frame construction has been shown to require a minimum energy and is hence an interesting process for gravitational physics, where examinations of localization vs energy cost have a long history. General covariance, however, precludes dependence on external frames. In order to enact an extraction protocol in a generally covariant theory, dependence on the external reference frame must first be removed and the states made relational. Here, we examine the implementation of an extraction protocol for Gaussian states, whose center of mass and relational degrees of freedom are entangled, in a relational toy model where translation invariance stands in for full diffeomorphism invariance. Constructing fully relational states and the corresponding extraction/localization can, in principle, be done in two ways. External frame position information can be removed through 𝐺-twirling over translations or one can spontaneously break the translation symmetry via the gradient of an auxiliary field, or 𝑍-model. We determine the energetics of quantum information localization after the states have been made fully relational via both the 𝐺-twirl and 𝑍-model. We also show one can obtain the 𝐺-twirl construction from a 𝑍-model as a limit of positive operator valued measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nondegenerate Noise-Resilient Superconducting Qubit

We propose a superconducting qubit based on engineering the first and second harmonics of the Josephson energy and phase relation 𝐸 𝐽⁢1 ⁢cos 𝜑 and 𝐸 𝐽⁢2 ⁢cos 2⁢𝜑. By constructing a circuit such that 𝐸 𝐽⁢2 is negative and |𝐸 𝐽⁢1 | ≪ |𝐸 𝐽⁢2 |, we create a periodic potential with two nondegenerate minima. The qubit, which we dub “harmonium,” is formed from the lowest-energy states of each minimum. Bit-flip protection of the qubit arises due to the localization of each qubit state to their respective minima, while phase-flip protection can be understood by considering the circuit within the Born-Oppenheimer approximation. We demonstrate with time-domain simulations that single- and two-qubit gates can be performed in approximately 100 ns. Finally, we compute the qubit coherence times using numerical diagonalization of the complete circuit in conjunction with state-of-the-art noise models. We estimate out-of-manifold heating times on the order of milliseconds, which can be treated as erasure errors using conventional dispersive readout. We estimate pure-dephasing times on the order of many tens of milliseconds, and bit-flip times on the order of seconds.

quantum error correction↗
DOE OSTI2024

Learning from metastable symmetric-tilt grain boundaries using physics-based descriptors

Grain boundaries (GBs) govern critical properties of polycrystalline materials. Although significant advancements have been made in characterizing minimum energy and ordered GBs, real GBs are seldom found in such well-defined states. This diversity of atomic arrangements in metastable states makes it challenging to establish structure-property relationships with physical insights. Here, to address this challenge, we use data-driven methods to explore these relationships and examine the underlying physics. In this study, we utilize a large atomistic database (~5000) of minimum energy and metastable states of symmetric-tilt copper GBs, combined with physically motivated local atomic environment (LAE) descriptors [strain functional descriptors (SFDs)], to predict GB properties and gain physical insights. Our regression models exhibit robust predictive capabilities using only 19 descriptors, generalizing to atomic environments in nanocrystals. A significant highlight of our work is the integration of an unsupervised method with SFDs to elucidate LAEs at GBs and their role in determining properties. The model, trained on these minimum energy and metastable GBs using SFDs, predicts the properties of unseen nanocrystals with good accuracy. Our research underscores the role of a physics-based representation of LAEs and the efficacy of data-driven methods in establishing GB structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Cage Reports on Its Contents: Spectroscopic Signatures of Cryo-Cooled K + - and Ba 2+ -Benzocryptand Complexes

UV photofragment spectroscopy and IR–UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K + or Ba 2+ imbedded inside it (labeled as K + -BzCrypt, Ba 2+ -BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH 2 -CH 2 -O-CH 2 -CH 2 -O-CH 2 -CH 2 -N chains provides a UV chromophore whose S 0 –S 1 transition we probe. K + -BzCrypt and Ba 2+ -BzCrypt have their S 0 –S 1 origin transitions at 35,925 and 36,446 cm –1 , respectively, blue-shifted by 174 and 695 cm –1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR–UV dip spectroscopy. The alkyl CH stretch region (2800–3000 cm –1 ) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K + -BzCrypt and Ba 2+ -BzCrypt, identifying two conformations (A and B) that lie below all others in energy. Here, we extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH 2 groups adjacent to these atoms. The localized CH 2 scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH 2 wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm –1 . As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH 2 symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗