Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “qualitative methods”

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,009 records · Page 56

Phase transitions and dimensional cross-over in layered confined solids

The nature of solid phases and cross-over of order–disorder phase transitions from two-dimensional (2D) layers to three-dimensional (3D) bulk in confined atomic systems remain largely unexplained. To this end, we consider noble gases and aluminum confined between graphene sheets at different pressures and temperatures. Using crystal structure search methods and molecular dynamics based on machine-learned potentials with quantum-mechanical accuracy, we identify structures of multilayer confined solids that deviate from simple close packing. Upon heating, we find that confined 2D monolayers melt according to the two-step continuous Kosterlitz–Thouless–Halperin–Nelson–Young theory. However, multilayer solids transition continuously into an intermediate layered-hexatic phase before melting discontinuously into an isotropic liquid. This intermediate phase persists at least up to 12 layers studied here. This change can be qualitatively understood based on the cross-over from 2D topological defects toward 3D ones during melting as the number of layers increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Convergent close-coupling calculations of electron scattering on LiH

Cross sections for electron scattering on the ground state of LiH are calculated using the molecular convergent close-coupling method. The fixed-nuclei approximation is utilized, and calculations are performed at the mean internuclear separation of 3.06 𝑎 0 for projectile energies ranging from 0.1 to 500 eV. Here, a model-potential approach is utilized to treat the LiH molecule as a two-electron system, and good agreement is found with the literature for various properties of the structure model. Cross sections are presented for elastic scattering, excitation to several electronic states, and ionization. Comparison is made with existing results, and qualitative agreement is found with small close-coupling calculations at low projectile energies. Good agreement is found with first-order results at high projectile energies.

74 ATOMIC AND MOLECULAR PHYSICS↗

Deformation and burst of a liquid droplet freely suspended in a linear shear field

A theoretical method is presented for predicting the deformation and the conditions for breakup of a liquid droplet freely suspended in a general linear shear field. This is achieved by expanding the solution to the creeping-flow equations in powers of the deformation parameter epsilon and using linear stability theory to determine the onset of bursting. When compared with numerical solutions and with the available experimental data, the theoretical results are generally found to be of acceptable accuracy although, in some cases, the agreement is only qualitative.

Barthes-Biesel, D.↗

A review of candidate multilayer insulation systems for potential use on wet-launched LH2 tankage for the space exploration initiative lunar missions

