Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytical models”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Probing the Mechanical Performance of Micro-architected Porous Structures Through In Situ Characterization and Analysis

Micropores play critical roles in both natural and man-made materials. Such pores take on a variety of shapes and sizes ranging from spherical to irregular sphere-like voids with diameters spanning from the nanometer to millimeter scales. When porous structures are mechanically loaded, the pores direct the stress around their free surfaces, altering the material’s mechanical response relative to fully dense materials. In this study, for the first time, we create micro-architected porous samples using nanolithography and investigate the role of pore morphology by conducting a series of in situ micropillar compression tests in scanning electron microscopy. The findings demonstrate that porosity is the primary factor influencing the mechanical response of these micro-architected materials, as often seen on a macroscopic level. Additionally, we observed that pore geometry had a significant impact on Young’s modulus, yield stress, and strain energy density as a secondary parameter. Then, the extracted Young’s modulus was compared to macroscopic empirical models and determined the analytical models sufficiently described the impact of porosity in the microscopic scale but failed to capture the impact of second-order parameters. Furthermore, these results suggest how porous materials can be tailored to achieve desired mechanical properties based on the engineering applications of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Characterizing Physical Properties of Streambed Interface Sediments Using In Situ Complex Electrical Conductivity Measurements

Streambed sediment physical properties such as surface area, are difficult to quantify in situ but exert a high-level control on a wide range of biogeochemical processes and sorption of contaminants. We introduce the use of complex electrical conductivity (CC) methods (also known as spectral-induced polarization (SIP)) that measure both real and imaginary conductivity to noninvasively and efficiently characterize shallow streambed sediments. We explore the method through synthetic modeling, laboratory, and field measurements to demonstrate the sensitivity of imaginary conductivity to sediment surface area, controlled in part by fine-grained iron oxides produced by anoxic groundwater discharge. Laboratory measurements verify expected relationships between CC parameters and sediment properties. Synthetic modeling using a 1D analytical model illustrates the influence of water layer depth and conductivity on the field CC measurements made at the streambed-stream water interface. Specifically, the inverted sediment imaginary conductivity is less impacted by uncertainty in the water layer depth and conductivity relative to the real conductivity and phase shift. Field CC measurements along a landfill-impacted river reveal discrete streambed zones with enhanced bulk surface area generally corresponding to anoxic groundwater discharge zones with high concentrations of fine-grained iron oxide precipitates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Comparison of steady-state analytical wake models implemented in wind farm analysis software

A common set of mathematical wind turbine wake models are implemented in a few, well-adopted computational tools for wind farm wake modelling. Although the referenced mathematical formulations are common, implementation details may lead to differences in results. This study presents a systematic comparison of the implementation of mathematical wake models in open source, Python-based wind turbine wake modelling software, and a set of the models are directly compared. Despite aligning only the mathematical model parameters and retaining the default computational model parameters, good agreement is found across most of the model implementations, and additional agreement is expected upon further parameters alignment.

17 WIND ENERGY↗
DOE OSTI2020

Evaporation Induced Entrainment of Contaminants from Evaporating and Burning Liquids.

This interim report details model development, theory, and a literature review focusing on the evaporation induced entrainment (sub-boiling) of contaminated liquids. Entrainment from a variety of sources is the topic of DOE Handbook 3010, and this report deals more broadly with fire related airborne sources of contaminants in hazardous operations. Relatively few studies have examined sub-boiling behavior in the past, however, it can be a phenomenon that presents a fire related risk under hazardous operations. Molecular dynamics simulations are used to infer the gaseous evolution of coordinated complexes, and a model for a water/plutonium/nitrate system is deduced from the simulation results by evaluating the statistical trends of the results. Questions remain as to the chemical reactivity and longevity of entrained species. A generalized computer model capability and simple analytical model assumptions are developed for predicting the results of these and other (boiling and solid entrainment) scenarios. Verification related predictions using these models are illustrated.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Testing of an Optical Fiber--Based Gamma Thermometer in the High Flux Isotope Reactor Gamma Irradiation Facility

