Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anisotropic pressure”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Isothermal magnetostatic atmospheres. III - Similarity solutions with current proportional to the magnetic potential squared

A family of isothermal magnetostatic atmospheres with one ignorable coordinate corresponding to a uniform gravitational field in a plane geometry is considered. It is assumed that the current (J) is proportional to the square of the magnetostatic potential and falls off exponentially with distance. Results are presented for the contributions of the anisotropic J x B force (where B is the magnetic field induction), the gravitational force, and the gas pressure gradient to the force balance.

Webb, G. M.↗
NASA NTRS2012

Integrated Structural/Acoustic Modeling of Heterogeneous Panels

A model for the dynamic response of heterogeneous media is presented. A given medium is discretized into a number of subvolumes, each of which may contain an elastic anisotropic material, void, or fluid, and time-dependent boundary conditions are applied to simulate impact or incident pressure waves. The full time-dependent displacement and stress response throughout the medium is then determined via an explicit solution procedure. The model is applied to simulate the coupled structural/acoustic response of foam core sandwich panels as well as aluminum panels with foam inserts. Emphasis is placed on the acoustic absorption performance of the panels versus weight and the effects of the arrangement of the materials and incident wave frequency.

Bednarcyk, Brett, A.↗
DOE OSTI2020

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Velocity-changing collisional effects in nonlinear atomic spectroscopy and photon echo decay in gases

A general theory of atomic dipole coherence under the influence of collisional phase changes, inelastic effects and optically active atom velocity changes, including those due to anisotropic interactions is presented. Velocity change effects are obtained in closed form. Line shapes appear as convolutions of standard pressure broadening contours with velocity-change contours. Width and shift parameters for the He-broadened Na D lines at 2 m bar pressure, 380 K are calculated, as are He-induced photon echo decay rates for these lines. Overall agreement with xperiment is reasonably good.

Herman, R. M.↗
NASA NTRS1990

Three dimensional LDV flow measurements and theoretical investigation in a radial inflow turbine scroll

A two-color LDV system was used in the measurement of three orthogonal velocity components at 758 points located throughout the scroll and the unvaned portion of the nozzle of a radial inflow turbine scroll. The cold flow experimental results are presented for the velocity field at the scroll tongue. In addition, a total pressure loss of 3.5 percent for the scroll is revealed from the velocity measurements combined with the static pressure readings. Moreover, the measurement of the three normal stresses of the turbulence has showed that the flow is anisotropic. Furthermore, the mean velocity components are compared with a numerical solution of the potential flow field using the finite element technique. The theoretical prediction of the exit flow angle variation agrees well with the experimental results. This variation leads to a higher scroll pattern factor which can be avoided by controlling the scroll cross sectional area distribution.

Malak, Malak Fouad↗
NASA NTRSDate not supplied

Evaluation of Wall-Modeled Les for Flow Over A Multi-Element Airfoil

We perform wall-modeled large-eddy simulations (WMLES) of turbulent flow over a 30P30N multielement airfoil at various angles of attack to assess its ability to accurately predict lift using an unstructured node-based finite-volume methodology. The flow conditions are based on the experiments reported in Klausmeyer & Lin [NASA /TM-112858]. The oncoming flow Mach number is 0.2, and the stowed chord-based Reynolds number (𝑅𝑒𝑐) is 9 million. The angle of attack (𝛼) is varied between 8◦ and 23◦ with the expected stall angle of around 21◦. The simulated geometry is periodic in the span, thus any three-dimensional effects present in the experiment are ignored. We first analyze the predictions with baseline isotropic and anisotropic grids containing over 16 grid points per boundary-layer thickness. The near-wall streamwise/nominal spacing ratio is 0.4 at the leading/trailing edges of the slat, main element and flap, and 1 elsewhere, and the first wall-normal/nominal spacing ratio is 0.375 over most of the airfoil for the anisotropic grid, while these are unity for the isotropic grid. It is found that the predicted lift and wall pressures, near-wall velocity profiles, and the location of onset of resolved turbulence are sensitive to the near-wall grid anisotropy. Further sensitivities were explored at lower angles of attack by independently varying the streamwise and wall-normal grid anisotropies. The effects of both off-wall and near-wall grid refinement were also quantified. For the grid resolutions used here, our results suggest that the first wall-normal spacing, which coincides with the WMLES exchange location, has a dominant effect on the predictions for this flow configuration.

