Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPERIMENTAL DATA”

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 811 records · Page 45

PAB3D: Its History in the Use of Turbulence Models in the Simulation of Jet and Nozzle Flows

This is a review paper for PAB3D s history in the implementation of turbulence models for simulating jet and nozzle flows. We describe different turbulence models used in the simulation of subsonic and supersonic jet and nozzle flows. The time-averaged simulations use modified linear or nonlinear two-equation models to account for supersonic flow as well as high temperature mixing. Two multiscale-type turbulence models are used for unsteady flow simulations. These models require modifications to the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k-epsilon) model with a RANS/LES transition function, dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) formulation. All of these models are implemented in the three-dimensional Navier-Stokes code PAB3D. This paper discusses computational methods, code implementation, computed results for a wide range of nozzle configurations at various operating conditions, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions.

Abdol-Hamid, Khaled S.↗

Spatial Pattern of Cell Damage in Tissue from Heavy Ions

A new Monte Carlo algorithm was developed that can model passage of heavy ions in a tissue, and their action on the cellular matrix for 2- or 3-dimensional cases. The build-up of secondaries such as projectile fragments, target fragments, other light fragments, and delta-rays was simulated. Cells were modeled as a cell culture monolayer in one example, where the data were taken directly from microscopy (2-d cell matrix). A simple model of tissue was given as abstract spheres with close approximation to real cell geometries (3-d cell matrix), as well as a realistic model of tissue was proposed based on microscopy images. Image segmentation was used to identify cells in an irradiated cell culture monolayer, or slices of tissue. The cells were then inserted into the model box pixel by pixel. In the case of cell monolayers (2-d), the image size may exceed the modeled box size. Such image was is moved with respect to the box in order to sample as many cells as possible. In the case of the simple tissue (3-d), the tissue box is modeled with periodic boundary conditions, which extrapolate the technique to macroscopic volumes of tissue. For real tissue, specific spatial patterns for cell apoptosis and necrosis are expected. The cell patterns were modeled based on action cross sections for apoptosis and necrosis estimated based on BNL data, and other experimental data.

Ponomarev, Artem L.↗

TPSAS-NF1676L-13635-DND

Objectives of AePW: Assess state-of-the-art Computational Aeroelasticity (CAe) methods as practical tools for the prediction of static and dynamic aeroelastic phenomena and responses on relevant geometries. - Perform comparative computational studies on selected test cases - Identify errors & uncertainties in computational aeroelastic methods - Identify gaps in existing aeroelastic databases - Provide roadmap of path forward (additional existing data sets, new experimental data sets, and analytical methods developments).

Walt Silva↗

A Case Study of Multimodal, Multi-institutional Data Management for the Combinatorial Materials Science Community

Although the convergence of high-performance computing, automation, and machine learning has significantly altered the materials design timeline, transformative advances in functional materials and acceleration of their design will require addressing the deficiencies that currently exist in materials informatics, particularly a lack of standardized experimental data management. The challenges associated with experimental data management are especially true for combinatorial materials science, where advancements in automation of experimental workflows have produced datasets that are often too large and too complex for human reasoning. The data management challenge is further compounded by the multimodal and multi-institutional nature of these datasets, as they tend to be distributed across multiple institutions and can vary substantially in format, size, and content. Furthermore, modern materials engineering requires the tuning of not only composition but also of phase and microstructure to elucidate processing–structure–property–performance relationships. To adequately map a materials design space from such datasets, an ideal materials data infrastructure would contain data and metadata describing (i) synthesis and processing conditions, (ii) characterization results, and (iii) property and performance measurements. In this work, we present a case study for the low-barrier development of such a dashboard that enables standardized organization, analysis, and visualization of a large data lake consisting of combinatorial datasets of synthesis and processing conditions, X-ray diffraction patterns, and materials property measurements generated at several different institutions. While this dashboard was developed specifically for data-driven thermoelectric materials discovery, we envision the adaptation of this prototype to other materials applications, and, more ambitiously, future integration into an all-encompassing materials data management infrastructure.

