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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measurements of the $^{5}D^{°}_{4}$–$^{5}P_{3}$ transition of singly ionized atomic iodine using intermodulated laser induced fluorescence

Iodine has been an element of recent interest for commercial use as fuel in electrostatic propulsion systems. A lingering problem when investigating ionized iodine using non-perturbative, laser-based techniques is determining the spectral width, i.e., the species temperature, of iodine. To this end, the hyperfine structure must be well understood to develop a spatially resolved diagnostic technique capable of ion flow and temperature measurements. Previous work investigated the lineshape of the transition between the $^{5}D^{°}_{4}$ and $^{5}P_{3}$ states of singly-ionized atomic iodine (I II) with laser induced fluorescence (LIF), but the hyperfine structure of the transition was unresolved in those measurements. In this work, an intermodulated LIF technique is used to measure an enhanced lineshape of the same I II transition. Here, a linear least squares fitting algorithm is used to fit the transition lineshape, where hyperfine transition locations and theorized relative amplitudes are constrained by theory. A lineshape model that incorporates hyperfine transition amplitude enhancement introduced from an intermodulated laser technique is implemented into the fitting function, as well as a nonlinear laser saturation effect. We report converged hyperfine coupling coefficients for these I II states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

NASA Technology Evaluation for Environmental Risk Mitigation Remediation Technology Collaboration Development

NASA is committed to finding solutions to agency cleanup problems that are better, cheaper, and more effective than the status quo. Unfortunately, some potential solutions involve innovative technologies for which NASA remediation managers may not have a high level of understanding or confidence. Since 2004, NASA's Stennis Space Center (SSC) in Mississippi has been pumping groundwater contaminated with trichloroethylene (TCE) and other halogenated volatile organic compounds (HVOC) from their cleanup location designated "Area G" through extraction wells to an aboveground treatment system. Over time, however, the effectiveness of this treatment strategy has diminished and an alternative approach is needed. In 2012, professionals from NASA's Principal Center for Technology Evaluation for Environmental Risk Mitigation (TEERM) introduced SSC managers to an innovative technology for enhancing the performance of SSC's existing pump and treat system. The technology, generally referred to as in situ chemical oxidation (ISCO), involves slowly and continuously injecting a strong but safe chemical oxidant into the groundwater. Treatment is enhanced by a "surfactant-type effect" which causes residual contamination from saturated soil to be released into the dissolved-phase where it can be readily oxidized. Any dissolved-phase contamination that was not oxidized can be collected by the extraction well network and treated aboveground. SSC was not familiar with the technology so to increase their confidence, TEERM identified a contractor who was willing to demonstrate their product and process at a significantly reduced price. An initial, small-scale demonstration of ISCO began at sse in March 2012 and completed in August 2012. This successful demonstration was followed by three larger-scale ISCO demonstrations between August and December 2012. The contractor's innovative Continuous Injection System (CIS) incorporated "green" and sustainable technologies and practices. A slow injection rate was maintained autonomously by the CIS, eliminating the need for multiple mobilizations of personnel and powered equipment. The CIS was calibrated to deliver only as much reagent as the formation would accept without "short circuiting", minimizing material waste. Public utility water pressure was used to mix, dilute, and inject the reagent. NASA personnel were trained to operate and maintain the system and remote monitoring and injection control capabilities were developed, further reducing the need for contractor mobilizations to the site. An integrated solar photovoltaic panel was used to power the control valves and web monitoring telemetry. SSC provided hydrogen peroxide for the process using their existing supplier for the aboveground pump and treat systems. The only consumable required from the contractor was their proprietary activator to catalyze SSC's peroxide.

Romeo, James↗
NASA NTRS2011

Validating Phasing and Geometry of Large Focal Plane Arrays

