Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Saturation problem”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Controls on Highly Siderophile Element Concentrations in Martian Basalt: Sulfide Saturation and Under-Saturation

Highly siderophile elements (HSE; Re, Au and the platinum group elements) in shergottites exhibit a wide range from very high, similar to the terrestrial mantle, to very low, similar to sulfide saturated mid ocean ridge basalt (e.g., [1]). This large range has been difficult to explain without good constraints on sulfide saturation or under-saturation [2]. A new model for prediction of sulfide saturation places new constraints on this problem [3]. Shergottite data: For primitive shergottites, pressure and temperature estimates are between 1.2-1.5 GPa, and 1350-1470 C [4]. The range of oxygen fugacities is from FMQ-2 to IW, where the amount of Fe2O3 is low and thus does not have a significant effect on the S saturation values. Finally, the bulk compositions of shergottites have been reported in many recent studies (e.g., [5]). All of this information will be used to test whether shergottites are sulfide saturated [3]. Modeling values and results: The database for HSE partition coefficients has been growing with many new data for silicates and oxides [6-8] to complement a large sulfide database [9- 11]. Combining these data with simple batch melting models allows HSE contents of mantle melts to be estimated for sulfide-bearing vs. sulfide-free mantle. Combining such models with fractional crystallization modeling (e.g., [12]) allows HSE contents of more evolved liquids to be modeled. Most primitive shergottites have high HSE contents (and low S contents) that can be explained by sulfide under-saturated melting of the mantle. An exception is Dhofar 019 which has high S contents and very low HSE contents suggesting sulfide saturation. Most evolved basaltic shergottites have lower S contents than saturation, and intermediate HSE contents that can be explained by olivine, pyroxene, and chromite fractionation. An exception is EET A79001 lithology B, which has very low HSE contents and S contents higher than sulfide saturation values . evidence for sulfide saturation during late fractional crystallization. These results show that shergottite HSE contents are controlled by silicates, oxides, and sulfides. In addition, the mantle producing the most primitive shergottites did not contain near chondritic relative ratios of the HSEs like the terrestrial mantle, and did not experience a late chondritic veneer.

Righter, Kevin↗
DOE OSTI2022

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Static actuator-sharing algorithm for concurrent control of multiple plasma properties

Simultaneous regulation of multiple properties in next-generation tokamaks like ITER and fusion pilot plant may require the integration of different plasma control algorithms. Such integration requires the conversion of individual controller commands into physical actuator requests while accounting for the coupling between different plasma properties. This work proposes a tokamak and scenario-agnostic actuator-sharing algorithm (ASA) to perform the above-mentioned command-request conversion and, hence, integrate multiple plasma controllers. The proposed algorithm implicitly solves a quadratic programming (QP) problem formulated to account for the saturation limits and the relation between the controller commands and physical actuator requests. Since the constraints arising in the QP program are linear, the proposed ASA is highly computationally efficient and can be implemented in the tokamak plasma control system in real time. Furthermore, the proposed algorithm is designed to handle real-time changes in the control objectives and actuators’ availability. Nonlinear simulations carried out using the Control Oriented Transport SIMulator illustrate the effectiveness of the proposed algorithm in achieving multiple control objectives simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Two-phase turbine staging

Two-phase turbines would be useful in Rankine cycles using low-quality or saturated-liquid working fluids. The chief problem is low efficiency due to slip losses in nozzles and due to liquid impact and friction losses in rotors. To raise efficiency, a staging method was tested in which the first stage operated at high speed with low relative liquid velocity, and the second stage recovered the first-stage exit energy. The theoretical efficiency for an experimental turbine operating with Refrigerant-22 was raised from 55 percent to 66 percent, and the measured efficiency was raised from 52 percent to 57 percent.

Elliott, D. G.↗
NASA NTRS1991

Adaptive control with an expert system based supervisory level

