Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state”

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,441 records · Page 80

Design, fabrication, and testing of a water-cooled tungsten-foam high-heat flux plasma-facing module

Here, we present a comprehensive study detailing the design, fabrication and testing a new modular “finger” component, featuring a hexagonal tungsten foam armor tile bonded to the top surface of the water-cooled finger unit. The mechanical, thermal, and particle transport properties of the tungsten foam armor were determined for a wide range of porosity and cell size, and then combined in an equivalent solid material with “effective” properties suitable for large-scale simulations. A series of 3-D multiphysics simulations for coupled thermal and mechanical behavior of the “finger” module are used to optimize the design before fabrication. The simulation results show that, for the expected experimental conditions, a safe range of heat flux and foam porosity is possible for both water-cooling and helium-cooling. Based on these multiphysics simulations, an optimum finger cooling module has been fabricated at Ultramet, Inc. A series of experiments were designed to study the thermomechanical performance of the “finger” unit under steady-state and transient thermal loading at surface heat flux up to 16 MW/m 2 in the High Energy Flux Test facility (HEFTY) at UCLA. Testing conditions included low-power, steady-state plasma exposure at 6 MW/m 2 , followed by severe high heat flux cycles at 16 MW/m 2 . Results of these tests are explained through comparisons with multiphysics simulations. To eliminate the effects of thermomechanical damage carried out on a module on its thermo-mechanical response on a subsequent test, three different yet identically designed and fabricated finger modules were tested. The results show that the current design is capable of cooling the walls of fusion devices at surface heat flux up to 6 MW/m 2 without failure, and that the failure heat flux limit is below 16 MW/m 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Generic Primary-Control Model for Grid-Forming Inverters: Towards Interoperable Operation & Control

This paper outlines an architectural vision centered around the notion of interoperability to integrate grid-forming inverter-based resources in large-scale grids. With the underlying principle of interoperability guiding developments, we focus on modeling the characteristics of droop, virtual synchronous machine, and virtual oscillator controls. Emphasis is placed on these three controllers since they are leading grid-forming control candidates and are likely to be commonplace as primary-control schemes in future systems. We show that these controllers can each be considered as instantiations of a more generic model and that all these controllers exhibit similar droop-like relations between pertinent terminal variables in steady state. This commonality between controllers gives interoperability among them such that automatic synchronization, power sharing, and voltage regulation can be achieved. Simulation results validate the models and demonstrate how the steady-state droop characteristics of these control methods can be aligned with the aid of the developed modeling paradigm.

control↗

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗

Entropy encoding applied to a video DM bit stream

This paper investigates bandwidth compression of delta modulated video signals using a variety of novel techniques. The video signal is encoded at 16 MHz with an adaptive delta modulator (ADM) using the Song video mode algorithm. One method of achieving bandwidth compression is through the detection and reduction of steady state patterns in the DM bit stream. Further reduction is obtained by deleting the horizontal and vertical drive information. Finally, a modified ADM algorithm is introduced to increase the frequency of steady state patterns. When the modified ADM is combined with entropy coding a 30-50% data compression can be realized.

Barba, J.↗

Was Venus wet? Deuterium reconsidered

The ratio of deuterium to hydrogen on Venus has been accepted as proof of a wetter, more earth-like part on that planet. However, the present-day water abundance and the nonthermal hydrogen escape flux on Venus imply that hydrogen is in a steady state and that a hydrogen source, most likely cometary infall, is present. An alternative interpretation of the D/H ratio is offered, in which the measured value is consistent with a steady-state evolution over the age of the solar system. No past water excess is required to explain the isotopic data.

Grinspoon, David H.↗

Transient response of joint dominated space structures - A new linearization technique

The linearization method presented for calculating the transient responses of nonlinear systems due to initial disturbances is an extension of the 'describing function' approach, in which the system's steady-state response is calculated by representing such nonlinear elements as space structure joints with impedances that are functions of response amplitude. It is shown that, for the transient case, the steady-state impedances can be averaged over the range of responses in order to furnish equivalent values of stiffness and damping; these, for a given set of initial displacements, may be treated as constant during calculations of system response.

