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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High temperature oxidation behaviors of bulk SiC with low partial pressures of air and water vapor in argon

Oxidation behaviors of bulk SiC at 1000–1400 °C in Ar-20 vol% H 2 O-10 vol% air were investigated. The thickness of the passivation oxide layer showed a parabolic increase with time. Dense and uniform layers formed below 1200 °C, whereas pores and cracks were observed at higher temperatures. An analytical model was proposed to estimate the oxide layer growth and SiC mass loss rate. The SiC layer should be effective in retaining the fission products after oxidation by excessive water and air. This work is the first detailed study of SiC layer oxidation during loss of coolant accidents (LOCA) for HTGRs.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Unusual size effects from tilted twin boundaries in nano-twinned metals

Tilted twin boundaries (TBs), whose plane normals are rotated at an angle from the parent grain axis, naturally occur in columnar-grained, nano-twinned (NT) metals. Here, using a combination of atomistic simulations and analytical modeling, we reveal that NT metals with the ideal, non-tilted TBs exhibit continuously increasing strength with decreasing twin thickness, and hence, no inverse twin thickness size effect on strength. In contrast, NT metals with tilted TBs exhibit an inverse size effect, and the critical twin thickness, below which strength decreases, increases as the TB-tilt angle increases. The analysis also identifies a critical value of TB tilt, for which strength becomes independent of twin thickness and is the weakest. The transition arises from a change in dislocation activity prevailing mostly on planes inclined to TBs to planes parallel to the TBs. Overall, these findings reveal a profound influence of TB tilt angle that could redirect the analysis and engineering of nano-twin structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On electrolyte wetting through lithium-ion battery separators

The wettability by electrolyte is a critical characteristic of lithium-ion battery separators since electrolyte absorption is essential for ionic transport. Most importantly, fast absorption of electrolyte facilitates the electrolyte filling and wetting processes during the battery assembly. However, theoretical or experimental investigation of separator wettability has rarely been reported, and thus very little quantitative data can be found in the existing literature. Here, innovative experimental techniques and associated analytical models have been developed to fully characterize the electrolyte wetting behavior through battery separators. Both in-plane and through-plane imbibition processes are investigated. The results show that, electrolyte needs significantly longer time to wet separators than it does for commercial electrodes. The wetting rate depends significantly on the pore size, porosity, and micro-structure of the separators. Electrolyte in-plane wetting is only observed in the separators manufactured by the wet process. Increasing separator porosity and pore size is shown to be substantially increasing through-plane wetting rate. In addition, ceramic coating is shown to increase in-plane wetting but retard through-plane wetting.

25 ENERGY STORAGE↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2025

Impact of molten salt inflow on the temperature distribution in thermal energy storage tanks at startup for central receiver concentrating solar power plants

Concentrating Solar Power (CSP) systems with molten salt thermal energy storage (TES) tanks are one of the most promising, renewable-based energy conversion technologies for larger-scale power generation. The TES tank is one of the most critical components in CSP plants due to its high-temperature operation (up to 565 °C), daily thermal cycling, and intermittent solar radiation conditions. The plant startup is one of the most challenging operation conditions that could lead to damaging thermal gradients due to low salt inventory levels. In this study an analytical model for the sparger ring was developed and integrated with a detailed computational fluid dynamics model of a commercial-scaled molten salt tank. The integrated model allows an accurate representation of the tank operation to evaluate the effect of molten salt inflow on the mixing process. The tank filling process during plant startup was analyzed considering sparger rings with variations in design features, including inlet orifice configurations, number of orifices, direction of the inlets, and orifice diameter. The results demonstrated that higher temperature gradients are obtained in the tank floor during the plant startup. A sparger ring configuration with a predetermined orifice inlet inclination (30°, 45° and 60°) leads to significant temperature differences in the floor, between 57 °C and 62 °C, but better homogeneity in the temperature of the salt inventory. Lower salt inflow velocities result in a more homogeneous floor temperature, with maximum temperature differences under 38 °C. The sparger ring configuration with 52 orifices of 1-in. diameter and vertical flow showed better homogeneity in the temperature differences as a function of the salt level and lower temperature gradients in the tank floor. Because large temperature gradients in the tank's floor have been identified as one of the main factors contributing to tank failures, assessing various sparger ring design features is fundamental to determining proper inflow conditions that lead to low-temperature gradients and reducing failure susceptibility.

