Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependent behavior”

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,387 records · Page 77

Procedures for characterizing an alloy and predicting cyclic life with the total strain version of Strainrange Partitioning

Procedures are presented for characterizing an alloy and predicting cyclic life for isothermal and thermomechanical fatigue conditions by using the total strain version of strainrange partitioning (TS-SRP). Numerical examples are given. Two independent alloy characteristics are deemed important: failure behavior, as reflected by the inelastic strainrange versus cyclic life relations; and flow behavior, as indicated by the cyclic stress-strain-time response (i.e., the constitutive behavior). Failure behavior is characterized by conducting creep-fatigue tests in the strain regime, wherein the testing times are reasonably short and the inelastic strains are large enough to be determined accurately. At large strainranges, stress-hold, strain-limited tests are preferred because a high rate of creep damage per cycle is inherent in this type of test. At small strainranges, strain-hold cycles are more appropriate. Flow behavior is characterized by conducting tests wherein the specimen is usually cycled far short of failure and the wave shape is appropriate for the duty cycle of interest. In characterizing an alloy pure fatigue, or PP, failure tests are conducted first. Then depending on the needs of the analyst a series of creep-fatigue tests are conducted. As many of the three generic SRP cycles are featured as are required to characterize the influence of creep on fatigue life (i.e., CP, PC, and CC cycles, respectively, for tensile creep only, compressive creep only, and both tensile and compressive creep). Any mean stress effects on life also must be determined and accounted for when determining the SRP inelastic strainrange versus life relations for cycles featuring creep. This is particularly true for small strainranges. The life relations thus are established for a theoretical zero mean stress condition.

Saltsman, James F.↗

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

Simulation of high-speed (orbital) releases of electron attachment materials in the ionosphere

The dynamics of ionospheric plasma irregularities produced by the release of electron attachment materials at orbital velocities across the geomagnetic field is studied. A two-dimensional electrostatic fluid model which includes electron attachment and mutual neutralization chemistry, self-consistent electric fields, and three-species transport is developed. Numerical simulations are performed to study the behavior at early and at late times after the release. At early times, of the order of or less than the attachment material neutral collision time, the negative ion cloud produced by the release may structure owing to the shear in the E x B velocity within the cloud. At high altitudes the cloud may bifurcate and form vortices on the back. At lower altitudes where ion-neutral collisional effects dominate, this structuring is suppressed. At late times, after a plasma depletion has formed due to neutralization chemistry, the cloud structures by the E x B interchange instability. Depending on the release altitude, the depletion structures by the collisional or inertial limit of this instability.

Scales, W. A.↗

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I - The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

From Well Log to Formation Model: A Novel Laboratory Calibrated Methodology with Demonstration

This work demonstrates how the characterization and modeling of both elastic and creep properties are essential to describe zones in layered rock formations as either low-stress targets for stimulation or high-stress barriers to fracture growth. Prediction of fracture height is critical for designing stimulation operations in oil and gas wells. Ideally, fractures are placed in target zones which will produce hydrocarbons and should not propagate into zones expected to be unproductive or to produce unwanted fluids such as water which in turn must be treated and/or disposed. The essential task in designing stimulation plans is predicting which zones have low horizontal stresses and which will be high-stress barriers to fracture growth. Despite this importance, there are gaps in current knowledge and a complete workflow from laboratory characterization to a finite element model which includes time dependent rock deformation is required. While the research and methodology presented here also have application to CO 2or hydrogen storage, wastewater injection, and geothermal applications, the focus will be on hydrocarbon extraction. This thesis presents the results of a characterization-to-prediction workflow for the Caney shale, which is an emerging hydrocarbon resource in Oklahoma, USA. It begins with an investigation to enable critical evaluation of the Caney zonation into nominally “brittle” and “ductile” zones based on properties observed from well logs. It shows none of the zones are consistently “brittle” or “ductile” mechanical behavior based on the variety of definitions of these terms. However, the nominally ductile zones are weaker and more prone to creep. A laboratory investigation of samples including strength, elastic, and creep properties, is then used in a finite element model of stress evolution. The model includes both elastic deformation and viscoplastic creep. Results predict the least creep-prone layers to have the lowest horizontal stresses, therefore comprising hydraulic fracturing targets. The most creep-prone layers attain a horizontal stress similar to the vertical stress and therefore are predicted to be high stress barriers to hydraulic fracture stimulation. In addition to defining stimulation target intervals, the model shows how as tectonic strain rate increases, there is a transition from creep-dominated stresses to stresses dominated by elasticity.