36 MATERIALS SCIENCE↗

Transonic Drag Prediction on a DLR-F6 Transport Configuration Using Unstructured Grid Solvers

A second international AIAA Drag Prediction Workshop (DPW-II) was organized and held in Orlando Florida on June 21-22, 2003. The primary purpose was to inves- tigate the code-to-code uncertainty. address the sensitivity of the drag prediction to grid size and quantify the uncertainty in predicting nacelle/pylon drag increments at a transonic cruise condition. This paper presents an in-depth analysis of the DPW-II computational results from three state-of-the-art unstructured grid Navier-Stokes flow solvers exercised on similar families of tetrahedral grids. The flow solvers are USM3D - a tetrahedral cell-centered upwind solver. FUN3D - a tetrahedral node-centered upwind solver, and NSU3D - a general element node-centered central-differenced solver. For the wingbody, the total drag predicted for a constant-lift transonic cruise condition showed a decrease in code-to-code variation with grid refinement as expected. For the same flight condition, the wing/body/nacelle/pylon total drag and the nacelle/pylon drag increment predicted showed an increase in code-to-code variation with grid refinement. Although the range in total drag for the wingbody fine grids was only 5 counts, a code-to-code comparison of surface pressures and surface restricted streamlines indicated that the three solvers were not all converging to the same flow solutions- different shock locations and separation patterns were evident. Similarly, the wing/body/nacelle/pylon solutions did not appear to be converging to the same flow solutions. Overall, grid refinement did not consistently improve the correlation with experimental data for either the wingbody or the wing/body/nacelle pylon configuration. Although the absolute values of total drag predicted by two of the solvers for the medium and fine grids did not compare well with the experiment, the incremental drag predictions were within plus or minus 3 counts of the experimental data. The correlation with experimental incremental drag was not significantly changed by specifying transition. Although the sources of code-to-code variation in force and moment predictions for the three unstructured grid codes have not yet been identified, the current study reinforces the necessity of applying multiple codes to the same application to assess uncertainty.

Lee-Rausch, E. M.↗

Optimal experimental design for prediction based on push-forward probability measures

Incorporating experimental data is essential for increasing the credibility of simulation-aided decision making and design. This paper presents a method which uses a computational model to guide the optimal acquisition of experimental data to produce data-informed predictions of quantities of interest (QoI). Many strategies for optimal experimental design (OED) select data that maximize some utility that measures the reduction in uncertainty of uncertain model parameters, for example the expected information gain between prior and posterior distributions of these parameters. In this paper, we seek to maximize the expected information gained from the pushforward of an initial (prior) density to the push-forward of the updated (posterior) density through the parameter-to-prediction map. The formulation presented is based upon the solution of a specific class of stochastic inverse problems which seeks a probability density that is consistent with the model and the data in the sense that the push-forward of this density through the parameter-to-observable map matches a target density on the observable data. While this stochastic inverse problem forms the mathematical basis for our approach, we develop a one-step algorithm, focused on push-forward probability measures, that leverages inference-for-prediction to bypass constructing the solution to the stochastic inverse problem. A number of numerical results are presented to demonstrate the utility of this optimal experimental design for prediction and facilitate comparison of our approach with traditional OED.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comprehensive Analysis of Streaming and Shutdown Dose Rate Experiments at JET with ORNL Fusion Neutronics Workflows