Computational Fluid Dynamics↗
DOE OSTI2021

Thermal conductivity tensor of β -1,3,5,7-tetranitro-1,3,5,7-tetrazoctane ( β -HMX) as a function of pressure and temperature

Here, we have used reverse non-equilibrium molecular dynamics (RNEMD) simulations to determine the full thermal conductivity tensor for the monoclinic high explosive crystal β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX). In order to do so for the monoclinic crystal, four directions for heat propagation are used. Effects of the temperature and pressure are investigated between 200 and 500 K and 0 and 5 GPa, respectively, which approximately covers the range where the β polymorph is stable. Simulations are carried out with the Smith–Bharadwaj non-reactive empirical potential [Smith and Bharadwaj, J. Phys. Chem. B 103, 3570 (1999)], which is known to reproduce well the thermo-elastic properties of HMX. Our results indicate that the thermal conductivity, κ, is highly anisotropic, with 36% difference between the two extreme values at 300 K and 0 GPa. A simple function is used to interpolate κ in the pressure-temperature regime considered in this study, which can be used in continuum models. The results from RNEMD simulations compare well with available experimental results from the literature and allow the determination of κ for any direction and temperature and pressure within or around the fitting interval.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Numerical computation of aerodynamics and heat transfer in a turbine cascade and a turn-around duct using advanced turbulence models

The objective of this research is to develop turbulence models to predict the flow and heat transfer fields dominated by the curvature effect such as those encountered in turbine cascades and turn-around ducts. A Navier-Stokes code has been developed using an explicit Runge-Kutta method with a two layer k-epsilon/ARSM (Algebraic Reynolds Stress Model), Chien's Low Reynolds Number (LRN) k-epsilon model and Coakley's LRN q-omega model. The near wall pressure strain correlation term was included in the ARSM. The formulation is applied to Favre-averaged N-S equations and no thin-layer approximations are made in either the mean flow or turbulence transport equations. Anisotropic scaling of artificial dissipation terms was used. Locally variable timestep was also used to improve convergence. Detailed comparisons were made between computations and data measured in a turbine cascade by Arts et al. at Von Karman Institute. The surface pressure distributions and wake profiles were predicted well by all the models. The blade heat transfer is predicted well by k-epsilon/ARSM model, as well as the k-epsilon model. It's found that the onset of boundary layer transition on both surfaces is highly dependent upon the level of local freestream turbulence intensity, which is strongly influenced by the streamline curvature. Detailed computation of the flow in the turn around duct has been carried out and validated against the data by Monson as well as Sandborn. The computed results at various streamwise locations both on the concave and convex sides are compared with flow and turbulence data including the separation zone on the inner well. The k-epsilon/ARSM model yielded relatively better results than the two-equation turbulence models. A detailed assessment of the turbulence models has been made with regard to their applicability to curved flows.

Lakshminarayana, B.↗
DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗
DOE OSTI2021

Decoupling Lattice and Magnetic Instabilities in Frustrated CuMnO2

The AMnO 2 delafossites (A = Na, Cu) are model frustrated antiferromagnets, with triangular layers of Mn 3+ spins. At low temperatures (T N = 65 K), a C2/m → P$\bar1$ transition is found in CuMnO 2 , which breaks frustration and establishes magnetic order. In contrast to this clean transition, A = Na only shows short-range distortions at T N . Here, we report a systematic crystallographic, spectroscopic, and theoretical investigation of CuMnO 2 . We show that, even in stoichiometric samples, nonzero anisotropic Cu displacements coexist with magnetic order. Using X-ray/neutron diffraction and Raman scattering, we show that high pressures act to decouple these degrees of freedom. This manifests as an isostuctural phase transition at ~10 GPa, with a reversible collapse of the c-axis. This is shown to be the high-pressure analogue of the c-axis negative thermal expansion seen at ambient pressure. Density functional theory (DFT) simulations confirm that dynamical instabilities of the Cu + cations and edge-shared MnO 6 layers are intertwined at ambient pressure. However, high pressure selectively activates the former, before an eventual predicted reemergence of magnetism at the highest pressures. Our results show that the lattice dynamics and local structure of CuMnO 2 are quantitatively different from nonmagnetic Cu delafossites and raise questions about the role of intrinsic inhomogeneity in frustrated antiferromagnets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Nontrivial nanostructure, stress relaxation mechanisms, and crystallography for pressure-induced Si-I → Si-II phase transformation

