SEARCH · Engineering Papers

Results for “dimensional analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Description of the Three-Dimensional Large-Scale Forcing Data from the 3D Constrained Variational Analysis (VARANAL3D)

This technical report introduces a Three-Dimensional Constrained Variational Analysis (3DCVA) (Tang and Zhang 2015) and its product of three-dimensional large-scale forcing data to drive single-column models (SCM), cloud-resolving models (CRM), and large-eddy simulation (LES) models, and to evaluate model results. The 3DCVA algorithm is an extension of the original 1D constrained variational analysis (1DCVA) (Zhang and Lin 1997, Zhang et al. 2001). The three-dimensional structure of the forcing data allows studies of spatial variation of the large-scale forcing fields and tests of physical parameterizations across scales. In the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility, the 3D forcing data are assigned the datastream name varanal3d. In this technical report, 3DCVA will be used to refer to the algorithm, while VARANAL3D will be used to refer to the data product.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Description of the Three-Dimensional Large-Scale Forcing Data from the 3D Constrained Variational Analysis (VARANAL3D)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Calibration of cohesive parameters for fracture test simulations of sub-sized bend specimens of AM steels

In this work, crack extensions in sub-sized side-grooved bend specimens of AM steels are examined with changing cohesive parameters. Two-dimensional FE analyses are first conducted to simulate crack extensions by nodal release. The simulation results indicate the near-tip maximum opening stress and separation work rate increase and then slightly decrease for the bend specimens. The hydrogen-charged bend specimens have lower near-tip maximum opening stresses and average separation work rates than those of the uncharged bend specimens. The changing near-tip maximum opening stresses and the separation work rates with increasing crack extension are determined as references for the changing cohesive strengths and energies, respectively. Two-dimensional FE analyses with different changing cohesive parameters are then performed. The simulation results with calibrated cohesive parameters can match well with the test results. The hydrogen-charged specimens have the lower changing cohesive strength vs crack growth curves compared with those of the uncharged specimens. The hydrogen-charged specimens have either lower or similar changing cohesive energy vs crack growth curves compared with those of the uncharged specimens. Finally, a three-dimensional FE analysis is conducted to simulate crack extension by nodal release in one bend specimen to examine the stress triaxialities in front of the growing crack fronts at different locations from the center plane to the planes near the side groove front with increasing crack extension to identify possible physical reasons for the changing cohesive parameters.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Evolution of global and local deformation in additively manufactured octet truss lattice structures

Additively manufactured lattice truss structures, often referred to as architected cellular materials, present significant advantages over conventional structures due to their unique characteristics such as high strength-to-weight ratios and surface area-to-volume ratios. These geometrically complex structures, however, come with concomitant challenges for qualification and inspection. In this study, compression testing interrupted with micro-computed tomography inspection was conducted to monitor the evolution of global and local deformation throughout the loading process of 304L stainless steel octet truss lattice structures. Both two- and three-dimensional image analysis techniques were leveraged to characterize geometric heterogeneities resulting from the laser powder bed fusion manufacturing process as well as track the structure throughout deformation. Variations from model-predicted behavior resulting from these heterogeneities are considered relative to the predicted and actual responses of the structures during compression to better understand, model, and predict the octet truss lattice structure compression response.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗

DOE OSTI2022

AIMSim : An accessible cheminformatics platform for similarity operations on chemicals datasets

The recent advances in deep learning, generative modeling, and statistical learning have ushered in a renewed interest in traditional cheminformatics tools and methods. Quantifying molecular similarity is essential in molecular generative modeling, exploratory molecular synthesis campaigns, and drug-discovery applications to assess how new molecules differ from existing ones. Further, most tools target advanced users and lack general implementations accessible to the larger community. In this work, we introduce Artificial Intelligence Molecular Similarity (AIMSim), an accessible cheminformatics platform for performing similarity operations on collections of molecules called molecular datasets. AIMSim provides a unified platform to perform similarity-based tasks on molecular datasets, such as diversity quantification, outlier and novelty analysis, clustering, dimensionality reduction, and inter-molecular comparisons. AIMSim implements all major binary similarity metrics and molecular fingerprints and is provided as a Python package that includes support for command-line use as well as a Graphical User Interface for code-free utilization with fully interactive plots.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A Decision-Relevant Factor-Fixing Framework: Application to Uncertainty Analysis of a High-Dimensional Water Quality Model