Benge, Margaret↗

A Synthetic Hydrogel Composite with the Mechanical Behavior and Durability of Cartilage

Here, we report the first hydrogel with the strength and modulus of cartilage in both tension and compression, and the first to exhibit cartilage-equivalent tensile fatigue strength at 100 000 cycles. These properties are achieved by infiltrating a bacterial cellulose (BC) nanofiber network with a poly(vinyl alcohol) (PVA)–poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt) (PAMPS) double network hydrogel. The BC provides tensile strength in a manner analogous to collagen in cartilage, while the PAMPS provides a fixed negative charge and osmotic restoring force similar to the role of aggrecan in cartilage. The hydrogel has the same aggregate modulus and permeability as cartilage, resulting in the same time-dependent deformation under confined compression. The hydrogel is not cytotoxic, has a coefficient of friction 45% lower than cartilage, and is 4.4 times more wear-resistant than a PVA hydrogel. The properties of this hydrogel make it an excellent candidate material for replacement of damaged cartilage.

36 MATERIALS SCIENCE↗

A comparison of the reduced and approximate systems for the time dependent computation of the polar wind and multiconstituent stellar winds

The paper derives the 'reduced' system of equations commonly used to describe the time evolution of the polar wind and multiconstituent stellar winds from the equations for a multispecies plasma with known temperature profiles by assuming that the electron thermal speed approaches infinity. The reduced system is proved to have unbounded growth near the sonic point of the protons for many of the standard parameter cases. For the same parameter cases, the unmodified system exhibits growth in some of the Fourier modes, but this growth is bounded. An alternate system (the 'approximate' system) in which the electron thermal speed is slowed down is introduced. The approximate system retains the mathematical behavior of the unmodified system and can be shown to accurately describe the smooth solutions of the unmodified system. Other advantages of the approximate system over the reduced system are discussed.

Browning, G. L.↗

Bounds on the Rate-Dependent Plastic Flow of Tantalum up to 75 GPa

In this study we report on the continued development of thermodynamics-based analysis of shock waves propagation with the objective of extracting information related to materials strength at high strain rates and pressures. Building on previous results reported for peak stresses of 10 GPa and 25 GPa, we present a series of three-step gas-gun shock experiments designed to explore the pressure and strain rate dependence of plastic flow in polycrystalline tantalum. These experiments at nominal peak stresses of 50 GPa and 75 GPa show the irreversible deformation before pullback to be almost entirely confined to the shock loading, with negligible plastic relaxation on the post-shock plateau. We also add a reverse-ballistics shot at 25 GPa, which was designed to reveal the pullback response with negligible interference from free-surface effects. General thermodynamic considerations allow us to place bounds on the plastic behavior even for parts of the curve that change far too rapidly for the velocimetric time resolution of (conservatively) ~ 5 ns. To analyze the data, we found it necessary to substantially improve the interpolation/extrapolation scheme in order to improve its robustness, flexibility and range of applicability. We describe the new scheme based on splines, as well as an extension of free-surface corrections to the post-shock rarefaction waves. Reanalysis with the new scheme produces results essentially within the error bars previously reported, showing that the known systematic errors associated with free-surface effects are relatively inconsequential for determining thermodynamic paths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the flow-stress model of Preston, Tonks, and Wallace - Observations, parameter constraints, and implementations (Rev. 1)

A detailed discussion about the PTW flow-stress model has been presented, including theoretical considerations, model calibrations to experimental data, and model implementations in LANL ASC hydro-codes. The main findings are: In order to reproduce the intended behavior of the PTW model at high strain rates, the following conditions on the model parameters need to be valid: y 1 ≥ s 0 and y 2 ≥ β ; When the conditions above are not satisfied, PTW and PTW_Mod1 have slightly different behavior in the phonon-drag regime, and should be treated as different models. PTW exhibits a discontinuity above the critical strain rate. PTW_Mod1 does not show such a discontinuity but implies a different strain-rate dependence at high strain rates. Examples of this behavior include older calibrations for stainless steel and depleted uranium; In the thermal-activation regime, some flexibility was introduced in the original PTW model to allow for different rate-dependent behaviors at intermediate strain rates. Example: copper. However, many material calibrations do not require this extra flexibility, and the complexity and ambiguity of the intended behavior can be partially blamed for the inconsistent model implementations across ASC codes; The limit for the material constant p going to zero is well defined and leads to a simplified flow-stress equation, discontinuity free; Modern calibrations of the PTW parameters are done using Impala. Constraints on the PTW parameters are enforced; Implementation of the PTW model in LANL ASC hydro-codes suffered of several shortcom ings. At the time of writing, FLAG and xRAGE contain consistent and correct implemen tations of the PTW and PTW_Mod1 models. Corresponding changes to Pagosa have been suggested.

