Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

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 37 records · Page 2

Time dependent model for non-inductive ECH X-I current ramp-up for SHPD tokamak facility

Non-inductive (NI) plasma current start-up and ramp-up is an important research topic for spherical tokamak (ST) based reactors and fusion pilot plant (FPP). For a compact FPP, the OH flux availability is highly restricted due to its compact geometry. Efficient fundamental extraordinary mode (X-I) electron cyclotron heating (ECH) current start-up and ramp-up regime was identified for a reactor-like high toroidal magnetic field range which has more than a hundred times higher current drive efficiency compared to more conventional ECH methods for the relevant start-up temperature range. High current drive efficiency is possible due to the strong X-I fundamental ECH interaction only with unidirectional passing electrons constrained by the wave accessibility conditions. Here, we extend the X-I electron cyclotron current drive (ECCD) investigation to a time dependent model to simulate the non-inductive current ramp-up to 10 MA for the Sustained High-Power Density (SHPD) facility. As the X-I ECCD driven current I EC rises, due to the back EMF driven negative current, the net plasma current Ip rises more slowly with the current resistive time scale. For tokamak confinement time (both L-mode and H-Mode) which tends to rise with I p , a positive feedback results and even with constant applied ECH power, T e0 , I EC , and I p can continue to rise to very high values. However, in a realistic situation, the T e0 rise should saturate due to a number of factors such as enhanced core radiation and increased power loss at high temperature. To simulate this effect, we adopt a maximum T e0 model which would limit the temperature rise to certain T e0 . With this model, we investigated the current ramp-up for various ECH power levels and the maximum T e0 of 15, 20, and 25 keV. We find that while the power required is reduced with increasing T e0 limit due to increased current drive efficiency, the time to reach 10 MA tends to go up due to the reduced plasma resistivity for the higher T e0 limit. We also find for a given Te0 limit, the time to reach 10 MA tends to be reduced by increasing the applied ECH power by over driving the current ramp-up where I EC is driven at significantly higher level than 10 MA. While the current ramp-up time may not be an issue for the steady-state reactor systems, if it is desirable to minimize the current ramp-up time, it is prudent to have a sufficient ECH power for current over-drive and have some T e0 limiting tools such as impurity seeding for enhanced radiation. In conclusion, a well-controlled NI ECH start-up also has a potential of improving the tokamak start-up reliability and avoid run-away electrons while the NI off-axis current drive could enhance MHD stability and plasma performance improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent models of single- and double-layer mantle convection

Time-dependent numerical models of two-dimensional convection which simulate mantle convection are presented. Starting from a fluid initially at rest with a purely conductive temperature profile, it is found that for a Rayleigh number of ten million, the model experiences a transient period of double-layer convection which lasts on the order of hundreds of millions of years when scaled to the earth's mantle. It is suggested that transient periods of double-layer convection in the earth's mantle may have existed in the past, whether or not such a period exists today.

Boss, A. P.↗

A three-dimensional time-dependent model of the polar wind

A time-dependent three-dimensional multiion model of the polar wind was developed, which covers the altitude range of from 120 to 9000 km and takes into account supersonic ion outflow, shock formation, and ion energization during plasma expansion events. The model was used to study the temporal response of global polar wind to changing magnetospheric conditions, for the winter solstice and for solar-minimum conditions in the northern polar region. Graphs illustrating temporal changes with changes in T(e), T(i), and T(n) along the dawn, the trough, and the dusk convection trajectories and in the O(+), O, and H densities along the same convection trajectories are presented together with conntours of the H(+) and the O(+) densities along the three convection trajectories.

Schunk, R. W.↗

Strain Rate Dependent Modeling of Polymer Matrix Composites

A research program is in progress to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Strain rate dependent inelastic constitutive equations have been developed to model the polymer matrix, and have been incorporated into a micromechanics approach to analyze polymer matrix composites. The Hashin failure criterion has been implemented within the micromechanics results to predict ply failure strengths. The deformation model has been implemented within LS-DYNA, a commercially available transient dynamic finite element code. The deformation response and ply failure stresses for the representative polymer matrix composite AS4/PEEK have been predicted for a variety of fiber orientations and strain rates. The predicted results compare favorably to experimentally obtained values.

Goldberg, Robert K.↗

A comparison of wind observations of the upper thermosphere from the Dynamics Explorer satellite with the predictions of a global time-dependent model

Experimental and theoretical investigations of high-latitude circulation in the upper thermosphere are reported. Vector-wind plots constructed by combining remotely sensed meridional-wind and in situ zonal-wind data obtained during four south-polar and three north-polar passes of the NASA Dynamics Explorer satellite in October and December, 1981, are presented and analyzed, taking the geomagnetic activity level and local solar time into account. The results are then compared with the predictions of 3D time-dependent global models of thermospheric neutral winds (Fuller-Rowell and Rees, 1980, 1981, 1983) adjusted to account for both solar-UV/EUV heating and ionization effects and quiet (Kp from 1 to 2) or moderately disturbed (Kp from 3 to 4) geomagnetic conditions (models Q and MD). An MD model incorporating a self-consistent description of the high-latitude ionosphere and a Q model excluding high-latitude effects are found to give the most accurate predictions for the respective geomagnetic conditions.