Crystallographic theory based on energy minimization suggests austenite-twinned martensite interfaces with specific orientation, which are confirmed experimentally for various materials. Pressure-induced phase transformation (PT) from semiconducting Si-I to metallic Si-II, due to very large and anisotropic transformation strain, may challenge this theory. Here, unexpected nanostructure evolution during Si-I → Si-II PT is revealed by combining molecular dynamics (MD), crystallographic theory, generalized for strained crystals, and in situ real-time Laue X-ray diffraction (XRD). Twinned Si-II, consisting of two martensitic variants, and unexpected nanobands, consisting of alternating strongly deformed and rotated residual Si-I and third variant of Si-II, form {111} interface with Si-I and produce almost self-accommodated nanostructure despite the large transformation volumetric strain of -0.237. The interfacial bands arrest the {111} interfaces, leading to repeating nucleation-growth-arrest process and to growth by propagating {110} interface, which (as well as {111} interface) do not appear in traditional crystallographic theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

New discretization and solution techniques for incompressible viscous flow problems

Several topics arising in the finite element solution of the incompressible Navier-Stokes equations are considered. Specifically, the question of choosing finite element velocity/pressure spaces is addressed, particularly from the viewpoint of achieving stable discretizations leading to convergent pressure approximations. The role of artificial viscosity in viscous flow calculations is studied, emphasizing work by several researchers for the anisotropic case. The last section treats the problem of solving the nonlinear systems of equations which arise from the discretization. Time marching methods and classical iterative techniques, as well as some modifications are mentioned.

Gunzburger, M. D.↗
NASA NTRS1989

Interpretation and modeling of the averaged equations for a fluid-solid flow

A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.

Shen, Hayley H.↗
NASA NTRS1994

Anisotropic Alfven-ballooning modes in Earth's magnetosphere

We have carried out a theoretical analysis of the stability and parallel structure of coupled shear Alfven and slow magnetosonic waves in Earth's inner magnetopause (i.e., at equatorial distances between about five and ten Earth radii) including effects of finite anisotropic Grad-Shafranov equation yields an approximate self-consistent magnetohydrodynamic (MHD) equilibrium. This MHD equilibrium is used in the numerical solution of a set of eigenmode equations which describe the field line eigenfrequency, linear stability, and parallel eigenmode structure. We call these modes anisotropic Alfven-ballooning modes. The main results are: (1) The field line eigenfrequency can be significantly lowered by finite pressure effects. (2) The parallel mode structure of the transverse wave components is fairly insensitive to changes in the plasma pressure, but the compressional magnetic component can become highly peaked near the magnetic equator as a result of increased pressure, especially when P(sub perpendicular to) is greater than P(sub parallel) (here P(sub perpendicular to) and P(sub parallel) are the perpendicular and parallel plasma pressure). (3) For the isotropic (P(sub parallel) = P(sub perpendicular to) = P) case ballooning instability can occur when the ratio of the plasma presure to the magnetic pressure, beta = P/(B squared/8 pi), exceeds a critical value beta(sup B)(sub 0) is approximately equal to 3.5 at the equator. (4) Compared to the isotropic case the critical beta value is lowered by anisotropy, either due to decreased field line bending stabilization when P(sub parallel) is greater than P(sub perpendicular to) or due to increased ballooning-mirror destabilization when P(sub perpendicular to) is greater than P(sub parallel). (5) We use a beta-delta stability diagram to display the regions of instability with respect to the equatorial values of the parameters bar beta and delta, where bar beta = (1/3)(beta(sub parallel) + 2 beta(sub perpendicular to)) is an average beta value and delta = 1 - P(sub parallel)/P(sub perpendicular to) is a measure of the plasma anisotropy. The diagram is divided into regions corresponding to the firehose, mirror and ballooning instabilities. It appears that observed values of the plasma pressure are below the critical value for the isotropic ballooning instability but it may be possible to approach a ballooning-mirror instability when P(sub perpendicular to)/P(sub parallel) is greater than or approximately 2.