The Kepler Mission is designed to survey our region of the Milky Way galaxy to discover hundreds of Earth-sized and smaller planets in or near the habitable zone. The Kepler photometer is an array of 42 CCDs (charge-coupled devices) in the focal plane of a 95-cm Schmidt camera onboard the Kepler spacecraft. Each 50x25-mm CCD has 2,200 x 1,024 pixels. The CCDs accumulate photons and are read out every six seconds to prevent saturation. The data is integrated for 30 minutes, and then the pixel data is transferred to onboard storage. The data is subsequently encoded and transmitted to the ground. During End-to-End Information System (EEIS) testing of the Kepler Mission System (KMS), there was a need to verify that the pixels requested by the science team operationally were correctly collected, encoded, compressed, stored, and transmitted by the FS, and subsequently received, decoded, uncompressed, and displayed by the Ground Segment (GS) without the outputs of any CCD modules being flipped, mirrored, or otherwise corrupted during the extensive FS and GS processing. This would normally be done by projecting an image on the focal plane array (FPA), collecting the data in a flight-like way, and making a comparison between the original data and the data reconstructed by the science data system. Projecting a focused image onto the FPA through the telescope would normally involve using a collimator suspended over the telescope opening. There were several problems with this approach: the collimation equipment is elaborate and expensive; as conceived, it could only illuminate a limited section of the FPA (.25 percent) during a given test; the telescope cover would have to be deployed during testing to allow the image to be projected into the telescope; the equipment was bulky and difficult to situate in temperature-controlled environments; and given all the above, test setup, execution, and repeatability were significant concerns. Instead of using this complicated approach of projecting an optical image on the FPA, the Kepler project developed a method using known defect features in the CCDs to verify proper collection and reassembly of the pixels, thereby avoiding the costs and risks of the optical projection approach. The CCDs composing the Kepler FPA, as all CCDs, had minor defects. At ambient temperature, some pixels look far brighter than they should. These ghot h pixels have a higher rate of charge leakage than the others due to manufacturing variations. They are usually stable over time, and appear at temperatures above 5 oC. The hot pixels on the Kepler FPA were mapped before photometer assembly during module testing. Selected hot pixels were used as target gstars h for the purposes of EEIS testing. gDead h pixels are permanently off, producing a permanently black pixel. These can also be used if there is some illumination of the FPA. During EEIS testing, Dark Current Full Frame Images (FFIs) taken at room temperature were used to create the hot pixel maps for all 84 Kepler photometer CCD channels. Data from two separate nights were used to create two hot pixel maps per channel, which were cross-correlated to remove cosmic ray events which appear to be hot pixels. These hot pixel maps obtained during EEIS testing were compared to the maps made during module testing to verify that the end-to-end data flow was correct.

Standley, Shaun P.↗
DOE OSTI2020

A new soil mixing layer model for simulating conservative solute loss from initially saturated soil to surface runoff

Surface water pollution due to solute loss from soil to surface runoff is a serious threat to the environment, and it is of great significance to predict the solute loss. While soil mixing layer theory has been widely used to predict the solute loss, it has been observed that the conventional soil mixing layer models with constant mixing layer depth underestimate solute concentration in surface runoff, and thus underestimate the solute loss after surface runoff starts. This study introduced a new soil mixing layer model with a new concept that the mixing layer depth increases downward with time after surface runoff occurs. This new concept was implemented by adding a mixing front between the soil mixing layer and its underlying soil layer, which is conceptually similar to the wetting front of soil infiltration given by the Green-Ampt equation. This is equivalent to treating the soil mixing layer as a deforming control volume. A mechanistic explanation for the new model was presented. This study also developed a closed form solution of the new model. The new model was evaluated by simulating five laboratory experiments conducted by two groups of researchers. The first two experiments have ponding water, and consider different drainage and runoff conditions at both early and late runoff times. The other three experiments do not have ponding water, and consider different infiltration conditions at the bottom of soil boxes. Simulation results of the five experiments show that the new model resolves the underestimation problem of conventional soil mixing layer models. The new model is subject to a total of six limits such as that the new model only considers conservative solute. Three limits are related to the assumptions used to derive the closed form solution, and they may be alleviated in a future research as briefly discussed at the end of the paper.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

A Report on Stochastic Fairness Queueing (SFQ) Experiments