Adaptive control is presently one of the methods available which may be used to control plants with poorly modelled dynamics or time varying dynamics. Although many variations of adaptive controllers exist, a common characteristic of all adaptive control schemes, is that input/output measurements from the plant are used to adjust a control law in an on-line fashion. Ideally the adjustment mechanism of the adaptive controller is able to learn enough about the dynamics of the plant from input/output measurements to effectively control the plant. In practice, problems such as measurement noise, controller saturation, and incorrect model order, to name a few, may prevent proper adjustment of the controller and poor performance or instability result. In this work we set out to avoid the inadequacies of procedurally implemented safety nets, by introducing a two level control scheme in which an expert system based 'supervisor' at the upper level provides all the safety net functions for an adaptive controller at the lower level. The expert system is based on a shell called IPEX, (Interactive Process EXpert), that we developed specifically for the diagnosis and treatment of dynamic systems. Some of the more important functions that the IPEX system provides are: (1) temporal reasoning; (2) planning of diagnostic activities; and (3) interactive diagnosis. Also, because knowledge and control logic are separate, the incorporation of new diagnostic and treatment knowledge is relatively simple. We note that the flexibility available in the system to express diagnostic and treatment knowledge, allows much greater functionality than could ever be reasonably expected from procedural implementations of safety nets. The remainder of this chapter is divided into three sections. In section 1.1 we give a detailed review of the literature in the area of supervisory systems for adaptive controllers. In particular, we describe the evolution of safety nets from simple ad hoc techniques, up to the use of expert systems for more advanced supervision capabilities.

Sullivan, Gerald A.↗
NASA NTRS1980

Multistage depressed collector with efficiency of 90 to 94 percent for operation of a dual-mode traveling wave tube in the linear region

An axisymmetric, multistage, depressed collector of fixed geometric design was evaluated in conjunction with an octave bandwidth, dual mode traveling wave tube (TWT). The TWT was operated over a wide range of conditions to simulate different applications. The collector performance was optimized (within the constraint of fixed geometric design) over the range of TWT operating conditions covered. For operation of the TWT in the linear, low distortion range, 90 percent and greater collector efficiencies were obtained leading to TWT overall efficiencies of 20 to 35 percent, as compared with 2 to 5 percent with an undepressed collector. With collectors of this efficiency and minimized beam interception losses, it becomes practical to design dual mode TWT's such that the low mode can represent operation well below saturation. Consequently, the required pulse up in beam current can be reduced or eliminated, and this mitigates beam control and dual mode TWT circuit design problems. For operation of the dual mode TWT at saturation, average collector efficiencies in excess of 85 percent were obtained for both the low and high modes across an octave bandwidth, leading to a three to fourfold increase in the TWT overall efficiency.

Ramins, P.↗
NASA NTRS2003

Recent Advances in Video Meteor Photometry

One of the most common (and obvious) problems with video meteor data involves the saturation of the output signal produced by bright meteors, resulting in the elimination of such meteors from photometric determinations. It is important to realize that a "bright" meteor recorded by intensified meteor camera is not what would be considered "bright" by a visual observer - indeed, many Generation II or III camera systems are saturated by meteors with a visual magnitude of 3, barely even noticeable to the untrained eye. As the relatively small fields of view (approx.30 ) of the camera systems captures at best modest numbers of meteors, even during storm peaks, the loss of meteors brighter than +3 renders the determination of shower population indices from video observations even more difficult. Considerable effort has been devoted by the authors to the study of the meteor camera systems employed during the Marshall Space Flight Center s Leonid ground-based campaigns, and a calibration scheme has been devised which can extend the useful dynamic range of such systems by approximately 4 magnitudes. The calibration setup involves only simple equipment, available to amateur and professional, and it is hoped that use of this technique will make for better meteor photometry, and move video meteor analysis beyond the realm of simple counts.

Swift, Wesley R.↗
DOE OSTI2022

Machine learning and atomic layer deposition: Predicting saturation times from reactor growth profiles using artificial neural networks