This report describes the design, thermal modeling, and gamma irradiation testing of an optical fiber–based gamma thermometer (OFBGT), which was irradiated in the High Flux Isotope Reactor (HFIR) Gamma Irradiation Facility (GIF). OFBGTs are a promising technology for application in nuclear reactors because they can provide a distributed measurement of gamma ray heating rate, unlike thermocouple-based gamma thermometers, which are fixed in-core sensors that can be used in boiling water reactors to calibrate local power range monitors. OFBGTs measure gamma ray heating rate by measuring the temperature difference between a pair of optical fibers; one fiber is in thermal contact with a heat sink (usually the reactor coolant), and the other is in thermal contact with a thermally isolated mass. The device can be calibrated with a heating wire within the thermal mass. The OFBGT that was designed and fabricated at Oak Ridge National Laboratory can measure distributed gamma ray heating rate over an effective measurement length of 61 cm, and the outer diameter of the sensor is 12.7 mm, giving the prototypical sensor design a relatively small footprint. The sensor housing is backfilled with Ar to ensure a well-predicted thermal response that is not affected by humidity or chemical interactions during operation. For calibration, the sensor design uses a Ni–Cr wire, which can be supplied with currents from 0 to 1 A to capture the wide range of potential gamma ray heating rates expected in HFIR’s spent fuel elements. The OFBGT was thermally modeled analytically and numerically; both models account for temperature-dependent thermal conductivities of the materials and show good agreement. The thermal response of the sensor inside spent HFIR fuel elements was simulated for times up to 1 year after discharge of the fuel element. Out-of-pile open-air tests indicated that the steady-state response of the sensor matches modeled results within experimental uncertainty. Calibration tests were performed using electrical heating in the HFIR spent fuel pool, above the fuel elements, to establish a relationship between the difference in spectral shift measured by optical fibers located inside and outside the OFBGT and the applied electrical heating. Subsequently, the OFBGT was placed within the the spent fuel element from HFIR cycle 501 to measure the spatial profile of the gamma heating rates. Results showed good agreement between the theoretical and measured gamma dose rate profiles, with maximum deviations of ~10% or less.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Techno-Economic Modelling of Tidal Energy Converter Arrays in the Tacoma Narrows

Hydrokinetic tidal energy converter (TEC) technology is yet to become cost competitive with other renewable energy sources. Understanding the interaction between energy production and the costs incurred harvesting that energy may unlock the economic potential of this technology. Although hydrodynamic simulation of TEC arrays has matured over time, including demonstration of how small and large arrays affect the resource, integration of cost modelling is often limited. The advanced ocean energy array techno-economic modelling tool ‘DTOcean’ enables designers to calculate and improve the levelised cost of energy (LCOE) of an array through parametric simulation of the energy extraction, design of the electrical network, moorings and foundations, and simulation of the installation and lifetime operations and maintenance of the array. This work presents a verification of DTOcean’s ability to simulate the techno-economic performance of TEC arrays by reproducing the hypothetical RM1 reference model, a semi-analytical model of a TEC array based in the Tacoma Narrows of Washington state, U.S.A. It is demonstrated that DTOcean can produce a reasonable estimate to the LCOE predicted by the reference model, giving (in Euro cents per kiloWatt hour) 36.69 ¢/kWh against the reference model’s 34.612 ¢/kWh for 10 TECs, while for 50 TECs, DTOcean calculated 20.34 ¢/kWh compared to 17.34 ¢/kWh for the reference model.

Topper, Mathew B. R. (ORCID:0000000347324347)↗
DOE OSTI2026

Introduction to and comparison of deep learning and optimization approaches to analytical wake modeling of a tilted wind turbine

This paper introduces innovative optimization and deep learning techniques to enhance the prediction of complex wake dynamics in the downstream wind velocity of tilted wind turbines. Traditional methods for calibrating the Bastankhah wake model often lead to increased errors in wind velocity distribution due to overfitting of the local wake characteristics. To address this issue, we propose an additional global optimization step to reduce errors in wind velocity predictions with respect to various wake parameters. Despite this improvement, the Bastankhah model's axisymmetric Gaussian wake shape limits its accuracy for complex wake structures. Therefore, we also propose a deep learning approach, which demonstrates promising results by accurately modeling complex wake shapes across a broader range of tilt angles with minimal computational cost. The deep learning approach achieves near-identical predictions to high-fidelity large-eddy simulations, representing a promising advancement in wake modeling.

17 WIND ENERGY↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗
DOE OSTI2020

Improved Frequency-Domain Steady-State Modeling of the Dual-Active-Bridge Converter Considering Finite ZVS Transition Time Effects

This article presents an improved analytical modeling (IAM) approach for the dual-active-bridge (DAB) converter in the frequency-domain. Specifically, finite transition times are incorporated into the improved modeling that are significant in DAB converters operated at high frequency. The IAM approach is first motivated by comparing traditional frequency-domain modeling against realistic simulation results in LTspice. Next, the IAM approach is mathematically developed in an iterative fashion, where the conventional frequency-domain approach is used as the initialization of the operational behavior of the converter. Subsequent iterations of the IAM algorithm incorporate finite rise- and falltimes of device output capacitances based on detailed switching behavior analyses, which may considerably perturb the desired modulation parameters (i.e., phase shift and zero states). While the initially proposed IAM approach demonstrates advantages over ideal modeling, it is still shown to exhibit errors when applied to device transitions with significant nonlinearity in parasitic output capacitance. To account for this, a three-slope approximation to the switching transition is then motivated and developed, which enables significantly improved steady-state prediction accuracy. Simulation and experimental results confirm unacceptable inaccuracies associated with the conventional frequency-domain modeling across a wide range of operating conditions, alongside enhanced accuracy provided by the proposed modeling approach.

