Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical model”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Surface Energy Modulated Self‐Assembly of Percolative Gallium Nanodroplet Films for Stretchable Electronics

Advancements in soft robotics and wearable technologies have led to an increased need for flexible and stretchable electronic systems. Liquid metals like Gallium (Ga) and its alloys offer high potential in the field of stretchable electronics and photonics due to their inherent high conductivity and malleability. However, their high surface tension and native oxide formation hinder uniform wetting and continuous film formation to get conducting films. Here, we report a fabrication strategy to realize electrically conductive Ga thin films with effective thicknesses around 100 nm on polymeric and elastomeric substrates. By enhancing the surface energy of the elastomer, we create uniform Ga nanodroplets capable of sustaining stretchability up to 200% strain. Our devices exhibit electronic conduction primarily through tunneling and display strain tunable piezoresistive characteristics, which are explained using Kirchhoff's resistive network and Simmon's analytical tunneling current model. To showcase the device's capability in fields of stretchable electronics, we demonstrate several functional prototypes, including flexible devices integrated with Light Emitting Diodes (LEDs), curvature and pressure sensors, and soft robotic grippers capable of sensing object curvature. The demonstrated devices highlight the potential of this approach for scalable integration of stretchable electronics in soft robotic and wearable platforms.

42 ENGINEERING↗
DOE OSTI2023

Identifying hydrologic signatures associated with streamflow depletion caused by groundwater pumping

Abstract Groundwater pumping can reduce streamflow in nearby waterways (‘streamflow depletion’), a process which must be accounted for in integrated management of surface and groundwater resources. However, causal identification of streamflow depletion from hydrographs alone is challenging because pumping impacts are masked by other drivers of hydrologic variability. To identify potential indicators of streamflow depletion, we used synthetic hydrographs and an analytical streamflow depletion model to assess potential pumping impacts on specific hydrograph characteristics (‘hydrologic signatures’) for 215 streamgages spanning the conterminous United States (CONUS). We found that streamflow depletion commonly impacts signatures associated with seasonal and annual low flows and low flow recessions. The largest impacts occurred during dry years, suggesting streamflow depletion may be evident in dry years even where impacts are unmeasurable in wet years. Random forest models indicated that streamflow depletion could significantly impact Annual, Summer, and Fall signatures in most streams. Our finding that multiple hydrologic signatures are consistently responsive to streamflow depletion across CONUS suggests that the underlying hydrological processes linking pumping to streamflow reductions are consistent across diverse settings, information that will aid in identifying indicators of streamflow depletion from streamflow hydrographs.

Lapides, Dana A.↗
DOE OSTI2022

Ultrafast Early Warning of Heart Attacks through Plasmon‐Enhanced Raman Spectroscopy using Collapsible Nanofingers and Machine Learning

Abstract As the leading cause of death, heart attacks result in millions of deaths annually, with no end in sight. Early intervention is the only strategy for rescuing lives threatened by heart disease. However, the detection time of the fastest heart‐attack detection system is >15 min, which is too long considering the rapid passage of life. In this study, a machine learning (ML)‐driven system with a simple process, low‐cost, short detection time (only 10 s), and high precision is developed. By utilizing a functionalized nanofinger structure, even a trace amount of biomarker leaked before a heart attack can be captured. Additionally, enhanced Raman profiles are constructed for predictive analytics. Five ML models are developed to harness the useful characteristics of each Raman spectrum and provide early warnings of heart attacks with >98% accuracy. Through the strategic combination of nanofingers and ML algorithms, the proposed warning system accurately provides alerts on silent heart‐attack attempts seconds ahead of actual attacks.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Sparsity Applications for Gradient‐Based Optimization of Wind Farms

