Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forcing fields”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Chemical Compensation Challenges in Processing Antiferroelectric PbZrO 3 Thin Films

The electric field-induced antipolar-to-polar phase transition in antiferroelectric materials is accompanied by large changes in the structural and functional response, making them attractive for many applications ranging from pulsed energy capacitors to high-strain-high (blocking-)force actuators, solid-state heating and cooling systems, and optoelectronic devices. PbZrO 3 , as an end member of the PZT, (x)PbZr 1–x Ti x O 3 , solid solution, has been the subject of many studies. However, processing of perovskite PbZrO 3 is extremely challenging due to phase instabilities induced by Pb loss at the temperatures necessary for perovskite crystallization. Here, we discuss the challenges associated with Pb loss, as well as its compensation through bulk Pb overstoichiometry and interfacial PbO additions, in chemical solution processed, highly oriented PbZrO 3 thin films. For both 042 o - and 001 o -oriented (o-orthorhombic distortion) PbZrO 3 thin films, the crystallization interfaces are the most important contributors to Pb loss and off-stoichiometric outcomes, and no single approach was sufficient to address these challenges. Pb-rich and Pb-deficient nonperovskite phases including ZrO x were only observable through microscopic characterization, and X-ray diffraction alone could not rule out the presence of such secondary phases. “Ideal” macroscopic antiferroelectric polarization-electric field double hysteresis curves were observed despite (at times) the large presence of secondary phases. However, reduced saturation polarization and increased phase transition electric fields were observed concomitantly with the presence of secondary phase(s) and tentatively assigned to the voltage drop across the ZrO x nanocrystals. Based on these results, it is imperative that the presence of secondary phases always be addressed (beyond X-ray diffraction spectra) in order to correctly evaluate the properties of antiferroelectric films.

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

Multimodal characterization of Te inclusions in Cd 1−x Zn x Te and Cd 1−x Zn x Te 1−y Se y for gamma and X-ray detectors

While CdZnTe (CZT) and CdZnTeSe (CZTS) semiconductors have emerged as compounds for room-temperature gamma and X-ray detection materials, they continue to be constrained by the formation of Te-inclusion defects generated during the growth and post-growth phases of the material, which adversely affect the detector performance. We demonstrate the utility of multimodal microscopic imaging and analysis for the characterization of the optical and electronic properties of Te inclusions in CZT and CZTS crystals at both micron and nanometer length scales. Having first identified regions with micron-scale Te inclusions using confocal Raman microscopy techniques, optically coupled infrared scattering near-field optical microscopic mapping was performed to map the distribution of these inclusions with nanometer spatial resolution and correlate the presence of Te inclusions in the matrix with other properties. Kelvin probe force microscopy was then utilized to characterize the variations of the work function associated with the presence of Te inclusions. Here, we observe an increase of ~ 240 mV in the work function associated with Te inclusions compared to the bulk CZT/CZTS crystals. Additionally, we observe that individual bulk grains in CZT can exhibit slight potential variations. Our findings develop a portrait of the charge trapping mechanisms in CZT and CZTS that act to degrade detector performance, while the demonstration of these combined microscopy techniques provides a new analytical tool that can be utilized for further optimization of the detector performance for these semiconducting compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Photoinduced patterning of oxygen vacancies to promote the ferroelectric phase of Hf 0.5 Zr 0.5 O 2

Photoinduced reductions in the oxygen vacancy concentration were leveraged to increase the ferroelectric phase fraction of Hf 0.5 Zr 0.5 O 2 thinfilms. Modest ( ~ 2 — 77 pJ=cm 2 ) laser doses of visible light (488 nm, 2.54 eV) spatially patterned the concentration of oxygen vacancies as monitored by photoluminescence imaging. Local, tip-based, near-field, nanoFTIR measurements showed that the photoinduced oxygen vacancy concentration reduction promoted formation of the ferroelectric phase (space group Pca2 1 ), resulting in an increase in the piezoelectric response measured by piezoresponse force microscopy. Photoinduced vacancy tailoring provides, therefore, a spatially prescriptive, postsynthesis, and low-entry method to modify phase in HfO 2 -based materials.

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

Bolometer tomography on Wendelstein 7-X for study of radiation asymmetry