Rees, D.↗

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗

Time-dependent modeling of field-aligned current-generated ion transients in the polar wind

The time evolution of field-aligned current-generated transient features in the high-latitude ionosphere is investigated. Ionospheric return currents generate significant downward heavy ion flows in the topside ionosphere with peak values well exceeding 10 to the 8th sq cm/s. When the return current ceases, the polar ionosphere rapidly returns to its previous equilibrium state. During the recovery phase of the return current event, an upward propagating heavy ion transient is formed, which is mainly characterized by a relatively short O(+) upwelling event. On the other hand, the H(+) escape flux remains relatively constant (within 10-20 percent) during field-aligned current events. It is also found that upward currents generate a transient heavy ion outflow, which exceeds the ambient H(+) escape flux by a factor of 3 to 5, depending on the duration and strength of the field-aligned current event.

Gombosi, T. I.↗

Model-dependence of minimal-twist OPEs in d > 2 holographic CFTs

Following recent work on heavy-light correlators in higher-dimensional conformal field theories (CFTs) with a large central charge C T , we clarify the properties of stress tensor composite primary operators of minimal twist, [ T m ], using arguments in both CFT and gravity. We provide an efficient proof that the three-point coupling $\left\langle {\mathcal{O}}_L{\mathcal{O}}_L\left[{T}^m\right]\right\rangle$, where ${\mathcal{O}}_L$ is any light primary operator, is independent of the purely gravitational action. Next, we consider corrections to this coupling due to additional interactions in AdS effective field theory and the corresponding dual CFT. When the CFT contains a non-zero three-point coupling $\left\langle TT{\mathcal{O}}_L\right\rangle$, the three-point coupling $\left\langle {\mathcal{O}}_L{\mathcal{O}}_L\left[{T}^2\right]\right\rangle$ is modified at large C T if $\left\langle TT{\mathcal{O}}_L\right\rangle \sim \sqrt{C_T}$. This scaling is obeyed by the dilaton, by Kaluza-Klein modes of prototypical supergravity compactifications, and by scalars in stress tensor multiplets of supersymmetric CFTs. Quartic derivative interactions involving the graviton and the light probe field dual to ${\mathcal{O}}_L $ can also modify the minimal-twist couplings; these local interactions may be generated by integrating out a spin- ℓ ≥ 2 bulk field at tree level, or any spin ℓ at loop level. These results show how the minimal-twist OPE coefficients can depend on the higher-spin gap scale, even perturbatively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Temperature and Composition Dependence Modeling of Viscosity and Electrical Conductivity of Low-Activity Waste Glass Melts

The development of models that accurately relate the properties of a glass melt to its temperature and composition is important for glass formulation, melter control, and modeling the melt flow, refractory corrosion, and production rate. Using a database consisting of more than 4,000 data points measured between 900 °C and 1250 °C for over 600 unique low-activity waste glass compositions, we developed models for the melt viscosity and electrical conductivity. Models based on the Gaussian process regression approach outperformed models based on the Vogel–Fulcher–Tammann equation according to four standard metrics and yielded reliable prediction intervals. The models found primarily linear effects between properties and individual components, except for the effect of the Na 2 O mass fraction on the electrical conductivity. The effects were found to be consistent with current theories on physical processes involved with those properties.

36 MATERIALS SCIENCE↗

Temperature-dependent model of helium bubble expansion and bursting in tungsten

Continuum-scale models that can reliably predict the behavior of helium in tungsten are of interest to the fusion community due to the projected impact of these materials on fusion reactor operation. In this work, we perform molecular dynamics simulations of spherical helium bubbles of various sizes in tungsten at different temperatures and depths with the goal of determining a mathematical model of the pressure and density at which the bubbles initially expand or burst as a function of depth, size, temperature, and surface orientation. The bubbles are small enough that their loop-punching pressures cannot be accurately predicted with continuum mechanics, and their expansion behavior is important, as it appears to cause many of the features observed on helium-irradiated tungsten surfaces. We vary the temperature, bubble size, bubble depth, and surface orientation in each case, recording the bubble pressure and density that result in bubble expansion. An exponential function with three adjustable parameters is found to fit the results well; the parameters that best fit our results are provided.

36 MATERIALS SCIENCE↗

OC7 phase I: Toward practical sea-state-dependent modeling of hydrodynamic viscous drag and damping

