Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipation rate”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Leakage From Coexisting Geologic Forcing and Injection‐Induced Pressurization: A Semi‐Analytical Solution for Multilayered Aquifers With Multiple Wells

Abstract Abnormal fluid pressures (above or below hydrostatic pressure) can develop and persist in sedimentary basins. The common occurrence of abnormal pressures may cause challenges for project permitting of geological carbon sequestration (GCS), particularly in reservoirs with pre‐injection overpressure. The leaky wells that may exist in some sedimentary basins can provide flow paths between deep brine aquifers and shallower freshwater aquifers. Pre‐injection relative overpressures can cause brine leakage through leaky wells even before any injection occurs. The tendency for flow through leaky wells is coupled with the process of pressure dissipation that occurs through aquitards. Specifically, with non‐zero permeability, aquitards can dissipate pressure over large areal extents and thereby reduce leakage rates through leaky wells. This study presents development of a semi‐analytical solution for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. The geologic forcing that causes abnormal pressures in the multilayered system can coexist with any number of injection, extraction, and leaky wells that also affect fluid pressure. The semi‐analytical model is applied to explore how leakage through leaky wells varies as functions of pressurization rate, along with aquitard and leaky well properties in an overpressured multilayered system. The results show that although injection‐induced pressures can dissipate rapidly through suitably permeable aquitards, coexisting geologic forcing may create sustained rates of brine leakage into freshwater aquifers through leaky wells. In GCS, a very low‐permeability aquitard with high capillary entry pressure to free‐phase CO 2 is desired to serve as the caprock to prevent leakage of CO 2 from the storage reservoir. Nevertheless, the results from this study show that the brine leakage impact to shallow freshwater aquifers through leaky wells might substantially decrease with increasing aquitard permeability values, as long as small aquitard permeability and high capillary entry pressure serve to prevent CO 2 leakage.

Cihan, A.↗
NASA NTRS1971

The fracture energy and some mechanical properties of a polyurethane elastomer.

The energy required to form a unit of new surface in the fracture of a polyurethane elastomer is determined. The rate sensitivity of the material has been reduced by swelling it in toluene. This paper primarily describes the experimental work of measuring the lower limit of the fracture energy. With this value and the creep compliance as a basis, the rate dependence of fracture energy for the unswollen material has been determined. It is thus shown that the dependence of the fracture energy on the rate of crack propagation can be explained by energy dissipation around the tip of the crack. Good agreement between the theoretically and experimentally determined relationships for the rate-sensitive fracture energy is demonstrated.

Mueller, H. K.↗
NASA NTRS1969

Automatic calorimetry system monitors RF power

Calorimetry system monitors the average power dissipated in a high power RF transmitter. Sensors measure the change in temperature and the flow rate of the coolant, while a multiplier computes the power dissipated in the RF load.

Harness, B. W.↗
NASA NTRS2010

Multi-Scale Modeling of Global of Magnetospheric Dynamics

To understand the role of magnetic reconnection in global evolution of magnetosphere and to place spacecraft observations into global context it is essential to perform global simulations with physically motivated model of dissipation that is capable to reproduce reconnection rates predicted by kinetic models. In our efforts to bridge the gap between small scale kinetic modeling and global simulations we introduced an approach that allows to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. We utilized the high resolution global MHD code BATSRUS and incorporate primary mechanism controlling the dissipation in the vicinity of reconnection sites in terms of kinetic corrections to induction and energy equations. One of the key elements of the multiscale modeling of magnetic reconnection is identification of reconnection sites and boundaries of surrounding diffusion regions where non-MHD corrections are required. Reconnection site search in the equatorial plane implemented in our previous studies is extended to cusp and magnetopause reconnection, as well as for magnetotail reconnection in realistic asymmetric configurations. The role of feedback between the non-ideal effects in diffusion regions and global magnetosphere structure and dynamics will be discussed.

Kuznetsova, M. M.↗
NASA NTRS2023

Origin and Evolution of Enceladus’s Tidal Dissipation