The algorithm for bolometer tomography at Wendelstein 7-X (W7-X) has been recently improved using a novel regularization functional, based on relative gradient smoothing (RGS) of the sought radiation profile. It has been validated using radiation patterns provided by 3D modeling under real plasma conditions as phantoms and then applied to bolometer measurements performed during the first divertor operation phase of W7-X. The following results are presented: 1) edge-localized 2D radiation patterns with clearly resolved magnetic island radiation structures, 2) an up-down asymmetry in the impurity radiation that is not captured by the 3D edge plasma transport modeling, 3) reversal of the asymmetry with reversed magnetic field direction. Further analysis reveals a poloidal variation of the emissivity in the outer confined plasma region with a field-direction dependent asymmetry, also supported by the soft X-ray measurements. This asymmetry is considered to be related to asymmetric impurity distributions, driven by the pronounced ion-impurity friction force at the plasma edge where the collisionality of the W7-X plasma is sufficiently high to develop impurity asymmetry as predicted by neoclassical theory of parallel impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Wafer-scale waveguide sidewall roughness scattering loss characterization by image processing

Photonic integrated circuits (PICs) are vital for developing affordable, high-performance optoelectronic devices that can be manufactured at an industrial scale, driving innovation and efficiency in various applications. Optical loss of modes in thin film waveguides and devices is a critical measure of their performance. Thin film growth, lithography, masking, and etching processes are imperfect processes that introduce significant sidewall and top-surface roughness and cause dominating optical losses in waveguides and photonic structures. This roughness, as perturbations couple light from guided to far-field radiation modes, leads to scattering losses that can be estimated from theoretical models. Typically, with UV-based lithography, sidewall roughness is significantly larger than wafer-top surface roughness. Atomic force microscopy (AFM) imaging measurement gives a 3D and high-resolution roughness profile, but the measurement is inconvenient, costly, and unscalable for large-scale PICs and at wafer-scale. Here, we evaluate the sidewall roughness profile based on 2D high-resolution scanning electron microscope (SEM) imaging. We characterized the loss on two homemade nitride and oxide films on 3-inch silicon wafers with 12 waveguide devices on each and correlated the scattering loss estimated from a 2D image-based sidewall profile and theoretical Payne model. The lowest loss of guided fundamental transverse electric (TE 0 ) mode is found at 0.075 dB/cm at 633 nm across 24 devices, a record at visible wavelength. Our work shows 100% success (edge continuity span exceeding 95% of image width/height) in edge detection in image processing of all images to estimate autocorrelation function and optical mode loss. These demonstrations offer valuable insights into waveguide sidewall roughness and a comparison of experimental and 2D SEM image processing based loss estimations with applications in loss characterization at wafer-scale PICs.

Khurana, Mohit (ORCID:0000000207898672)↗
DOE OSTI2023

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Wind and Structural Loads on Parabolic Trough Solar Collectors at Nevada Solar One

Wind loading is a main contributor to structural design costs of Concentrating Solar Power (CSP) collectors, such as heliostats and parabolic troughs. These structures must resist the mechanical forces generated by turbulent wind. At the same time, the reflector surfaces must exhibit the necessary rigidity to maintain their optimal optical performance in windy conditions. Over two years, NREL conducted comprehensive field measurements of the atmospheric turbulent wind conditions and the resulting structural wind loads on parabolic troughs at the Nevada Solar One (NSO) plant. The measurement set-up included meteorological masts and structural load sensors on four trough rows. Additionally, we commissioned a lidar scanning the horizontal plane over the trough field. This data set catalogs the high-resolution data set characterizing the complex flow field and resulting structural loads on parabolic trough collectors. By providing this first-of-its-kind data set to the CSP community, we aim to enhance the community's understanding of wind-loading experienced by CSP collector structures. This data set will also help design next-generation solar collectors and photovoltaic trackers.

Array↗
DOE OSTI2024

Field demonstration of a tracer method to track simulated exhaled air trajectories and mixing in three connected rooms with upper-room GUV