The storage of cryogenic propellants such as liquid hydrogen (LH2) and liquid oxygen (LO2) for the future Space Exploration Initiative (SEI) will require lightweight, high performance thermal protection systems (TPS's). For the near-term lunar missions, the major weight element for most of the TPS's will be multilayer insulation (MLI) and/or the special structures/systems required to accommodate the MLI. Methods of applying MLI to LH2 tankage to avoid condensation or freezing of condensible gases such as nitrogen or oxygen while in the atmosphere are discussed. Because relatively thick layers of MLI will be required for storage times of a month or more, the transient performance from ground-hold to space-hold of the systems will become important in optimizing the TPS's for many of the missions. The ground-hold performance of several candidate systems are given as well as a qualitative assessment of the transient performance effects.

Knoll, Richard H.↗

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser pre-sintering for denudation reduction in the laser powder bed fusion additive manufacturing of Ti-6Al-4V alloy

As a fundamentally different process for creation of complex parts, laser powder bed fusion (LPBF) additive manufacturing exhibits complex behavior at the melt pool that has a direct influence on the macroscopic properties of the final material. Powder denudation, a previously studied deleterious phenomenon, is the displacement of powder due to entrainment in a jet of rapidly vaporized metal created by the process laser. To minimize the effect of denudation and inspired by methods in adjacent fields, a laser pre-sintering (LPS) method was studied in an effort to reduce denudation. LPS was found to be most effective in the 1–4 J/mm3 volumetric energy density processing regimes, minimizing the common experimentally observed problems with the LPBF process. Scanning electron microscopy of the LPS-treated powder bed revealed insight into the accelerating effect morphological irregularities along the powder particle surface have in the pre-sintering process. To further understand the process, an estimation of the sintered neck sizes necessary to eliminate denudation in LPS is made, qualitatively affirming the low energy densities needed to prevent denudation. Here, the results presented in this work represent an initial step towards understanding and implementing LPS as a quality-enhancing step in LPBF.

36 MATERIALS SCIENCE↗

Taylor wave solution for a general equation of state

This document describes a solution procedure for calculating the Taylor wave behind an unsupported Chapman–Jouguet (CJ) detonation in planar, cylindrical, and spherical geometries given a general equation of state. The resulting semi-analytic solution can be utilized to examine new equation of state models for detonation products and during the verification of hydrodynamic codes. The governing partial differential equations are reduced to ordinary differential equations in both characteristic and self-similar forms. The first-order systems corresponding to each geometry are amenable to solution numerically using commonly available methods. A difficulty arises at the CJ point in radial coordinates where the similarity equations become singular. Two separate strategies are proposed to integrate the first-order system. The first one uses an asymptotic approximation near the CJ point that can be used to perturb the boundary conditions. The second one applies a change of variables which removes the singularity at the expense of an additional equation to be integrated. A test problem is provided for the Davis products equation of state to illustrate the qualitative features of the Taylor wave in each geometric configuration and compared with a Lagrangian hydrodynamics research code. A Python code listing gives an implementation using the SciPy library to assists users in generating the results.

97 MATHEMATICS AND COMPUTING↗

Computed Tomography Scanning and Petrophysical Measurements of the Lively Grove #1 Well Core

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core from the Lively Grove #1 Well (API 121892494700), drilled near Marissa, Washington County, Illinois. Core from the well was obtained as part of the Illinois Storage Corridor’s Carbon Storage Assurance Facility Enterprise (CarbonSAFE) project (DE FE0031892). The primary impetus of this work was to capture a detailed digital representation of the core from the Lively Grove #1 Well. The collaboration between the U.S. Department of Energy’s (DOE) NETL and the Illinois State Geological Survey (ISGS) enables other research entities to access information about this potential carbon storage location and Cambro-Ordovician Storage Complex formations of the Illinois Basin. The resultant datasets are presented in this report and can be accessed from NETL’s Energy Data eXchange (EDX) online system. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. Fractures, discontinuities, and millimeter-scale features were readily detectable with the medical CT scanner-acquired images. Imaging with the NETL medical CT scanner was performed on entire cores. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), gamma density, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for potential future analysis. Higher-resolution industrial and micro-CT images were acquired from selected zones along the depth of the core to visualize the structure in higher detail. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; with the resulting macro- and micro-descriptions relevant to many subsurface energy related examinations traditionally performed at NETL.

54 ENVIRONMENTAL SCIENCES↗

Computed Tomography Scanning and Petrophysical Measurements of Eastern Williston Basin Twin Buttes and Hagel Formations

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) in Morgantown, West Virginia, were used to characterize core from two wells that represent coal resources across North Dakota. These include the MC23080C Well in Mercer County and the 23-B001 Well in Oliver County. The primary impetus of this work was to capture a detailed digital representation of the core from the MC23080C and 23-B001 Wells. The collaboration between the NETL and the Energy and Environment Research Center (EERC) enables other research entities to access information about this potential carbon ore, rare earth, and critical mineral resource plays in the Williston Basin. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. Fractures, discontinuities, and millimeter-scale features were readily detectable with the medical CT scanner acquired images. Imaging with the NETL medical CT scanner was performed on entire cores. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), gamma density, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for potential future analysis. Higher-resolution industrial and micro-CT images were acquired from selected zones along the depth of the core to visualize the structure in higher detail. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core, with the resulting macro- and micro-descriptions relevant to many subsurface energy-related examinations traditionally performed at NETL.

01 COAL, LIGNITE, AND PEAT↗

On the instabilities of supersonic mixing layers - A high-Mach-number asymptotic theory