Enceladus possesses a subsurface ocean beneath a conductive ice shell. Based on shell thickness models, the estimated total conductive heat loss from Enceladus is 25–40 GW; the measured heat output from the South Polar Terrain (SPT) is 4–19 GW. The present-day SPT heat flux is of order 100 mWm -21 −, comparable to estimated paleo-heat fluxes for other regions of Enceladus. These regions have nominal ages of about 2 Ga, but the estimates are uncertain because the impactor flux in the Saturnian system may not resemble that elsewhere. Enceladus’s measured rate of orbital expansion implies a low dissipation factor Q p for Saturn, with Q p ≈3×10 -3 (neglecting the role of Dione). This value implies that Enceladus’s present-day equilibrium tidal heat production (roughly 50 GW, but with large uncertainties) is in approximate balance with its heat loss. If Q p is constant, Enceladus cannot be older than 1.5 Gyr (because otherwise it would have migrated more than is permissible). However, Saturn’s dissipation may be better described by the “resonance-locking” theory, in which case Enceladus’s orbit may have only evolved outwards by about 35% over the age of the Solar System. In the constant-Q p scenario, any ancient tidal heating events would have been too energetic to be consistent with the observations. Because resonance-locking makes capture into earlier mean-motion orbital resonances less likely, the inferred ancient heating episodes probably took place when the current orbital resonance was already established. In the resonance-locking scenario, tidal heating did not change significantly over time, allowing for a long-lived ocean and a relatively stable ice shell. If so, Enceladus is an attractive target for future exploration from a habitability standpoint.

Satellites of Saturn↗
NASA NTRS1988

Slim accretion disks around black holes

Black hole accretion disk models at moderately high accretion rates have been constructed, taking correctly into account dissipation and heat transport in the transonic part of the flow. The present models differ from those of Shakura and Sunyaev (1973) in many important respects. The surface density at a given radius relates to the accretion rate through an S-shaped function.

Szuszkiewicz, E.↗
DOE OSTI2021

Singularities in nearly uniform one-dimensional condensates due to quantum diffusion

Dissipative systems often exhibit wavelength-dependent loss rates. One prominent example is Rydberg polaritons formed by electromagnetically induced transparency, which have long been a leading candidate for studying the physics of interacting photons and also hold promise as a platform for quantum information. In this system, dissipation is in the form of quantum diffusion, i.e., proportional to k 2 (k being the wavevector) and vanishing at long wavelengths as k→0. Here, we show that one-dimensional condensates subject to this type of loss are unstable to long-wavelength density fluctuations in an unusual manner: after a prolonged period in which the condensate appears to relax to a uniform state, local depleted regions quickly form and spread ballistically throughout the system. We connect this behavior to the leading-order equation for the nearly uniform condensate—a dispersive analog to the Kardar-Parisi-Zhang equation—which develops singularities in finite time. Furthermore, we show that the wavefronts of the depleted regions are described by purely dissipative solitons within a pair of hydrodynamic equations, with no counterpart in lossless condensates. We close by discussing conditions under which such singularities and the resulting solitons can be physically realized.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗
DOE OSTI2024

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗
NASA NTRS1986

Viscous nongeostrophic baroclinic instability

Calculations have been performed of the (linear) stability of a baroclinic flow to three-dimensional perturbations. Both the simple Eady basic state and the rotating Hadley cell of Antar and Fowlis are considered. The independent influences of the Richardson (Ri), thermal Rossby (baroclinicity), Ekman, and Prandtl numbers are examined, as well as the influences of the angle of orientation of the horizontal wave vector and the wavelength. It is shown that if the wavelength is allowed to vary freely, disturbances of the Eady type are preferred (i.e., have greatest growth rate) unless Ri and Ekman numbers are small enough and the thermal Rossby number is large enough. In the latter case, disturbances whose angles of orientation are almost symmetric and whose wavelengths are mesoscale are preferred. If, on the other hand, the wavelength is fixed at a mesoscale size, only the symmetric and almost symmetric modes have growth. By allowing the wave vector orientation to deviate from purely symmetric, it is noted that the region of instability (i.e., critical Ri) is increased, the extent of which is greater for longer wavelength. For Prandtl number = 1, permitting the angle to be nonsymmetric demonstrates the existence of two maxima in growth rate at opposite angles of orientation and with very different energetics. For Prandtl number far enough from one and for large enough dissipation, only one of these two modes has positive growth rates. Growing oscillatory modes were found for some cases.