Current experimental fusion systems and conceptual designs of fusion pilot plants (FPPs) are growing in complexity and size. Several radiation metrics are crucial to the safe operation of fusion machines, including neutron flux streaming through openings and the shutdown dose rate (SDDR). Most current designs of advanced experimental fusion systems—and the most probable candidates for FPPs—are based on the tokamak concept, which is prone to neutron streaming through the myriad openings needed for diagnostic and support systems. SDDR is caused by decay gamma rays from radionuclides that become activated by neutrons during the operation of a fusion system that use deuterium-deuterium (DD), tritium-tritium, or deuterium-tritium plasma. Because computational tools have become essential for determining these radiation metrics, they must be validated against reliable and applicable experimental data. Experiments at the Joint European Torus (JET) provide a unique source of experimental data for validating computational tools and nuclear data used to determine SDDR and neutron fluxes in streaming-dominated geometries. Here, this paper presents the comprehensive analysis of the high-performance DD JET SDDR, and streaming experiments performed using Oak Ridge National Laboratory (ORNL) fusion workflows. The computational results were compared with experimental results that consist of online SDDR measurements with ionization chambers and neutron fluence streaming measurements using thermoluminescent detectors. The ratio of calculated-to-experimental SDDR values ranges from 0.6 to 2.5, and the streaming results range from 0.5 to 8.0. Future work will include analyzing the JET 2021 DTE2 campaign alongside the integration of the Shift Monte Carlo transport code into all ORNL fusion neutronics workflows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experiment and analysis of directional effects on radiant heat transfer.

Critical examination of the validity of the commonly used simplified models for predicting local and overall radiant heat interchange among real surfaces. The spectral local irradiation was measured at wavelengths of 3.08 and 4.51 microns in a configuration consisting of three plane parallel surfaces of finite extent. The test surfaces were gold with mechanical roughness and temperature varying from 0.02 to 7.1 microns and 77 to 760 K, respectively. Comparisons between experimental data and predictions for six different models ranging from the simple diffuse to the most detailed diffuse-plus-specular-directional-property model with specular component calculated according to Beckmann are presented. It is shown that the calculation of the overall irradiation using appropriate constant-property models agrees well with the experimental data to within combined experimental and analytical uncertainty. In general, it is concluded that for accurate prediction of local irradiation the directional characteristics of the surfaces must be considered.

Toor, J. S.↗

Verification of an Analytical Method for Measuring Crystal Nucleation Rates in Glasses from DTA Data

A recently proposed analytical (DTA) method for estimating the nucleation rates in glasses has been evaluated by comparing experimental data with numerically computed nucleation rates for a model lithium disilicate glass. The time and temperature dependent nucleation rates were predicted using the model and compared with those values from an analysis of numerically calculated DTA curves. The validity of the numerical approach was demonstrated earlier by a comparison with experimental data. The excellent agreement between the nucleation rates from the model calculations and fiom the computer generated DTA data demonstrates the validity of the proposed analytical DTA method.

Ranasinghe, K. S.↗

Development and preliminary results of 270 GHz microwave forward scattering diagnostic system on the experimental advanced superconducting tokamak (EAST)

To measure localized (intermediate, high) poloidal wavenumber electron density fluctuations, a tangential millimeter-wave collective scattering system has been designed and successfully developed on the experimental advanced superconducting tokamak (EAST). This innovative system employs a 270 GHz mm-wave probe beam, emitted from the mid-plane of Port K and directed towards Port P after passing through a combination of two in-vessel mirrors. Here these two ports are located on the EAST device’s midplane and are 110° apart. The downward scattered signals pass through two in-vessel mirrors, and exit through the P window. The received multi-channel separated scattering signals are arranged along the poloidal direction. This optical arrangement makes the monitored turbulent wave number is poloidal dominant and up to 40 cm −1 . The tangential microwave scattering scheme enables this diagnostic to monitor local density fluctuations, with typically a radial spatial resolution △ R = 5 cm, and a wavenumber resolution of △ k = 0.4 cm −1 . The scattering system underwent comprehensive laboratory testing in 2023, with installation finalized in 2024. The system features steerable launch and receiver optics, allowing remote control to position the scattering volume from normalized radius ρ –0 to the pedestal region on a shot-by-shot basis. Effective experimental data, characterized by distinct off-center spectral peaks, have been successfully obtained in neutral beam-heated plasmas on EAST. The beam tracing forward modeling has been employed for experimental data interpretation. The advanced millimeter-wave scattering system provides fluctuation measurement capability from ion temperature gradient to electron temperature gradient scale on low field side. Preliminary experimental density fluctuation data have been successfully obtained in neutral beam-heated plasmas on EAST. The localized measurement will be used for kinetic turbulence transport numerical simulation validation.