SRI International (SRI) has developed an improved queueing algorithm, known as Stochastic Fairness Queueing (SFQ), for best-effort traffic (i.e., traffic that does not require any guaranteed service). SFQ is a probablistic variant of strict fair queueing where instead of a single queue being allocated per flow, a fixed number of queues are used and a hash function maps the IP source and destination to a particular queue. A seed to the hash function is also perturbed occasionally to help distribute the flows amongst different queues when more than one flow maps to the same queue during the lifetime of the flow. SFQ provides 'fair' access by trying to ensure that each flow from source to destination host obtains equal access to the available bandwidth. This report covers a series of experiments performed on DARTnet evaluating the behavior and performance of SFQ against a FIFO queueing discipline. These experiments were designed to show SFQ's advantages and performance, and include tests demonstrating: Fair utilization of available resources; Starvation prevention; Graceful degradation under overload conditions; and Resource usage. In general, the experiments do show that SFQ is better than FIFO queueing at allocating bandwidth equally among a set of flows. SFQ also prevents a stream from dominating the available bandwidth, which seems to be a tendency with FIFO queueing (i.e., if a flow demands more than its share of the available bandwidth, with FIFO queueing that stream receives a disproportionate amount when compared to flows demanding less than their share). Furthermore, SFQ seems to reward 'nice' users of the network by providing a lower variance in delay and more throughput when their resource demand is less than their available share. Both SFQ and FIFO queueing seem to degrade fairly well as the network becomes saturated and to recover well as the network becomes less congested. Not unexpectedly, FIFO queueing is a little more efficient than SFQ-the delays are less and the throughput slightly higher because SFQ requires more processing. However, the performance difference between the two queueing disciplines is relatively small. However, the experiments do point out some interesting behavior. FIFO queueing can behave better than SFQ with seed perturbation. We recommend further evaluation of the hash function and the seed perturbation technique. There are probably weaknesses in their current selection that cause this unexpected behavior. SFQ also seems to possess good scaling properties. To verify this, more experiments with a larger number of streams from more hosts need to be executed and examined, including the staggered introduction of streams. Staggering the streams may prove important, because graphs in the degradation experiment revealed some unexpected increases and decreases in throughput, which should be examined. This may again be due to the interaction of the hash function with the seed perturbation but it may also be related to some other unknown problem.

Denny, Barbara A.↗
NASA NTRS2022

Coronal Polarimetry: Determining the Magnetic Origins of Coronal Mass Ejections

Synopsis: The mechanism for the release of stored magnetic energy in solar eruptions remains a major unsolved problem of Heliophysics. Choosing between triggers requires knowledge of pre-eruptive magnetic fields (B). Although linear polarization in visible/infrared (VIR) coronal emission lines reveals intriguing clues about coronal mass ejection (CME) precursor topology, small telescope apertures limit current capability for measuring vector field. Current coronal observations are insufficient to diagnose 3D coronal magnetic fields in CME precursors. Coronal cavities are the ideal candidates for CME precursor studies. BLOS in cavities is a direct measure of stored magnetic free energy, and the presence and location of topological X-points (reconnection locations) and O-points (circulation of B about axis) distinguish predictions of the flux rope-torus-instability and sheared-arcade-breakout CME models. - Finding: Coronal cavities are ubiquitous throughout the solar cycle, erupt as CMEs, and have coronal polarimetric signatures that distinguish between model predictions. - Recommendation: Coronal cavities should be targeted in a comprehensive multiwavelength study of CME precursors & associated eruption-driving mechanisms. Large ground-based VIR telescopes can measure both coronal magnetic field topology and strength through a combination of the saturated Hanle and BLOS -sensitive Zeeman effects. - Finding: The 4m Daniel K. Inouye Solar Telescope (DKIST) and proposed 1.5m Coronal Solar Magnetism Observatory Large Coronagraph (COSMO-LC) make a major leap forward in VIR coronal sensitivity, enabling measurement of vector B. - Recommendation: A dedicated coronal synoptic telescope with a large field-of-view (i.e., COSMO-LC) is needed to identify the dominant CME driving mechanisms, by measuring coronal magnetic fields from precursor state into eruption. In the UV, a new opportunity has arisen to make use of the strong H I Lyman-α coronal line, which obtains a measurement of B that is independent of the Zeeman effect. - Finding: Small-telescope spectropolarimetric capability in the unsaturated Hanle regime provides an independent and complementary coronal magnetic diagnostic to large ground-based telescope measurements. - Recommendation: The proposed 12-cm Coronal Lyman-α Resonance Observatory (CLARO) spectropolarimetric coronagraph demonstrates a path forward for space-based observations of the coronal magnetic field that should be incorporated into future missions away from the Sun-Earth line (e.g., COMPLETE).