Miller, T. L.↗
DOE OSTI2025

Suppression of Composition g -modes in Chemically Equilibrating Warm Neutron Stars

We investigate the impact of chemical equilibration and the resulting bulk viscosity on nonradial oscillation modes of warm neutron stars at temperatures up to T ≈ 5 MeV, relevant for protoneutron stars and neutron star postmerger remnants. In this regime, the relaxation rate of weak interactions becomes comparable to the characteristic frequencies of composition g-modes in the core, resulting in resonant damping. To capture this effect, we introduce the dynamical sound speed, a complex, frequency-dependent generalization of the adiabatic sound speed that encodes both the restoring force and the dissipative effects of bulk compression. Using realistic weak reaction rates and three representative equations of state, we compute the complex frequencies of composition g-modes with finite-temperature profiles. We find that bulk viscous damping becomes increasingly significant with temperature and can completely suppress composition g-modes. In contrast, the f-mode remains largely unaffected by bulk viscosity due to its nearly divergence-free character. Our results highlight the sensitivity of g-mode behavior to thermal structure, weak reaction rates, and the equation of state, and establish the dynamical sound speed as a valuable descriptor characterizing oscillation properties in dissipative neutron star matter.

Neutron stars↗
DOE OSTI2025

Electrochemical Random-Access Memory: Progress, Perspectives, and Opportunities

Non-von Neumann computing using neuromorphic systems based on analogue synaptic and neuronal elements has emerged as a potential solution to tackle the growing need for more efficient data processing, but progress toward practical systems has been stymied due to a lack of materials and devices with the appropriate attributes. Recently, solid state electrochemical ion-insertion, also known as electrochemical random access memory (ECRAM) has emerged as a promising approach to realize the needed device characteristics. ECRAM is a three terminal device that operates by tuning electronic conductance in functional materials through solid-state electrochemical redox reactions. This mechanism can be considered as a gate-controlled bulk modulation of dopants and/or phases in the channel. Early work demonstrating that ECRAM can achieve nearly ideal analogue synaptic characteristics has sparked tremendous interest in this approach. More recently, the realization that electrochemical ion insertion can be used to tune the electronic properties of many types of materials including transition metal oxides, layered two-dimensional materials, organic and coordination polymers, and that the changes in conductance can span orders of magnitude has further attracted interest in ECRAM as the basis for analogue synaptic elements for inference accelerators as well as for dynamical devices that can emulate a wide range of neuronal characteristics for implementation in analogue spiking neural networks. At its core, ECRAM shares many fundamental aspects with rechargeable batteries, where ion insertion materials are used extensively for their ability to reversibly store charge and energy. Computing applications, however, present drastically different requirements: systems will require many millions of devices, scaled down to tens of nanometers, all while achieving reliable electronic-state tuning at scaled-up rates and endurances, and with minimal energy dissipation and noise. Further, in this review, we discuss the history, basic concepts, recent progress, as well as the challenges and opportunities for different types of ECRAM, broadly grouped by their primary mobile ionic charge carrier, including Li, protons, and oxygen vacancies.

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

Gravitomagnetic vorticity generation in black hole accretion discs: a potential spatial constraint on plasma flow stability

We calculate the vorticity generation rate in the accretion disc near a slowly rotating black hole in the low-velocity, weak-field limit of general relativity. Specifically, we find that the frame-dragging effect due to the black hole’s rotation – manifested through the gravitomagnetic field – can generate vorticity in a moving plasma in the accretion disc. The mechanism remains operational as long as the accretion disc has non-negligible vertical height and is independent of the exact thermodynamical profile of the disc. The enstrophy density generation rate, as a measure of turbulence and dissipation, is presented, which indicates that the frame-dragging effect can disrupt the stability of the disc away from the z = 0 plane.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Can oscillations grow in a sunspot umbra.