14 SOLAR ENERGY↗
DOE OSTI2020

Analytical solution for partial heating on the exterior of the pipe with application to measuring fluid flow rate

An analytical model using Green’s functions for partial external heating of a pipe is developed, which results in an exact mathematical solution for the radial and axial temperature distribution in the pipe wall. Partial heating consists of a constant heat flux function imposed over a small section of the exterior of a pipe and for a limited time duration. The solution comprises steady-state and transient parts, and an algebraic identity is used to decrease the number of summation terms in the slowly-converging steady- state part. Intrinsic verification principles are used to verify the solution. As an example application, this transient solution is applied toward the development of a simple, noninvasive method for in-field measurement of the flow rate in pipes. Furthermore, to simulate this application, a pulse of energy is imposed to the wall of the pipe, and the developed mathematical solution is used to find the flow rate of the fluid inside the pipe. An optimal experiment is designed to find the best measurement location and time.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Generalized Analysis of Dynamic Flow Fouling on Heat Transfer Surfaces

Several energy applications involve flow of untreated heat transfer fluids that contain trace amounts of mineral salts, such as calcium carbonate and calcium sulfate, which show reverse solubility with temperature. The gradual accretion of these salts onto the heat transfer surfaces causes fouling that degrades thermal performance over time. While fouling of bare surfaces has been extensively studied, fouling of nonwetting superhydrophobic or liquid-infused surfaces is less understood, especially under dynamic flow conditions. This study presents a systematic experimental analysis to elucidate dynamic flow fouling of calcium carbonate and calcium sulfate on smooth, superhydrophobic and liquid-infused surfaces in a generalized manner. Here, copper tubes with modified inner surface wettability are fabricated and subjected to forced convection at different flow rates and fouling salt concentrations. Fouling characteristics of the two salts on the different surfaces are examined through microstructural analysis, and quantified in terms of asymptotic fouling resistance that measures the total extent of fouling, and a threshold time, which correlates to the onset of significant fouling. Using Taguchi design of experiments, the study reports, for the first time, closed-form analytical relationships for the asymptotic fouling thermal resistance and threshold time on the Reynolds number, dimensionless concentration of fouling agent and dimensionless infused liquid viscosity that represents the different surface types in a unified manner. It is shown that the analytical model accurately predicts the fouling characteristics in up to about 97% of the data, based on which contour maps of optimum surface designs for minimizing fouling resistance are presented. The results provide information on designing surfaces for reduced fouling or for estimating the fouling characteristics for a given surface in applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Novel nonwetting solid-infused surfaces for superior fouling mitigation

Fouling is a ubiquitous issue in several environmental and energy applications. Here we introduce novel nonwetting solid-infused surfaces (SIS) with superior anti-fouling characteristics that are durable than conventional nonwetting surfaces in a dynamic flow environment. A systematic study is presented to elucidate the fouling mitigation performance of SIS in comparison to lubricant-infused surface (LIS) and conventional smooth surface. Copper tubes with SIS, LIS or smooth inner walls are fabricated and subjected to accelerated calcium sulfate fouling in a flow fouling experimental setup. Fouling on the various surface types is quantified in terms of asymptotic fouling resistance, and the fundamental morphological differences in the interactions of the foulant and the various surface types are analyzed. Based on a systematic sweep of the parameter combinations using design of experiments and Taguchi analysis, an analytical dependence of asymptotic fouling resistance on the governing parameters namely, Reynolds number, foulant concentration and temperature is derived. The analytical model is shown to predict the asymptotic fouling resistance to within 20% accuracy with a 95% confidence. In addition, for the first time, the effects of shear durability on the fouling mitigation performance of LIS vis-à-vis SIS are studied. It is shown that the novel nonwetting SIS offers a robust option for superior fouling mitigation over LIS in the long run.

42 ENGINEERING↗
DOE OSTI2022

Depth-resolved thermal conductivity and damage in swift heavy ion irradiated metal oxides

