Bayesian Spatially Varying Multi-regularization Image Deblurring
Presentation to be presented during the MSTS Science and Technology Work-in-Progress Seminar (Webex) scheduled for November 16, 2021.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Engineering Papers
Results for “regularity”
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.
Generalized Regularization of Constrained Optimal Control Problems.
Abstract not provided.
Heidrich, Casey↗
Construction of Differentially Private Empirical Distributions from a low-order Marginals Set through Solving Linear Equations with l2 Regularization.
Abstract not provided.
Eugenio, Evercita Cuevas↗
Assessing arsenic species in foods using regularized linear regression of the arsenic K-edge X-ray absorption near edge structure
Explore the source record for details and available documents.
99 GENERAL AND MISCELLANEOUS↗
Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows
The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.
high pressure reverse osmosis↗
Automated shaker placement and regularized input estimation for MIMO testing
Explore the source record for details and available documents.
Coletti, Keaton Patrick [Sandia National Laborator↗
Physiological identification of morphologically distinct afferent classes innervating the cristae ampullares of the squirrel monkey
1. Semicircular-canal afferents in the squirrel monkey were characterized by their resting discharge, discharge regularity, sensitivity to galvanic currents delivered to the ear (beta *), the gain (g2Hz), and phase lead (phi 2Hz) of their response to 2-Hz sinusoidal head rotations, and their antidromic conduction velocity. Discharge regularity was measured by a normalized coefficient of variation (CV*); the higher the CV*, the more irregular the discharge. g2Hz and phi 2Hz were expressed relative to angular head velocity. 2. These physiological measures were used in an attempt to discern the discharge properties of the three morphological classes of afferents innervating the crista. Presumed bouton (B) fibers were identified as slowly conducting afferents. Presumed calyx (C) fibers were recognized by their irregular discharge and low rotational gains. The remaining fibers were considered to be dimorphic (D) units. Single letters (B, C, and D) are used to emphasize that the classification is based on circumstantial evidence and may be wrong for individual fibers. Of the 125 identified fibers, 13 (10%) were B units, 36 (29%) were C units, and 76 (61%) were D units. 3. B units were regularly discharging D units ranged from regularly to irregularly discharging. C units were the most irregularly discharging afferents encountered. The mean resting discharge for the entire sample was 74 spikes/s. Resting rates were similar for regularly discharging B and D units and higher than those for irregularly discharging C and D units. 4. Except for their lower conduction velocities, the discharge properties of B units are indistinguishable from those of regularly discharging D units. Many of the discharge properties of B and D units vary with discharge regularity. There is a strong, positive relation when beta *, g2Hz, or phi 2Hz is plotted against CV*. For beta * or phi 2Hz, C units conform to the relation for B and D units. In contrast, values of g2Hz for C units are three to four times lower than predicted from the relation for the other two classes. 5. Internal (axon) diameters (dp) of peripheral vestibular-nerve fibers were estimated from central antidromic conduction velocities. Thick fibers (dp > or = 49 microns) were irregularly discharging, mostly C units. Medium-sized fibers (dp = 1.5-4 microns) included regular, intermediate, and irregular D units, as well as C units. Thin fibers (dp < or = 1.5 microns) were defined as B units.(ABSTRACT TRUNCATED AT 400 WORDS).
Non-NASA Center↗
Automatic Aircraft Collision Avoidance System and Method
The invention is a system and method of compressing a DTM to be used in an Auto-GCAS system using a semi-regular geometric compression algorithm. In general, the invention operates by first selecting the boundaries of the three dimensional map to be compressed and dividing the three dimensional map data into regular areas. Next, a type of free-edged, flat geometric surface is selected which will be used to approximate terrain data of the three dimensional map data. The flat geometric surface is used to approximate terrain data for each regular area. The approximations are checked to determine if they fall within selected tolerances. If the approximation for a specific regular area is within specified tolerance, the data is saved for that specific regular area. If the approximation for a specific area falls outside the specified tolerances, the regular area is divided and a flat geometric surface approximation is made for each of the divided areas. This process is recursively repeated until all of the regular areas are approximated by flat geometric surfaces. Finally, the compressed three dimensional map data is provided to the automatic ground collision system for an aircraft.
Skoog, Mark↗
Digital Map Rendering Method
The invention is a system and method of compressing terrain data to be used to render a three-dimensional map representation of the terrain data. In general, the invention operates by first selecting the boundaries of the terrain data and to be compressed and dividing the terrain data into regular geometric areas. Next, a type of free-edged, planar geometric surface is defined which is used to approximate the terrain data for each regular geometric area. The approximations are checked to determine if they fall within user selected tolerances. If the approximation for a specific regular geometric area is within the user specified tolerances, the data is saved for that specific regular geometric area. If the approximation for a specific regular geometric area falls outside the user specified tolerances, the regular geometric area is divided and a free-edged, planar geometric surface approximation is made for each of the divided areas. This process is recursively repeated until all of the regular geometric areas are approximated by free-edged, planar surfaces. Finally, the compressed terrain data is used to render a three-dimensional map.
Skoog, Mark↗
Error Propagation in Microwave Soil Moisture and Vegetation Optical Depth Retrievals
Satellite soil moisture and vegetation optical depth (VOD; related to the total vegetation water mass per unit area) are increasingly being used to study water relations in the soil-plant continuum across the globe. However, soil moisture and VOD are typically jointly estimated, where errors in the optimization approach can cause compensation between both variables and confound such studies. It is thus critical to quantify how satellite microwave measurement errors propagate into soil moisture and VOD. Such a study is especially important for VOD given limited investigations of whether VOD reflects in-situ plant physiology. Furthermore, despite new approaches that constrain (or regularize) VOD dynamics to reduce soil moisture errors, there is limited study of whether regularization reduces VOD errors without obscuring true vegetation temporal dynamics. Here, we find that, across the globe, VOD is less robust to measurement error (more difficult for optimization methods to find the true solution) than soil moisture in their joint estimation. However, a moderate degree of regularization (via time-constrained VOD) reduces errors in VOD to a greater degree than soil moisture and reduces spurious soil moisture-VOD coupling. Furthermore, despite constraining VOD time dynamics, regularized VOD variations on sub-weekly scales are both closer to simulated true VOD time series and have global VOD post-rainfall responses with reduced error signatures compared to VOD retrievals without regularization. Ultimately, we recommend moderately regularized VOD for use in large scale studies of soil-plant water relations because it suppresses noise and spurious soil moisture-VOD coupling without removing the physical signal.
Vegetation optical depth↗
Complement to thermodynamics of dyonic Taub-NUT-AdS spacetime
We examine the thermodynamics of Euclidean dyonic Taub-Nut/Bolt-AdS 4 black holes for a variety of horizon geometries to understand how gauge field regularity conditions influence the thermodynamic relations. We find several distinct features that distinguish the NUT-charged case from its dyonic Reissner-Nordstrom counterpart. For the Nut solution, the gauge field vanishes at the horizon and so regularity is ensured. For the Bolt solution we find that the norm of the gauge field is required to vanish at the horizon in order to satisfy both regularity and the first law of thermodynamics. This regularity condition yields a constraint on the electric and magnetic charges and so reduces cohomogeneity of the system; for spherical horizons, the regularity condition removing the Misner string singularity further reduces cohomogeneity. We observe that bolt solutions with increasing electric charge have positive heat capacity, but upon turning on the magnetic charge to make the solution dyonic, we find that the properties of the uncharged one are retained, having both positive and negative heat capacity. We also study the extremal Bolt solution, finding that Misner string disappears at the horizon in the zero temperature limit. We find that the extremal solution has finite-temperature-like behaviour, with the electric potential playing a role similar to temperature.
79 ASTRONOMY AND ASTROPHYSICS↗
ORNL_AISD_NiNb
This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.
36 MATERIALS SCIENCE↗
Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data
Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.
54 ENVIRONMENTAL SCIENCES↗
SPILC: An expert student advisor
The Lamar University Computer Science Department serves about 350 undergraduate C.S. majors, and 70 graduate majors. B.S. degrees are offered in Computer Science and Computer and Information Science, and an M.S. degree is offered in Computer Science. In addition, the Computer Science Department plays a strong service role, offering approximately sixteen service course sections per long semester. The department has eight regular full-time faculty members, including the Department Chairman and the Undergraduate Advisor, and from three to seven part-time faculty members. Due to the small number of regular faculty members and the resulting very heavy teaching loads, undergraduate advising has become a difficult problem for the department. There is a one week early registration period and a three-day regular registration period once each semester. The Undergraduate Advisor's regular teaching load of two classes, 6 - 8 semester hours, per semester, together with the large number of majors and small number of regular faculty, cause long queues and short tempers during these advising periods. The situation is aggravated by the fact that entering freshmen are rarely accompanied by adequate documentation containing the facts necessary for proper counselling. There has been no good method of obtaining necessary facts and documenting both the information provided by the student and the resulting advice offered by the counsellors.
Read, D. R.↗
Multilevel-in-width training for deep neural network regression
A common challenge in regression is that for many problems, the degrees of freedom required for a high-quality solution also allows for overfitting. Regularization is a class of strategies that seek to restrict the range of possible solutions so as to discourage overfitting while still enabling good solutions, and different regularization strategies impose different types of restrictions. In this paper, we present a multilevel regularization strategy that constructs and trains a hierarchy of neural networks, each of which has layers that are wider versions of the previous network's layers. Here, we draw intuition and techniques from the field of Algebraic Multigrid (AMG), traditionally used for solving linear and nonlinear systems of equations, and specifically adapt the Full Approximation Scheme (FAS) for nonlinear systems of equations to the problem of deep learning. Training through V-cycles then encourage the neural networks to build a hierarchical understanding of the problem. We refer to this approach as multilevel-in-width to distinguish from prior multilevel works which hierarchically alter the depth of neural networks. The resulting approach is a highly flexible framework that can be applied to a variety of layer types, which we demonstrate with both fully connected and convolutional layers. We experimentally show with PDE regression problems that our multilevel training approach is an effective regularizer, improving the generalize performance of the neural networks studied.
97 MATHEMATICS AND COMPUTING↗
A Particle Method for the Multispecies Landau Equation
Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.
Mathematics↗
A chemical kinetic analysis of knock propensity of methanol-to-gasoline fuel
Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. Here, in this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 % vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 % vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).
02 PETROLEUM↗
Janus Monolayers for Ultrafast and Directional Charge Transfer in Transition Metal Dichalcogenide Heterostructures
Charge transfer properties of van der Waals heterostructures formed by Janus and regular transition metal dichalcogenide monolayers are studied by time-resolved pump–probe measurements and photoluminescence spectroscopy. Measurements of electron and hole transfer in three heterostructures with atomic layer sequences of S–W–Se/S–W–S, Se–W–S/S–W–S, and S–W–Se/Se–W–Se reveal that charge transfer from regular to Janus monolayers is ultrafast regardless of the direction of the built-in electric field of the Janus monolayer (Janus field). However, the charge transfer from Janus to regular layers is directional and controlled by the Janus field. When the current direction is along the field, the charge transfer is ultrafast and efficient, while the field blocks the charge transfer with an opposite charge current direction. The transferred carriers form interlayer excitons with extended lifetimes compared to the intralayer excitons. The demonstrated ultrafast and directional charge transfer between Janus and regular TMD layers shows that the Janus structures can be used to make 2D heterostructures with efficient and directional charge transfer properties.
36 MATERIALS SCIENCE↗