Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Formulation, Implementation and Validation of a 1D Boundary Layer Inflow Scheme for the QUIC Modeling System

Recent studies have highlighted the importance of accurate meteorological conditions for urban transport and dispersion calculations. In this work, we present a novel scheme to compute the meteorological input in the Quick Urban & Industrial Complex () diagnostic urban wind solver to improve the characterization of upstream wind veer and shear in the Atmospheric Boundary Layer (ABL). The new formulation is based on a coupled set of Ordinary Differential Equations (ODEs) derived from the Reynolds Averaged Navier–Stokes (RANS) equations, and is fast to compute. Building upon recent progress in modeling the idealized ABL, we include effects from surface roughness, turbulent stress, Coriolis force, buoyancy and baroclinicity. We verify the performance of the new scheme with canonical Large Eddy Simulation (LES) tests with the GPU-accelerated FastEddy"Equation missing" solver in neutral, stable, unstable and baroclinic conditions with different surface roughness. Furthermore, we evaluate QUIC calculations with and without the new inflow scheme with real data from the Urban Threat Dispersion (UTD) field experiment, which includes Lidar-based wind measurements as well as concentration observations from multiple outdoor releases of a non-reactive tracer in downtown New York City. Compared to previous inflow capabilities that were limited to a constant wind direction with height, we show that the new scheme can model wind veer in the ABL and enhance the prediction of the surface cross-isobaric angle, improving evaluation statistics of simulated concentrations paired in time and space with UTD measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Flavoured warped axion

We consider a 5D extension of the DFSZ axion model that addresses both the axion quality and fermion mass hierarchy problems, and predicts flavour-dependent, off-diagonal axion-fermion couplings. The axion is part of a 5D complex scalar field charged under a U(1)PQ symmetry that is spontaneously broken in the bulk, and is insensitive to explicit PQ breaking on the UV boundary. Bulk Standard Model fermions interact with two Higgs doublets that can be localized on the UV boundary or propagate in the bulk to explain the fermion masses and mixings. When the Higgs doublets are localized on the UV boundary, they induce flavour diagonal couplings between the fermions and the axion. However, when the Higgs doublets propagate in the bulk, the overlap of the axion and fermion profiles generates flavour off-diagonal couplings. The effective scale of these off-diagonal couplings in both the quark and lepton sectors can be as small as 10 11 GeV, and therefore will be probed in future precision flavour experiments.

Global Symmetries↗
DOE OSTI2022

High Entropy Oxide Relaxor Ferroelectrics

Relaxor ferroelectrics are important in technological applications due to strong electromechanical response, energy storage capacity, electrocaloric effect, and pyroelectric energy conversion properties. Current efforts to discover and design materials in this class generally rely on substitutional doping as slight changes to local compositional order can significantly affect the Curie temperature, morphotropic phase boundary, and electromechanical responses. In this work, we demonstrate that moving to the strong limit of compositional complexity in an ABO 3 perovskite allows stabilization of relaxor responses that do not rely on a single narrow phase transition region. Entropy-assisted synthesis approaches are utilized to synthesize single-crystal Ba(Ti 0.2 Sn 0.2 Zr 0.2 Hf 0.2 Nb 0.2 )O 3 [Ba(5B)O] films. The high levels of configurational disorder present in this system are found to influence dielectric relaxation, phase transitions, nanopolar domain formation, and Curie temperature. Temperature-dependent dielectric, Raman spectroscopy, and second-harmonic generation measurements reveal multiple phase transitions, a high Curie temperature of 570 K, and the relaxor ferroelectric nature of Ba(5B)O films. The first-principles theory calculations are used to predict possible combinations of cations to design relaxor ferroelectrics and quantify the relative feasibility of synthesizing these highly disordered single-phase perovskite systems. Further, the ability to stabilize single-phase perovskites with various cations on the B-sites offers possibilities for designing high-performance relaxor ferroelectric materials for piezoelectric, pyroelectric, and electrocaloric applications.

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