We investigated thermal transport in swift heavy ion (SHI) irradiated insulating single crystalline oxide materials: yttrium aluminum garnet- Y 3 Al 5 O 12 (YAG), sapphire (Al 2 O 3 ), zinc oxide (ZnO) and magnesium oxide (MgO) irradiated by 167 MeV Xe ions at 10 12 – 10 14 ions/cm 2 fluences. Depth profiling of the thermal transport on nano- and micro- meter scales was assessed by time-domain thermoreflectance (TDTR) and modulated thermoreflectance (MTR) methods, respectively. This combination allowed us to isolate the conductivities of different sub-surface damage-regions characterized by their distinct microstructure evolution regimes. Thermal conductivity degradation in SHI irradiated YAG and Al 2 O 3 is attributed to formation of ion tracks and subsequent amorphization, while in ZnO and MgO it is mostly due to point defects. Additionally, notably lower conductivity when probed by very low penetrating thermal waves is consistent with surface hillock formation. An analytical model based on Klemens-Callaway method for thermal conductivity coupled with a simplified microstructure evolution capturing saturation in defect concentration was used to obtain depth dependent damage across the ion impacted region. The studies showed that YAG has the highest damage profile resulting in the less dependence of thermal conductivity with the depth, while MgO on the contrary has the strongest dependence. The presented work sheds new light on how SHI induced defects affect thermal transport degradation and recovery of oxide ceramics as promising candidates for next generation nuclear reactor applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Energy transfer for storage or recovery in capacitive deionization using a DC-DC converter

Energy recovery from capacitive deionization (CDI) has the potential to increase overall desalination efficiency. We here define the storage (one-way) and utilization (round-trip) efficiencies between a CDI cell and an energy storage device using a generic direct current/direct current (DC/DC) converter circuit. Presented is a closed-form analytical model for the case of a supercapacitor (SC) as the energy storage device and a buck-boost converter as the DC/DC converter. The model is benchmarked with a numerical model, showing good agreement. Also presented is a comparison among energy transfer methods wherein desalination productivity is fixed. For constant current operation, this condition requires higher currents to compensate for the inactive time introduced by the converter. For fixed productivity, the mean current of the converter circuit approaches the constant current value for higher initial voltages and lower capacitances in the supercapacitor. Finally, we show the effect of relevant parameters of the SC in the storage and utilization efficiency. The model predicts storage and utilization efficiencies of at least 90% for an initial voltage of 1 V or higher and reasonable CDI and SC parameters. Finally, we provide engineering operational parameters to maximize the efficiency of energy transfer and guidance in the selection of electronic components.

25 ENERGY STORAGE↗
DOE OSTI2020

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗
DOE OSTI2025

Glassy carbon formation from pyrolysis of polymeric coatings on fiber-optic sensors

Deploying fiber-optic sensors in nuclear reactors requires a detailed understanding of radiation effects on the fiber materials and the transmitted signals. Previous work has shown large wavelength shifts in the reflected spectra obtained from polymer-coated fiber-optic temperature sensors exposed to high neutron fluences. The sensor drift resulting from these wavelength shifts cannot be explained by radiation effects on fused silica glass. These shifts are hypothesized to be caused by the conversion of the polymeric fiber coating to a glassy carbon via radiolysis and/or pyrolysis and subsequent radiation-induced compaction. Here, thermal degradation of these polymeric coatings was studied to provide insight into the potential origins of the sensor drift phenomenon. Acrylate- and polyimide-coated fibers were heated under various temperatures (250–1300 °C) and environments (oxidative and inert), and the resulting coating products were characterized via mass-loss data, scanning electron microscope imaging, and Raman spectroscopy. Results suggest that the polymer decomposition product of both coating types, at least under inert conditions, is indeed a glassy carbon. Analytical models that account for radiation-induced glassy carbon coating compaction show significant compressive fiber strains and predicted wavelength shifts that agree well with experimental measurements, providing additional evidence that supports the hypothesized origins of the sensor drift.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of alloying on intrinsic ductility in WTaCrV high entropy alloys

Tungsten (W) exhibits desirable properties for extreme applications, such as the divertor in magnetic fusion reactors, but its practicality remains limited due to poor formability and insufficient irradiation resistance. In this work, we study the intrinsic ductility of body-centered cubic WTaCrV based high entropy alloys (HEAs), which are known to exhibit excellent irradiation resistance. The ductility evaluations are carried out using a criterion based on the competition between the critical stress intensity factors for emission (K Ie ) and cleavage (K Ic ) in the {110} slip planes and {110} crack planes, which are evaluated within the linear elastic fracture mechanics framework and computed using density functional theory. The results suggest that increasing the alloying concentrations of V and reducing the concentrations of W can significantly improve the ductility in these HEAs. The elastic anisotropy for these HEAs is analyzed using the Zener anisotropy ratio and its correlation with the concentration of W in the alloys is studied. Results indicate that these alloys tend to be fairly isotropic independently from the concentration of W in them. The computed data for the elastic constants of these HEAs is also compared against available experimental data. The results are in good agreement, validating the robustness and accuracy of the computational methods. Multiple phenomenological ductility metrics were also computed and analyzed against the analytical model. Some metrics, mainly the surrogate D parameter, show a good correlation with the Rice model. The potential of these empirical metrics to serve as surrogate screening models for optimizing the compositional space is also discussed.