Optimizing wind farms is essential for designing efficient energy systems, especially as farms grow larger and span multiple sites. However, this optimization becomes increasingly challenging due to the rising computational cost associated with more turbines. Gradient‐based optimization methods scale better than gradient‐free approaches for large problems, but the most computationally expensive component remains the calculation of gradients for the objective function and constraint Jacobians. To address this, we propose leveraging sparsity to accelerate gradient evaluations and reduce the size of the constraint Jacobian. Wind farms naturally exhibit sparsity—many turbines do not influence each other under certain wind directions. However, unlike traditional sparse problems with fixed patterns, wind farm sparsity is dynamic, requiring new strategies to handle changing interactions efficiently. This paper presents a study of sparsity in wind farm optimization and introduces several methods to exploit it. These strategies are tested on multiple farms using the analytic Cumulative Curl model, with gradients computed via automatic differentiation (AD). The same sparsity‐aware techniques are also applicable to finite difference (FD) methods, where they can yield even greater speedups due to the high cost of directional evaluations. Results show that sparse methods achieve up to a 10x speedup with less than ± 5% variance in optimized wake losses compared to traditional methods. These findings suggest that sparsity‐aware optimization not only maintains solution quality but also scales efficiently with farm size, enabling more comprehensive design exploration at reduced computational cost.

17 WIND ENERGY↗
DOE OSTI2021

First measurement of large area jet transverse momentum spectra in heavy-ion collisions