Factor Fixing (FF) is a common method for reducing the number of model parameters to lower computational cost. FF typically starts with distinguishing the insensitive parameters from the sensitive and pursues uncertainty quantification (UQ) on the resulting reduced-order model, fixing each insensitive parameter at a fixed value. There is a need, however, to expand such a common approach to consider the effects of decision choices in the FF-UQ procedure on metrics of interest. Therefore, to guide the use of FF and increase confidence in the resulting dimension-reduced model, we propose a new adaptive framework consisting of four principles: (a) re-parameterize the model first to reduce obvious non-identifiable parameter combinations, (b) focus on decision relevance especially with respect to errors in quantities of interest (QoI), (c) conduct adaptive evaluation and robustness assessment of errors in the QoI across FF choices as sample size increases, and (d) reconsider whether fixing is warranted. The framework is demonstrated on a spatially-distributed water quality model. The error in estimates of QoI caused by FF can be estimated using a Polynomial Chaos Expansion (PCE) surrogate model. Built with 70 model runs, the surrogate is computationally inexpensive to evaluate and can provide global sensitivity indices for free. For the selected catchment, just two factors may provide an acceptably accurate estimate of model uncertainty in the average annual load of Total Suspended Solids (TSS), suggesting that reducing the uncertainty in these two parameters is a priority for future work before undertaking further formal uncertainty quantification.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Reduced‐Order Probabilistic Emulation of Physics‐Based Ring Current Models: Application to RAM‐SCB Particle Flux

Abstract In this work, we address the computational challenge of large‐scale physics‐based simulation models for the ring current. Reduced computational cost allows for significantly faster than real‐time forecasting, enhancing our ability to predict and respond to dynamic changes in the ring current, valuable for space weather monitoring and mitigation efforts. Additionally, it can also be used for a comprehensive investigation of the system. Thus, we aim to create an emulator for the Ring current‐Atmosphere interactions Model with Self‐Consistent magnetic field (RAM‐SCB) particle flux that not only improves efficiency but also facilitates forecasting with reliable estimates of prediction uncertainties. The probabilistic emulator is built upon the methodology developed by Licata and Mehta (2023), https://doi.org/10.1029/2022sw003345 . A novel discrete sampling is used to identify 30 simulation periods over 20 years of solar and geomagnetic activity. Focusing on a subset of particle flux, we use Principal Component Analysis for dimensionality reduction and Long Short‐Term Memory (LSTM) neural networks to perform dynamic modeling. Hyperparameter space was explored extensively resulting in about 5% median symmetric accuracy across all data sets for one‐step dynamic prediction. Using a hierarchical ensemble of LSTMs, we have developed a reduced‐order probabilistic emulator (ROPE) tailored for time‐series forecasting of particle flux in the ring current. This ROPE offers accurate predictions of omnidirectional flux at a single energy with no pitch angle information, providing robust predictions on the test set with an error score below 11% and calibration scores under 8% with bias under 2% providing a significant speed up as compared to the full RAM‐SCB run.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Discrepancy in the Critical State Void Ratio of Poorly Graded Sand due to Shear Strain Localization

The critical state (CS) concept is a theoretical framework that models the constitutive behavior of soils, including sand and other granular materials. It supports the notion of a unique postfailure state, where the soil ultimately experiences continuous shearing with no change in the plastic volumetric strain. However, the published literature has frequently noted the nonconvergence of sand specimens with different initial densities to a unique CS in the compression plane due to many factors such as specimen fabric, particle morphology, breakage, and grain size distribution. This paper examines the CS for poorly graded (uniform) glass beads and 3 different types of silica sands using 50 conventional triaxial compression (CTC) experiments, 12 oedometer tests, and in situ synchrotron microcomputed tomography (SMT) scans for 10 CTC experiments. The results of the 50 CTC experiments revealed a diffused CS zone in the compression plane, which was further examined using the in situ SMT scans. A thorough three-dimensional image analysis of the SMT scans accurately quantified the evolution of the local void ratio (elocal) versus axial compression within zones of intensive shearing toward the center of the specimen. The evolution of the void ratio was also measured using the entire volume of the specimen (eglobal). At the CS, the elocal/eglobal ratio was assessed to be ~1.25 when a single shear band developed within the scanned specimens and ~1.1-1.15 for specimens that failed via external bulging that was internally manifested by the development of multiple shear bands. This finding suggests that the CS zone in the compression plane can be attributed to the common wrong consideration of eglobal evolution in lieu of elocal within the developing shear bands. Furthermore, the lack of shear band development in uniaxial compression has made the results of the oedometer test reliable in quantifying the CS parameters in the compression plane.

58 GEOSCIENCES↗