In this work, we explore the application of deep neural networks to the optimization of atomic layer deposition (ALD) processes. In particular, we focus on a one-shot optimization problem, where we try to predict the optimal dose time that leads to saturation everywhere in the reactor based on thickness values measured at different points of an ALD reactor after a single trial growth. In order to tackle this problem, we introduce a dataset designed to train neural networks to predict saturation times based on these inputs for a cross-flow ALD reactor. Here, we then explore the predictive ability of artificial neural networks of different depths and sizes using a separate testing dataset to evaluate their accuracies. The results obtained show that networks trained using stochastic gradient descent methods can accurately predict saturation times without requiring any additional information on the surface kinetics. This provides a viable approach to minimize the number of experiments required to optimize new ALD processes in a known reactor, and it highlights the way machine learning can be leveraged for thin film growth and manufacturing. While the datasets and training procedure depend on the reactor geometry, the trained neural networks provide a general surrogate model connecting thickness values and trial dose times with optimal saturation times that can be reused for different ALD processes within the same reactor.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Immortal rays: Rethinking random ray neutron transport on GPU architectures

The Random Ray Method (TRRM) is a recently developed adaptation of the Method of Characteristics for neutral particle transport simulations. TRRM has demonstrated excellent performance on 3D nuclear reactor benchmark problems using CPU-based compute systems. When porting to GPU-based systems, however, new performance challenges arise that are unique to processors targeting massive fine-grained parallelism. For smaller problems, or for large problems that are domain decomposed across many computational nodes, the problem size per node has insufficient parallelism to saturate GPU node resources, thus greatly limiting speedup. In this study, we report on a newly developed “immortal ray” variant of TRRM. Here, the immortal ray technique exposes significantly more fine-grained parallelism by fundamentally reformulating the numerical details of ray discretization, resulting in performance tradeoffs with significant overall benefit on GPUs. For very small 2D simulation problems we found the new immortal ray variant allowed for up to a 4.4x speedup when run on a single GPU. For larger 3D simulation problems we found the new variant improved strong scaling by 3x when run on the Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A General Decentralized Dynamic State Estimation With Synchronous Generator Magnetic Saturation

We report the saturation of the nonlinear magnetic circuit of synchronous generators is often neglected when performing Kalman filter-based dynamic state estimation (DSE), yielding significant estimation bias. This letter addresses this problem and proposes a generalized DSE framework to handle magnetic saturation. Moreover, this letter derives a state initialization procedure that improves the Kalman filter tracking speed. The framework is flexible in dealing with different saturation functions and generator models. Numerical results on the Texas 2000-bus system verify the effectiveness of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Safe microburst penetration techniques: A deterministic, nonlinear, optimal control approach

A relatively large amount of computer time was used for the calculation of a optimal trajectory, but it is subject to reduction with moderate effort. The Deterministic, Nonlinear, Optimal Control algorithm yielded excellent aircraft performance in trajectory tracking for the given microburst. It did so by varying the angle of attack to counteract the lift effects of microburst induced airspeed variations. Throttle saturation and aerodynamic stall limits were not a problem for the case considered, proving that the aircraft's performance capabilities were not violated by the given wind field. All closed loop control laws previously considered performed very poorly in comparison, and therefore do not come near to taking full advantage of aircraft performance.

Psiaki, Mark L.↗
NASA NTRS1992

The Zhamanshin impact feature: A new class of complex crater?