14 SOLAR ENERGY↗
DOE OSTI2020

The Universe at z > 10: predictions for JWST from the universemachine DR1

The James Webb Space Telescope (JWST) is expected to observe galaxies at z > 10 that are presently inaccessible. Here, we use a self-consistent empirical model, the universemachine, to generate mock galaxy catalogues and light-cones over the redshift range z = 0-15. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy–halo relationships, and galaxy–galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at z < 10. We predict that Cycle 1 JWST surveys will very likely detect galaxies with M* > 10 7 M⊙ and/or M 1500 < -17 out to at least z ~ 13.5. Number density uncertainties at z > 12 expand dramatically, so efforts to detect z > 12 galaxies will provide the most valuable constraints on galaxy formation models. Additionally, the faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits (M* < 10 7 M⊙ or M 1500 > -17). For reionization models, extrapolating luminosity functions with a constant faint-end slope from M 1500 = -17 down to M 1500 = -12 gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to z ~ 12. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and light-cones for common fields are available online.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗
DOE OSTI2020

Reexamining supercritical gas adsorption theories in nano-porous shales under geological conditions

Adsorption is an important phenomenon in surface chemistry, especially for nano-porous shales. In the shale gas-in-place, adsorbed gas could contribute up to 85%. However, adsorption is hard to quantitatively characterize due to different adsorption mechanisms, patterns, surfaces, and pore sizes. Moreover, key thermodynamic parameters, such as the enthalpy of adsorption, are challenging to determine due to uncertainties in adsorbed gas densities used in constructing absolute isotherms. In this paper, we revisit the Brunauer, Emmett, and Teller (BET) model, analytically simplify Ono-Kondo (OK) models for subsurface shales, and compare commonly used mono- and multilayer adsorption models (e.g., Langmuir, supercritical Dubinin-Radushkevich (SDR), supercritical BET (SBET), and simplified OK (OKs) models) with recently proposed pressure-dependent adsorption densities to develop a practical and reliable methodology that can be used in the supercritical state, typical for subsurface black shale conditions.Three independent data sets were used for nitrogen and methane adsorption isotherms at different temperatures. Furthermore, we demonstrate that adsorption predicted by the SDR model is comparable or lower than that predicted by the SBET model, but higher than the amounts predicted by Langmuir and OK models. The nitrogen BET method tends to underestimate the accessible SSA for methane. Measurement of isosteric heat of adsorption is suggested instead of the experimental fitting method, due to the significant difference between the calculated results by two commonly used methods. Experimental fitting and simulation methods are also briefly reviewed to guide future research on shale gas adsorption.

58 GEOSCIENCES↗
DOE OSTI2020

Prediction for Pressure Differential Across HEPA Filter Media Based on Media Characteristics and Particle Size Distribution

A new method for predicting the pressure drop across High Efficiency Particulate Air (HEPA) filter media is proposed based upon mass deposited onto the filter and known physical characteristics of the filter media. Detailed are the methods used in conjunction with current filter loading models to predict the pressure drop, as well as tests conducted to validate the methods. The benefit of a prediction model for practical use lies in the manufacturing and service life of nuclear grade HEPA filters. HEPA filters for use in nuclear facilities have a prescribed expiration date and a maximum allowable operating pressure drop. Therefore, the ability to predict the pressure drop across a filter and relate it to an expected length of service time can enable a reduction in wasted filters. This will allow for more informed decisions to be made based upon the dictated life cycle of the filters. Also, understanding and predicting how the pressure drop of HEPA filter media behaves as a function of physical characteristics and loaded mass can assist in future design and manufacturing of filter media. Currently, most existing pressure drop models are either computationally based or analytical methods relying on data gathered during tests. Neither are practical for prediction; the computational methods are difficult to implement, and the current analytical methods are more useful as tools for analysis. The current analytical model developed by Bergman et al. is based upon the pressure drop from each media fiber and modeling the deposited particles as newly formed fibers. Using Bergman's model and the media properties, air properties, and known aerosol particle size distribution, a pressure drop curve as a function of loaded mass can be created. This curve implies initial loading in the depth of the filter media with a transition to pure surface loading on the filter face, and a sensitivity to an evolving particle size distribution as mass continues to load onto the filter. To implement the predictive model, knowledge of the mean filter fiber diameter and porosity of the filter media is required. Traditionally, the mean fiber diameter is calculated as an effective diameter from prior media testing, however, in this study a Scanning Electron Microscope (SEM) technique is used to acquire this variable. Validation of the predictive model is provided by flat sheet media tests under a controlled environment with a measured particle size distribution of the challenge aerosol and shows reasonable preliminary agreement.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A two-dimensional analytical unit cell model for redox flow battery evaluation and optimization