Deep learning of electrochemical CO 2 conversion literature reveals research trends and directions

Large-scale and openly available material science databases are mainly composed of computer simulation results rather than experimental data. Some examples include the Materials Project, Open Quantum Materials Database, and Open Catalyst 2022. Unfortunately, building large-scale experimental databases remains challenging due to the difficulties in consolidating locally distributed datasets. In this work, focusing on the catalysis literature of CO 2 reduction reactions (CO 2 RRs), we present a machine learning (ML)-based protocol for selecting highly relevant papers and extracting important experimental data. First, we report a document embedding method (Doc2Vec) for collecting papers of greatest relevance to the specific target domain, which yielded 3154 CO 2 RR-related papers from six publishers. Next, we developed named entity recognition (NER) models to extract twelve entities related to material names (catalyst, electrolyte, etc.) and catalytic performance (Faradaic efficiency, current density, etc.). Further, among several tested models, the MatBERT-based approach achieved the highest accuracy, with an average F1-score of 90.4% and an F1-score of 95.2% in a boundary relaxation evaluation scheme. The accurate and accelerated NER-based data extraction from a large volume of catalysis literature enables temporal trend analyses of the CO 2 RR catalysts, products, and performances, revealing the potentially effective material space in CO 2 RRs. While this work demonstrates the effectiveness of our ML-based text mining methods for specifically CO 2 RR literature, the methods and approach are applicable to and may be used to accelerate the development of other catalytic chemical reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Laser cooling limits in fields with a polarisation gradient of atoms with different recoil energies

Based on the numerical solution of the quantum kinetic equation for the atomic density matrix, which makes it possible to accurately take into account the recoil effects in the interaction of atoms with field photons, we have studied the limits of laser cooling of atoms using closed optical transitions characterised by different recoil parameters (the ratio of the recoil energy to the natural linewidth). It is shown that for optical transitions with an insufficiently small recoil parameter, the polarisation effects, which lead to the possibility of sub-Doppler laser cooling, lose their efficiency and do not ensure an attainment of the temperature below the Doppler limit. The analysis performed allows one to outline the boundaries of the sub-Doppler theory of laser cooling of atoms. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Understanding aerosol properties in convective outflows during TRACER

Convective outflow boundaries commonly form across the southeastern U.S., often originating from storms initiated along sea-breeze fronts. As outflows spread, merge, and collide, they frequently trigger secondary convection in air masses already modified by primary storm outflow. This resulting air mass carries a distinctive aerosol population, which can influence the development of secondary convective cells. During the U.S. DOE TRACER field campaign in summer 2022, we investigated how outflows affect aerosol size, composition, and hygroscopicity across Greater Houston, Texas. Using a novel detection algorithm, we identified and verified 72 outflow boundaries from surface meteorological observations. In the relatively simple continental environment northwest of Houston, a common relationship emerged between CCN-inferred hygroscopicity (κ) and size-resolved aerosol concentrations. The κ of aerosols typically decreased in size ranges where aerosol concentrations increased following outflow passage and increased where concentrations decreased. In contrast, at the Atmospheric Radiation Measurement’s (ARM) fixed site in La Porte, Texas, the aerosol response depended on outflow direction, shaped by concentrated industrial and shipping activity to the north and south/southeast along the Houston Ship Channel. Concurrent size distribution and composition changes suggested case-specific variability in κ, with either reinforcing or offsetting effects. Outflows passing over industrial/shipping sectors often promoted less hygroscopic aerosol populations while cleaner residential outflows often favored more hygroscopic aerosol populations. These results highlight key aerosol changes following outflow passage, providing a foundation for improved understanding of aerosol-cloud interactions related to secondary convection initiation along or in the wake of outflow boundaries.