36 MATERIALS SCIENCE↗

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

59 BASIC BIOLOGICAL SCIENCES↗

Longitudinal Proximity Effects in Superconducting Transition-Edge Sensors

We have found experimentally that the critical current of a square thin-film superconducting transition-edge sensor (TES) depends exponentially upon the side length L and the square root of the temperature T, a behavior that has a natural theoretical explanation in terms of longitudinal proximity effects if the TES is regarded as a weak link between superconducting leads. As a consequence, the effective transition temperature T(sub c) of the TES is current-dependent and at fixed current scales as 1/L(sup 2). We also have found that the critical current can show clear Fraunhofer-like oscillations in an applied magnetic field, similar to those found in Josephson junctions. We have observed the longitudinal proximity effect in these devices over extraordinarily long lengths up to 290 micrometers, 1450 times the mean-free path.

Sadlier, John E.↗

Longitudinal Proximity Effects in Superconducting Transition-Edge Sensors

We have found experimentally that the critical current of a square thin-film superconducting transition-edge sensor (TES) depends exponentially upon the side length L and the square root of the temperature T, a behavior that has a natural theoretical explanation in terms of longitudinal proximity effects if the TES is regarded as a weak link between superconducting leads. As a consequence, the effective transition temperature T(sub c) of the TES is current-dependent and at fixed current scales as 1/L(sup 2). We also have found that the critical current can show clear Fraunhofer-like oscillations in an applied magnetic field, similar to those found in Josephson junctions. We have observed the longitudinal proximity effect in these devices over extraordinarily long lengths up to 290 micrometers, 1450 times the mean-free path.

Sadleir, John E.↗

Adsorption and intracellular uptake of mercuric mercury and methylmercury by methanotrophs and methylating bacteria

The cell surface adsorption and intracellular uptake of mercuric mercury Hg(II) and methylmercury (MeHg) are important in determining the fate and transformation of Hg in the environment. However, current information is limited about their interactions with two important groups of microorganisms, i.e., methanotrophs and Hg(II)-methylating bacteria, in aquatic systems. This study investigated the adsorption and uptake dynamics of Hg(II) and MeHg by three strains of methanotrophs, Methylomonas sp. strain EFPC3, Methylosinus trichosporium OB3b, and Methylococcus capsulatus Bath, and two Hg(II)-methylating bacteria, Pseudodesulfovibrio mercurii ND132 and Geobacter sulfurreducens PCA. Distinctive behaviors of these microorganisms towards Hg(II) and MeHg adsorption and intracellular uptake were observed. The methanotrophs took up 55–80% of inorganic Hg(II) inside cells after 24 h incubation, lower than methylating bacteria (>90%). Approximately 80–95% of MeHg was rapidly taken up by all the tested methanotrophs within 24 h. In contrast, after the same time, G. sulfurreducens PCA adsorbed 70% but took up <20% of MeHg, while P. mercurii ND132 adsorbed <20% but took up negligible amounts of MeHg. These results suggest that microbial surface adsorption and intracellular uptake of Hg(II) and MeHg depend on the specific types of microbes and appear to be related to microbial physiology that requires further detailed investigation. Despite being incapable of methylating Hg(II), methanotrophs play important roles in immobilizing both Hg(II) and MeHg, potentially influencing their bioavailability and trophic transfer. Furthermore, methanotrophs are not only important sinks for methane but also for Hg(II) and MeHg and can influence the global cycling of C and Hg.

59 BASIC BIOLOGICAL SCIENCES↗

Direct quantification of ion composition and mobility in organic mixed ionic-electronic conductors

Ion transport in organic mixed ionic-electronic conductors (OMIECs) is crucial due to its direct impact on device response time and operating mechanisms but is often assessed indirectly or necessitates extra assumptions. Operando x-ray fluorescence (XRF) is a powerful, direct probe for elemental characterization of bulk OMIECs and was used to directly quantify ion composition and mobility in a model OMIEC, poly(3,4-ethylenedioxythiophene)-poly(styrene sulfonate) (PEDOT:PSS), during device operation. The first cycle revealed slow electrowetting and cation-proton exchange. Subsequent cycles showed rapid response with minor cation fluctuation (~5%). Comparison with optical-tracked electrochromic fronts revealed mesoscale structure–dependent proton transport. The calculated effective ion mobility demonstrated thickness-dependent behavior, emphasizing an interfacial ion transport pathway with a higher mobile ion density. The decoupling of interfacial effects on bulk ion mobility and the decoupling of cation and proton migration elucidate ion transport in conventional and emerging OMIEC-based devices and has broader implications for other ionic conductors writ large.