The stability of a family of tanh mixing layers is studied at large Mach numbers using perturbation methods. It is found that the eigenfunction develops a multilayered structure, and the eigenvalue is obtained by solving a simplified version of the Rayleigh equation (with homogeneous boundary conditions) in one of these layers which lies in either of the external streams. This analysis leads to a simple hypersonic similarity law which explains how spatial and temporal phase speeds and growth rates scale with Mach number and temperature ratio. Comparisons are made with numerical results, and it is found that this similarity law provides a good qualitative guide for the behavior of the instability at high Mach numbers. In addition to this asymptotic theory, some fully numerical results are also presented (with no limitation on the Mach number) in order to explain the origin of the hypersonic modes (through mode splitting) and to discuss the role of oblique modes over a very wide range of Mach number and temperature ratio.

Balsa, Thomas F.↗

A systems approach to a DSMC calculation of a control jet interaction experiment

This study deals with the development of a methodology for numerically simulating the interaction of a reaction control system (RCS) jet with a low-density external flow. A European Space Agency (ESA) experiment was chosen as a test case since it provided experimental data that could validate some of the numerical results. The initial approach was to focus on several subproblems having direct relevance to the full interaction problem. This enabled different numerical methods to be separately investigated and validated for each part of the interaction problem. In this manner, the best methodology for solving the full interaction problem was developed. The subproblems considered in this study included typical RCS nozzle and plume flows, a flat plate at zero incidence, and the flow past the experimental test model without the control jet firing. Once these calculations were completed, a simulation was performed with a control jet operating at a relatively low density. The results from this latter simulation provided qualitative insight into the complex interaction process.

Tartabini, Paul V.↗

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗

Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms

This numerical study compares the wave field generated by the spectral wave action balance code, SNL-SWAN, to the linear-wave boundary-element method (BEM) code, WAMIT. The objective of this study is to assess the performance of SNL-SWAN for modeling wave field effects produced by individual wave energy converters (WECs) and wave farms comprising multiple WECs by comparing results from SNL-SWAN with those produced by the BEM code WAMIT. BEM codes better model the physics of wave-body interactions and thus simulate a more accurate near-field wave field than spectral codes. In SNL-SWAN, the wave field’s energy extraction is modeled parametrically based on the WEC’s power curve. The comparison between SNL-SWAN and WAMIT is made over a range of incident wave conditions, including short-, medium-, and long-wavelength waves with various amounts of directional spreading, and for three WEC archetypes: a point absorber (PA), a pitching flap (PF) terminator, and a hinged raft (HR) attenuator. Individual WECs and wave farms of five WECs in various configuration were studied with qualitative comparisons made of wave height and spectra at specific locations, and quantitative comparisons of the wave fields over circular arcs around the WECs as a function of radial distance. Results from this numerical study demonstrate that in the near-field, the difference between SNL-SWAN and WAMIT is relatively large (between 20% and 50%), but in the far-field from the array the differences are minimal (between 1% and 5%). The resultant wave field generated by the two different numerical approaches is highly dependent on parameters such as: directional wave spreading, wave reflection or scattering, and the WEC’s power curve.

16 TIDAL AND WAVE POWER↗

Swept shock/boundary-layer interactions: Scaling laws, flowfield structure, and experimental methods

A general review is given of several decades of research on the scaling laws and flowfield structures of swept shock wave/turbulent boundary layer interactions. Attention is further restricted to the experimental study and physical understanding of the steady-state aspects of these flows. The interaction produced by a sharp, upright fin mounted on a flat plate is taken as an archetype. An overall framework of quasiconical symmetry describing such interactions is first developed. Boundary-layer separation, the interaction footprint, Mach number scaling, and Reynolds number scaling are then considered, followed by a discussion of the quasiconical similarity of interactions produced by geometrically-dissimilar shock generators. The detailed structure of these interaction flowfields is next reviewed, and is illustrated by both qualitative visualizations and quantitative flow images in the quasiconical framework. Finally, the experimental techniques used to investigate such flows are reviewed, with emphasis on modern non-intrusive optical flow diagnostics.

Settles, Gary S.↗

A helical actuator driven by biased SMA: design, model, and experiment