high-k collective scattering↗

Prediction of fatigue-crack growth in a high-strength aluminum alloy under variable-amplitude loading

The present paper is concerned with the application of an analytical crack-closure model to study crack growth under various load histories. The model was based on a crack-tip plasticity concept like the Dugdale model, but modified to leave plastically deformed material in the wake of the advancing crack tip. The effect of material thickness on plasticity was accounted for by using a constraint factor on tensile yielding at the crack tip. The model was used to correlate crack-growth rates under constant-amplitude loading, and to predict crack growth under variable-amplitude loading on a high-strength aluminum alloy (7475-T7351) sheet material. The experimental data were obtained from Zhang et al. Predicted crack-growth lives agreed well with experimental data. For ten crack-growth tests subjected to various variable-amplitude load histories, the ratio of predicted-to-experimental lives ranged from 0.54 to 1.19. The mean value of the ratio of predicted-to-experimental lives was 0.95 and the standard error was 0.2 using a constraint factor of 1.9 in the model. Crack-opening stresses calculated from the model were significantly different from those determined by Zhang et al. using a striation-based experimental method.

Newman, J. C., Jr.↗

Prediction of fatigue-crack growth in a high-strength aluminum alloy under variable-amplitude loading

The present paper is concerned with the application of an analytical crack-closure model to study crack growth under various load histories. The model was based on a concept like the Dugdale model, but modified to leave plastically deformed material in the wake of the advancing crack tip. The thickness effect was accounted for by using a 'constraint' factor on tensile yielding at the crack tip. The model was used to correlate crack-growth rates under constant-amplitude loading, and to predict crack growth under variable-amplitude loading on a high-strength aluminum alloy (7475-T7351) sheet material. The experimental data was obtained from Zhang et al. Predicted crack-growth lives agreed well with experimental data. For ten crack-growth tests subjected to various variable-amplitude load histories, the ratio of predicted-to-experimental lives ranged from 0.54 to 1.19. The mean value of predicted-to-experimental lives was 0.95 with a standard error of 0.2 for a constraint factor of 1.9.

Newman, J. C., Jr.↗

High Temperature Erosion Modeling in Particle Based CSP Systems