42 ENGINEERING↗

Fermilab Booster loss modelling and rebalancing using Bayesian methods

Fermilab Booster is being upgraded for the PIP-II project to support 20Hz ramp rate at higher intensities. Loss trip limits determine the achievable peak power. To meet PIP-II requirements, losses need to be halved as compared to current levels. Losses primarily occur at injection and transition crossing, with both gradually increasing and threshold-like intensity-dependent behaviors. The existing simulation models are not yet good enough for quantitative loss predictions. In practice, it will be necessary to tune up the Booster using iterative methods and operator intuition. In this paper we present an effort to systematically model Booster losses using active learning (Bayesian exploration) techniques, and subsequently to rebalance them for higher trip limit margins. We first created several sets of spatially and temporally isolated orbit and optics knobs, and trained Gaussian process models for each beam loss monitor as well as beam current. This is a complex task due to safety and timing requirements – we discuss mitigations such as uncertainty constraints and approximate fitting. Once models are stable, we perform large-scale single and multi-objective tuning using scalarized objectives made up of critical beam loss locations. Our results demonstrate significant rebalancing of losses, increasing trip margins, as well as an overall improvement in beam transmission efficiency. We are exploring how to combine existing simulations with experimental data and automate the collection procedure so that more advanced surrogate models can be created over time.

Kuklev, Nikita [Fermilab]↗

Cusp region particle precipitation and ion convection for northward interplanetary magnetic field

Data from Atmosphere Explorer D for periods of strong northward interplanetary magnetic field show the following characteristic behavior in the dayside magnetospheric cusp region: energy-time spectrograms of suprathermal positive ion fluxes exhibit a characteristic 'V' pattern as the spacecraft moves toward higher latitudes; that is, with the peak in the energy spectrum falling in energy and then rising again. Convection velocities follow this pattern closely with strong east-west flows (with antisunward components) occurring in the equatorward half of the 'V' and significant sunward flows occurring in the poleward half of the 'V'. These patterns can be understood qualitatively in terms of a model of ionospheric electric potential produced by the known dependence of Birkeland current densities on magnetic activity.

Burch, J. L.↗

Excitation of slosh waves associated with low frequency impulsive reverse gravity acceleration of geyser initiation

The objective of the cryogenic fluid management of the spacecraft propulsion system is to develop the technology necessary for acquisition or positioning of liquid and vapor within a tank in reduced gravity to enable liquid outflow or vapor venting. The requirement to settle or to position liquid fuel over the outlet end of the spacecraft propellant tank prior to main engine restart poses a microgravity fluid behavior problem. The purpose of the present study is to investigate the stability of the most efficient technique for propellant resettling through the minimization of propellant usage and weight penalties. In this study slosh wave excitation induced by the resettling flow field activated by 0.1 Hz low frequency, impulsive reverse gravity acceleration, during reorientation with the initiation of geyser for liquid fill levels of 30, 50, 65, 70 and 80 percent, have been studied. Characteristics of slosh waves with various frequencies excited are discussed. Slosh wave excitation shift the fluid mass distribution in the container which impose time dependent variations in spacecraft moment of inertia. This information is important for the spacecraft control during the course of liquid reorientation.

Hung, R. J.↗

Thermal expansion and phonon anharmonicity of cuprite studied by inelastic neutron scattering and ab initio calculations

Inelastic neutron scattering measurements were performed with a time-of-flight chopper spectrometer to observe phonons in all parts of the Brillouin zone of a single crystal of cuprite Cu 2 O. We reduced the experimental data to phonon dispersions in the high-symmetry directions, and changes between 10 and 300 K are reported. In this paper, we show ab initio quasiharmonic (QH) and anharmonic (AH) calculations of phonon dispersions. We performed all AH calculations with a temperature-dependent effective potential method. Both QH and AH calculations account for the small negative thermal expansion of cuprite at low temperatures. However, the measured temperature-dependent phonon behavior was predicted more accurately with the AH calculations than the QH ones. Nevertheless, at 300 K, the cubic AH used in this paper did not entirely account for the experimental phonon dispersions in cuprite.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