Jet production in lead-lead (PbPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV is studied with the CMS detector at the LHC, using PbPb and pp data samples corresponding to integrated luminosities of 404 μb$^{−1}$ and 27.4 pb$^{−1}$, respectively. Jets with different areas are reconstructed using the anti-k$_{T}$ algorithm by varying the distance parameter R. The measurements are performed using jets with transverse momenta (p$_{T}$) greater than 200 GeV and in a pseudorapidity range of |η| < 2. To reveal the medium modification of the jet spectra in PbPb collisions, the properly normalized ratio of spectra from PbPb and pp data is used to extract jet nuclear modification factors as functions of the PbPb collision centrality, p$_{T}$ and, for the first time, as a function of R up to 1.0. For the most central collisions, a strong suppression is observed for high-p$_{T}$ jets reconstructed with all distance parameters, implying that a significant amount of jet energy is scattered to large angles. The dependence of jet suppression on R is expected to be sensitive to both the jet energy loss mechanism and the medium response, and so the data are compared to several modern event generators and analytic calculations. The models considered do not fully reproduce the data.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Direct Fabrication of Vertically Stacked Double Barrier Tunnel Junctions Based on Graphene and h-BN

Direct manufacturing of two-dimensional material-based double barrier (DB) tunnel junctions, based on a lithography-free approach was developed. Graphene/h-BN/Graphene/h-BN/Graphene devices were deposited on Si/SiO2 substrates by employing a plasma enhanced chemical vapor deposition technique in a sequential manner. DB tunneling junctions with varying barrier widths (by varying the thickness of the second graphene layer) were studied. Samples were characterized using Raman, Atomic Force Microscopy and X-ray photoemission spectroscopy. The I–V characteristics of tunneling current showed resonant tunneling behavior at room temperature with a negative differential conductance. The behavior could be explained with quantum mechanical double barrier tunneling model in which analytic solutions to Schrödinger’s equation were obtained in each region of the system. Resonances in transmission probability coefficient for varying barrier widths were evaluated and compared with the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On the competition between mixing rate and uniformity in a coaxial hydrodynamic focusing mixer

Fast microfluidic mixers for use with line-of-sight integrating detection schemes pose unique challenges. Such detectors typically cannot discriminate signal from slow moving (e.g. near internal walls) and fast-moving portions of the fluid stream. This convolves reaction rate dynamics with fluid flow residence time dynamics. Further, the small cross sections of typical three-dimensional hydrodynamic focusing devices lead to lower detection signals. The current study focuses on achieving both small time scales of mixing and homogenous residence times. This is achieved by injecting sample through a center capillary and hydrodynamically focusing using a sheath flow within a tapered second capillary. The current design also features a third, larger coaxial capillary. The mixed stream flows into the large cross-section of this third capillary to decelerate and expand the stream by up to 14-fold to improve line-of-sight signal strength of reaction products. Hydrodynamic focusing, mixing, and expansion are studied using analytical and numerical models and also studied experimentally using a fluorescein-iodide quenching reaction. The experimentally validated models are used to explore trade-offs between mixing rate and uniformity. For the first time, this work presents detailed analysis of the Lagrangian time history of species transport during mixing inside coaxial capillaries to measure mixing nonuniformity. The mixing region enables order 100 μs mixing times and residence time widths of the same order (140 μs).

42 ENGINEERING↗
DOE OSTI2020

Geometry-independent microstructure optimization for electron beam powder bed fusion additive manufacturing

The ability to manipulate the formation of microstructure is one of the potential advantages of additive manufacturing. Yet, the additive manufacturing process is riddled with complex interactions between processing parameters, geometry, and spot melting sequence which makes the task of controlling microstructure challenging. This paper proposes a method to solve the microstructure texture optimization problem by leveraging established optimization techniques and increasingly fast and accurate analytical heat transfer models. Respectively, we apply two different optimization algorithms, genetic algorithm and concurrent greedy line search, in order to induce arbitrarily shaped microstructure regions and we use a heat transfer model to evaluate the fitness of solutions in an iterative approach that is far more efficient than pure trial and error. The objective of the optimization algorithm is to control local thermal conditions – hence microstructure formation – during the melting process by manipulating the order and the dwell time of all the points of the infill pattern sequence. We then analyze the optimized sequence to determine important path characteristics that affect the microstructure formation. Through examples, we show that relative order and frequency, as measured by weighted average and variance of spot order, are important factors that influence the final characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Numerical simulation of falling film sensible heat transfer over round horizontal tubes

Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.

42 ENGINEERING↗
DOE OSTI2026

Optimization of well design and CO 2 injection strategy for risk reduction in Class VI geological carbon sequestration wells

The safety and durability of Class VI wells are critical for geological carbon sequestration (GCS). However, current GCS operations face unique challenges: unlike traditional Class II wells, Class VI CO 2 injection wells operate at rates up to 100 times higher, dramatically increasing the risk of wellbore leakage and structural compromise due to severe temperature drops and associated mechanical stresses. Despite existing guidelines on material selection, there remains a substantial gap in understanding how rapid CO 2 injection rates, low surface temperatures, and variable reservoir conditions interact to threaten long-term well integrity. This study presents a comprehensive, original workflow integrating advanced analytical and numerical models for both well flow and well integrity analysis. By systematically simulating a wide range of field-relevant scenarios—including variations in injection rate, CO 2 temperature, and reservoir pressure—this work provides the first cross-validated assessment of cooling effects on wellbore. The results reveal that extreme temperature drops, up to 60 °C, can occur under high injection rates, particularly in depleted reservoirs, significantly increasing the risk of cement failure. Building on these insights, the study proposes innovative, practical well design and operational strategies, including ductile cement formulations, pre-stressing techniques, advanced insulation coatings, and proactive management of injection rates. The safety of Class VI well extends beyond simply using CO 2 resistant materials. Cement materials should possess optimal thermo-hydraulic-mechanical-chemical properties for effective performance. This work provides a scientific basis for optimizing Class VI well designs, with direct benefits for minimizing environmental risk, lowering operational costs, and enhancing the long-term reliability of GCS.

25 ENERGY STORAGE↗
DOE OSTI2020

Performance and wake flow characterization of a 1:8.7-scale reference USDOE MHKF1 hydrokinetic turbine to establish a verification and validation test database

As hydrokinetic turbine technologies continue to advance towards commercialization, public datasets on the performance characteristics for these devices and their flow field effects are invaluable to advance our understanding of these technologies and to validate analytical and numerical models. The Applied Research Laboratory at The Pennsylvania State University (ARL Penn State) collaborated with Sandia National Laboratories and the University of California at Davis to design, fabricate (at a 1:8.7 scale), and experimentally test a novel hydrokinetic turbine rotor design to provide an open platform and dataset for further study and development. The water tunnel test of this three-bladed, horizontal-axis rotor recorded power production, blade loading, the near-wake flow, cavitation effects, and noise generation. These state-of-the-art measurements demonstrate much of the complex physics associated with the flow through an unducted, horizontal-axis turbine, and they elucidate the performance characteristics and flow field effects at an unprecedented fidelity, accuracy and resolution. Measurements of powering coefficients (power, torque and thrust) as a function of tip-speed-ratio were performed. The dataset also includes unsteady measurements of driveshaft loading, blade strain, tower pressures, and radiated noise. Detailed flow mapping using laser Doppler velocimetry, and planar and stereo particle image velocimetry includes measurements of mean velocity and Reynolds stresses. Although the wake measurements are limited to less than half a diameter, they reveal the complex flow patterns in the near-wake structure of the rotor. The full database, available at the United States Department of Energy’s marine and hydrokinetic data repository, includes tunnel and model Computer Aided Design geometry files and inflow data sufficient for a “Model-the-Test” computational Verification and Validation study.

13 HYDRO ENERGY↗
DOE OSTI2024

A simple model for short-range ordering kinetics in multi-principal element alloys

Short-range ordering (SRO) in multi-principal element alloys influences material properties such as strength and corrosion. While some degree of SRO is expected at equilibrium, predicting the kinetics of its formation is challenging. We present a simplified isothermal concentration-wave (CW) model to estimate an effective relaxation time of SRO formation. Estimates from the CW model agree to within a factor of five with relaxation times obtained from kinetic Monte Carlo (kMC) simulations when above the highest ordering instability temperature. Further, the advantage of the CW model is that it only requires mobility and thermodynamic parameters, which are readily obtained from alloy mobility databases and Metropolis Monte Carlo simulations, respectively. The simple parameterization of the CW model and its analytical nature makes it an attractive tool for the design of processing conditions to promote or suppress SRO in multicomponent alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An adjoint method for determining the sensitivity of island size to magnetic field variations

An adjoint method to calculate the gradient of island width in stellarators is presented and applied to a set of magnetic field configurations. The underlying method for calculation of the island width is that of Cary & Hanson ( Phys. Fluids B, vol. 3, issue 4, 1991, pp. 1006–1014) (with a minor modification), and requires that the residue of the island centre be small. Therefore, the gradient of the residue is calculated in addition. Both the island width and the gradient calculations are verified using an analytical magnetic field configuration introduced by Reiman & Greenside ( Comput. Phys. Commun. , vol. 43, issue 1, 1986, pp. 157–167). The method is also applied to the calculation of the shape gradient of the width of a magnetic island in a National Compact Stellarator Experiment (NCSX) vacuum configuration with respect to positions on a coil. A gradient-based optimization is applied to a magnetic field configuration studied by Hanson & Cary ( Phys. Fluids , vol. 27, issue 4, 1984, pp. 767–769) to minimize stochasticity by adding perturbations to a pair of helical coils. Although only vacuum magnetic fields and an analytical magnetic field model are considered in this work, the adjoint calculation of the island width gradient could also be applied to a magnetohydrodynamic (MHD) equilibrium if the derivative of the magnetic field, with respect to the equilibrium parameters, is known. Using the island width gradient calculation presented here, more general gradient-based optimization methods can be applied to design stellarators with small magnetic islands. Moreover, the sensitivity of the island size may itself be optimized to ensure that coil tolerances, with respect to island size, are kept as high as possible.

Physics↗
DOE OSTI2025

Decoupling Carrier Dynamics and Energy Transport in Ultrafast Near-Field Nanoscopy

Ultrafast near-field optical nanoscopy has emerged as a powerful platform to characterize low-dimensional materials. While analytical and numerical models have been established to account for photoexcited carrier dynamics, quantitative evaluation of the associated pulsed laser heating remains elusive. Here, we decouple the photocarrier density and temperature increase in near-field nanoscopy by integrating the two-temperature model (TTM) with finite-difference time-domain (FDTD) simulations. These results reveal that the electron–phonon coupling in a silicon film after femtosecond laser excitation is most pronounced within approximately 3 ps–substantially shorter than the photocarrier decay time scale at tens of picoseconds. Moreover, the coupled TTM-FDTD method indicates that ultrafast laser heating can cause up to a 14% variation in the near-field signal at a 220 μJ/cm 2 pump pulse fluence. Our numerical results are further validated by transient experiments, highlighting the potential of this method for investigations of carrier and thermal phenomena in emerging nanomaterials and nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Emergent biaxiality in chiral hybrid liquid crystals

Biaxial nematic liquid crystals are fascinating systems sometimes referred to as the Higgs boson of soft matter because of experimental observation challenges. Here we describe unexpected states of matter that feature biaxial orientational order of colloidal supercritical fluids and gases formed by sparse rodlike particles. Colloidal rods with perpendicular surface boundary conditions exhibit a strong biaxial symmetry breaking when doped into conventional chiral nematic fluids. Minimization of free energy prompts these particles to orient perpendicular to the local molecular director and the helical axis, thereby imparting biaxiality on the hybrid molecular-colloidal system. The ensuing phase diagram features colloidal gas and liquid and supercritical colloidal fluid states with long-range biaxial orientational symmetry, as supported by analytical and numerical modeling at all hierarchical levels of ordering. Unlike for nonchiral hybrid systems, dispersions in chiral nematic hosts display biaxial orientational order at vanishing colloid volume fractions, promising both technological and fundamental research utility.

Colloids↗
DOE OSTI2020

Measurements of enhanced performance in an indirect drive inertial confinement fusion experiment when reducing the contact area of the capsule support

Experimental results from indirectly driven inertial confinement fusion experiments testing the performance gained from using an alternate capsule tent support are reported. The polar tent describes an alternate geometry for the thin membrane used to support the Deuterium–Tritium (DT) filled capsule. In this work, the contact area is reduced by 23 times by locating the tent support close to the poles of the capsule. The polar tent experiments are repeats of previous 3 shock 1.63 MJ, 400 TW high foot experiments and use a 165 μm thick silicon doped carbon hydrogen plastic (CH) shell. Using the polar tent support, we report a DT neutron yield of 1.07 ×1016, 76% higher than the expected Y DT ∝V 7.7 scaling. This is, at the time of writing, the highest neutron yield to date from a CH shell implosion. Furthermore, we find that the inferred pressure when using the polar tent is significantly above the model based on analytic scaling even when accounting for tent effects. Analysis of x-ray and neutron images shows the reduction of lobes produced by nominal tent features. The reduction of these features in the polar tent experiments leads to decreased low mode (P2 and P4) asymmetry compared to the nominal tent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nonlinear bubble competition of the multimode ablative Rayleigh–Taylor instability and applications to inertial confinement fusion

The self-similar nonlinear evolution of the multimode ablative Rayleigh–Taylor instability (RTI) and the ablation-generated vorticity effect are studied for a range of initial conditions. We show that, unlike classical RTI, the nonlinear multimode bubble-front evolution remains in the bubble competition regime due to ablation-generated vorticity, which accelerates the bubbles, thereby preventing a transition into the bubble-merger regime. We develop an analytical bubble competition model to describe the linear and nonlinear stages of ablative RTI. Here, we show that vorticity inside the multimode bubbles is most significant at small scales with large initial perturbation. Since these small scales persist in the bubble competition regime, the self-similar growth coefficient ab can be enhanced by up to 30% relative to ablative bubble competition without vorticity effects. We use the ablative bubble competition model to explain the hydrodynamic stability boundary observed in OMEGA low-adiabat implosion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