Chan, Anthony A.↗
NASA NTRS2002

Generalized Reliability Methodology Applied to Brittle Anisotropic Single Crystals

A generalized reliability model was developed for use in the design of structural components made from brittle, homogeneous anisotropic materials such as single crystals. The model is based on the Weibull distribution and incorporates a variable strength distribution and any equivalent stress failure criteria. In addition to the reliability model, an energy based failure criterion for elastically anisotropic materials was formulated. The model is different from typical Weibull-based models in that it accounts for strength anisotropy arising from fracture toughness anisotropy and thereby allows for strength and reliability predictions of brittle, anisotropic single crystals subjected to multiaxial stresses. The model is also applicable to elastically isotropic materials exhibiting strength anisotropy due to an anisotropic distribution of flaws. In order to develop and experimentally verify the model, the uniaxial and biaxial strengths of a single crystal nickel aluminide were measured. The uniaxial strengths of the <100> and <110> crystal directions were measured in three and four-point flexure. The biaxial strength was measured by subjecting <100> plates to a uniform pressure in a test apparatus that was developed and experimentally verified. The biaxial strengths of the single crystal plates were estimated by extending and verifying the displacement solution for a circular, anisotropic plate to the case of a variable radius and thickness. The best correlation between the experimental strength data and the model predictions occurred when an anisotropic stress analysis was combined with the normal stress criterion and the strength parameters associated with the <110> crystal direction.

Salem, Jonathan A.↗
DOE OSTI2021

Realization of an Ideal Cairo Tessellation in Nickel Diazenide NiN 2 : High-Pressure Route to Pentagonal 2D Materials

Most of the studied two-dimensional (2D) materials are based on highly symmetric hexagonal structural motifs. In contrast, lower-symmetry structures may have exciting anisotropic properties leading to various applications in nanoelectronics. In this work we report the synthesis of nickel diazenide NiN 2 which possesses atomic-thick layers comprised of Ni 2 N 3 pentagons forming Cairo-type tessellation. The layers of NiN 2 are weakly bonded with the calculated exfoliation energy of 0.72 J/m 2 , which is just slightly larger than that of graphene. The compound crystallizes in the space group of the ideal Cairo tiling (P4/mbm) and possesses significant anisotropy of elastic properties. The single-layer NiN 2 is a direct-band-gap semiconductor, while the bulk material is metallic. Furthermore, this indicates the promise of NiN 2 to be a precursor of a pentagonal 2D material with a tunable direct band gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Two-dimensional burn velocity and analysis of burn products of nitromethane at high pressure

We have studied the laser-initiated deflagration of nitromethane (CH 3 NO 2 ) under high static pressure in the diamond anvil cell. Time-resolved images of the deflagrations have been obtained using intensified CCDs (ICCD). In contrast with previous work, we rely on spontaneous emission from the reaction, rather than changes in the speckle pattern produced by artificial illumination of the sample. Furthermore, as opposed to the 1D records obtained previously with streak cameras, ICCDs permit imaging of the burn in two dimensions providing the ability to directly observe anisotropic deflagration behavior. Here, we report several examples of this behavior and discuss its possible origins. We have also investigated the products of the reaction using Raman spectroscopy. At pressures below 25 GPa, the burn product is observed to be opaque and has a Raman spectrum consistent with a carbonaceous soot. At pressures above 25 GPa, the burn product is observed to be transparent (consistent with earlier reports) and its Raman spectrum reveals the presence of molecular N 2 and a nitrogen-carbonate species. The latter species can be recovered to ambient pressure and has a Raman spectrum consistent with ammonium carbonate or ammonium bicarbonate. To aid interpretation of the Raman spectra, measurements have also been made on the isotopologues CH 3 15 NO 2 and 13 CH 3 NO 2 . This work establishes the importance of advanced time-resolved imaging to reveal the details of deflagration under high pressure and further advances our understanding of the corresponding chemistry.

Deflagration↗