We conducted tracer gas experiments in three connected nursery rooms to track simulated exhaled air trajectories and mixing. We emulated exhaled air with a pulse release of ethanol and measured its concentration with 2s time resolution with spatially distributed metal oxide sensors in both the upper and occupied levels of the room. We found that the overhead cooling supply air enhanced vertical air mixing within the rooms. When a room had higher supply airflow than the others, the tracer gas mixed quickly in that room because the supply air dominated the room mixing, thereby isolating it from air mixing with the other rooms. When the forced air system was not operating, the tracer gas resided longer in the upper room before descending. The tracer method also depicts air trajectories from different release locations by detecting the release point and affecting nearby sensors. These experiments support the potential of this method to be used in the field for understanding air trajectories in rooms with upper-room GUV.

Um, Chai Yoon↗
DOE OSTI2026

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗
DOE OSTI2024

Nuclear-matter saturation and symmetry energy within Δ -full chiral effective field theory

Nuclear saturation and the symmetry energy are key properties of low-energy nuclear physics that depend on fine details of the nuclear interaction. The equation of state around saturation is also an important anchor for extrapolations to higher densities and studies of neutron stars. Here we develop a unified statistical framework that uses realistic nuclear forces to link the theoretical modeling of finite nuclei and infinite nuclear matter. We construct fast and accurate emulators for nuclear-matter observables and employ an iterative history-matching approach to explore and reduce the enormous parameter domain of Δ -full chiral interactions. We perform rigorous uncertainty quantification and find that model calibration including O 16 observables gives saturation predictions that are more precise than those that only use few-body data. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Microstructure-Based Degradation Modeling in Solid Oxide Cells

The presentation is intended for a graduate class in Georgia Southern University taught by Dr. Hayri Sezer on electrochemistry, SOFC and Li ion battery. It summaries NETL's recent works on microstructure-based modeling of degradation in solid oxide cells. It covers the degradation mechanism of Ni coarsening, Ni migration and the build-up of oxygen partial pressure inside the electrolyte. The phase-field models developed for Ni coarsening and migration are reviewed and the results are compared to experiments in literature to examine different diffusion mechanisms (Ni self-diffusion, Ni(OH)2 diffusion in pores and Ni(OH) surface diffusion) and driving forces (Ni(OH)x concentration and Ni-YSZ contact angle change). It also demonstrates the developed model for oxygen partial pressure and shows that the predicted oxygen partial pressure can explain the Ni oxidation in the hydrogen electrode and crack growth on the electrolyte-oxygen electrode interface under certain operating conditions.

nickel coarsening↗
DOE OSTI2021

Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows

The guiding center and gyrokinetic theory of magnetized particle motion is extended to the regime of large electric field gradients perpendicular to the magnetic field. A gradient in the electric field directly modifies the oscillation frequency and causes the Larmor orbits to deform from circular to elliptical trajectories. In order to retain a good adiabatic invariant, there can only be strong dependence on a single coordinate at lowest order, so that resonances do not generate chaotic motion that destroys the invariant. When the gradient across magnetic flux surfaces is dominant, the guiding center drift velocity becomes anisotropic in response to external forces and additional curvature drifts must be included. The electric polarization density remains gyrotropic, but both the polarization and magnetization are modified by the change in gyrofrequency. The theory can be applied to shear flows that are even stronger than those observed in the edge transport barrier of a high-performance tokamak (H-mode) pedestal, even if the toroidal field is as small as or even smaller than the poloidal field. Yet, the theory retains a mathematical form that is similar to the standard case and can readily be implemented within existing simulation tools.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Database and deep-learning scalability of anharmonic phonon properties by automated brute-force first-principles calculations

Understanding the anharmonic phonon properties of crystal compounds—such as phonon lifetimes and thermal conductivities—is essential for investigating and optimizing their thermal transport behaviors. These properties also impact optical, electronic, and magnetic characteristics through interactions between phonons and other quasiparticles and fields. In this study, we develop an automated first-principles workflow to calculate anharmonic phonon properties and build a comprehensive database encompassing more than 6500 inorganic compounds. Utilizing this dataset, we train a graph neural network model to predict thermal conductivity values and spectra from structural parameters, demonstrating a scaling law in which prediction accuracy improves with increasing training data size. High-throughput screening with the model enables the identification of materials exhibiting extreme thermal conductivities—both high and low. The resulting database offers valuable insights into the anharmonic behavior of phonons, thereby accelerating the design and development of advanced functional materials.