Umbral flashes and running penumbral waves have been attributed by Moore (1972) to overstable oscillations in the umbra. His numerical results were derived by inserting physical conditions at two particular depths beneath the umbral surface. Seven variables must be specified at each point. We have extended Moore's analysis to examine the depth-dependence of overstable oscillations in a recently computed umbral model. Electrical conductivity is evaluated taking full account of partial ionization and magnetic fields. In the surface layers, within 250 km of the top of the umbral convection zone, the conductivity is so low that Joule dissipation is more rapid than the growth rate of oscillations. In these layers, Moore's results are therefore not applicable. At greater depths, oscillations can grow and we agree with Moore that both umbral flashes and penumbral waves may be due to overstable oscillations. However, we suggest that both phenomena can arise at the same depth in the spot, and not in two layers, as Moore suggests.

Mullan, D. J.↗
NASA NTRS1985

Development and life-testing of 10 kW class thermal arcjet engines

A facility to develop and test thermal arcjet engines over extended periods of time has been constructed and is described in this paper. It consists of a large vacuum tank, high capacity vacuum pumps, a 100 kW power supply and a large ammonia propellant storage and delivery system. The facility is instrumented to measure electrical power dissipated in the engine, propellant mass flow rate and developed thrust. Pressures and temperatures up to 2400 K can also be measured. The entire facility is computer-controlled and can be operated unattended for many weeks. Two 30 kW thermal arcjet engines that have been designed, built and are being tested in this facility are also described.

Pivirotto, T. J.↗
NASA NTRS1988

The self-consistent dynamic pole tide in non-global oceans

The dynamic pole tide is determined by solving Laplace tide equations which take into account the presence of continents in oceans, oceanic self-gravitation and loading, and mantle elasticity. Dynamical effects are found to be only mild. It is shown that the dynamical pole tide contributes about one day more to the Chandler period than a static pole tide would, and dissipates wobble energy at a very weak rate. It is noted that, depending on the wobble period predicted for an oceanless elastic earth, mantle anelasticity at low frequencies may nevertheless contribute negligibly to the Chandler period.

Dickman, S. R.↗
NASA NTRS1989

An investigation of cell centered and cell vertex multigrid schemes for the Navier-Stokes equations

Two efficient and robust finite-volume multigrid schemes for solving the Navier-Stokes equations are investigated. These schemes employ either a cell centered or a cell vertex discretization technique. An explicit Runge-Kutta algorithm is used to advance the solution in time. Acceleration techniques are applied to obtain faster steady-state convergence. Accuracy and convergence of the schemes are examined. Computational results for transonic airfoil flows are essentially the same, even for a coarse mesh. Both schemes exhibit good convergence rates for a broad range of artificial dissipation coefficients.

Radespiel, R.↗
NASA NTRS1991

Nitrogen airglow sources - Comparison of Triton, Titan, and earth

The individual contributions of direct solar excitation, photoelectron excitation, and magnetospheric electron excitation of Triton and Titan airglow observed by the Voyager Ultraviolet Spectrometer (UVS) are quantified. The principal spectral features of Triton's airglow are shown to be consistent with precipitation of magnetospheric electrons with power dissipation about 500 million W. Solar excitation rates of the dominant N2 and N(+) emission features are factors of 2-7 weaker than magnetospheric electron excitation. On Titan, the calculated disk center and bright limb N(+) 1085 A intensities due to solar excitation agree with observed values, while the 970 A feature is mostly N21 c5 band emission. The calculated LBH intensity by photoelectrons suggests that magnetospheric electrons play a minor role in Titan's UV airglow. On earth, solar/photoelectron excitation explains the observed N(+) 1085 A and LBH intensites and accounts for only 40 percent of the N(+) 916 A intensity.

Strobel, Darrell F.↗