Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Simulation”

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,153 records · Page 64

Adaptive relay autotuning under static and non-static disturbances with application to friction stir welding

When a traditional relay feedback test is biased, asymmetry is induced, erroneous system identification occurs, and the test may stop oscillating. Adaptive methods which restore symmetry to the relay test are developed, and additional methods are presented to detect and recover from a non-oscillation condition. This is demonstrated both in simulation and in a process with significant transients and disturbances - friction stir welding (FSW). In testing, it is shown that the adaptive method is able to restore symmetry even in the presence of severe and changing disturbances, including after the weld plunge where temperature gradients change rapidly. Finally, these improvements enable system identification of systems such as FSW in a more robust and accurate manner.

42 ENGINEERING↗

In-situ determination of strain during transient burst testing and the temperature dependence of Zircaloy-4 claddings

Understanding fuel system behavior during postulated loss-of-coolant accidents is pertinent for continued safe and efficient operation of light water reactors, particularly as higher burnups are being pursued and safety margins re-evaluated. Conventional mechanical models for the incumbent Zr alloys typically rely on the assumption that steady-state creep is the dominant fuel cladding response during transient accident conditions. To investigate this assumption, simulated accident burst testing was performed on Zircaloy-4 claddings with balloon behavior measured in-situ. Here, two distinct loading conditions were utilized during burst testing: (1) constant-gas-inventory where pressure was allowed to increase with temperature and (2) constant pressure. In-situ strains and strain rates were measured via 2-dimensional digital image correlation techniques and synchronized with temperature to determine deformation dependencies. The temperature dependence of strain rate was characterized by a two segment Arrhenius relationship, with a distinct transition between the high and low temperature/strain regimes. The average activation energy of the lower temperature/strain regime was 328 ± 25 kJ/mol, in agreement with the ~320 kJ/mol used for conventional LOCA models. However, the higher temperature/strain segment, which encompassed most of ballooning, showed increased activation energies as well as a dependence on whether the burst region was in view. For tests that burst away from the camera view, the average high temperature/strain segment activation energy was 635 ± 150 kJ/mol. For samples where the rupture opening formed in view, the average activation energy was 1015 ± 179 kJ/mol. This observed shift in temperature dependence indicates a transition in deformation mechanism at the end of life, possibly to time independent failure mechanisms, which has not yet been visualized in the literature for Zr alloys. Parameters at the transition points were analyzed to determine thresholds for this change in behavior, which occurred at an average hoop strain of 6.9 ± 2.1 %.

36 MATERIALS SCIENCE↗

Why Does Atmospheric Radiative Heating Weaken Midlatitude Cyclones?

Abstract Recent work has indicated that atmospheric radiative heating reduces the kinetic energy of large‐scale eddies in the midlatitudes. However, a physical mechanism that connects radiation to the midlatitude eddy kinetic energy is still uncertain. Using a high‐resolution general circulation model we perform an experiment in which the radiative cooling profile at each model time step is overwritten with the climatological mean, computed from a control simulation. This approach separates the mean and transient effects of radiative heating on the extratropical circulation. We find that, when radiative heating is fixed, the globally‐averaged eddy kinetic energy is enhanced by ∼6%. We show that thermal radiation dampens temperature anomalies near the surface and tropopause in low‐pressure systems, destroying eddy available potential energy and eddy kinetic energy. We identify this as a possible mechanism by which atmospheric radiative heating weakens midlatitude cyclones.

Mischell, Eric↗

Control over emergence and alignment of transient blisters in thermo-responsive gels using hierarchically patterned substrates

Surfaces with tunable microscale textures are vital in a large variety of technological applications, including heat transfer, antifouling and adhesion. To facilitate such broad-scale use, there is a need to create surfaces that undergo reconfigurable changes in topology and thus, enable switchable functionality. To date, there is a relative dearth of methods for engineering surfaces that can be actuated to change topography over a range of length scales, and hence, form tunable hierarchically structured layers. Combining modeling and experiments, we design a geometrically patterned, thermo-responsive poly (N-isopropylacrylamide) gel film that undergoes controllable hierarchical changes in topology with changes in temperature. At the bottom, the film is covalently bound to a solid, curved substrate; at the top, the film encompasses longitudinal rectangular ridges that are oriented perpendicular to the underlying cylindrical curves. At temperatures below lower critical solution temperature (LCST), the swollen gel exhibits 3D variations in polymer density and thickness defined by the gel's top and bottom topography. As the temperature rises above LCST, the interplay between the upper ridges and lower curves in the gel drives non-uniform, directional solvent transport, the nucleation and propagation of a phase-separated higher-density skin layer, and the resulting pressure buildup within the film. These different, interacting kinetic processes lead to an instability, which produces transient microscopic blisters in the film. Through simulations, we show how tuning the width of the ridges modifies the propagation of a skin layer and creates localized pressure build-up points, which enables control over the emergence, distribution, and alignment of the microscopic blisters. Futhermore, we provide a simple argument to predict the size of such microscopic features. Experiments confirm our predictions and further highlight how our computational model enables the rational design of topographical transitions in these tunable films. The development of actuatable, hierarchically structured films provides new routes for achieving switchable functionality in actuators, drug release systems and adhesives.