Ohnishi, Masato [University of Tokyo (Japan); Inst↗
DOE OSTI2025

Phenome‐to‐genome insights for evaluating root system architecture in field studies of maize

Abstract Understanding the genetic basis of root system architecture (RSA) in crops requires innovative approaches that enable both high‐throughput and precise phenotyping in field conditions. In this study, we evaluated multiple phenotyping and analytical frameworks for quantifying RSA in mature, field‐grown maize in three field experiments. We used forward and reverse genetic approaches to evaluate >1700 maize root crowns, including a diversity panel, a biparental mapping population, and maize mutant and wild‐type alleles at two known RSA genes,DEEPER ROOTING 1(DRO1) andRootless1(Rt1). We show the utility of increasing the dimensionality of traditional two‐dimensional (2D) techniques, referred to as the “2D multi‐view” method, to improve the capture of whole root system information for mapping genetic variation influencing RSA. Comparison of univariate and multivariate genome‐wide association study (GWAS) approaches revealed that multivariate traits were effective at dissecting complex RSA phenotypes and identifying pleiotropic quantitative trait loci (QTLs). Overall, three‐dimensional (3D) root models generated from X‐ray computed tomography and digital phenotyping captured a larger proportion of RSA trait variations compared to other methods of root phenotyping, as evidenced by both genome‐wide and single‐gene analyses. Among the individual root traits, root pulling force emerged as a highly heritable estimate of RSA that identified the largest number of shared QTLs with 3D phenotypes. Our study shows that integrating complementary phenotyping technologies helps to provide a more comprehensive understanding of the genetic architecture of RSA in field‐grown maize.

Genetics & Heredity↗
DOE OSTI2020

Deep learning interfacial momentum closures in coarse-mesh CFD two-phase flow simulation using validation data

Multiphase flow phenomena have been widely observed in the industrial applications while it remains a challenging yet unsolved problems. Three-dimensional computational fluid dynamics (CFD) approaches resolve the flow fields on a finer special and temporal scales which can complement the dedicated experimental study. However, closures have to be introduced to reflect the underlying physics in multiphase flow. Among them, the interfacial forces, including drag, lift, turbulent dispersion and wall lubrication forces, play in important role on the bubble’s distribution and migration in liquid-vapor two-phase flow. Development of those closures traditionally rely on the experimental data and analytical derivation with simplified assumptions which usually cannot deliver a universal solution across wide range of flow conditions. In this paper, a data-driven approach, named as Feature Similarity Measurement (FSM), is developed and applied to improve the simulation capability of two-phase flow with coarse-mesh CFD approach. Interfacial momentum transfer in adiabatic bubbly flow serves as the focus of the present study. Both a mature and a simplified set of interfacial closures are taken as the low fidelity data. Experimental data and fine mesh CFD simulations results are adopted as high-fidelity data. Qualitative and quantitative analysis are performed in this paper which reveals that FSM can substantially improve the prediction of coarse mesh CFD model regardless of the choice of interfacial closures and it provides scalability and consistency across discontinuous flow regimes. Furthermore, it demonstrates that data-driven method can aid the multiphase flow modeling by exploring the connections between local physical features and simulation errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Mean-field dynamo as a quantum-like modulational instability

Presented here is a novel formulation of the mean-field dynamo as a modulational instability of magnetohydrodynamic (MHD) turbulence. This formulation, termed mean-field wave kinetics (MFWK), is based on the Weyl symbol calculus and allows describing the interaction between the mean fields (magnetic field and fluid velocity) and turbulence without requiring scale separation that is commonly assumed in the literature. The turbulence is described by the Wigner–Moyal equation for the spectrum of the two-point correlation matrix (Wigner matrix) of magnetic-field and velocity fluctuations and depicts the turbulence as an effective plasma of quantum-like particles that interact via the mean fields. Eddy–eddy interactions, which serve as ‘collisions’ in this effective plasma, are modelled within the standard minimal tau approximation to aid comparison with existing theories. Using MFWK, the non-local electromotive force is calculated for generic turbulence from first principles, modulo the limitations of MFWK. This result is then used to study, both analytically and numerically, the modulational modes of MHD turbulence, which appear as linear instabilities of the said effective quantum-like plasma of fluctuations. The standard α 2 -dynamo and other known results are reproduced as special cases. A new dynamo effect is predicted that is driven by correlations between the turbulent flow velocity and the turbulent current.

astrophysical plasmas↗