magnetic energy↗
DOE OSTI2025

Deep Learning for Subsurface Flow: A Comparative Study of U‐Net, Fourier Neural Operators, and Transformers in Underground Hydrogen Storage

Subsurface flow research is essential for the sustainable management of natural resources and the environment. Deep learning (DL) has significantly advanced this field by developing efficient and accurate surrogate models to replace computationally expensive physics‐based simulations. These surrogate models are commonly used to predict the spatiotemporal evolution of state variables, such as gas saturation and reservoir pressure, in heterogeneous geological formations. Despite the various DL models applied to this task, there is a lack of studies systematically comparing their performance. This absence of comparative analysis leads to somewhat arbitrary DL model selection in subsurface flow research, resulting in suboptimal performance and potentially inaccurate predictions. To bridge this gap, we conduct a systematic comparison study of three popular DL architectures—U‐Net, Fourier Neural Operators (FNO), and Segmentation Transformer (SETR)—in surrogate modeling of underground hydrogen storage (UHS). We focus on UHS due to its promise of enhancing clean energy resilience and its cyclic operational conditions that represent common scenarios in various subsurface applications. We evaluate the models based on accuracy, training cost, and inference speed. The comparison shows that U‐Net achieves the highest accuracy, followed by SETR and FNO. Despite its lower accuracy, FNO has the highest inference speed. SETR offers competitive accuracy with the least training memory usage, demonstrating the potential of transformers in learning subsurface flow. Our results provide guidance for selecting DL models for surrogate modeling in a wide range of subsurface flow problems.

42 ENGINEERING↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
NASA NTRS1996

The infrared spectrum of the Galactic center and the composition of interstellar dust

We have obtained 5-8 micrometers spectra of the Galactic center from the Kuiper Airborne Observatory at resolving powers of approximately 50, approximately 150, and approximately 300. These spectra show absorption features at 5.5, 5.8, 6.1, and 6.8 micrometers. Together with previously observed features in the 3 micrometers region, these features are compared with laboratory spectra of candidate materials. The 3.0 and 6.1 micrometers features are due to the OH stretching and bending variations of H2O and are well fitted by water of hydration in silicates (e.g., talc). The 3.0 micrometer band is equally well fitted by ice mixtures containing 30% H2O, but such mixtures do not provide a good fit to the observed 6.1 micrometer band. The 3.4 and 6.8 micrometers features are identified with the CH stretching and deformation modes in CH2 and CH3 groups in saturated aliphatic hydrocarbons. The 6.1 micrometer band shows a short wavelength shoulder centered on 5.8 micrometer, attributed to carbonyl (C double bond O) groups in this interstellar hydrocarbon dust component. Finally, the narrow 5.5 micrometer feature is also attributed to carbonyl groups, but in the form of metal carbonyls [e.g., Fe(CO)4]. We have derived column densities and abundances along the line of sight toward the Galactic center for the various identified dust components. This analysis shows that hydrocarbon grains contain only 0.08 of the elemental abundance of C and contribute only a relatively minor fraction (0.1) of the total dust volume. Most of the interstellar dust volume is made up of silicates (approximately 0.6). Small graphite grains, responsible for the 2200 angstroms bump, account for 0.07 of the total dust volume. The remaining one-quarter of the interstellar dust volume consists of a material(s) without strong IR absorption features. Likely candidates include large graphite grains, diamonds, or amorphous carbon grains, which all have weak or no IR active modes. Finally, various models for the origin of the hydrocarbon dust component of the interstellar dust are discussed. All of them face some problems in explaining the observations, in particular, the absence of the spectroscopic signature of hydrocrbon grains in sources associated with molecular clouds.