The record of 10-km-scale impact events of Quaternary age includes only two 'proven' impact structures: the Zhamanshin Impact Feature (ZIF) and the Bosumtwi Impact Crater (BIC). What makes these impact landforms interesting from the standpoint of recent Earth history is their almost total lack of morphologic similarity, in spite of similar absolute ages and dimensions. The BIC resembles pristine complex craters on the Moon to first order (i.e., 'U'-shaped topographic cross section with preserved rim), while the ZIF displays virtually none of the typical morphologic elements of a 13- to 14-km-diameter complex crater. Indeed, this apparent lack of a craterlike surficial topographic expression initially led Soviet geologists to conclude that the structure was only 5.5 to 6 km in diameter and at least 4.5 Ma in age. However, more recent drilling and geophysical observations at the ZIF have indicated that its pre-erosional diameter is at least 13.5 km, and that its age is most probably 0.87 Ma. Why the present topographic expression of a 13.5-km complex impact crater less than 1 m.y. old most closely resembles heavily degraded Mesozoic shield craters such as Lappajarvi is a question of considerable debate. Hypotheses for the lack of a clearly defined craterlike form at the ZIF include a highly oblique impact, a low-strength 'cometary' projectile, weak or water-saturated target materials, and anomalous erosion patterns. The problem remains unresolved because typical erosion rates within the arid sedimentary platform environment of central Kazakhstan in which the ZIF is located are typically low; it would require at least a factor of 10 greater erosion at the ZIF in order to degrade the near-rim ejecta typical of a 13.5-km complex crater by hundreds of meters in only 0.87 Ma, and to partially infill an inner cavity with 27 cu km (an equivalent uniform thickness of infill of 166 m). Our analysis of the degree of erosion and infill at the ZIF calls for rates in the 0.19 to 0.38 mm/yr range over the lifetime of the landform, which are a factor of 10 to 20 in excess of typical rates for the Kazakhstan semidesert.

Garvin, J. B.↗
NASA NTRS1998

Uncertainty in Calibration, Detection and Estimation of Metal Concentrations in Engine Plumes Using OPAD

Improvements in uncertainties in the values of radiant intensity (I) can be accomplished mainly by improvements in the calibration process and in minimizing the difference between the background and engine plume radiance. For engine tests in which the plume is extremely bright, the difference in luminance between the calibration lamp and the engine plume radiance can be so large as to cause relatively large uncertainties in the values of R. This is due to the small aperture necessary on the receiving optics to avoid saturating the instrument. However, this is not a problem with the SSME engine since the liquid oxygen/hydrogen combustion is not as bright as some other fuels. Applying the instrumentation to other type engine tests may require a much brighter calibration lamp.

Hopkins, Randall C.↗
NASA NTRS2002

An Optical Index of Phytoplankton Photoacclimation and Its Relation to Light-Saturated Photosynthesis in the Sea

In relation to understanding ocean biology at the global scale, one of NASA's primary foci has been measurements of near-surface concentrations of phytoplankton chlorophyll. Chlorophyll is an important light-absorbing pigment in phytoplankton. The absorbed light energy is used to fix carbon in the process of photosynthesis. Photosynthesis, in turn, is critical to the growth of phytoplankton and the function of entire marine ecosystems. Thus, the use of satellite surface chlorophyll data to estimate primary production in the ocean has been a key focus of much biological oceanography research. One of the major challenges in this research is to develop relationships that allow a given chlorophyll concentration (a standing stock) to be interpreted in terms of carbon fixation (a rate). This problem centers on the description of the light-saturated photosynthetic rate, Pbmax. In this paper, we describe how optical measurements of light attenuation provide information on particulate organic carbon (POC) concentrations. We then show how the ratio of POC to chlorophyll (Theta) provides critical information on variability in Pbmax. We then test this relationship between Theta and Pbmax using field data from a variety of open ocean ecosystems.

Behrenfeld, Michael J.↗
NASA NTRS1986

Models of the plasmaspheric thermal plasma distribution

Current understanding of the thermal plasma in the atmosphere and its coupling to the ionosphere is reviewed. Existing models appear adequate to explain the gross behavior of the cold thermal plasma, but there remain some vexing problems. Notably, (1) why does the density in flux tubes appear to saturate at lower values than are predicted theoretically, (2) what causes the sunset peak in measured Te, and (3) why does the equatorial plasmapause signature differ in latitude from the ionosphere signatures. The more difficult problem of what happens during the early stages of refilling after a magnetic storm, when the high altitude plasma is likely to be supersonic and collisionless, has received much attention, but the results are not definite. A number of papers have dealt with the interaction of supersonic counterstreaming fluxes and there are now models that can handle the transition from supersonic to subsonic flows although the transition from a collisionless to a collision-dominated plasma remains difficult to deal with.

Richards, P. G.↗
DOE OSTI2022

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