Experimental and numerical study of heat transfer characteristics of single-phase free-surface fan jet impingement with automatic transmission fluid
Not Available
Feng, Xuhui (ORCID:0000000261559356)↗
SEARCH · Engineering Papers
Results for “numerical modelling”
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.
Effect of time and spatial domains on monolayer 2D material interface thermal conductance measurement using ns ET-Raman
Raman-based approaches to characterizing thermal properties of nanoscale and microscale materials have proven to be powerful experimental methods. Here, this work explores the potential of energy-transport state-resolved Raman (ET-Raman) for the precise determination of interfacial thermal conductance (G) between monolayer 2D materials and their substrates. ET-Raman employs a continuous wave laser and nanosecond (ns) pulsed laser to induce both steady state and transient heating while simultaneously collecting 2D material temperature information from the Raman signal. A 3D numerical model built using the finite volume method simulates the two differentiated thermal transport processes. The simulated temperature rises are compared with experimental Raman temperature information to extract G. In this study, we investigate the temporal effects of the pulse frequency and pulse width of the ns laser on precise characterization of G. By running multiple simulations for a range of pulsed laser frequencies and pulse widths, we determine that ET-Raman is most sensitive to precise interfacial thermal conductance measurements with short pulse widths (~below 20 ns) and with pulsing frequencies that ensure sufficient cooling time between consecutive pulses (~below 500 kHz). We also conclude that substrate thickness and laser spot size play critical roles in heat distribution and subsequent sensitivity to meaningful G measurements. As both laser spot size and SiO 2 substrate thickness decrease, detection of interfacial thermal transport becomes more difficult. These results offer essential guidance in the design of future Raman-based approaches to thermal measurements using pulsed lasers.
2D materials↗
Numerical simulation of evaporating wavy falling liquid films in laminar gas streams
Numerical simulations are performed to investigate the interfacial heat and mass transfer from evaporating wavy falling liquid films in interaction with laminar gas streams. The OpenFOAM solver has been used to conduct the simulations where the liquid-gas interface is resolved using the Volume of Fluid method of solving for three phases, i.e., liquid, vapor, and air. The configuration considered is a falling liquid (water) film on a heated vertical plate with a confined laminar moist air (gas) flow that is either (a) co-current or (b) counter-current to the downward liquid flow. The evaporation at the liquid-gas interface is driven by the interfacial gradient of the vapor mass fraction. Interfacial waves are triggered using a monochromatic forcing disturbance that leads to sinusoidal or solitary waves forming at the liquid-gas interface under respective forcing frequencies. Further, the numerical model is validated with the available experimental data. The results show nearly a 15% enhancement in time-averaged Sherwood number (Sh) due to film waviness (sinusoidal or solitary) at the lower volumetric gas flow rate, Q g = +50 (co-current) and Q g = -50 (counter-current). This enhancement in the Sh for both the waves further increases by 11% with Q g = +800 and 196% with Q g = -800. A closer examination of the mass transfer process over a wave demonstrates that with Q g = +50, the concentration of the gas side streamlines at the trough locations of the wave leads to higher values of Sh at these locations. However, with Q g = +800, although the overall Sh increases, vortices appear at the wave trough locations, leading to a corresponding decrease in the local Sh values. Correlations are proposed for predicting Sh under co-current and counter-current gas flow effects.
36 MATERIALS SCIENCE↗
Material analysis of metal hydrides for bulk hydrogen storage
The United States Department of Energy (DOE) Hydrogen Storage Engineering Center of Excellence (HSECoE), 2009–2014, was established to evaluate the viability of material-based hydrogen storage for light duty vehicles. The HSECoE concluded that there were no metal hydrides that simultaneously met the DOE targets for gravimetric and volumetric capacity, and rate of onboard refueling. However, it was noted by the DOE that potential applications existed for stationary hydrogen storage systems. An application of particular interest is using a metal-hydride to store a hydrogen supply for fuel cell backup power to a 5 MWe data center. A generalized numerical model that coupled mass and energy conservation, chemical kinetics and thermodynamics was developed to simulate the discharge process for a given metal-hydride. The model was applied to an interstitial and a complex metal-hydride. In conclusion, it was found that the interstitial metal-hydride, which is discharged by a pressure swing is preferrable in a practical sense to the complex metal-hydride, which requires a temperature swing discharge process.
08 HYDROGEN↗
Optimization of a Packed Particle Magnetocaloric Refrigerator: A Combined Experimental and Theoretical Study
Magnetocaloric refrigeration has strong potential towards achieving high efficiency hydrogen liquefaction. Optimising parameters such as particle size and operating cycles can have a significant impact on liquefaction performance. Here, this work reports on a room temperature magnetocaloric refrigeration prototype designed with helium as the heat transfer fluid and packed particle Gd 0.8875 Ce 0.1025 Si 0.84 Cr 0.19 as the refrigerant. The temperature spans and cooling power with three different particle sizes at two different operating cycles were measured. A maximum temperature span of 16.7 K and a maximum cooling power of 8.6 W·kg –1 were obtained. Reduction of performance in the faster cycle was observed in larger particles due to insufficient thermal diffusion time and its implications for cryogenic temperatures are discussed. A numerical model was developed based on the experimentally observed variation of the temperature span and cooling power with particle size and cycle duration. This model was then used to study the competing effects of viscous dissipation and thermal diffusion for different particle sizes and cycle durations.
08 HYDROGEN↗
Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times
Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.
42 ENGINEERING↗
3D Simulations of gas injection on callovo-oxfordian claystone assuming spatial heterogeneity and anisotropy
A series of gas injection tests on Callovo-Oxfordian (COx) claystone from the Bure underground research laboratory (URL) in France were carried out at the British Geological Survey (BGS). The tests were performed using a triaxial apparatus specifically designed to capture small volumetric strains induced by the injected gas flow and consequent material dilatancy. The long-duration experiments were monitored throughout. Measurements also included pressure, stresses (axial and radial stresses prescribed for each test stage), rate of gas inflow, gas outflow volume as well as pore-pressures observed at various points of the sample. A coupled hydro-gas-mechanical 3D numerical model has been developed to simulate the tests. Initial permeability is assumed heterogeneous throughout the specimen and embedded fractures are incorporated in the formulation. Gas pressure-induced deformations during the test lead to variations of permeability due to changes in matrix porosity and, especially, fracture aperture as well as fracture orientation due to material anisotropy. A programme of sensitivity analyses involving the variation of different aspects and parameters of the model contributes to a better understanding of the phenomena and highlights its complexity. The model is able to reproduce the observed behaviour of the tests.
58 GEOSCIENCES↗
An adaptive Hessian approximated stochastic gradient MCMC method
Bayesian approaches have been successfully integrated into training deep neural networks. One popular family is stochastic gradient Markov chain Monte Carlo methods (SG-MCMC), which have gained increasing interest due to their ability to handle large datasets and the potential to avoid overfitting. Although standard SG-MCMC methods have shown great performance in a variety of problems, they may be inefficient when the random variables in the target posterior densities have scale differences or are highly correlated. Here, we present an adaptive Hessian approximated stochastic gradient MCMC method to incorporate local geometric information while sampling from the posterior. The idea is to apply stochastic approximation (SA) to sequentially update a preconditioning matrix at each iteration. The preconditioner possesses second-order information and can guide the random walk of a sampler efficiently. Instead of computing and saving the full Hessian of the log posterior, we use limited memory of the samples and their stochastic gradients to approximate the inverse Hessian-vector multiplication in the updating formula. Moreover, by smoothly optimizing the preconditioning matrix via SA, our proposed algorithm can asymptotically converge to the target distribution with a controllable bias under mild conditions. To reduce the training and testing computational burden, we adopt a magnitude-based weight pruning method to enforce the sparsity of the network. Our method is user-friendly and demonstrates better learning results compared to standard SG-MCMC updating rules. The approximation of inverse Hessian alleviates storage and computational complexities for large dimensional models. Numerical experiments are performed on several problems, including sampling from 2D correlated distribution, synthetic regression problems, and learning the numerical solutions of heterogeneous elliptic PDE. The numerical results demonstrate great improvement in both the convergence rate and accuracy.
97 MATHEMATICS AND COMPUTING↗
Stability-based optimization of ACRT for the growth of CZT by the traveling heater method
Numerical modeling is employed to study the effects of several accelerated crucible rotation technique (ACRT) rotation schedules on the growth of cadmium zinc telluride (CZT) by the traveling heater method (THM). In conjunction with these analyses, a measure is developed to estimate the amplitude of disturbances associated with cellular interface growth morphologies that develop over time and across the melt-solid interface. Furthermore, this disturbance amplitude measure is used to assess the ability of ACRT to stabilize the growth interface and thereby minimize inclusion formation processes. Notably, our analyses reveal that interfacial instability is not well correlated with melt undercooling, explaining why classical ACRT approaches to reduce undercooling by mixing may not stabilize growth. Computations show that interface stability is best achieved via rotation schedules that accentuate outward Ekman flows during spin-up and, during spin-down, eliminate Taylor-Görtler flows and minimize inward Ekman flows.
36 MATERIALS SCIENCE↗
Three-dimensional mapping of mineral in intact shark centra with energy dispersive x-ray diffraction
The centra of shark vertebrae consist of cartilage mineralized by a bioapatite similar to bone's carbonated hy-droxyapatite, and, without a repair mechanism analogous to remodeling in bone, these structures still survive millions of cycles of high-strain loading. The main structures of the centrum are an hourglass-shaped double cone and the intermedialia which supports the cones. Little is known about the nanostructure of shark centra, spe-cifically the relationship between bioapatite and cartilage fibers, and this study uses energy dispersive diffraction (EDD) with polychromatic synchrotron x-radiation to study the spatial organization of the mineral phase and its crystallographic texture. The unique energy-sensitive detector array at beamline 6-BM-B, the Advanced Photon Source, enables EDD to quantify the texture within each sampling volume with one exposure while constructing 3D maps via specimen translation across the sampling volume. Herein this study maps a centrum from two shark orders, a carcharhiniform and a lamniform, with different intermedialia structures. In the blue shark (Prionace glauca, Carcharhiniformes), the bioapatite's c-axes are oriented laterally within the centrum's cone walls but axially within the wide wedges of the intermedialia; the former is interpreted to resist lateral deformation, the latter to support axial loads. In the shortfin mako (Isurus oxyrinchus, Lamniformes), there is some tendency for c-axis variation with position, but the situation is unclear because one dimension of the sampling volume is consid-erably larger than the thickness and spacing of the intermedialia's radially-oriented lamellae. Because elastic modulus in collagen plus bioapatite mineralized tissues varies significantly with both volume fraction of bio-apatite and crystallographic texture, the present 3D EDD-derived maps should inform future 3D numerical models of shark centra under applied load.
59 BASIC BIOLOGICAL SCIENCES↗
A hybrid numerical and machine learning framework for evaluating the performance of a 780 cm 2 aqueous organic redox flow battery
Aqueous organic redox flow battery (AORFB) is a promising cost-competitive technology for large-scale energy storage. Among existing work, the dihydroxyphenazine (DHP)-based AORFB has demonstrated high energy density and low capacity degradation in 10 cm2 cells during lab tests. However, its commercial-scale performance in more complex environments remains unknown, posing a barrier for commercialization. To address this gap, this work presents a comprehensive performance evaluation of a 780 cm 2 DHP-based AORFB by combining physics-based numerical model, machine learning (ML)-based surrogate models, and ML-derived sensitivity quantification. Specifically, we first select 12 key battery parameters that include 10 physicochemical quantities and 2 operation quantities, then select 6 performance metrics that include energy efficiency (EE), discharging capacity, charging energy, and power losses due to concentration, activation, and ohmic over-potentials. With such selection, 12800 combinations of the 12 parameters are subsequently generated using the Latin Hypercube Sampling method. These combinations, together with 38 pre-defined State of Charge, are then integrated to a validated AORFB model developed in COMSOL to compute the performance metrics. With both input parameters and performance metrics, 60 deep neural network (DNN) surrogate models are then trained to approximate the relationship between the 10 physicochemical quantities and 6 performance metrics at each flow rate and current density. Sensitivity scores are then calculated based on the DNN models. Two additional sensitivity analysis tools, i.e., MARS, and SHAP, are also used to cross-validate the sensitivity scores from the DNN. The results demonstrate that 1) the standard potential ranks the first in controlling EE and charging energy, 2) the membrane conductivity is most critical for power loss and EE, and 3) specific area and reaction rate control activation power loss.
25 ENERGY STORAGE↗
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↗
Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact
Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.
42 ENGINEERING↗
Mapping 3D grain and precipitate structure during in situ mechanical testing of open-cell metal foam using micro-computed tomography and high-energy X-ray diffraction microscopy
Open-cell metal foams are ultra-low-density cellular metals with complex hierarchical structures that span bulk, cell, ligament, and sub-ligament scales and give rise to desirable properties such as high strength-to-weight ratio and excellent energy absorption. Although literature suggests that intrinsic material structures at sub-ligament length scales (e.g., grains and precipitates) play an important role in mechanical behavior of open-cell metal foams, there are very few experimental measurements of such structures in three dimensions and for meaningful volumes of foam. This study seeks to map and track the three-dimensional (3D) grain and precipitate structures of an intact volume of open-cell aluminum foam by advancing microstructural characterization techniques that leverage X-ray micro-computed tomography (μCT) and far-field high-energy X-ray diffraction microscopy (FF-HEDM). A 6%-relative-density aluminum foam sample was mechanically tested in compression while μCT and FF-HEDM measurements were collected at interrupted loading states at beamline 1-ID of the Advanced Photon Source. Further, a new scanning strategy and reconstruction algorithm were established to enable characterization of a foam volume with diameter approximately four times wider than the nominal width of the X-ray beam. The result is a set of maps that detail both the 3D grain and precipitate structures throughout the foam volume at successive strain steps. A novel grain tracking procedure was developed to track individual grains within the foam volume by accounting for the large rigid-body motions that individual ligaments can undergo during mechanical loading. The ability to track grains and precipitate structures in three dimensions throughout large bulk deformation of ultra-low-density polycrystalline materials enables new possibilities for validating numerical models and investigating local failure mechanisms. Furthermore, the methods and procedures developed in this study could be applied to other ultra-low-density structures, such as additively manufactured lattices.
36 MATERIALS SCIENCE↗
Machine tool cross beam design, fabrication, and testing using metal big area additive manufacturing
This paper describes the application of metal Big Area Additive Manufacturing (mBAAM) to the fabrication of a machine tool cross beam. The replacement of a traditional box design weldment with a new design printed by wire arc additive manufacturing using the MedUSA system at Oak Ridge National Laboratory (ORNL) is detailed. This requires a new design strategy based on the unique mBAAM capabilities. The intent of the new design is to reduce mass, while maintaining the dynamic stiffness. To compare the two designs, the natural frequencies and mode shapes are measured using impact testing and predicted using finite element analysis. It is confirmed that the printed structure dynamics agreed with the numerical model predictions, which demonstrates that it is feasible to model a large-scale mBAAM part and understand its behavior prior to printing. Another notable outcome of this study is that the significant residual stress and distortion in the print indicate that knowledge gaps remain for widespread implementation of mBAAM.
42 ENGINEERING↗
Modeling snowflake divertors in MAST-U tokamak using UEDGE code
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Numerical analysis of dynamic load following response in a natural circulation molten salt power reactor system
The Molten Salt Reactor (MSR) concept is a rapidly evolving Generation IV design that has recently attracted favorable attention due to the potential for reducing waste generation, realizing passive safety features, and seizing on the opportunity for cost effective economics. An investigation into the power transient behavior of an autonomous load-following, closed-loop, natural circulation MSR system is important to quantifying operational and safety performance under dynamic conditions. This paper presents the results of a STAR-CCM+ and a comparatively simple asymmetric, one-dimensional, numerical model to solve the compound dynamic MSR power behavior subject to flow and temperature reactivity feedback only. Results show that reactor power is affected by fuel salt flow velocity (global) and temperatures (local bulk volume averages) in a coupled, time-delayed manner that results in a unique compound dynamic, closed-loop power feedback mechanism. The 1-D simulation approach opens the possibility of performing inexpensive computations to evaluate time-dependent reactor performance relative to thermo-physical fuel salt limitations. Results also potentially motivate a convincing conclusion that natural circulation MSRs provide a leap in safety and reliability.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Mid-height seismic isolation of equipment in nuclear power plants: Numerical simulations and design recommendations
Seismic isolation can mitigate earthquake demands on safety-class equipment and is being considered for application to some advanced nuclear reactors. For tall, slender vessels that could represent an advanced reactor, a steam generator or a heat exchanger, mid-height seismic isolation has been shown to be beneficial and practical. Here, this paper focuses on numerical modeling of a mid-height isolated tall, slender vessel with two primary goals 1) investigate and quantify the benefits of mid-height seismic isolation for a range of support structure stiffnesses, isolation systems, and seismic inputs, and 2) provide analysis and design recommendations for isolation of safety-class equipment. Results of response-history analysis were in good agreement with experimental measurements and demonstrated that mid-height isolation can substantially reduce seismic demands on tall, slender vessels for a range of support structure stiffnesses, isolation systems, and seismic inputs. Importantly, the reductions in horizontal spectral accelerations in the mid-height isolated vessel from the non-isolated condition were not affected by the stiffness of the support structure. Recommendations are made for analysis and design of isolated equipment and testing of isolators used for equipment protection in nuclear facilities, which are also applicable to other industries.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