DOE OSTI2020

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving the linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the resistive DCON code [A. H. Glasser, Z. R. Wang, and J.-K. Park, Phys. Plasmas 23, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Additionally, the solution of the ideal outer region is critical to the asymptotic matching and is challenging in toroidal geometry due to the singular power series solution at the resonant surfaces. Thus, a systematic verification of the outer region Δ' matrix is made by reproducing the well-known analytical Δ' result in Furth et al. [Phys. Fluids 16, 1054–1063 (1073)] and by making a quantitative benchmark with the PEST3 code [A. Pletzer and R. L. Dewar, J. Plasma Phys. 45, 427–451 (1991)]. Finally, the reconstructed numerical solution of the resistive plasma response from the toroidal matching matrix is presented. Compared with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Two-photon polymerization printed lattices as support structures in multi-shell ICF targets: Platform development and initial assessment

Recent interest in fielding direct drive multi-shell targets on the NIF [K. Molvig et al., Phys. Rev. Lett. 116, 255003 (2016) and S. X. Hu et al., Phys. Rev. E 100, 063204 (2019)] has highlighted the need for a low density structure to support the inner shell(s) and to avoid energy loss in the acceleration and collision process. We have developed a two-shell platform to evaluate the use of low density two-photon polymerization (2PP) printed lattices as a support structure between the shells. Here, 2PP structures are an attractive option because they can be produced at densities as low as 5 mg/cc, which is ideal for multi-shell targets, and their 3D structure can be tailored to the user's needs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

DOE OSTI2022

Experimental observations of multiple modes of filamentary structures in the magnetized dusty plasma experiment (MDPX) device

Due to the small charge-to-mass ratio of dust particles, it is often necessary to use large magnetic fields of B ≥ 1 T, in order to observe the influence of magnetic forces in laboratory dusty plasmas. However, when experiments are performed at high magnetic fields in capacitively coupled, radio frequency discharges used for these dusty plasma experiments, the plasma is often observed to form filamentary structures between the electrodes that are aligned to the external magnetic field which disrupt the uniformity of the plasma and adversely impact some of our dusty plasma experiments. Recent experiments performed in the Magnetized Dusty Plasma Experiment device seek to identify and characterize these filamentary structures. This paper discusses the morphology of several distinct filamentary modes that are formed in low temperature argon plasmas. Three different spatial modes of filaments that appear to be strongly dependent on the magnetic field and neutral pressure are identified. There is strong evidence that each spatial mode has a threshold condition that is dependent on the ion Hall parameter. The criteria for the formation of the filaments are shown to be consistent with predictions of recent numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Ambient-temperature liquid jet targets for high-repetition-rate HED discovery science

High-power lasers can generate energetic particle beams and astrophysically relevant pressure and temperature states in the high-energy-density (HED) regime. Recently-commissioned high-repetition-rate (HRR) laser drivers are capable of producing these conditions at rates exceeding 1 Hz. However, experimental output from these systems is often limited by the difficulty of designing targets that match these repetition rates. To overcome this challenge, we have developed tungsten microfluidic nozzles, which produce a continuously replenishing jet that operates at flow speeds of approximately 10 m/s and can sustain shot frequencies up to 1 kHz. The ambient-temperature planar liquid jets produced by these nozzles can have thicknesses ranging from hundreds of nanometers to tens of micrometers. In this work, we illustrate the operational principle of the microfluidic nozzle and describe its implementation in a vacuum environment. Further, we provide evidence of successful laser-driven ion acceleration using this target and discuss the prospect of optimizing the ion acceleration performance through an in situ jet thickness scan. Future applications for the jet throughout HED science include shock compression and studies of strongly heated nonequilibrium plasmas. When fielded in concert with HRR-compatible laser, diagnostic, and active feedback technology, this target will facilitate advanced automated studies in HRR HED science, including machine learning-based optimization and high-dimensional statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

A windowed mean trajectory approximation for condensed phase dynamics

We propose a trajectory-based quasi-classical method for approximating dynamics in condensed phase systems. Building upon the previously developed optimized mean trajectory approximation that has been used to compute linear and nonlinear spectra, we borrow some ideas from filtering trajectory methods to obtain a novel semiclassical method for the dynamical propagation of density matrices. This new approximation is tested rigorously against standard multistate electronic models, spin-boson models, and models of the Fenna–Matthews–Olson complex. For dissipative systems, the current method is significantly better or as good as many other semiclassical methods available, especially at low temperatures and for off-diagonal density matrix elements, whereas for scattering models, the current method bears similar limitations as mean-field propagation schemes. All results are tested against the numerically exact hierarchical equations of motion method. In conclusion, the new method shows excellent agreement across various parameter regimes with numerically exact results, highlighting the robustness and accuracy of our approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Three-Dimensional Heat Flux and Thermal Analysis of Angled Tungsten Samples on DIII-D

ITER-grade tungsten and dispersoid-strengthened tungsten samples with the top surface angled at ~15° towards the incident plasma flux were exposed to 9 H-mode discharges with edge-localized modes (ELMs) in the lower divertor of DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Surface damage included cracking and flaking of material on the two samples farthest away from the plasma strike point, and significant melting of the two samples closest to the strike point. Heat flux and thermal analysis tools new to DIII-D have been applied to better understand this material response and to help optimize the exposure conditions for future experiments. SMITER field-line tracing simulations based on IRTV data and EFIT equilibria estimate an average inter-ELM perpendicular heat flux, 𝑞⊥,𝑖nter−𝐸LM , on the angled surfaces of 10.1 – 19.6 MW/m² for a majority of the 9 discharges, increasing to 15.6 – 24.5 MW/m² for the single, higher-power shot where samples melted. Fast camera data showed shallow intra-ELM melting and re-solidification, which transitioned to bulk inter-ELM melting with melt motion in the 𝐽⃗ 𝑥 𝐵⃗ direction. About 50% of the protruding volume of the most affected sample was displaced via melt-motion. SIERRA thermal modeling software was able to reproduce an onset time of melting consistent with fast camera data and final sample conditions, within < 200 ms. Maximum surface temperatures of 3122 K and 2787 K are estimated for the samples farthest away from the strike point, while the closest samples achieve melting at 4067 ms and 4750 ms into the ~5000 ms plasma exposure. A +10% increase in both the SMITER 𝑞⊥,𝑖nter−𝐸LM calculations and the estimated ELM heat loads 𝑞⊥, 𝐸LM was required to achieve this result, which is within the uncertainty of the diagnostic data but likely accounts for non-ideal geometry effects plus other physics uncertainties not included in this first iteration of modeling. This work provided valuable estimates of the 3D temperature evolution to help better understand the observed surface morphology and internal recrystallization of samples, which are discussed in detail in a complementary manuscript [1]. Benchmarking efforts with more diagnosed DIII-D experiments are underway to further refine the SMITER and SIERRA models for DiMES. Future use of these tools will enable researchers to precisely target heat flux exposure conditions in DIII-D to test, but not exceed, the thermomechanical limitations of novel plasma-facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Magnetic properties of γ-Fe 2 O 3 nanoparticles in a porous SiO 2 shell for drug delivery

A method is presented for synthesizing core–shell nanoparticles with a magnetic core and a porous shell suitable for drug delivery and other medical applications. The core contains multiple γ -Fe 2 O 3 nanoparticles (~15 nm) enclosed in a SiO 2 (~100–200 nm) matrix using either methyl (denoted TMOS- γ -Fe 2 O 3 ) or ethyl (TEOS- γ -Fe 2 O 3 ) template groups. Low-temperature Mössbauer spectroscopy showed that the magnetic nanoparticles have the maghemite structure, γ -Fe 2 O 3 , with all the vacancies in the octahedral sites. Saturation magnetization measurements revealed that the density of γ -Fe 2 O 3 was greater in the TMOS- γ -Fe 2 O 3 nanoparticles than TEOS- γ -Fe 2 O 3 nanoparticles, presumably because of the smaller methyl group. Magnetization measurements showed that the blocking temperature is around room temperature for the TMOS- γ -Fe 2 O 3 and around 250 K for the TEOS- γ -Fe 2 O 3 . Three dimensional topography analysis shows clearly that the magnetic nanoparticles are not only at the surface but have penetrated deep in the silica to form the core–shell structure.

Physics↗

DOE OSTI2020

Time and space resolved diagnostics for plasma thermal-chemical instability of fuel oxidation in nanosecond plasma discharges

Here, an instability in a nanosecond pulsed dielectric barrier discharge plasma occurring in methane–oxygen–argon mixtures is experimentally observed and measured by 1D time-resolved in situ electric field measurements. This instability, which seems to be created by the positive feedback between plasma kinetics and plasma-assisted low temperature fuel oxidation, is studied using electric field induced second harmonic generation and direct ICCD imaging. The rapid formation of streamers from an originally uniform discharge appears to be caused by the chemical kinetics of plasma-assisted low temperature methane oxidation, which may be resulting in a new type of plasma instability: a thermal-chemical instability. The results also revealed that the occurrence of this possible thermal-chemical instability in a reactive flow drastically changes the plasma properties by forming multiple secondary discharges and possibly leads to micron-sized non-uniform electric distributions. Single shot uncalibrated measurements of the electric field of the micron sized streamers appears to show much greater strengths than the average electric field. Furthermore, one-dimensional data analysis shows the positive feedback loop between the streamers and the low temperature plasma assisted oxidation chemistry in the plasma thermal-chemical instability. The present finding advances the understanding plasma instability growth and provides a new way to control plasma uniformity in plasma-assisted combustion and plasma fuel reforming.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