NASA Discipline Exobiology↗
NASA NTRS1996

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗
NASA NTRS2016

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗
NASA NTRS1995

Nonequilibrium, large-amplitude MHD fluctuations in the solar wind

Compressible MHD simulations in one dimension with three-dimensional vectors are used to investigate a number of processes relevant to problems in interplanetary physics. The simulations indicate that a large-amplitude nonequilibrium (e.g., linearly polarized) Alfvenic wave, which always starts with small relative fluctuations in the magnitude B of the magnetic field, typically evolves to flatten the magnetic profile in most regions. Under a wide variety of conditions B and the density rho become anticorrelated on average. If the mean magnetic field is allowed to decrease in time, the point where the transverse magnetic fluctuation amplitude delta B(sub T) is greater than the mean field B(sub 0) is not special, and large values of delta B(sub T)/B(sub 0) do not cause the compressive thermal energy to increase remarkably or the wave energy to dissipate at an unusually high rate. Nor does the 'backscatter' of the waves that occurs when the sound speed is less than the Alfven speed result, in itself, in substantial energy dissipation, but rather primarily in a phase change between the magnetic and velocity fields. For isolated wave packets the backscatter does not occur for any of the parameters examined; an initial radiation of acoustic waves away from the packet establishes a stable traveling structure. Thus these simulations, although greatly idealized compared to reality, suggest a picture in which the interplanetary fluctuations should have small deltaB and increasingly quasi-pressure balanced compressive fluctuations, as observed, and in which the dissipation and 'saturation' at delta B(sub T)/B(sub 0) approximately = 1 required by some theories of wave acceleration of the solar wind do not occur. The simulations also provide simple ways to understand the processes of nonlinear steepening and backscattering of Alfven waves and demonstrate the existence of previously unreported types of quasi-steady MHD states.

Roberts, D. Aaron↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
NASA NTRS2017

Integrated Demand Management: Minimizing Unanticipated Excessive Departure Delay while Ensuring Fairness from a Traffic Management Initiative

This paper introduces NASA's Integrated Demand Management (IDM) concept and presents the results from an early proof-of-concept evaluation and an exploratory experiment. An initial development of the concept was focused on integrating two systems - i.e. the FAA's newly deployed Traffic Flow Management System (TFMS) tool called the Collaborative Trajectory Options Program (CTOP) and the Time-Based Flow Management (TBFM) system with Extended Metering (XM) capabilities to manage projected heavy traffic demand into a capacity-constrained airport. A human-in-the-loop (HITL) simulation experiment was conducted to demonstrate the feasibility of the initial development of the concept by adapting it to an arrival traffic problem at Newark Liberty International Airport (EWR) during clear weather conditions. In this study, the CTOP was utilized to strategically plan the arrival traffic demand by controlling take-off times of both short- and long-haul flights (long-hauls specify aircraft outside TBFM regions and short-hauls specify aircraft within TBFM regions) in a way that results in equitable delays among the groups. Such strategic planning allows less airborne delay to occur within TBFM by feeding manageable long-haul traffic demand while reserving sufficient slots in the overhead streams for the short-haul departures. The manageable traffic demand indicates the TBFM scheduler assigns no more airborne delay than its assigned airspace is capable of absorbing. TBFM then uses its time-based metering capabilities to deliver the desirable throughput by tactically rescheduling the TBFM entered long-haul flights and short-haul departures. Additional research was also performed to explore use of Required Time of Arrival (RTA) capabilities as a potential control mechanism for the airborne flights to improve arrival traffic delivery accuracy of scheduled long-haul traffic demand. The study results show that both short- and long-haul flights received similar ground delays. In addition, there was a noticeable reduction in the total amount of excessive unanticipated last-minute ground delays, i.e. delays that are frequently imposed on the short-haul flight in current day operations due to saturation in the overhead stream, commonly referred to as 'double penalty'. Furthermore, the concept achieved the target throughput while minimizing the expected cost associated with overall delays in arrival traffic. Assessment of the RTA capabilities showed that there was indeed improvement of the scheduled entry times into TBFM regions by using RTA capabilities. However, with respect to reduction in delays incurred within TBFM, there was no observable benefit of improving the precision of long-haul flights entry times.