42 ENGINEERING↗

Gravitational perturbations from oscillons and transients after inflation

We study the scalar and tensor perturbations generated by the fragmentation of the inflaton condensate into oscillons or transients after inflation, using nonlinear classical lattice simulations. Without including the backreaction of metric perturbations, we find that the magnitude of scalar metric perturbations never exceeds a few ×10 -3 , whereas the maximal strength of the gravitational wave signal today is $O$(10 -9 ) for standard post-inflationary expansion histories. We provide parameter scalings for the α-attractor models of inflation, which can be easily applied to other models. We also discuss the likelihood of primordial black hole formation, as well as conditions under which the gravitational wave signal can be at observationally interesting frequencies and amplitudes. Finally, we provide an upper bound on the frequency of the peak of the gravitational wave signal, which applies to all preheating scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Model Identification of Inverter Nonlinear Control Dynamics

To comprehensively predict in silico potential disruptive dynamic performance of converter-based generation and develop solutions, the application of electromagnetic transient programs is desirable. However, these accurate simulations highly depend on inverter models of actual field implementations, which are often unavailable because they are proprietary. To address this problem, we propose a time-domain model identification method that can be used to determine the control models of inverters in continuous time from multiple experimental studies. The focus is on identifying the dynamic representation of controls because they dominate the dynamic response of renewable assets. To illustrate the applicability of the technique, we consider a subclass of inverter controllers for autonomous grid applications that rely on nonlinear oscillators.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Improving power grid transient stability and transfer capability using HVDC emergency power controls

In the future, with high levels of renewable energy sources penetration, the transmission corridors will have to be flexible to operate under faster transients and variable power flows. Fast power ramping and overloading capabilities of HVDC system can be utilized to ensure transient stability of the interconnected power grid. In this paper, HVDC auxiliary controllers have been studied that can provide the power grid with additional synchronizing torques through fast DC power changes in response to arising system conditions. The effectiveness of the auxiliary controllers in improving transient stability is demonstrated by performing dynamic simulations on the detailed model of the Western Interconnection system (WECC).

Sundaresh, Lakshmi↗

Atomic-scale polarization switching in wurtzite ferroelectrics

Ferroelectric wurtzites have the potential to revolutionize modern microelectronics because they are easily integrated with multiple mainstream semiconductor platforms. However, the electric fields required to reverse their polarization direction and unlock electronic and optical functions need substantial reduction for operational compatibility with complementary metal-oxide semiconductor (CMOS) electronics. To understand this process, we observed and quantified real-time polarization switching of a representative ferroelectric wurtzite (Al 0.94 B 0.06 N) at the atomic scale with scanning transmission electron microscopy. The analysis revealed a polarization reversal model in which puckered aluminum/boron nitride rings in the wurtzite basal planes gradually flatten and adopt a transient nonpolar geometry. Independent first-principles simulations reveal the details and energetics of the reversal process through an antipolar phase. Furthermore, this model and local mechanistic understanding are a critical initial step for property engineering efforts in this emerging material class.

36 MATERIALS SCIENCE↗

Inferring small-sample-reactivity worth from oscillation experiment with auto-corrected local flux

To estimate reactivity changes as a function of time from the signal of detectors count rates, one usually uses point kinetic equations. However, when detectors are located inside the core, it is often necessary to compute a correction factor to account for the flux-shape deformation (i.e. the change of detector efficiency that may have been induced by neutron absorbers movements). Various simulation methods have been developed depending on transients observed. They range from the resolution of static equations up to full-time-dependent neutron transport in 3D. However, those methods require high computational-cost calculations of local reaction rates that rely on nuclear data evaluations. Therefore, a 'measured' reactivity is never determined through pure measurements. Taking advantage of the known oscillation movement of a small sample at the center of the MINERVE reactor, the innovative algorithm presented in this paper derives the correction factor needed to get no time drift on reactivity steps. Derived correction factors are in good agreement with those obtained independently with TRIPOLI-4 through sensitivity analysis. Low sensitivity to the sample composition and moderate sensitivity to the detectors position reinforce the reliability of those results. This application paves the way to almost pure nuclear data reactivity measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Impact of Integration Scheme on Performance of Anisotropic Plasticity Models

Given the prevalent role of metals in a variety of industries, schemes to integrate corresponding constitutive models in finite element applications have long been studied. A number of formulations have been developed to accomplish this task; each with their own advantages and costs. Often the focus has been on ensuring the accuracy and numerical stability of these algorithms to enable robust integration. While important, emphasis on these performance metrics may often come at the cost of computational expense potentially neglecting the needs of individual problems. In the current work, the performance of two of the most common integration methods for anisotropic plasticity -- the convex cutting plane (CCP) and closest point projection (CPP) -- across a variety of metrics is assessed; including accuracy and cost. A variety of problems are considered ranging from single elements to large representative simulations including both implicit quasistatic and explicit transient dynamic type responses. The relative performance of each scheme in the different instances is presented with an eye towards guidance on when the different algorithms may be beneficial.

36 MATERIALS SCIENCE↗