1. Introduction. Wear and erosion damage of materials from solid particle and surface interactions is a major issue in various industries. Although more common in mining and oil and gas production, erosion is becoming a critical issue in renewable technologies as well such as particle based concentrated solar thermal power (CSP) systems. In particle based CSP systems, solid particles are used to absorb solar energy and as thermal storage. However, these particles may cause significant amount of wear to system components while through the system at high temperatures. This damage can be costly, and therefore, requires a greater understanding of solid particle erosion in CSP systems. Throughout the years, models and tools have been developed to predict and control erosion in industries such as oil and gas production. However, these erosion models and erosion prediction tools have been mainly developed based on erosion data for much higher velocities and lower temperatures, i.e. the operating conditions for which the existing models have been built are not pertinent to those expected in CSP system. It is known that erosion depends on many factors such as material properties, erodent particle properties, and particle impact speed and angle. However, the effect of temperature on erosion is not vastly investigated. In this work, a temperature-based correlation is introduced that will be used to modify the existing erosion models to predict erosion rates at velocities and temperatures relevant to Gen3 CSP systems. The existing models are first validated against erosion experiments run at low temperatures and low velocities. Subsequently, a temperature correction term is developed that can extend the existing models to high temperatures, based on the available experimental data showing the effect of temperature on impact erosion. 2. Erosion Models. Continuing sequence of impacts from solid particles on surfaces would result in loss of material due to mechanical interaction between solid surface and particles. Erosion equations have been developed to predict erosion under different conditions and erosion mechanism, including cutting and deformation erosion. These models are mainly either mechanistic, empirical, and semi-mechanistic models, when the latter combine the theory of the erosion mechanism and particle motion with the available experimental data. One of the first empirical correlation was introduced by American Petroleum Institute (API) Recommended Practice (RP) 14E [1]. This correlation was very conservative, and several improved empirical and semi-mechanistic equations were developed in the years following. More recently, the models introduced by Erosion/Corrosion Research Center (E/CRC) at the University of Tulsa and Arabnejad. et al. [2] are frequently used in the literature and industry, as they account for many parameters affecting erosion including particle impact speed and angle, material density and hardness, and particle size and shape. However, both these models were developed based on data at relatively high velocities and also do not take into account any mechanical changes in material associated with thermal cycling at high temperatures. 3. Erosion Prediction at High Temperatures. In this work, the accuracy of Arabnejad et al. model and E/CRC model is investigated to predict impact erosion at conditions relevant to Gen3 CSP systems. Two significant deviations expected in Gen3 systems compared to operating conditions used to build these models are particle impact velocities and system temperature. As a first step, the performance of these models was validated at low velocities and temperatures. Impact erosion experiments were conducted on SS316 coupons using HSP 40/70 ceramic particles. Assuming particle rate of 1 kg/s/m, 0.0254 m of particle curtain thickness, and particle velocity of 1.5 m/s in the system, an overall erosion of 0.59 mm/year was calculated. The experimental results were subsequently compared to those from computational simulations and erosion of 0.47 mm/year was obtained using the Arabnejad et al. model. The results indicated that the models work well at low velocities and low temperature conditions. To account for temperature effect, a mathematical correlation was developed using data published by DUCOM [3] for Inconel 600 eroded by alumina particles at high velocities. The correlation was applied to both Arabnejad et. al model and the E/CRC model. The prediction results from these modified models were within 20% of this experimental data. Predictions of of erosion by the modified E/CRC model at three different temperatures are made. Similar to the calculation at low temperature, a typical CSP system with 1 kg/s/m of particle flow rate per unit length of the particle curtain and a curtain thickness of 0.0254 m is assumed. The annual thickness loss calculations were run for different particle impact velocities, assuming a uniform particle impact area equal to the cross-sectional area of the curtain (particle-particle interactions and dispersion of particles are not considered in obtaining the results). It is observed that erosion increases exponentially as temperature increases. Furthermore, it is also noted that, erosion changes non-linearly with impact velocity. Currently, experiments are also being conducted to measure erosion of SS316 at 800 ºC with HSP 40/70 particles at low impact velocities. We expect to use the results from high temperature testing to further improve the temperature correlation function. Similar models are also being developed for abrasion erosion resulting from particle sliding along the surfaces as well as attrition from particle to particle and particle to surface interactions. References. [1] Institute, A. P. (1991). API Recommended Practice for Design and Installation of Offshore Production Platform Piping System, API RP 14E. [2] Arabnejad, H., Mansouri, A., Shirazi, S. A., and McLaury, B. S. (2015a). Development of mechanistic erosion equation for solid particles. Wear, 332–333, 1044–1050. http://doi.org/10.1016/j.wear.2015.01.031. [3] https://ducom.com/high-temperature-erosion-evaluating-sample-wear/

14 SOLAR ENERGY↗

Modeling the dynamics of phospholipids in the fluid phase of liposomes