Cell performance optimization is important for improving the overall system efficiency of a redox flow battery. To gain better insights into key controlling factors of system efficiency, this work first proposed a theoretical model for a unit flow battery cell by extending a two-dimensional analytic model to a full battery cell. Such a model is then used for cell performance optimization after validating it with experimental and numerical modeling data. With the model results, the activation, equilibrium, and pump energy losses are identified as the dominant sources of battery energy losses. Further, a guideline for reducing these sources is also proposed. Following the guideline, the mass transport coefficient is shown as a key control factor of the equilibrium energy loss and Coulombic efficiency (CE). Approaches are then proposed to improve CE and the overall system efficiency. The mass transport is also revealed as the mechanism of pump rate optimization where an optimal pump rate significantly reduces the equilibrium energy loss. With both low equilibrium and pump energy losses, an optimal electrode porosity or specific area design can further improve a battery's system efficiency based on an optimal porosity predicted by the present model. The model also demonstrates distinct behaviors and overestimation in the system efficiency when reduced to a zero-dimensional model with neglected mass transport resistance. With the new model, the guideline, and new insights, this work provides a reliable and efficient tool for the evaluation and optimization of redox flow battery design in practical applications.

25 ENERGY STORAGE↗
DOE OSTI2023

Vortex Structure of Head Bubble in Convective Cloud Starting Plume

Abstract A developing cumulus cloud (Cu) was modeled, and dynamic, thermodynamic, and microphysical properties of an ascending head bubble reproducing the upper part of a developing Cu were investigated. The data for analysis are taken from 10-m-resolution LES of trade wind Cu under BOMEX conditions. The detection of a rising bubble is carried out using wavelet filtering of the velocity fields and microphysical fields, while a low-frequency signal of the filtering is associated with the convective-scale structure of cloud. We substantiate and discuss the representation of the bubble as a vortex ring, and estimate the parameters of this vortex ring. The simplest Hill’s vortex was chosen as a model of a vortex ring inside cloud. Analytical approximations of the radial profiles of the vertical velocity and of conservative quantities (such as total water mixing ratio and liquid water potential temperature inside and outside the bubble) are obtained. The spatial structure of these quantities is investigated using analytical expressions. Analytical models for spatial distributions of liquid water content (LWC) and adiabatic fraction (AF) are also designed and analyzed. The results demonstrate the existence of a cloud core with high values of LWC and AF up to the height of 1800 m. The horizontally averaged value of the adiabatic fraction, calculated analytically using the Hill’s vortex concept, is evaluated as 0.39, which is the typical AF value in the upper parts of such Cu. The vertical profiles of different important quantities characterizing cloud structure are presented. The analysis performed in this study allows us to conclude that a rising vortex ring plays the dominating role in formation of the thermodynamic and microphysical structure of developing Cu. Significance Statement 1) Dynamic and thermodynamic fields of a developing cumulus cloud simulated by high-resolution LES with spectral bin microphysics are separated into convective and turbulent components by means of the wavelet technique. 2) The analysis of convective component of the cloud revealed the existence of a vortex ring at the developing stage of the cloud and evaluate its parameters. 3) The analysis performed in this study allows us to conclude that a rising vortex ring plays the important role in formation of the thermodynamic and microphysical structure of developing Cu. 4) The study provides a novel insight into the cloud–environment interaction. 5) The approximating equations describing the vortex ring can be usefully applied for developing new schemes of convective parameterization.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

More Than 30 000-fold Field Enhancement of Terahertz Nanoresonators Enabled by Rapid Inverse Design

The rapid development of 6G communications using terahertz (THz) electromagnetic waves has created a demand for highly sensitive THz nanoresonators capable of detecting these waves. Among the potential candidates, THz nanogap loop arrays show promising characteristics but require significant computational resources for accurate simulation. This requirement arises because their unit cells are 10 times smaller than millimeter wavelengths, with nanogap regions that are 1 000 000 times smaller. To address this challenge, we propose a rapid inverse design method using physics-informed machine learning, employing double deep Q-learning with an analytical model of the THz nanogap loop array. In ~39 h on a middle-level personal computer, our approach identifies the optimal structure through 200 000 iterations, achieving an experimental electric field enhancement of 32 000 at 0.2 THz, 300% stronger than prior results. In conclusion, our analytical model-based approach significantly reduces the amount of computational resources required, offering a practical alternative to numerical simulation-based inverse design for THz nanodevices.

21 terahertz time-domain spectroscopy↗