Here, this article presents a collaborative research campaign under the OC7 project on refining the engineering modeling approach for hydrodynamic viscous drag and damping on floating wind platforms, focusing on the adjustment of hydrodynamic drag and damping coefficients for different sea states. The participant simulation results show significant improvements over the previous OC6 project in predicting the low-frequency resonance motion under nonoperational conditions. The improvements are mainly due to enhanced modeling, including the adoption of wave stretching, and directly tuning the coefficients to measured platform motion in waves instead of free decay. For accurate predictions of mean- and slow-drift motion, the better performing models use a decreasing column splash zone drag coefficient and increasing surge damping/drag with increasing wave height. The model tuning for heave and pitch resonance shows less consistency. Generally, both quadratic drag and additional heave or pitch damping are needed for accurate predictions. Alternatively, a quadratic drag formulation with velocity filtering for the rectangular pontoons leads to improved predictions without additional damping. This model is also potentially more predictive, requiring minimal adjustment to its parameters for different conditions.

17 WIND ENERGY↗

Revised Nitrogen Reduction Scaling Relations from Potential-Dependent Modeling of Chemical and Electrochemical Steps

The electrochemical nitrogen reduction reaction (NRR) is a promising route to enable carbon-free ammonia production. However, this reaction is limited by the poor activity and selectivity of current catalysts. The rational design of superior NRR electrocatalysts requires a detailed mechanistic understanding of current material limitations to inform how these might be overcome. The current understanding of how scaling limits NRR on metal catalysts is predicated on a simplified reaction pathway that considers only proton-coupled electron transfer (PCET) steps. Here, we apply grand-canonical density functional theory to investigate a more comprehensive NRR mechanism that includes both electrochemical and chemical steps on 30 metal surfaces in solvent under an applied potential. We applied ..phi..max, a grand-canonical adaptation of the Gmax thermodynamic descriptor, to evaluate trends in catalyst activity. This approach produces a ..phi..max "volcano" diagram for NRR activity scaling on metals that qualitatively differs from the scaling relations identified when only PCET steps are considered. NH3* desorption was found to limit the NRR activity for materials at the top of the volcano and truncate the volcano's peak at increasingly reducing potentials. These revised scaling relations may inform the rational design of superior NRR electrocatalysts. This approach is transferable to study other materials and reaction chemistries where both electrochemical and chemical steps are modeled under an applied potential.

electrocatalysis↗

Spatially dependent modeling and simulation of runaway electron mitigation in DIII-D

New simulations with the Kinetic Orbit Runaway electron (RE) Code (KORC) show RE deconfinement losses to the wall during plasma scrape off are the primary current dissipation mechanism in DIII-D experiments with high-Z impurity injection, and not collisional slowing down. The majority of simulations also exhibit an increase in the RE beam energy due to acceleration by the induced toroidal electric field, even while the RE beam current is decreasing. In this study, KORC integrates RE orbits using the relativistic guiding center equations of motion and incorporates time-sequenced, experimental reconstructions of the magnetic and electric fields and line integrated electron density to construct spatiotemporal models of electron and partially ionized impurity transport in the companion plasma. Comparisons of experimental current evolution and KORC results demonstrate the importance of including Coulomb collisions with partially ionized impurity physics, initial RE energy, pitch angle, and spatial distributions, and spatiotemporal electron and partially ionized impurity transport. The research presented here provides an initial quantification of the efficacy of RE mitigation via injected impurities and identification of the critical role played by loss of confinement due to plasma scrape off on the inner wall as compared to the relatively slow collisional damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A three-dimensional, time-dependent model of Mobile Bay

A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.

Pitts, F. H.↗

A time dependent model of the Jovian current sheet

Analysis of vector helium magnetometer measurements during the Pioneer 10 outbound pass through the Jovian magnetosphere reveals that the average location of the near-equatorial current sheet is a surface whose shape is determined by the velocity at which the magnetic perturbations produced by a rotating tilted dipole propagate to large radial distances. A model is presented which is a surface encircling the planet and rotating rigidly with it. Based on a linear fit to the times of the current sheet crossings, a formula is derived for the surface representing the average position of the near-equatorial current sheet in the Jovian magnetosphere.

Kivelson, M. G.↗

A one-dimensional time-dependent model of the magnetized ionosphere of Venus

The behavior and time evolution of the large-scale magnetic fields and plasma of the dayside Venus ionosphere are studied using a one-dimensional model. The coupled continuity, momentum, and Maxwell's equations are solved simultaneously for O(+), O2(+), and H(+), and the magnetic field. The calculated magnetic field profiles are in good agreement with Pioneer Venus orbiter magnetometer observations. The magnetic field structure is quasi-steady for slow changes of the solar wind dynamic pressure. The peak at 165 km is maintained by downward convection from higher altitudes. The time scale for the decay of the field by the pure one-dimensional vertical diffusion/convection process is several hours unless the flux is resupplied from the top of the ionosphere.

Shinagawa, H.↗

Time-dependent models of radiatively driven stellar winds. I - Nonlinear evolution of instabilities for a pure absorption model

Numerical radiation-hydrodynamics simulations of the nonlinear evolution of instabilities in radiatively driven stellar winds have been performed. The results show a strong tendency for the unstable flow to form rather sharp rarefactions in which the highest speed material has very low density. The qualitative features of the model agree well with the reqirements of displaced narrow absorption components in UV lines.

Owocki, Stanley P.↗