Development of SAM Code Capabilities for Safety Analysis of GCR Air-ingress Events

Air-ingress following a depressurized loss-of-forced-cooling (DLOFC) event is a challenging, multiphysics safety scenario for High-Temperature Gas-Cooled Reactors (HTGRs), involving coupled gas composition transport, buoyancy-driven flow redistribution, graphite oxidation, and structural heat-up. Despite its importance — air ingress is a key scenario identified in the PIRT process for the HTGRs — existing system-level safety codes have lacked the integrated capability to simulate the complete event sequence with high confidence. This report documents the development, validation, and demonstration of three new capabilities in the SAM code to address this gap: (1) a multi-component gas mixture flow model with binary diffusion to track the helium-air composition and its effect on system density and flow; (2) a 0-D graphite oxidation model based on the Roes correlation, including oxygen consumption and exothermic heat release; and (3) an isentropic critical flow model for accurate representation of primary system depressurization through a break. These capabilities are validated against two benchmark experiments. The NSTF heavy-gas ingress experiment validates the multi-component flow model: SAM correctly reproduces the rapid buoyancydriven flow stagnation and subsequent natural circulation recovery driven by composition-dependent density changes. The NACOK graphite oxidation experiment validates the oxidation model: SAM predicts a bottom-level graphite weight loss of 25%, in close agreement with the measured 24%, and reproduces the strong axial nonuniformity and block-geometry dependence of oxidation, at a level comparable to the SPECTRA and TINTE codes. The validated capabilities are then exercised together in an integrated, reactor-scale simulation of a DLOFC air-ingress transient in a simplified HTR-PM pebble-bed reactor. In a single calculation spanning approximately 8 days, SAM reproduces the complete accident sequence: rapid depressurization, densityand diffusion-driven air ingress over ˜15 hours, onset of buoyancy-driven natural circulation, exothermic graphite oxidation with a peak fuel temperature at ˜62 hours, and eventual passive cooldown. These results demonstrate that SAM now provides the nuclear community with a preliminarily validated, modern systemlevel tool for HTGR air-ingress safety analysis, filling a recognized capability gap. Future extensions to broaden species tracking, improve oxidation chemistry, and refine the reactor model are discussed.

Yang, Gang↗

Fast differential thermal analysis.

A study is conducted of the limits of time-dependent differential thermal analysis (DTA). Applications of DTA to the hysteresis phenomena of glasses are also reported. The computation of DTA sample temperatures is discussed, giving attention to feasible heating rates, transient effects to be measured, and the simulation of DTA responses to changes in thermal diffusivity similar to hysteresis phenomena in glasses. The reproducibility of temperature recording is checked with zone-refined organic chemicals.

Wunderlich, B.↗

Kinetic approach to the evaporation and condensation problem

In the paper, the Boltzmann equation governing the evaporation and condensation phenomena is solved by the Monte Carlo method. Based on the kinetic theory of gas the role of the non-equilibrium Knudsen layer and the growth of the hydrodynamic region outside the layer as time proceeds are simulated. Results show two possible types of transient developments in the vapor phase. The effects of the molecular absorption coefficient of the phase surface are examined. Except in the case of very strong evaporation the kinematic effects of binary collisions among vapor molecules on the mass flux rate are not serious. The limiting case of the quasi-steady evaporation and the maximal value of the evaporation rate are obtained.

Murakami, M.↗

Engine dynamic analysis with general nonlinear finite-element codes. I Overall approach and development of bearing damper element

NASA-sponsored research on engine dynamic simulation using general finite element nonlinear time transient computer codes available on the open market is reviewed. The approach taken was to develop software packages to model engine components which are not typically found on dynamical structures and are therefore not already computer codes. The software package developed for squeeze-film bearing dampers is outlined, and the results of a parametric study of damper pressure for a variety of specified circular orbits are presented for both long-bearing and short-bearing solutions. The data from a four-degree-of-freedom rotor-damper-stator model under conditions of small rotor unbalance through large rotor unbalance are also given.

Adams, M. L.↗

Nonlinear structural analysis of a turbine airfoil using the Walker viscoplastic material model for B1900 + Hf

A viscoplastic material model for the high temperature turbine airfoil material B1900 + Hf was developed and was demonstrated in a three dimensional finite element analysis of a typical turbine airfoil. The demonstration problem is a simulated flight cycle and includes the appropriate transient thermal and mechanical loads typically experienced by these components. The Walker viscoplastic material model was shown to be efficient, stable and easily used. The demonstration is summarized and the performance of the material model is evaluated.

Meyer, T. G.↗

Systems Improved Numerical Fluids Analysis Code

Systems Improved Numerical Fluids Analysis Code, SINFAC, consists of additional routines added to April, 1983, version of SINDA. Additional routines provide for mathematical modeling of active heat-transfer loops. Simulates steady-state and pseudo-transient operations of 16 different components of heat-transfer loops, including radiators, evaporators, condensers, mechanical pumps, reservoirs, and many types of valves and fittings. Program contains property-analysis routine used to compute thermodynamic properties of 20 different refrigerants. Source code written in FORTRAN 77.

Costello, F. A.↗