A helical actuator driven by biased shape memory alloy (SMA) patterns embedded into a soft composite ribbon base is presented in this work. Instead of common U-shape SMA wires, a single SMA wire is woven into planar patterns, which enable helical deformation of the composite ribbon from an initially flat geometry. An analytical static model is established for accurate and rapid prediction of the helical reconfiguration arising from the shape memory effect of woven SMA patterns, followed by validation of the static model using the finite element method (FEM). The finite element results are compared with the analytical solutions given by this static model, which show a high agreement. Parametric study of the influences of eight independent design variables on the dependent helical parameters, such as combined curvatures, pitches, and helical angles, is completed. It is found that the helically deformed geometry is mainly dominated by diameters, biased positions, inclined angles, and numbers of skewed segments of the SMA wire. Fabrication and in-situ experimental test of a prototype of such helical actuators qualitatively demonstrate its dramatic three-dimensional (3D) spiral reconfiguration from a two-dimensional (2D) flat ribbon. In conclusion, such SMA patterns will allow more diverse designs of soft actuators for a wider range of robotic applications.

36 MATERIALS SCIENCE↗

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗

A multi-dimensional dataset for two-phase instability in low pressure natural circulation based on direct transient local measurement

Thermal hydraulic instability can challenge the control and the safety of nuclear reactors. This experimental research captures and investigates the instability in a low-pressure (1–10 atmospheres) natural circulation loop with an adiabatic chimney installed above the heated section. Direct transient local measurement of two-phase flow is performed at multiple axial locations along the flow channel with thermocouples and conductivity probes traversed by customized high-precision mechanisms. Based on the periodicity of oscillations, a post-processing method is proposed to align the data captured non-simultaneously to reconstruct the instability behavior in one representative period. The averaging operation over multiple oscillation cycles is performed to reduce the random errors. Cyclic phenomena are clearly demonstrated by the post-processed data, and rapid void fraction change induced by the flashing is observed to be significant in the driving mechanism for low pressure oscillations. Parametric effects of inlet subcooling, heat flux, inlet flow restriction and system pressure are studied within the experiment range, and qualitative results are explained and compared with existing literature. The ratio between the oscillation period and the travelling time of the fluid particle through the test section is compared with the theoretical value for the density wave oscillation. The new dataset contains directly measured transient multi-dimensional two-phase data along the flow channel for natural circulation under density wave oscillations. The generated dataset is valuable for the validation of system analysis codes and CFD codes in the simulation of natural circulation phenomena and dynamic two-phase flow.

42 ENGINEERING↗

Microstructure of irradiated AGR TRISO particle buffer layers as measured by X-ray computed tomography

Shrinkage of the initially low-density buffer layer in tristructural isotropic (TRISO) coated fuel particles during irradiation is a well-known phenomenon with potential implications for fission product and actinide transport, as well as potential pyrocarbon fracture that in some cases has been observed to impact particle performance. During post-irradiation examination, the buffer layer's structure is commonly determined using 2D microscopy of a particle cross section or a series of particle cross sections at staggered depths. Although these methods provide a general idea of the irradiated buffer microstructure, they do not provide a full picture of the TRISO coating microstructure. By contrast, x-ray computed tomography (XCT) provides full 3D imaging of the TRISO particle. Particles from various compacts irradiated for the Advanced Gas Reactor Fuel Development and Qualification Program were imaged using XCT. Some of these particles were selected because of abnormal fission product inventories (e.g., low 137 Cs), whereas others were randomly selected from the center of the fission product inventory distribution. Qualitative and quantitative analysis techniques were applied to the randomly selected particles as representatives of typical TRISO behavior to study the post-irradiation structure of the buffer layer. These results showed that, while separation of the buffer and inner pyrolytic carbon layers was a common behavior during the radiation-induced shrinkage of the buffer, a portion of the original buffer/inner pyrolytic carbon interface remained intact throughout irradiation in nearly all cases. Furthermore, these results also indicated clear trends in the degree of buffer densification with irradiation temperature and fluence.

36 MATERIALS SCIENCE↗