Collaborative Trajectory Options Program (CTOP)↗
DOE OSTI2021

Proceedings of RIKEN BNL Research Center Workshop: Small-x Physics in the EIC Era [Slides]

Understanding the high energy limit of hadronic and nuclear collisions is at the forefront of nuclear and particle physics. When boosted to ultrahigh energies, all hadrons and nuclei eventually transform into a universal form of matter called the Color Glass Condensate (CGC). The CGC is characterized by the high density (saturation) of small-x gluons which leads to distinct experimental signatures. Tantalizing hints of the CGC have been observed at HERA, RHIC and the LHC, but the prospects for the discovery of the CGC are more promising at future experiments such as the Electron-Ion Collider (EIC) at BNL Indeed, according to the National Science Academy report published in 2019, one of the three major goals of the EIC is to address the nature of the gluon saturation. (The other two are the mass and spin structure of the nucleons.) With this in mind, the small-x community is gearing up to meet the challenges of the EIC era. Over the past several years, there has been impressive progress in the next-to-leading order (NLO) calculations in the CGC framework of various observables such as single hadron production in proton-nucleus collisions, inclusive and exclusive dijet and trijet production in Deep Inelastic Scattering (DIS), jet-plus-photon production in DIS, etc. We expect that NLO calculations will be the standard tool to confront future experimental data at the EIC. Another emerging trend of the community is the interplay between small-x physics and spin physics. The RHIC result for the gluon helicity ΔG has underscored the necessity to understand the longitudinal spin structure of the proton at small-x. There have been theoretical indications that a significant fraction of spin and orbital angular momentum is stored in the small-x region. As for the transversely polarized proton, a surprising new connection between the gluon Sivers function at small-x and the QCD Odderon has been pointed out and its implications at the EIC has been discussed. In view of these developments, we think it is timely to organize a dedicated workshop on small-x physics to summarize the present status of the field and to discuss future directions. A major focus of this workshop will be to identify outstanding problems that could significantly benefit from collaborative efforts amongst scientists working on formal, phenomenological, and computational aspects of small-x physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Mitigation of Silver Ion Loss from Solution by Polymer Coating of Metal Surfaces, Part IV

Ionic silver (Ag+) biocide is a leading candidate to provide residual microbial control in spacecraft potable water systems, but suffers from rapid concentration loss due to interactions with the metallic storage containers and tubing. One approach to resolve this problem is the coating of metal surfaces with an inert barrier film. In previous reports, we have described our investigations addressing Ag+ loss mitigation and adhesion performance of parylene barrier coatings on coupons of several metal alloys and 316Ltubing under static immersion. In such conditions, parylene-C and -AF4 coatings have shown excellent Ag+loss mitigation and mixed long-term adhesion performance, depending on parylene species and substrate surface chemistry/structure. In Part IV of this series, we report on our work to investigate the performance of parylene-C barrier coatings, under more challenging and realistic conditions, in order to evaluate potential suitability for use. The resilience and associated Ag+ loss with Parylene-C coating on½” (12.7 mm)316L tubing and fittings under medium-term immersion and repeated fitting dis/reassembly were investigated. Potential mechanical challenges to barrier coatings in spacecraft potable water systems were investigated, with two main focuses being identified: liquid flow/pressure cycling in tubing and cyclic operation of bellows in positive-expulsion storage tanks. Two corresponding testbeds for experimental characterization of coating performance were developed. In the Flow/Pressure Testbed System, Ag+ solution will be pumped through internally coated 316L tubing, with optional pressure cycling. In the Bellows Tank Testbed System, an internally coated edge-welded 316L bellows containing Ag+ solution will be cyclically extended and compressed, analogous to the operation of bellows tanks used for potable water storage and delivery on the International Space Station. The design and operation of these testbeds are described. Finally, the Ag+ adsorption and saturation behaviors of Kalrez and Viton, for potential use as seal materials in fittings and other connections, were characterized in limited experiments.