Foelsche, G. A.↗

Radiation-Driven Warping: Nonisothermal Disks - 2

Recent work by Pringle and by Maloney, Begelman, & Pringle has shown that geometrically thin, optically thick, accretion disks are unstable to warping driven by radiation torque from the central source. This work was confined to isothermal (i.e., surface density Sigma varies as R(sup -3/2) disks. In this paper we generalize the study of radiation-driven warping to include general power-law surface density distributions, Sigma varies as R(sup -delta).We consider the range from Delta = 3/2 (the isothermal case) to Delta = -3/2, which corresponds to a radiation-pressure-supported disk; this spans the range of surface density distributions likely to be found in real astrophysical disks. In all cases there are an infinite number of zero-crossing solutions (i.e., solutions that cross the equator), which are the physically relevant modes if the outer boundary of the disk is required to lie in a specified plane. However, unlike the isothermal disk, which is the degenerate case, the frequency eigenvalues for Delta does not equal 3/2 are all distinct. In all cases the location of the zero moves outward from the steady state (pure precession) value with increasing growth rate; thus, there is a critical minimum size for unstable disks. Modes with zeros at smaller radii are damped. The critical radius and the steady state precession rate depend only weakly on Delta. An additional analytic solution has been found for Delta = 1. The case Delta = 1 divides the solutions into two qualitatively different regimes. For Delta greater than or equal to 1, the fastest growing modes have maximum warp amplitude, close to the disk outer edge, and the ratio of Beta(sub max) to the warp amplitude at the disk inner edge, Beta(sub o), is much greater than 1. For Delta less than 1, Beta(sub max/Beta(sub o) approximately equals 1, and the warp maximum steadily approaches the origin as Delta decreases. This implies that nonlinear effects must be important if the warp extends to the disk inner edge for Delta greater than or equal to 1, but for Delta less than 1 nonlinearity will be important only if the warp amplitude is large at the origin. Because of this qualitative difference in the shapes of the warps, the effects of shadowing of the central source by the warp will also be very different in the two regimes of Delta. This has important implications for radiation-driven warping in X-ray binaries, for which the value of Delta characterizing the disk is likely to be less than unity. In real accretion disks the outer boundary condition is likely to be different from the zero-crossing condition that we have assumed. In accretion disks around massive black holes in active galactic nuclei, the disk will probably become optically thin before the outer disk boundary is reached, whereas in X-ray binaries there will be an outer disk region (outside the circularization radius) in which the inflow velocity is zero but angular momentum is still transported. We show that in both these cases the solutions are similar to the zero-crossing eigenfunctions.

Maloney, Philip R.↗

Thermodynamic Modeling of a Solid Oxide Fuel Cell to Couple with an Existing Gas Turbine Engine Model

The Controls and Dynamics Technology Branch at NASA Glenn Research Center are interested in combining a solid oxide fuel cell (SOFC) to operate in conjunction with a gas turbine engine. A detailed engine model currently exists in the Matlab/Simulink environment. The idea is to incorporate a SOFC model within the turbine engine simulation and observe the hybrid system's performance. The fuel cell will be heated to its appropriate operating condition by the engine s combustor. Once the fuel cell is operating at its steady-state temperature, the gas burner will back down slowly until the engine is fully operating on the hot gases exhausted from the SOFC. The SOFC code is based on a steady-state model developed by The U.S. Department of Energy (DOE). In its current form, the DOE SOFC model exists in Microsoft Excel and uses Visual Basics to create an I-V (current-voltage) profile. For the project's application, the main issue with this model is that the gas path flow and fuel flow temperatures are used as input parameters instead of outputs. The objective is to create a SOFC model based on the DOE model that inputs the fuel cells flow rates and outputs temperature of the flow streams; therefore, creating a temperature profile as a function of fuel flow rate. This will be done by applying the First Law of Thermodynamics for a flow system to the fuel cell. Validation of this model will be done in two procedures. First, for a given flow rate the exit stream temperature will be calculated and compared to DOE SOFC temperature as a point comparison. Next, an I-V curve and temperature curve will be generated where the I-V curve will be compared with the DOE SOFC I-V curve. Matching I-V curves will suggest validation of the temperature curve because voltage is a function of temperature. Once the temperature profile is created and validated, the model will then be placed into the turbine engine simulation for system analysis.

Brinson, Thomas E.↗

Variable gap thermal conductivity apparatus and method

An apparatus and a method for determining the thermal conductivity of a fluid specimen are provided. The apparatus and the method include determining thermal conductivity using a quasi-steady state variable gap axial flow technique. The fluid specimen is heated on one side by a heat source with a known power output and cooled on the other side. After reaching steady state, a resulting temperature drop through the fluid specimen exists. This temperature drop, the known fluid specimen thickness (or gap distance), and the known power output are used to calculate the thermal resistance of the fluid specimen. The thermal conductivity of the fluid specimen is then determined using a curve fit of thermal resistance with respect to gap distance.

Gallagher, Ryan C.↗

Photoinduced phase separation in the lead halides is a polaronic effect

We present a perspective on recent observations of the photoinduced phase separation of halides in multi-component lead-halide perovskites. The spontaneous phase separation of an initial homogeneous solid solution under steady-state illumination conditions is found experimentally to be reversible, stochastic, weakly dependent on morphology, yet strongly dependent on composition and thermodynamic state. Regions enriched in a specific halide species that form upon phase separation are self-limiting in size, pinned to specific compositions, and grow in number in proportion to the steady-state carrier concentration until saturation. These empirical observations of robustness rule out explanations based on specific defect structures and point to the local modulation of an existing miscibility phase transition in the presence of excess charge carriers. A model for rationalizing existing observations based on the coupling between composition, strain, and charge density fluctuations through the formation of polarons is reviewed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relative permeability as a stationary process: Energy fluctuations in immiscible displacement

Relative permeability is commonly used to model immiscible fluid flow through porous materials. In this work, we derive the relative permeability relationship from conservation of energy, assuming that the system to be non-ergodic at large length scales and relying on averaging in both space and time to homogenize the behavior. Explicit criteria are obtained to define stationary conditions: (1) there can be no net change for extensive measures of the system state over the time averaging interval; (2) the net energy inputs into the system are zero, meaning that the net rate of work done on the system must balance with the heat removed; and (3) there is no net work performed due to the contribution of internal energy fluctuations. Results are then evaluated based on direct numerical simulation. Dynamic connectivity is observed during steady-state flow, which is quantitatively assessed based the Euler characteristic. We show that even during steady-state flow at low capillary number ([Formula: see text]), typical flow processes will explore multiple connectivity states. The residence time for each connectivity state is captured based on the time-and-space average. The distribution for energy fluctuations is shown to be multi-modal and non-Gaussian when terms are considered independently. However, we demonstrate that their sum is zero. Given an appropriate choice of the thermodynamic driving force, we show that the conventional relative permeability relationship is sufficient to model the energy dissipation in systems with complex pore-scale dynamics that routinely alter the structure of fluid connected pathways.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A new solution to the accretion problem.

The Bondi-Hoyle theory of line-accretion is briefly reviewed, and it is shown that there are an infinity of steady-state solutions satisfying the requirements originally imposed as boundary conditions. Questions are raised as to the applicability of these solutions to physical reality. The existence of an entirely different type of steady-state solution that implies slow velocity in the accretion stream beyond the neutral point is demonstrated. Accurate numerical values of this solution are obtained for a number of cases. The braking action given by this solution, for any selected cut-off distance, is far stronger than in the Bondi-Hoyle solution. Consideration of the physical state in the accretion stream shows that for this slow solution the appropriate cutoff range would be only of the order of 100 AU rather than of interstellar-distance order.

Lyttleton, R. A.↗

Investigation of space shuttle vehicle 140C configuration orbiter (model 16-0) wheel well pressure loads in the Rockwell International 7.75 x 11 foot wind tunnel (OA143)

Experimental aerodynamic investigations were conducted on a sting mounted .0405-scale representation of the 140C outer mold line space shuttle orbiter configuration in the Rockwell International 7.75 x 11.00 foot low speed wind tunnel. The primary test objectives were to define the orbiter wheel well pressure loading and its effects on landing gear thermal insulation and to investigate the pressure environment experienced by both the horizontal flight nose probe and air vent door probes. Steady state and dynamic pressure values were recorded in the orbiter nose gear well, left main landing gear well, horizontal flight nose probe, and both left and right air vent door probe. All steady state pressure levels were measured by Statham differential pressure transducers while dynamic pressure levels were recorded by Kulite high frequency response pressure sensors.

Mennell, R. C.↗

Characteristics of the NASA Lewis bumpy-torus plasma generated with positive applied potentials

Experimental observations were made during steady-state operation of a bumpy-torus plasma at input powers up to 150 kW in deuterium and helium gas and with positive potentials applied to the midplane electrodes. In this steady-state ion heating method a modified Penning discharge is operated such that the plasma is acted upon by a combination of strong electric and magnetic fields. Experimental investigation of a deuterium plasma revealed electron temperatures from 14 to 140 eV and ion kinetic temperatures from 160 to 1785 eV. At least two distinct modes of operation exist. Experimental data shows that the average ion residence time in the plasma is virtually independent of the magnetic field strength. Data was taken when all 12 anode rings were at high voltage, and in other symmetric configurations in which the toroidal plasma was generated by applying positive potentials to six anode rings, three anode rings, and a single anode ring.

Roth, J. R.↗

Characteristics of the NASA Lewis Bumpy Torus plasma generated with positive applied potentials

Experimental observations have been made during steady-state operation of the NASA Lewis Bumpy Torus experiment at input powers up to 150 kW in deuterium and helium gas, and with positive potentials applied to the midplane electrodes. This steady-state ion heating method utilizes a modified Penning discharge operated in a bumpy torus confinement geometry such that the plasma is acted upon by a combination of strong electric and magnetic fields. Experimental investigation of a deuterium plasma revealed electron temperatures from 14 to 140 eV and ion kinetic temperatures from 160 to 1785 eV. At least two distinct modes of operation exist, each of which is associated with a characteristic range of background pressure and electron temperature. Experimental data show that the average ion residence time in the plasma is virtually independent of the magnetic field strength.

Roth, J. R.↗

Analysis of pressure distortion testing

The development of a distortion methodology, method D, was documented, and its application to steady state and unsteady data was demonstrated. Three methodologies based upon DIDENT, a NASA-LeRC distortion methodology based upon the parallel compressor model, were investigated by applying them to a set of steady state data. The best formulation was then applied to an independent data set. The good correlation achieved with this data set showed that method E, one of the above methodologies, is a viable concept. Unsteady data were analyzed by using the method E methodology. This analysis pointed out that the method E sensitivities are functions of pressure defect level as well as corrected speed and pattern.

Koch, K. E.↗

The low-latitude circulation of Mars

A steady-state two-layer atmospheric model is used in an attempt to simulate the circulation of the equatorial regions of Mars, with consideration of both thermal and mechanical orographic effects. A steady-state solution is sought with the aim of producing mean seasonal circulations upon which tides can be superimposed, and the relative importance of orographic forcing at low latitudes is evaluated. It is assumed that the atmospheres over a topographic feature and over flat terrain possess the same radiative and convective properties, simple dissipation laws are introduced which represent a surface drag effect, and six first-order difference equations are obtained and then reduced to two linear complex second-order coupled difference equations. Analysis of the response of the equinoctial basic wind field to a Gaussian orographic feature situated at the equator and zero deg longitude reveals a number of similarities and differences between the Martian and terrestrial responses. Most of the similarities are explained in terms of the dynamics of a rotating atmosphere, while the major differences are at least partially explained by the different optical properties, lower boundaries, and smaller mass of the Martian atmosphere.

Webster, P. J.↗