Here, we present the derivation of a new model to describe neutron spin echo spectroscopy and quasi-elastic neutron scattering data on liposomes. We compare the new model with existing approaches and benchmark it with experimental data. The analysis indicates the importance of including all major contributions in the modeling of the intermediate scattering function. Simultaneous analysis of the experimental data on lipids with full contrast and tail contrast matched samples reveals highly confined lipid tail motion. A comparison of their dynamics demonstrates the statistical independence of tail-motion and height–height correlation of the membrane. A more detailed analysis indicates that the lipid tail relaxation is confined to a potential with cylindrical symmetry, in addition to the undulation and diffusive motion of the liposome. Despite substantial differences in the chemistry of the fatty acid tails, the observation indicates a universal behavior. The analysis of partially deuterated systems confirms the strong contribution of the lipid tail to the intermediate scattering function. Within the time range from 5 to 100 ns, the intermediate scattering function can be described by the height–height correlation function. The existence of the fast-localized tail motion and the contribution of slow translational diffusion of liposomes determine the intermediate scattering function for t < 5 ns and t > 100 ns, respectively. Taking into account the limited time window lowers the bending moduli by a factor of 1.3 (DOPC) to 2 (DMPC) compared to the full range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Simulation of the High Temperature Test Facility (HTTF) Core Using the 2D Ring Model with SAM

A two-dimensional ring model is developed with SAM to model the core of the High Temperature Test Facility (HTTF) at the system level. The ring model simplifies the complex structure of the HTTF core by converting the hexagonal rows of heaters and flow channels into layers of concentric annular rings. The ring model is first compared against a 3D-1D model where the solid structures are fully resolved in three dimensions while the fluid structures are modeled as one-dimensional flows. Comparison between the 3D-1D and the ring models shows that the latter can predict major parameters reasonably well under steady-state normal operating condition, but the heater temperatures are underpredicted. Adjustment is made to the effective thermal conductivity of the ceramic core of the ring model to improve the heater temperature predictions. Furthermore, the ring model is also used to simulate a transient Pressurized Conduction Cooldown (PCC) condition and is benchmarked with the experimental data from the HTTF Test PG-27. Good agreement is obtained between the experimental data and the predictions by the ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Zero power experiment for molten salt reactor - code verification and validation request

Molten Salt Reactors (MSRs) are one of the six perspective Generation IV nuclear technologies that should fulfil sustainability, economics, safety, and proliferation resistance goals. In recent years these concepts have attracted increased interest in research as well as in industry. In the context of developing new, innovative technology, a zero-power experiment would be an ideal option for improving the understanding of MSR systems under low-risk and low-cost conditions. The given research was focused on cross-verification of the nuclear codes for MSR modelling and simulations as required for the development of the first design of a potential facility and a draft of an experimental programme. The analysis included the evaluation of lattice and Monte Carlo codes. The study showed that all considered nuclear tools were in good agreement with each other. Although this does not mean the results will match the real experimental data, which is essential for the regulatory approval of a novel technology. Currently, experimental data is almost absent for MSR technology. Hence, a zero power experiment is needed to assure we can rely on modelling and simulation studies for further development of MSR technology and deliver real-world experiments to support the licensing of future MSRs. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An approximate calculation of the strong interaction on a transonic airfoil

This paper describes an interactive computational scheme that calculates the transonic flow field over an airfoil. The inviscid flow, the turbulent boundary layer in the weak interaction region, and the elliptic flow in the strong interaction region are computed separately, and the solutions are matched iteratively in a direct, second-order interaction scheme. The pressure distribution, obtained from a transonic, small-perturbation solution, is prescribed on the boundary layer and the Navier-Stokes equations and the ensuing displacement thicknesses are used to modify the inviscid flow. The best overall results are obtained by using a turbulence model that eliminates the eddy viscosity in the strong interaction region. Computed pressure distributions are compared with experimental data and with Navier-Stokes solutions. In weak interactions, the agreement with experimental data is very good, as it is also in strong interactions except for the pressures in the separated-flow region. There the agreement is qualitative only, and the pressure recovery is too high. The results in both cases are superior or equal to those of the Navier-Stokes solutions, and the present method is significantly faster.

Seginer, A.↗