Anisotropy↗
DOE OSTI2019

The small-angle neutron scattering extension in MCNPX and the SANS cross section for nanodiamonds

Detonation nanodiamonds provide an intriguing prospect for the development of a very cold neutron (VCN) source or as a high albedo material for VCN reflection helping to control leakage of cold and very cold neutrons in a moderator. The investigation of source and reflector designs necessitates the development of new simulation tools for existing Monte Carlo codes that include small angle neutron scattering (SANS). The SANS extension in MCNPX 2.7 provides an analytical model of hard-sphere scattering with variable scattering particle sizes, size distributions, and packing fractions in order to supplement currently existing scattering kernels as well as the capability to read an experimentally determined small angle neutron scattering profile from a data file. The SANS extension is discussed and compared to nanodiamond SANS measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Characterization of the new Ultracold Neutron beamline at the LANL UCN facility

The neutron electric dipole moment (nEDM) experiment that is currently being developed at Los Alamos National Laboratory (LANL) will use ultracold neutrons (UCN) and Ramsey’s method of separated oscillatory fields to search for a nEDM. Here, in this paper, we present measurements of UCN storage and UCN transport performed during the commissioning of a new beamline at the LANL UCN source and demonstrate a sufficient number of stored polarized UCN to achieve a statistical uncertainty of δd n = 2 x 10 -27 e cm in 5 calendar years of running. We also present an analytical model describing data that provides a simple parameterization of the input UCN energy spectrum on the new beamline.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Measurement of the Higgs boson mass in the H → ZZ*→ 4ℓ decay channel using 139 fb –1 of s = 13 TeV pp collisions recorded by the ATLAS detector at the LHC

The mass of the Higgs boson is measured in the H→ZZ * →4ℓ decay channel. The analysis uses proton–proton collision data from the Large Hadron Collider at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector between 2015 and 2018, corresponding to an integrated luminosity of 139 fb –1 . The measured value of the Higgs boson mass is 124.99 ± 0.18(stat.) ± 0.04(syst.) GeV. In final states with muons, this measurement benefits from an improved momentum-scale calibration relative to that adopted in previous publications. The measurement also employs an analytic model that takes into account the invariant-mass resolution of the four-lepton system on a per-event basis and the output of a deep neural network discriminating signal from background events. This measurement is combined with the corresponding measurement using 7 and 8 TeV pp collision data, resulting in a Higgs boson mass of 124.94 ± 0.17(stat.) ± 0.03(syst.) GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The exit-wave power-cepstrum transform for scanning nanobeam electron diffraction: robust strain mapping at subnanometer resolution and subpicometer precision

Scanning nanobeam electron diffraction (NBED) with fast pixelated detectors is a valuable technique for rapid, spatially resolved mapping of lattice structure over a wide range of length scales. However, intensity variations caused by dynamical diffraction and sample mistilts can hinder the measurement of diffracted disk centers as necessary for quantification. Robust data processing techniques are needed to provide accurate and precise measurements for complex samples and non-ideal conditions. Here we present an approach to address these challenges using a transform, called the exit wave power cepstrum (EWPC), inspired by cepstral analysis in audio signal processing. The EWPC transforms NBED patterns into real-space patterns with sharp peaks corresponding to inter-atomic spacings. We describe a simple analytical model for interpretation of these patterns that cleanly decouples lattice information from the intensity variations in NBED patterns caused by tilt and thickness. With tracking the inter-atomic spacing peaks in EWPC patterns, strain mapping is demonstrated for two practical applications: mapping of ferroelectric domains in epitaxially strained PbTiO 3 films and mapping of strain profiles in arbitrarily oriented core-shell Pt-Co nanoparticle fuel-cell catalysts. The EWPC transform enables lattice structure measurement at sub-pm precision and sub-nm resolution that is robust to small sample mistilts and random orientations.

42 ENGINEERING↗