Thompson, Seth A. [Texas A & M Univ., College Stat↗
DOE OSTI2024

Simulations of tokamak edge plasma turbulent fluctuations based on a minimal 3D model

In this study, a new simulation model for tokamak boundary plasma, SOLT3D, is implemented in the BOUT++ framework. The simulation model includes a set of dynamic equations describing collisional boundary plasma and neutral gas in the tokamak scrape-off layer and divertor region. The model is verified against standard linear plasma instabilities and available nonlinear results. For L-mode like conditions, SOLT3D reproduces characteristics of boundary plasma turbulent fluctuations that are similar to published experimental data, in terms of the amplitude and spatial dependence of the fluctuations. It also reproduces realistic plasma fluxes on material surfaces and expected Bohm-like effective radial transport. Plasma fluctuations inferred from the simulations imply inevitably a significant level of intrinsic 'noise' for axisymmetric tokamak plasma transport modeling, introducing errors on the order of unity. In particular, the toroidally averaged atomic rates below 5–10 eV are strongly modified by turbulent plasma fluctuations, which should significantly affect the standard axisymmetric modeling of the tokamak edge plasma and divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of finite-element models using full-field experimental data: Levelling finite-element analysis data through a digital image correlation engine

Full-field data from digital image correlation (DIC) provide rich information for finite-element analysis (FEA) validation. However, there are several inherent inconsistencies between FEA and DIC data that must be rectified before meaningful, quantitative comparisons can be made, including strain formulations, coordinate systems, data locations, strain calculation algorithms, spatial resolutions and data filtering. As such, in this paper, we investigate two full-field validation approaches: (1) the direct interpolation approach, which addresses the first three inconsistencies by interpolating the quantity of interest from one mesh to the other, and (2) the proposed DIC-levelling approach, which addresses all six inconsistencies simultaneously by processing the FEA data through a stereo-DIC simulator to ‘level’ the FEA data to the DIC data in a regularisation sense. Synthetic ‘experimental’ DIC data were generated based on a reference FEA of an exemplar test specimen. The direct interpolation approach was applied, and significant strain errorswere computed, even though therewas no model form error, because the filtering effect of theDIC enginewas neglected. In contrast, the levelling approach provided accurate validation results, with no strain error when no model form error was present. Next, model form error was purposefully introduced via a mismatch of boundary conditions. With the direct interpolation approach, the mismatch in boundary conditions was completely obfuscated, while with the levelling approach, it was clearly observed. Finally, the ‘experimental’ DIC datawere purposefully misaligned slightly fromthe FEA data. Both validation techniques suffered from the misalignment, thus motivating continued efforts to develop a robust alignment process. In conclusion, direct interpolation is insufficient, and the proposed levelling approach is required to ensure that the FEA and the DIC data have the same spatial resolution and data filtering. Only after the FEA data have been ‘levelled’ to the DIC data can meaningful, quantitative error maps be computed.

42 ENGINEERING↗
DOE OSTI2021

Simulating Boundary Region Plasma Instabilities

The boundary regions of low-temperature plasmas are known to be susceptible to kinetic instabilities, which can affect the energies and fluxes of particles directed at the material boundary. For example, both the ion acoustic instability as well as an instability near the electron plasma frequency have been observed. Particle-in-cell (PIC) simulation is a tool that, alongside experiments, can capture the effects these instabilities have on the particle distribution functions. Ultimately, simulations can determine under what conditions these effects are significant by comparing to theoretical predictions and explore conditions unamenable to experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗
DOE OSTI2020

Dirichlet absorbing boundary conditions for classical and peridynamic diffusion-type models

Diffusion-type problems in (nearly) unbounded domains play important roles in various fields of fluid dynamics, biology, and materials science. The aim of this paper is to construct accurate absorbing boundary conditions (ABCs) suitable for classical (local) as well as nonlocal peridynamic (PD) diffusion models. The main focus of the present study is on the PD diffusion formulation. The majority of the PD diffusion models proposed so far are applied to bounded domains only. In this study, we propose an effective way to handle unbounded domains both with PD and classical diffusion models. For the former, we employ a meshfree discretization, whereas for the latter the finite element method (FEM) is employed. The proposed ABCs are time-dependent and Dirichlet-type, making the approach easy to implement in the available models. The performance of the approach, in terms of accuracy and stability, is illustrated by numerical examples in 1D, 2D, and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A comparative study of two-way and offline coupled WRF v3.4 and CMAQ v5.0.2 over the contiguous US: performance evaluation and impacts of chemistry–meteorology feedbacks on air quality

The two-way coupled Weather Research and Forecasting and Community Multiscale Air Quality (WRF-CMAQ) model has been developed to more realistically represent the atmosphere by accounting for complex chemistry–meteorology feedbacks. In this study, we present a comparative analysis of two-way (with consideration of both aerosol direct and indirect effects) and offline coupled WRF v3.4 and CMAQ v5.0.2 over the contiguous US. Long-term (5 years from 2008 to 2012) simulations using WRF-CMAQ with both offline and two-way coupling modes are carried out with anthropogenic emissions based on multiple years of the U.S. National Emission Inventory and chemical initial and boundary conditions derived from an advanced Earth system model (i.e., a modified version of the Community Earth System Model/Community Atmospheric Model). The comprehensive model evaluations show that both two-way WRF-CMAQ and WRF-only simulations perform well for major meteorological variables such as temperature at 2 m, relative humidity at 2 m, wind speed at 10 m, precipitation (except for against the National Climatic Data Center data), and shortwave and longwave radiation. Both two-way and offline CMAQ also show good performance for ozone (O3) and fine particulate matter (PM 2.5 ). Due to the consideration of aerosol direct and indirect effects, two-way WRF-CMAQ shows improved performance over offline coupled WRF and CMAQ in terms of spatiotemporal distributions and statistics, especially for radiation, cloud forcing, O 3 , sulfate, nitrate, ammonium, elemental carbon, tropospheric O 3 residual, and column nitrogen dioxide (NO 2 ). For example, the mean biases have been reduced by more than 10 W m –2 for shortwave radiation and cloud radiative forcing and by more than 2 ppb for max 8 h O 3 . However, relatively large biases still exist for cloud predictions, some PM 2.5 species, and PM 10 that warrant follow-up studies to better understand those issues. The impacts of chemistry–meteorological feedbacks are found to play important roles in affecting regional air quality in the US by reducing domain-average concentrations of carbon monoxide (CO), O 3 , nitrogen oxide (NO x ), volatile organic compounds (VOCs), and PM 2.5 by 3.1 % (up to 27.8 %), 4.2 % (up to 16.2 %), 6.6 % (up to 50.9 %), 5.8 % (up to 46.6 %), and 8.6 % (up to 49.1 %), respectively, mainly due to reduced radiation, temperature, and wind speed. The overall performance of the two-way coupled WRF-CMAQ model achieved in this work is generally good or satisfactory and the improved performance for two-way coupled WRF-CMAQ should be considered along with other factors in developing future model applications to inform policy making.

58 GEOSCIENCES↗
DOE OSTI2021

Anisotropic elastoplastic response of double-porosity media

Here we present a continuum framework to simulate fluid flow through anisotropic elastoplastic media with double porosity. Two effective stress measures σ' and δ" emerge from the thermodynamic formulation, which are energy-conjugate to the elastic and plastic components of strain, respectively. Both effective stress measures can be expressed as a combination of the total Cauchy stress σ and the average pore pressure $\bar{p}$ in the two pore scales. In the effective stress for elasticity, $\bar{p}$ is scaled with a rank-2 Biot tensor, whereas the effective stress for plasticity follows the Terzaghi form in which $\bar{p}$ is scaled by the Kronecker delta. The Biot tensor and storage coefficients are derived as functions of elasticity parameters and porosities. A mixed finite element formulation is introduced to discretize the domain and solve initial boundary-value problems. A stabilization scheme is employed on equal-order interpolation for both displacement and pressure fields throughout the entire range of drainage responses. Numerical simulations reproduce the hydromechanical response of Opalinus shale in one-dimensional consolidation tests throughout the range of primary and secondary consolidation under different external loads. Numerical simulations of the consolidation of a rectangular domain subjected to a strip load demonstrate the efficacy of the proposed stabilization scheme, as well as illustrate the impacts of stress history, mass transfer, and different pore systems on the hydromechanical response.

42 ENGINEERING↗
DOE OSTI2025

AI‐Driven Defect Engineering for Advanced Thermoelectric Materials

Thermoelectric materials offer a promising pathway to directly convert waste heat to electricity. However, achieving high performance remains challenging due to intrinsic trade-offs between electrical conductivity, the Seebeck coefficient, and thermal conductivity, which are further complicated by the presence of defects. This review explores how artificial intelligence (AI) and machine learning (ML) are transforming thermoelectric materials design. Advanced ML approaches including deep neural networks, graph-based models, and transformer architectures, integrated with high-throughput simulations and growing databases, effectively capture structure-property relationships in a complex multiscale defect space and overcome the “curse of dimensionality”. This review discusses AI-enhanced defect engineering strategies such as composition optimization, entropy and dislocation engineering, and grain boundary design, along with emerging inverse design techniques for generating materials with targeted properties. Finally, it outlines future opportunities in novel physics mechanisms and sustainability, highlighting the critical role of AI in accelerating the discovery of thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

RNAseq-based transcriptome assembly of Clostridium acetobutylicum for functional genome annotation and discovery

Accurate genome annotations are essential in modern biology and biotechnology, yet they are still largely based on genome sequencing and comparative analyses. We show that the Clostridium acetobutylicum genome annotation can be markedly improved by integrating bioinformatic predictions with RNA sequencing (RNAseq) data. Samples were acquired under butanol, butyrate, and unstressed treatments across various growth conditions. Analysis of an initial assembly revealed errors due to background signals and limitations of assembly algorithms. Hurdles for RNAseq transcriptome mapping include optimizing library complexity and sequencing depth, yet most studies report low sequencing depth and ignore the effect of ribosomal RNA abundance. An integrative analysis was developed to combine motif predictions, single-nucleotide resolution sequencing depth, and library complexity to resolve difficulties in assembly curation. This minimized false positive error and determined gene boundaries, in some cases, to the exact base-pair of prior studies. This will be the first strand-specific transcriptome assembly in a Clostridium organism.

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

The repulsive Casimir-type forces in binary mixture of weakly interacting Bose–Einstein condensates gases

Highlights: • The equations of the state of a binary mixture weakly interacting Bose gases are established based on the Lagrangian of the system. • The formulae of the Casimir-type energy and the resulting Casimir-type force are produced from the quantum part of the grand canonical energy. • The results for the Casimir-type effect in the one-field case can be re-obtained by using the results of the two-field case. Within the framework of the quantum field theory (QFT) in the one-loop approximation, we investigate the Casimir-type effect in a mixture of two-component weakly interacting Bose-Einstein condensates (BECs) gases confined between two parallel plates at zero temperature. The Casimir-type forces in the miscible state of BECs for different boundary conditions (BCs) are discovered and compared to each other. The repulsive nature of the Casimir-type force in BECs is defined for the first time. Our calculations are finished in both grand canonical ensemble (GCE) and canonical ensemble (CE).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Corrosion of 316L & 316H stainless steel in molten LiF-NaF-KF (FLiNaK)

The corrosion behavior of 316L and 316H stainless steels (low and high carbon, respectively) in molten LiF-NaF-KF (FLiNaK) salt at 700 °C for 500 h has been investigated. Tests were performed in 316L stainless steel, pyrolytic boron nitride, and glassy carbon capsule materials. The capsule material had a notable effect on the corrosion behavior of the two stainless steels. Further, the modes of corrosion for the two types of stainless steel in various capsule materials were classified into three categories, namely, uniform surface recession, selective grain boundary Cr dissolution, and the dissolution of Cr in the bulk grains. By most metrics, the 316H stainless steel exhibited slightly superior corrosion resistance, while by others the two stainless steels were comparable. These differences in corrosion behavior are attributed to the differences in grain size and the degree of sensitization.

36 MATERIALS SCIENCE↗