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.
PSH Transient Simulation Modeling presentation from the WPTO FY14 - FY16 Peer Review. Transient effects are an important consideration when designing a PSH system, yet numerical techniques for hydraulic transient analysis still need improvements for adjustable-speed (AS) reversible pump-turbine applications.
During loss of coolant accident (LOCA) and reactivity-initiated accident (RIA), nuclear fuel rods experience high heating rates that change the microstructure and properties of zirconium cladding materials. The aim was to determine how different fast heating rates affect the kinetics of transformation and whether its consequences on the material yield stresses in the dual-phase region could be predicted. Phase fraction and texture of cold-rolled Zircaloy-4 were thus measured at fast heating rates similar to those experienced during LOCA and RIA by in-situ high-energy synchrotron x-ray diffraction; the yield stress was also measured at different heating rates to the same dual-phase temperature by electrical resistance method. Crystal plasticity finite element model was used to simulate the flow stresses with the inputs of the measured phase fraction and texture. The model is calibrated using the flow stress at the heating rate of 10 deg. Cs{sup -1} and used to predict the flow stress at the heating rate of 50 deg. Cs{sup -1} at the same temperature. A significant shift of the β transus temperature was identified for fast heating rates in the range of 10-100 deg. Cs{sup -1}. There is great similarity in texture at different heating rates. In the dual-phase temperature region, the a texture is almost identical to the starting texture and the β texture is nearly random. Yield stresses of cold-rolled Zircaloy-4 decrease at faster heating rate at the same temperature of 920 deg. C. A good agreement was found between the experimental and predicted flow stress at different heating rates, indicating this methodology has potential to be used to predict flow behavior in a fast transient regime. (authors)
A method using a Green's function is developed for computing transient temperatures in a semitransparent layer by using the two-flux method coupled with the transient energy equation. Each boundary of the layer is exposed to a hot or cold radiative environment, and is heated or cooled by convection. The layer refractive index is larger than one, and the effect of internal reflections is included with the boundaries assumed diffuse. The analysis accounts for internal emission, absorption, heat conduction, and isotropic scattering. Spectrally dependent radiative properties are included, and transient results are given to illustrate two-band spectral behavior with optically thin and thick bands. Transient results using the present Green's function method are verified for a gray layer by comparison with a finite difference solution of the exact radiative transfer equations; excellent agreement is obtained. The present method requires only moderate computing times and incorporates isotropic scattering without additional complexity. Typical temperature distributions are given to illustrate application of the method by examining the effect of strong radiative heating on one side of a layer with convective cooling on the other side, and the interaction of strong convective heating with radiative cooling from the layer interior.
Trapping effects that include significant delays in permeation saturation, abrupt changes in permeation rate associated with temperature changes, and larger than expected inventories of hydrogen isotopes in the material, were seen in implantation-driven permeation experiments using 25-μm and 50-μm thick tungsten foils at temperatures of 638-825 K. Computer models that simulate permeation transients reproduce the steady-state permeation and reemission behavior of these experiments with expected values of material parameters. However, the transient time characteristics were not successfully simulated without the assumption of traps of substantial trap energy and concentration. An analytical model based on the assumptions of thermodynamic equilibrium between trapped hydrogen atoms and a comparatively low mobile atom concentration successfully accounts for the observed behavior. Using steady-state and transient permeation data from experiments at different temperatures, the effective trap binding energy may be inferred. We analyze a tungsten coated divertor plate design representative of those proposed for ITER and ARIES and consider the implications for tritium permeation and retention if the same trapping we observed was present in that tungsten. Inventory increases of several orders of magnitude may result.
Previous work of Anand and al. (1969) on the effect of aerodynamic torques on the steady-state behavior of a slowly spinning satellite without damping is extended to consideration of the effects of external torques and of a nutational damper on the transient and state motions of a satellite. Both the transient and steady-state behaviors are analyzed, and the attitude errors relative to a fixed star are established. The results indicate that, in the case of a particular satellite to be used for charting X-ray emitting sources in the celestial sphere, the effect of a nutational damper on the accuracy of the charting computations is significant.
The purpose of this experiment is to localize celestial sources of X-rays to a precision of 5 to 20 arcsec, so that they may be identified with optical or radio counterparts. We also study their angular size and structure, to elucidate the X-ray emission mechanism. Four grids, wound with nickel alloy wires, are aligned to produce a modulated X-ray signal when the spacecraft rotation carries the collimator field of view across an X-ray source. An image dissector star camera and rate integrating gyroscopes are used to deduce the celestial position from which a peak X-ray signal arises. On-orbit performance and thermal stability have been very good. Based on the first six weeks of quick look data, we have the following results: the location of an X-ray burst source within a globular cluster; the location of 4U1608-52, a unique object which shows bursting, transient, and steady X-ray behavior and is also an optical source; the location of the X-ray transient H1705-25, which provided Anglo-Australian astronomers the basis for an optical identification with a nova-like object; the location of the X-ray transient H1743-32 in the galactic center; and a precise measurement of the angular extent of the Crab Nebula.
The macrothermomechanical behavior of materials subjected to rapid thermal or mechanical loading such as occurs in most heat treatments is described. The equations are developed for Lagrangian, Eulerian, and intermediary kinematic descriptions and are independent of the constitutive laws and the equation of state; they can be solved numerically for a specified material and boundary conditions. The coupled transport effects between dissipation and energy are included. The conventional linearized stability approach indicates the numerical procedure to be stable, with certain restriction on the time step size.
The transient current systems between tethered plasmas in a large magnetoplasma are investigated experimentally for extrapolation to electrodynamic tethers in space. The studies measure the perturbed magnetic fields and the current density associated with pulsed currents to electrodes in three-dimensional space and time. The electrodes excite electron whistlers because they produce fields that dominantly couple to electrons, allowing pulsed currents to propagate and disperse as whistler wave packets. The wave packets evolve into force-free, flux-ropelike field configurations, and a whistler 'wedge' is formed in the plasma due to 'eddy' currents caused by insulated tethers with dc currents. Substantial radiation into the whistler mode happens with moving VLF antennas as well as tethers, and the wave spread within the ray cone is the most significant characteristic event. The wave spread widens the current channel, incites current closure, and is also associated with a 'phantom loop' phenomenon.
The present computational study of transient heat transfer and fluid flow in a circular pool of n-decane which is undergoing central radiative heating from above gives attention to the volumetric absorption of the radiation incident on the pool surface. The first part of this study notes that buoyancy influences the number and recirculation rates of the subsurface vortices by stabilizing hot subsurface fluid above the colder core fluid; this affects the liquid surface temperature profile and in turn governs the velocity profile that is due to thermocapillarity. In the second part, the effects of gas-liquid phase coupling, variable density and thermophysical properties, and vaporization are considered.
Integration of advanced nuclear reactors to industrial processes in an integrated energy system (IES) has many potential advantages. The coupling of nuclear power to varied industrial processes introduces additional accident scenarios unique to the specific systems included in the IES. Similarly, noise in one system of the IES could potentially create negative effects in the other systems. Because of the limited amount of comparable nuclear IES operating experience, there is added value in scoping the inherent behavior of these systems in response to these transients. Here, this paper investigates the safety-related behavior of an open-source IES dynamic model of a fast spectrum molten salt reactor (MSR) powering a regenerative Rankine cycle for electricity production and a hybrid sulfur (HyS) cycle for hydrogen production. The simulations include a loss of cooling water to the Rankine cycle, a loss of sulfur in the HyS cycle, and the frequency characteristics of the entire IES across six orders of magnitude of frequency. The results show safe behavior in response to the two accident scenarios, wherein a disruption in heat transfer leads to increased salt temperatures and a decrease in reactor power. The frequency analysis shows that at higher frequencies, temperature changes propagating through an IES are increasingly damped to the point of negligibility.
Four different configurations of Multi-Effect-Desalination (MED): backward feed (BF), forward feed (FF), parallel feed (PF) and parallel/cross feed (PCF) are modeled in steady and dynamic operation to examine the behavior of these configurations under transient operation. Adding thermal vapor compressor (TVC) unit with the last effect of parallel/cross configuration is also studied and compared to the other configurations to show the advantage of this type of integration. The transient operation is achieved by applying ramp changes (disturbances) in the main input parameters that include steam flow rate at the heat source and cooling seawater temperature and flow rate at the heat sink. Steady state results indicate that parallel/cross feed has a relatively better performance characteristic than the other feed configurations regarding Gain Output Ratio (GOR) and specific heat consumption (SHC). Adding TVC unit to the parallel/cross feed achieves the lowest water production cost. Dynamic results show that MED-TVC with parallel/cross feed has the fastest response while backward feed and forward feed have the slower response to the applied disturbances. In the case of heat source disturbance, MED plants operating in the backward feed arrangement may be exposed to shut down due to the significant increase in brine level for the first effect. The MED-TVC process is susceptible to the heat sink disturbances compared to the PCF and other configurations. Also, changes in input parameters lead to the highest reduction in GOR for the MED-TVC process compared to the different configurations, especially for cooling seawater temperature. The highest reduction in brine level occurs in MED-TVC for a disturbance in the cooling seawater flow rate and may lead to MED-TVC evaporator drying condition. Hence, MED-TVC will require a reliable control system to avoid this type of plant operational disturbance.
A theoretical model for describing bias-dependent transient and steady-state behaviors of dark current in amorphous selenium (a-Se) avalanche detector structures has been developed. The analytical model considers bulk thermal generation current from mid-gap sates, transient carrier depletion, and carrier injection from the electrodes incorporating avalanche multiplication. The proposed physics-based dark current model is compared with the published experimental results on three potential a-Se avalanche detector structures. The steady-state dark current is the minimum for the structures that have effective blocking layers for both holes and electrons. The transient decay time to reach a plateau decreases considerably with increasing electric field.
The accuracy of current flexible multibody formalisms based on assumed modes is examined in the context of standard spacecraft motions involving structural components undergoing both slow and fast overall translational and rotational motions as well as small deformations. Limitations of current techniques in treating: (1) element-specific coupling behavior of large motion and small deformation, and (2) motion-induced structural stiffness variations, are noted. The roles of nonlinear and linear elastic structural theories in accurately predicting transient large-displacement dynamic behavior of flexible multibody systems are examined in detail. Coupling effects between deformation and overall motion are carefully scrutinized in the context of assumed-mode discretization techniques. Consistently linearized beam, plate, and shell formulations involving in-plane stretch variables are proposed and shown to yield very accurate simulation results and extremely fast modal convergence for most motions involving small strains. In some particular cases, however, in which membrane stiffness dominates bending stiffness, a nonlinear strain formulation is required in order to capture proper coupling between deformation and overall motion.
Dynamic hemiwicking behavior is observable in both nature and a wide range of industrial applications ranging from biomedical devices to thermal management. We present a semi-analytical modeling framework (without empirical fitting coefficients) to predict transient capillary-driven hemiwicking behavior of a liquid through a nano/microstructured surface, specifically a micropillar array. In our model framework, the liquid domain is discretized into micropillar unit cells to enable the time marching of the hemiwicking front. A simplified linear pressure drop is assumed along the hemiwicking length such that the local meniscus curvature, contact angle, and effective liquid height are determined at each time step in our transient model. This semi-analytical model is validated with experimental data from our own experiments and from published literature for different fluids. Our model predicts hemiwicking dynamics with <20% error over a broad range of micropillar geometries with height-to-pitch ratio ranging between ≈0.34 and 6.7 and diameter-to-pitch ratio in the range of ≈0.25–0.7 and without any fitting parameters. For lower diameter-to-pitch ratio data points related to sparse micropillar array arrangements, we suggest modifications to the semi-analytical model. This work sheds light on complex and dynamic solid–liquid–vapor interfacial interactions which could serve as a guide for the design of textured surfaces for wicking enhancement in multi-phase thermal and mass transport technologies and applications.
The purpose of this paper is to study the dynamic behavior of a three-wave journal bearing using a transient approach. The transient analysis permits the determination of the rotor behavior after the fractional frequency whirl appears. The journal trajectory is determined by solving a set of nonlinear equations of motion using the Runge-Katta method. The fluid film forces are computed by integrating the transient Reynolds equation at each time step location of the shaft with respect to the bearing. Because of the large values of the rotational speeds, turbulent effects were included in the computations. The influence of the temperature on the viscosity was also considered. Numerical results were compared to experimenta1 results obtained at the NASA Glenn Research Center. Comparisons of the theoretical results with experimental data were found to be in good agreement. The numerical and experimental results showed that the fluid film of a three-wave journal bearing having a diameter of 30 mm, a length of 27 mm, and a wave amplitude ratio greater than 0.15 is stable even at rotational speeds of 60,000 RPM. For lower wave amplitude ratios, the threshold speed at which the fluid film becomes unstable depends on the wave amplitude and on the supply pocket pressure. Even if the fluid film is unstable, the wave bearing maintains the whirl orbit inside the bearing clearance.
Dynamic simulation of superconducting magnets is critical for the design of quench protection systems to prevent potentially damaging temperatures and high voltage from developing after magnet quench. Modeling these scenarios is challenging due to the many multiscale phenomena which impact magnet behavior. These range from conductor scale effects of quench and interfilament coupling currents up to the behavior of the magnet in its powering and protection circuit. In addition, a strong coupling between electromagnetic and thermal domains is required to capture temperature and field dependent material properties and quench behavior. We present a finite element approach which integrates the various effects into the commercial software ANSYS by means of programming new element types. This is shown capable of simulating the strongly coupled transient electromagnetic, thermal, and circuit behavior of superconducting magnets required for quench protection studies. A benchmarking study is presented which shows close agreement between the new ANSYS elements and a COMSOL Multiphysics implementation developed at CERN for dump resistor and coupling loss induced quench based magnet protection of a Nb 3 Sn block dipole. Following this, the ANSYS implementation is shown reproducing strongly coupled quench back behavior observed during the test of a Nb 3 Sn superconducting undulator prototype at Lawrence Berkeley National Laboratory.
Computational methods are described for simulating the formation and the material behavior of ice in prevailing transient environments. The methodology developed at the NASA Lewis Research Center was adopted. A three dimensional finite-element heat transfer analyzer was used to predict the thickness of ice formed under prevailing environmental conditions. A multi-factor interaction model for simulating the material behavior of time-variant ice layers is presented. The model, used in conjunction with laminated composite mechanics, updates the material properties of an ice block as its thickness increases with time. A sample case of ice formation in a body of water was used to demonstrate the methodology. The results showed that the formation and the material behavior of ice can be computationally simulated using the available composites technology.
This paper describes highlights of an ongoing validation effort conducted to assess the viability of applying a set of analytic point source transient free molecule equations to model behavior ranging from molecular effusion to rocket plumes. The validation effort includes encouraging comparisons to both steady and transient studies involving experimental data and direct simulation Monte Carlo results. Finally, this model is applied to describe features of two exotic transient scenarios involving NASA Goddard Space Flight Center satellite programs.