SEARCH · Engineering Papers
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.
We report the synthesis of layered anisotropic semiconductor GeSe and GeSe 2 nanomaterials through low temperature and atmospheric pressure chemical vapor deposition using halide based precursors. The crystal phase is controlled by simply changing selenium vapor pressure.
A fabric-enriched continuum breakage mechanics (F-CBM) framework is developed to examine the relation between elastic and inelastic anisotropy in granular materials subjected to high pressure. For this purpose, emphasis is given to the influence of the anisotropic fabric on the plastic cap regime – that is, a macroscopic signature of the inelasticity of granular materials widely dependent on the particle crushability. Anisotropic behaviour prior to yielding is introduced through a symmetric second-order fabric tensor embedded in the expression of the elastic energy potential. It is shown that, thanks to its energy-based formulation, the F-CBM framework provides a platform to predict the rotation and distortion of the yield cap of granular materials as an outcome of anisotropic strain energy storage prior to grain crushing. Parametric analyses indicate that both the degree of anisotropy and the orientation of the fabric have an impact on the energy release resulting from grain crushing and contribute to the resulting stress–strain behaviour. It is shown that, despite the use of very few parameters, the proposed framework can accurately predict the anisotropic yielding and inelastic stress–strain response solely on the basis of the directional properties inferred from measurements of elastic anisotropy.
A closed form solution is presented that predicts the response of filament wound composite shells subjected to internal pressure. The material of the shell is assumed to be general cylindrically anisotropic. Based on the theory of cylindrical anisotropic elasticity coupled differential equations are developed using Lekhnitskii's stress function approach. Two composite systems, graphite/epoxy and glass/epoxy, are selected to demonstrate the influence of degree of material anisotropy and fiber orientations on the axial and induced twisting deformation. Detailed stress distributions of (45) off-axis unidirectional and (45/-45)s symmetric angle-ply fiber-reinforced laminated shells are shown to illustrate the effect of radius-to-thickness ratio and stacking sequence.
Development of an algorithm to perform the optimal sizing of buckling resistant, imperfect, anisotropic ring-stiffened cylinders subjected to axial compression, torsion, and internal pressure is presented. The enforcement of stability constraints is treated in a way that does not require any eigenvalue analysis. Assumption of criticality of these stability constraints during the optimal sizing of the cylinders produced designs that nevertheless satisfied all of the stress constraints as well as the stability constraints. Case studies performed using a combination of penalty function and feasible direction optimization methods indicate that the presence of the axisymmetric initial imperfection in the cylinder wall can significantly affect the optimal designs. Weight savings associated with the addition of two rings to the unstiffened cylinder and/or the addition of internal pressure is substantial when torsion makes up a significant fraction of the combined load state.
In situ x-ray synchrotron diffraction experiments were carried out on nickel-based high-strength superalloys under pressure to understand their deformation mechanism using a diamond anvil cell (DAC). Radial x-ray diffraction determines the room-temperature equations of state and yield strengths of binary Ni 3 Al alloy and 500 ppm boron-doped Ni 3 Al to pressures of 20 and 46 GPa, respectively. Crystallographic preferred orientations observed in these superalloys due to anisotropic stress field in DAC indicate the onset of plastic deformation. Inverse pole figure analysis reveals that the underlying deformation mechanisms change from an octahedral slip to a simultaneous activation of octahedral and cube slips upon doping with boron. The yield-strength values were found to increase with pressure and are comparable to those determined from axial diffraction experiments. The results indicate that the yield strength of Ni 3 Al:B is about 0.5 GPa higher (at pressures below 20 GPa) due to grain boundary strengthening by boron. It is shown that due to high elastic anisotropy of Ni 3 Al alloy, the yield-strength estimations from diffraction experiments strongly depend on the micromechanical model used to convert the measured elastic strains into stresses.
Controlling discretization error is a remaining challenge for computational fluid dynamics simulation. Grid adaptation is applied to reduce estimated discretization error in drag or pressure integral output functions. To enable application to high O(10(exp 7)) Reynolds number turbulent flows, a hybrid approach is utilized that freezes the near-wall boundary layer grids and adapts the grid away from the no slip boundaries. The hybrid approach is not applicable to problems with under resolved initial boundary layer grids, but is a powerful technique for problems with important off-body anisotropic features. Supersonic nozzle plume, turbulent flat plate, and shock-boundary layer interaction examples are presented with comparisons to experimental measurements of pressure and velocity. Adapted grids are produced that resolve off-body features in locations that are not known a priori.
CHESS is a direct geometry neutron spectrometer designed to detect and analyze weak signals intrinsic to small cross-sections (e.g., small mass, small magnetic moments or neutron absorbing materials). This instrument is optimized to enable unprecedented characterization of quantum materials, spin liquids, thermoelectric and battery materials, liquids, and soft matter. The ability to simultaneously measure dynamic processes over a wide energy range for very small samples will make CHESS the spectrometer of choice for the initial exploration of new materials. The broad dynamic range will also be well matched to measurements of relaxation processes and excitations in soft and biological matter. The 15 Hz repetition rate of STS enables use of multiple incident energies within a single source pulse, greatly expanding the information gained in a single experiment. An essential feature of CHESS is the capability for polarization analysis to separate nuclear from magnetic scattering or coherent from incoherent scattering in hydrogenous materials, and better understanding spin-anisotropic correlations. This instrument will employ advanced sample environments such as high-pressure cells, dilution refrigerators, high field cryo-magnets and polarization devices, as well as combinations of these, to solve problems at the forefront of materials research. CHESS will be one of the flagship spectrometers of the Second Target Station (STS), providing world leading capabilities.
The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.
Anisotropic porous media model for fluid motion in rectangular fuel cell cavities, analyzing pressure, velocity, stream functions and purge time
Four current research projects are summarized: (1) active control of rotor system dynamics; (2) attenuation of rotor vibration using controlled pressure hydrostatic bearings; (3) a new seal test facility for measuring isotropic and anisotropic linear rotordynamic characteristics; and (4) the use of rotordynamic instability thresholds to accurately measure bearing rotordynamic characteristics.
A recently developed K-tau model for near wall turbulent flows is applied to a variety of test cases. The turbulent flows considered include the incompressible flat plate boundary layer with adverse pressure gradients, incompressible flow past a backward facing step, and the supersonic flat plate boundary layer at zero pressure gradient. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.
Limit loads under plane stress and plane strain are found for a circumferentially reinforced elastoplastic ring subjected to interior pressure. These are used as bounds on an estimate of the failure pressure of a SiC/Ti test rig that is being fabricated and tested. The ring is to serve as a benchmark against which deformation and failure analysis methods can be assessed. An anisotropic perfect plasticity idealization of the SiC/Ti ring material is made and used in the limit load calculations. An estimate of the failure pressure of the NASA/PW benchmark test ring is given.
Limit loads under plane stress and plane strain are found for a circumferentially reinforced elastoplastic ring subjected to interior pressure. These are used as bounds on an estimate of the failure pressure of a SiC/Ti test rig that is being fabricated and tested. The ring is to serve as a benchmark against which deformation and failure analysis methods can be assessed. An anisotropic perfect plasticity idealization of the SiC/Ti ring material is made and used in the limit load calculations. An estimate of the failure pressure of the NASA/PW benchmark test ring is given.
A class of parametric instabilities of finite-amplitude, circularly polarized Alfven waves in a plasma with pressure anisotropy is studied by application of the CGL equations. A linear perturbation analysis is used to find the dispersion relation governing the instabilities, which is a fifth-order polynomial and is solved numerically. A large-amplitude, circularly polarized wave is unstable with respect to decay into three waves: one sound-like wave and two side-band Alfven-like waves. It is found that, in addition to the decay instability, two new instabilities that are absent in the framework of the MHD equations can occur, depending on the plasma parameters.
Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.
Recent seismic studies suggest an anisotropic Mantle Transition Zone (MTZ) in areas adjacent to subducted slabs. Wadsleyite is the main anisotropy contributor in the upper MTZ, therefore the interpretation of these seismic observations requires the knowledge of single-crystal elastic moduli (C ij s) and the deformation-induced lattice preferred orientation (LPO) of wadsleyite. Wadsleyite can host up to 3 wt% water in its crystal structure as point defects in the form of hydroxyl groups, however, the combined effect of water content, Fe content, and pressure on the C ij s of wadsleyite remains unclear. In this study, we measured the high-pressure single-crystal elasticity of a synthetic hydrous Fe-bearing wadsleyite (0.14(5) wt% water, Fe#=9.4, Fe 3+ /ΣFe=0.3) up to 18.2 (2) GPa. In combination with previous experimental data, we separated the effects of pressure, water, and Fe contents on the Cijs and intrinsic elastic anisotropy of wadsleyite. Overall, our results suggest that the intrinsic elastic anisotropy of wadsleyite decreases with pressure, water, and Fe contents. At 15 GPa, increasing the water content by 0.1 wt% or Fe# by 1 decreases the V P and V S anisotropy of wadsleyite by ~1.1-1.3%, and ~0.8-1.3% in average, respectively. Combined with the LPO determined in previous deformation experiments, we modeled the seismic anisotropy in the upper MTZ generated by a sub-vertical mantle flow near cold subducted slabs and a sub-horizontal mantle flow in the ambient mantle. In both scenarios, the LPO of wadsleyite leads to V SV (vertically polarized shear wave velocity) > V SH (horizontally polarized shear wave velocity). Our results suggest that wadsleyite may account for a weak anisotropic MTZ (<1%) on the global scale. Considering the fact that water decreases the elastic anisotropy but promotes LPO of wadsleyite, seismic anisotropy may not be a good water sensor.
Here, this work is motivated by the desire to devise an internal pressure test that can mimic a displacement-controlled loading scenario and demonstrate how to apply the multiaxial stress and strain data from the test to develop an elastic/plastic constitutive model for a thin-walled tubular component. This is achieved by conducting simultaneous measurements of tangential and axial strain during the pressure test and integrating these strain measures into a feedback loop with the pressure controller. It is shown how data from such a test can be used to develop a large mechanical property data set relevant to biaxial loading conditions. The data obtained have high confidence evidenced by their low variability and alignment with other literature studies. Additionally, data from these internal pressure tests combined with full-tube axial tensile tests allow for the derivation of the Hill anisotropic yield function. The developed Hill yield function is validated by comparing the plastic strain ratios from the full-tube tension tests and by comparing the predicted yield stress in the tangential direction with measured values from ring tension tests in a previous study.
Earth deformation signals caused by atmospheric pressure loading are detected in vertical position estimates at Global Positioning System (GPS) stations. Surface displacements due to changes in atmospheric pressure account for up to 24% of the total variance in the GPS height estimates. The detected loading signals are larger at higher latitudes where pressure variations are greatest; the largest effect is observed at Fairbanks, Alaska (latitude 65 deg), with a signal root mean square (RMS) of 5 mm. Out of 19 continuously operating GPS sites (with a mean of 281 daily solutions per site), 18 show a positive correlation between the GPS vertical estimates and the modeled loading displacements. Accounting for loading reduces the variance of the vertical station positions on 12 of the 19 sites investigated. Removing the modeled pressure loading from GPS determinations of baseline length for baselines longer than 6000 km reduces the variance on 73 of the 117 baselines investigated. The slight increase in variance for some of the sites and baselines is consistent with expected statistical fluctuations. The results from most stations are consistent with approximately 65% of the modeled pressure load being found in the GPS vertical position measurements. Removing an annual signal from both the measured heights and the modeled load time series leaves this value unchanged. The source of the remaining discrepancy between the modeled and observed loading signal may be the result of (1) anisotropic effects in the Earth's loading response, (2) errors in GPS estimates of tropospheric delay, (3) errors in the surface pressure data, or (4) annual signals in the time series of loading and station heights. In addition, we find that using site dependent coefficients, determined by fitting local pressure to the modeled radial displacements, reduces the variance of the measured station heights as well as or better than using the global convolution sum.