Smooth Mathematical Representation of the DER_A Aggregated Model
Not Available
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “smooth saturation function”
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.
Quartz dissolution. I - Negative crystal experiments and a rate law. II - Theory of rough and smooth surfaces
The range of the measured quartz dissolution rates, as a function of temperature and pOH, extent of saturation, and ionic strength, is extended to cover a wider range of solution chemistries, using the negative crystal methodology of Gratz et al. (1990) to measure the dissolution rate. A simple rate law describing the quartz dissolution kinetics above the point of zero charge of quartz is derived for ionic strengths above 0.003 m. Measurements were performed on some defective crystals, and the mathematics of step motion was developed for quartz dissolution and was compared with rough-face behavior using two different models.
Gratz, Andrew J.↗
A Preprocessing Tool for Enhanced Ion Mobility–Mass Spectrometry-Based Omics Workflows
The ability to improve the data quality of ion mobility–mass spectrometry (IM-MS) measurements is of great importance for enabling modular and efficient computational workflows and gaining better qualitative and quantitative insights from complex biological and environmental samples. We developed the PNNL PreProcessor, a standalone and user-friendly software housing various algorithmic implementations to generate new MS-files with enhanced signal quality and in the same instrument format. Different experimental approaches are supported for IM-MS based on Drift-Tube (DT) and Structures for Lossless Ion Manipulations (SLIM), including liquid chromatography (LC) and infusion analyses. The algorithms extend the dynamic range of the detection system, while reducing file sizes for faster and memory-efficient downstream processing. Specifically, multidimensional smoothing improves peak shapes of poorly defined low-abundance signals, and saturation repair reconstructs the intensity profile of high-abundance peaks from various analyte types. Further, other functionalities are data compression and interpolation, IM demultiplexing, noise filtering by low intensity threshold and spike removal, and exporting of acquisition metadata. Several advantages of the tool are illustrated, including an increase of 19.4% in lipid annotations and a two-times faster processing of LC-DT IM-MS data-independent acquisition spectra from a complex lipid extract of a standard human plasma sample. The software is freely available at https://omics.pnl.gov/software/pnnl-preprocessor.
59 BASIC BIOLOGICAL SCIENCES↗
Magnetic field generated resistivity maximum in graphite
In zero magnetic field, B, the electrical resistivity, rho(O,T) of highly oriented pyrolytic (polycrystalline) graphite drops smoothly with decreasing T, becoming constant below 4 K. However, in a fixed applied magnetic field B, the resistivity rho(B,T) goes through a maximum as a function of T, with larger maximum for larger B. The temperature of the maximum increases with B, but saturates to a constant value near 25 K (exact T depends on sample) at high B. In single crystal graphite a maximum in rho(B,T) as a function of T is also present, but has the effects of Landau level quantization superimposed. Several possible explanations for the rho(B,T) maximum are proposed, but a complete explanation awaits detailed calculations involving the energy band structure of graphite, and the particular scattering mechanisms involved.
Wollam, J. A.↗
Magnetic field generated resistivity maximum in graphite
In a zero magnetic field, the electrical resistivity of highly oriented pyrolytic (polycrystalline) graphite drops smoothly with decreasing temperature, becoming constant below 4 K. However, in a fixed applied magnetic field, the resistivity goes through a maximum as a function of temperature, with a larger maximum for a larger field intensity. The temperature of the maximum increases with field intensity, but saturates to a constant value near 25 K (the exact temperature depends on the sample) at high intensity. In single-crystal graphite, a maximum in resistivity as a function of temperature is also present, but has the effects of Landau level quantization superimposed. Several possible explanations for the resistivity maximum are proposed, but a complete explanation awaits detailed calculations involving the energy band structure of graphite, and the particular scattering mechanisms involved.
Woollam, J. A.↗
Geometry of conformal manifolds and the inversion formula
Families of conformal field theories are naturally endowed with a Riemannian geometry which is locally encoded by correlation functions of exactly marginal operators. We show that the curvature of such conformal manifolds can be computed using Euclidean and Lorentzian inversion formulae, which combine the operator content of the conformal field theory into an analytic function. Analogously, operators of fixed dimension define bundles over the conformal manifold whose curvatures can also be computed using inversion formulae. These results relate curvatures to integrated four-point correlation functions which are sensitive only to the behavior of the theory at separated points. We apply these inversion formulae to derive convergent sum rules expressing the curvature in terms of the spectrum of local operators and their three-point function coefficients. We further show that the curvature can smoothly diverge only if a conserved current appears in the spectrum, or if the theory develops a continuum. We verify our results explicitly in 2d examples. In particular, for 2d (2,2) superconformal field theories we derive a lower bound on the scalar curvature, which is saturated by free theories when the central charge is a multiple of three.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling
Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.
58 GEOSCIENCES↗
Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks
Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).
42 ENGINEERING↗
A reference open-source controller for fixed and floating offshore wind turbines
Abstract. This paper describes the development of a new reference controller framework for fixed and floating offshore wind turbines that greatly facilitates controller tuning and represents standard industry practices. The reference wind turbine controllers that are most commonly cited in the literature have been developed to work with specific reference wind turbines. Although these controllers have provided standard control functionalities, they are often not easy to modify for use on other turbines, so it has been challenging for researchers to run representative, fully dynamic simulations of other wind turbine designs. The Reference Open-Source Controller (ROSCO) has been developed to provide a modular reference wind turbine controller that represents industry standards and performs comparably to or better than existing reference controllers. The formulation of the ROSCO controller logic and tuning processes is presented in this paper. Control capabilities such as tip speed ratio tracking generator torque control, minimum pitch saturation, wind speed estimation, and a smoothing algorithm at near-rated operation are included to provide modern controller features. A floating offshore wind turbine feedback module is also included to facilitate growing research in the floating offshore arena. All of the standard controller implementations and control modules are automatically tuned such that a non-controls engineer or automated optimization routine can easily improve the controller performance. This article provides the framework and theoretical basis for the ROSCO controller modules and generic tuning processes. Simulations of the National Renewable Energy Laboratory (NREL) 5 MW reference wind turbine and International Energy Agency 15 MW reference turbine on the University of Maine semisubmersible platform are analyzed to demonstrate the controller's performance in both fixed and floating configurations, respectively. The simulation results demonstrate ROSCO's peak shaving routine to reduce maximum rotor thrusts by over 10 % compared to the NREL 5 MW reference wind turbine controller on the land-based turbine and to reduce maximum platform pitch angles by nearly 30 % when using the platform feedback routine instead of a more traditional low-bandwidth controller.
17 WIND ENERGY↗
A Revised Calibration Function and Results for the Phoenix Mission TECP Relative Humidity Sensor
The original calibration function of the R(sub H) sensor on the Phoenix Thermal and Electrical Conductivity Sensor (TECP) has been revised in order to extend the range of the valid calibration, and to improve accuracy. The original function returned non-physical R(sub H) values at the lowest temperatures. To resolve this, and because the original calibration was performed against a pair of hygrometers that measured frost point (T(sub f)), the revised calibration equation is also cast in terms of frost point. Because of the complexity of maintaining very low temperatures and high R(sub H) in the laboratory, no calibration data exists at T is greater than 203K. However, sensor response during the mission was smooth and continuous down to 181 K. Therefore we have opted to include flight data in the calibration data set; selection was limited to data acquired during periods when the atmosphere is known to have been saturated. T(sub f) remained below 210 K throughout the mission(P is greater than 0.75 Pa). R(sub H), conversely, ranged from 1 to well under 0.01 diurnally, due to approximately 50 K temperature variations. To first order, both vapor pressure and its variance are greater during daylight hours. Variance in overnight humidity is almost entirely explained by temperature, while atmospheric turbulence contributes substantial variance to daytime humidity. Likewise, data gathered with the TECP aloft reflect higher H2O abundances than at the surface, as well as greater variance. There is evidence for saturation of the atmosphere overnight throughout much of the mission. In virtually every overnight observation, once the atmosphere cooled to T(sub f), water vapor begins to decrease, and tracks air temperature. There is no evidence for substantial decreases in water vapor prior to saturation, as expected for adsorptive exchange. Likewise, there is no evidence of local control of vapor by phases such as perchlorate hydrates hydrated minerals. The daytime average H2O pressure does not change substantially over the course of the mission, although the H2O column abundance varies by a factor of 2. Column abundances calculated from TECP data are lower than orbital measurements if one assumes that H2O is uniformly mixed through a single scale height. These results argue that the vertical distribution of H2O begins to change well in advance of surface concentrations as northern autumn approaches.
Phoenix Thermal and Electrical Conductivity Sensor↗
Multi-energy reconstructions, central electron temperature measurements, and early detection of the birth and growth of runaway electrons using a versatile soft x-ray pinhole camera at MST
A multi-energy soft x-ray pinhole camera has been designed, built, and deployed at the Madison Symmetric Torus to aid the study of particle and thermal transport, as well as MHD stability physics. This novel imaging diagnostic technique employs a pixelated x-ray detector in which the lower energy threshold for photon detection can be adjusted independently on each pixel. The detector of choice is a PILATUS3 100 K with a 450 μm thick silicon sensor and nearly 100 000 pixels sensitive to photon energies between 1.6 and 30 keV. An ensemble of cubic spline smoothing functions has been applied to the line-integrated data for each time-frame and energy-range, obtaining a reduced standard-deviation when compared to that dominated by photon-noise. The multi-energy local emissivity profiles are obtained from a 1D matrix-based Abel-inversion procedure. Central values of Te can be obtained by modeling the slope of the continuum radiation from ratios of the inverted radial emissivity profiles over multiple energy ranges with no a priori assumptions of plasma profiles, magnetic field reconstruction constraints, high-density limitations, or need of shot-to-shot reproducibility. In tokamak plasmas, a novel application has recently been tested for early detection, 1D imaging, and study of the birth, exponential growth, and saturation of runaway electrons at energies comparable to 100 × Te,0; thus, early results are also presented.
Delgado-Aparicio, L. F. (ORCID:0000000167394380)↗
Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches
Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.
58 GEOSCIENCES↗
Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations
The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.
79 ASTRONOMY AND ASTROPHYSICS↗
TEM Analyses of Itokawa Regolith Grains and Lunar Soil Grains to Directly Determine Space Weathering Rates on Airless Bodies
Samples returned from the moon and Asteroid Itokawa by NASA's Apollo Missions and JAXA's Hayabusa Mission, respectively, provide a unique record of their interaction with the space environment. Space weathering effects result from micrometeorite impact activity and interactions with the solar wind. While the effects of solar wind interactions, ion implantation and solar flare particle track accumulation, have been studied extensively, the rate at which these effects accumulate in samples on airless bodies has not been conclusively determined. Results of numerical modeling and experimental simulations do not converge with observations from natural samples. We measured track densities and rim thicknesses of three olivine grains from Itokawa and multiple olivine and anorthite grains from lunar soils of varying exposure ages. Samples were prepared for analysis using a Leica EM UC6 ultramicrotome and an FEI Quanta 3D dual beam focused ion beam scanning electron microscope (FIB-SEM). Transmission electron microscope (TEM) analyses were performed on the JEOL 2500SE 200kV field emission STEM. The solar wind damaged rims on lunar anorthite grains are amorphous, lack inclusions, and are compositionally similar to the host grain. The rim width increases as a smooth function of exposure age until it levels off at approximately 180 nm after approximately 20 My (Fig. 1). While solar wind ion damage can only accumulate while the grain is in a direct line of sight to the Sun, solar flare particles can penetrate to mm-depths. To assess whether the track density accurately predicts surface exposure, we measured the rim width and track density in olivine and anorthite from the surface of rock 64455, which was never buried and has a surface exposure age of 2 My based on isotopic measurements. The rim width from 64455 (60-70nm) plots within error of the well-defined trend for solar wind amorphized rims in Fig. 1. Measured solar flare track densities are accurately reflecting the surface exposure of the grains. Track densities correlate with the amorphous rim thicknesses. While the space-weathered rims of anorthite grains are amorphous, the space-weathered rims on both Itokawa and lunar olivine grains show solar wind damaged rims that are not amorphous. Instead, the rims are nanocrystalline with high dislocation densities and sparse inclusions of nanophase Fe metal. The rim thicknesses on the olivine grains also correlate with track density. The Itokawa olivine grains have track densities that indicate surface exposures of approximately 10(exp 5) years. Longer exposures (up to approximately 10(exp 7) years) do not amorphize the rims, as evidenced by lunar soil olivines with high track densities (approximately 10(exp 11) cm(exp -2)). From the combined data, shown in Fig. 1, it is clear that olivine is damaged (but not amorphized) more rapidly by the solar wind compared to anorthite. The olivine damaged rim forms quickly (in approximately 10(exp 6) y) and saturates at approximately 120nm with longer exposure time. The anorthite damaged rims form more slowly, amorphize, and grow thicker than the olivine rims. This is in agreement with numerical modeling data which predicts that solar wind damaged rims on anorthite will be thicker than olivine. However, the models predict that both olivine and anorthite rims will amorphize and reach equilibrium widths in less than 10(exp 3) y, in contrast to what is observed for natural samples. Laboratory irradiation experiments, which show rapid formation of fully amorphous and blistered surfaces from simulated solar wind exposures are also in contrast to observations of natural samples. These results suggest that there is a flux dependence on the type and extent of irradiation damage that develops in olivine. This flux dependence suggests that great caution be used in extrapolating between high-flux laboratory experiments and the natural case, as demonstrated by. We constrain the space weathering rate through analysis of returned samples. Provided that the track densities and the solar wind damaged rim widths exhibited by the Itokawa grains are typical of the fine-grained regions of Itokawa, then the space weathering rate is on the order of 10(exp 5) y. Space weathering effects in lunar soils saturate within a few My of exposure while those in Itokawa regolith grains formed in approximately 10(exp 5) y. Olivine and anorthite respond differently to solar wind irradiation. The space weathering effects in olivine are particularly difficult to reconcile with laboratory irradiation studies and numerical models. Additional measurements, experiments, and modeling are required to resolve the discrepancies among the observations and calculations involving solar wind amorphization of different minerals on airless bodies.
Berger, Eve L.↗