mitigation↗
NASA NTRS2005

Experimental Studies of Phase Equilibria of Meteorites and Planetary Bodies

The primary theme of this project was the application of experimental petrology and geochemistry to a variety of problems in meteoritics and planetary geology. The studies were designed to help develop constraints on the histories of primitive meteorites and their components, the environments in which they formed and evolved, and to understand quantitatively the processes involved in the evolution of igneous rocks on the earth and other planetary bodies. We undertook several projects relating to the origin of CAIs and chondrules. Systematics in the thermodynamic properties of CAI-like liquids were investigated and used to elucidate speciation of multi-valent cations and sulfide capacity of silicate melts and to constrain redox conditions and the vapor pressures of volatile species over molten chondrules. We experimentally determined vanadium speciation in meteoritic pyroxenes and in pyroxenes crystallized from CAI-like melts under very reducing conditions. We also found that bulk oxygen isotope compositions of chondrules in the moderately unequilibrated LL chondrites are related to the relative timing of plagioclase crystallization. We completed an experimental study on the vaporization of beta-SiC and SiO2 (glass or cristobalite) in reducing gases and established the conditions under which these presolar grains could have survived in the solar nebula. We expanded our technique for determining the thermodynamic properties of minerals and liquids to iron-bearing systems. We determined activity-composition relationships in Pt-Fe, Pt-Cr and Pt-Fe-Cr alloys. Results were used to determine the thermodynamic properties of chromite-picrochromite spinels including the free energy of formation of end-member FeCr2O4. We also established a new approach for evaluating Pt-Fe saturation experiments. We calculated the T-fO2 relationships in equilibrated ordinary chondrites and thereby constrained the conditions of metamorphism in their parent bodies.

Stolper, Edward M.↗
DOE OSTI2026

Selective Sorbent Design: CaS Aerogel for Rapid Remediation of Aqueous Pb (II)

Heavy metals are a persistent environmental problem due to their high toxicity, even at very low concentrations (parts per billion, ppb). The removal of such diluted heavy metals is challenging because of the competition the counterions (Ca 2+ , Na + , Mg 2+ , etc.) present in natural water bodies. The design of sorbents capable of removing ions below the action limit (15 ppb for Pb 2+ ) requires a strong driving force for selective uptake and rapid removal. In this work, we report the synthesis of porous CaS aerogels (surface area = 143.6 m 2 /g) by oxidative assembly of CaS nanoparticles and describe their use in selective Pb 2+ ion remediation from water. Despite the presence of amorphous CaCO 3 (up to 50 wt %) in the gel network, the gels demonstrated a capacity of 17.1 mmol Pb/g aerogel (3543 mg/g), and this could be augmented to 22.5 mmol Pb/g aerogel (4593 mg/g) by modifying the synthesis to reduce CaCO 3 content to ca. 15 wt %. Moreover, the selectivity of CaS aerogels toward Pb 2+ ions is high, as evidenced by little-to-no change in the distribution constant (K d ∼ 10 4 ) in the presence of competing ions (1 M) such as Na + , Mg 2+ , and Ca 2+ . During remediation with low concentrations (100 ppb) of Pb 2+ with CaS aerogels, the level of Pb 2+ dropped to 5.4 ppb (below the 15 ppb EPA limit) within 1 h with a 95.4% removal efficiency. In contrast to the CO 2 supercritically dried aerogels, lower surface area ambient dried gels (xerogels) only remove 40% of the lead ions from a 100 ppb solution, saturating within 1 h. The efficiency and rapidity of selective Pb 2+ uptake using CdS aerogels arise from a combination of a strong thermodynamic driving force for cation exchange (K eq = 2.5 × 10 27 ) and chemisorption along with favorable kinetics associated with the high surface area porous architecture. These results show that formation of high surface area metal chalcogenide aerogels by oxidative assembly to form nanocrystalline architectures, as previously demonstrated for II−VI and IV−VI semiconductors, can be extended to the more highly ionic alkaline earth sulfides.

Aerogels↗