Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mathematics & Computing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,729 records · Page 96

Trigonometric polynomial approximation, K-functionals and generalized moduli of smoothness

Best approximation and approximation by families of linear polynomial operators (FLPO) in the spaces L{sub p}, 0<p⩽+∞, are investigated for periodic functions of an arbitrary number of variables in terms of the generalized modulus of smoothness generated by a periodic generator which, near the origin, is assumed to be close in a certain sense to some homogeneous function of positive order. Direct and inverse theorems (Jackson- and Bernstein-type estimates) are proved; conditions on the generators are obtained under which the approximation error by an FLPO is equivalent to an appropriate modulus of smoothness. These problems are solved by going over from the modulus to an equivalent K-functional. The general results obtained are applied to classical objects in the theory of approximation and smoothness. In particular, they are applied to the methods of approximation generated by Fejér, Riesz and Bochner-Riesz kernels, and also to the moduli of smoothness and K-functionals corresponding to the conventional, Weyl and Riesz derivatives and to the Laplace operator. Bibliography: 24 titles. (paper)

97 MATHEMATICS AND COMPUTING↗

Weighted inequalities for quasilinear integral operators on the semi-axis and applications to Lorentz spaces

A precise characterization of inequalities in weighted Lebesgue spaces with positive quasilinear integral operators of iterative type on the half-axis is given. All cases of positive integration parameters are treated, including the case of supremum. Applications to the solution of the well-known problem of the boundedness of the Hardy-Littlewood maximal operator in weighted Lorentz Γ-spaces are given. Bibliography: 41 titles.

97 MATHEMATICS AND COMPUTING↗

On some problems related to the Hilbert-Smith conjecture

The Hilbert-Smith conjecture claims that if a compact group G acts freely on a manifold, then it is a Lie group. For a finite-dimensional orbit space a reduction of the Hilbert-Smith conjecture to certain other problems in geometric topology is presented; in these the key problem is the existence of an essential sequence of lens spaces of increasing dimension. Bibliography: 52 titles. (paper)

97 MATHEMATICS AND COMPUTING↗

Necessary and sufficient conditions for the existence and uniqueness of a bounded solution of the equation ((dx(t))/(dt))=f(x(t)+h{sub 1}(t))+h{sub 2}(t)

Necessary and sufficient conditions for a bounded solution of the nonlinear scalar differential equation dx(t)/dt=f(x(t)+h{sub 1}(t))+h{sub 2}(t), t∈R, to exist and be unique are presented in the case when f(x) is a continuous function and the functions h{sub 1}(t) and h{sub 2}(t) are bounded and continuous. The case when h{sub 1}(t) and h{sub 2}(t) are almost periodic functions is also investigated. Bibliography: 31 titles. (paper)

97 MATHEMATICS AND COMPUTING↗

Collective Dispersal Leads to Variance in Fitness and Maintains Offspring Size Variation within Marine Populations

Variance in fitness is well known to influence the outcome of evolution but is rarely considered in the theory of marine reproductive strategies. In coastal environments, turbulent mesoscale eddies can collect larvae into packets, resulting in collective dispersal. Larvae in packets return to the coast or are lost offshore in groups, producing variance in fitness. Using a Markov process to calculate fixation probabilities for competing phenotypes, we examine the evolution of offspring size and spawning duration in species with benthic adults and pelagic offspring. The offspring size that provides mothers with the highest mean fitness also generates the greatest variance in fitness, but pairwise invasion plots show that bet-hedging strategies are not evolutionarily stable; maximizing expected fitness correctly predicts the unique evolutionarily stable strategy. Nonetheless, fixation can take a long time. We find that selection to increase spawning duration as a risk avoidance strategy to reduce the negative impacts of stochastic recruitment success can allow multiple offspring sizes to coexist in a population for extended periods. This has two important consequences for offspring size: (1) coexistence occurs over a broader range of sizes and is longer when spawning duration is longer because longer spawning durations reduce variation in fitness and increase the time to fixation, and (2) longer spawning durations can compensate for having a nonoptimal size and even allow less optimal sizes to reach fixation. Collective dispersal and longer spawning durations could effectively maintain offspring size variation even in the absence of good and bad years or locations. Finally, empirical comparisons of offspring size would therefore not always reflect environment-specific differences in the optimal size.

59 BASIC BIOLOGICAL SCIENCES↗

Theoretical and numerical studies of inverse source problem for the linear parabolic equation with sparse boundary measurements

We consider the inverse source problem in the parabolic equation, where the unknown source possesses the semi-discrete formulation. Theoretically, we prove that the flux data from any nonempty open subset of the boundary can uniquely determine the semi-discrete source. This means the observed area can be extremely small, and that is the reason we call it sparse boundary data. For the numerical reconstruction, we formulate the problem from the Bayesian sequential prediction perspective and conduct the numerical examples which estimate the space-time-dependent source state by state. To better demonstrate the method’s performance, we solve two common multiscale problems from two models with a long source sequence. The numerical results illustrate that the inversion is accurate and efficient.

97 MATHEMATICS AND COMPUTING↗

Approximate symmetries of guiding-centre motion

In a strong, inhomogeneous magnetic field, charged particle dynamics may be studied in the guiding-centre approximation, which is known to be Hamiltonian. When the magnetic field is quasisymmetric, the first-order guiding-centre (FGC) Hamiltonian structure admits a continuous symmetry, and therefore a conserved quantity in addition to the energy. Since the FGC system is only an approximation, it is also interesting to consider approximate symmetries of the guiding-centre Hamiltonian structure. We find that any approximate spatial symmetry coincides with quasisymmetry to leading order. For approximate phase-space symmetries, we derive weaker conditions than quasisymmetry. The latter include 'weak quasisymmetry' as a subcase, recently proposed by Rodríguez et al. Our results, however, show that weak quasisymmetry is necessarily non-spatial at first order. Finally, we demonstrate that if the magnetic field is constrained to satisfy magnetohydrostatic force balance then an approximate symmetry must agree with quasisymmetry to leading order.

97 MATHEMATICS AND COMPUTING↗

Sparsifying priors for Bayesian uncertainty quantification in model discovery

We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.

97 MATHEMATICS AND COMPUTING↗

Efficient molecular dynamics using geodesic integration and solvent–solute splitting

We present an approach to Langevin dynamics in the presence of holonomic constraints based on decomposition of the system into components representing geodesic flow, constrained impulse and constrained diffusion. We show that a particular ordering of the components results in an integrator that is an order of magnitude more accurate for configurational averages than existing alternatives. Moreover, by combining the geodesic integration method with a solvent–solute force splitting, we demonstrate that stepsizes of at least 8 fs can be used for solvated biomolecules with high sampling accuracy and without substantially altering diffusion rates, approximately increasing by a factor of two the efficiency of molecular dynamics sampling for such systems. The methods described in this article are easily implemented using the standard apparatus of modern simulation codes.

97 MATHEMATICS AND COMPUTING↗

Peak effect and dynamics of stripe- and pattern-forming systems on a periodic one-dimensional substrate

Here we examine the ordering, pinning, and dynamics of two-dimensional pattern-forming systems interacting with a periodic one-dimensional substrate. In the absence of the substrate, particles with competing long-range repulsion and short-range attraction form anisotropic crystal, stripe, and bubble states. When the system is tuned across the stripe transition in the presence of a substrate, we find that there is a peak effect in the critical depinning force when the stripes align and become commensurate with the substrate. Under an applied drive, the anisotropic crystal and stripe states can exhibit soliton depinning and plastic flow. When the stripes depin plastically, they dynamically reorder into a moving stripe state that is perpendicular to the substrate trough direction. We also find that when the substrate spacing is smaller than the widths of the bubbles or stripes, the system forms pinned stripe states that are perpendicular to the substrate trough direction. The system exhibits multiple reentrant pinning effects as a function of increasing attraction, with the anisotropic crystal and large bubble states experiencing weak pinning but the stripe and smaller bubble states showing stronger pinning. We map out the different dynamic phases as a function of filling, the strength of the attractive interaction term, the substrate strength, and the drive, and demonstrate that the different phases produce identifiable features in the transport curves and particle orderings.

97 MATHEMATICS AND COMPUTING↗

Autonomy Loops for Monitoring, Operational Data Analytics, Feedback, and Response in HPC Operations

Many High Performance Computing (HPC) facilities have developed and deployed frameworks in support of continuous monitoring and operational data analytics (MODA) to help improve efficiency and throughput. Because of the complexity and scale of systems and workflows and the need for low-latency response to address dynamic circumstances, automated feedback and response have the potential to be more effective than current human-in-the-loop approaches which are laborious and error prone. Progress has been limited, however, by factors such as the lack of infrastructure and feedback hooks, and successful deployment is often site- and case-specific. In this position paper we report on the outcomes and plans from a recent Dagstuhl Seminar, seeking to carve a path for community progress in the development of autonomous feedback loops for MODA, based on the established formalism of similar (MAPE-K) loops in autonomous computing and self-adaptive systems. By defining and developing such loops for significant cases experienced across HPC sites, we seek to extract commonalities and develop conventions that will facilitate interoperability and interchangeability with system hardware, software, and applications across different sites, and will motivate vendors and others to provide telemetry interfaces and feedback hooks to enable community development and pervasive deployment of MODA autonomy loops.

autonomy loops↗

DEEP NEURAL NETWORKS FOR NON-LINEAR MODEL-BASED ULTRASOUND RECONSTRUCTION

Ultrasound reflection tomography is widely used to image large complex specimens that are only accessible from a single side, such as well systems and nuclear power plant containment walls. Typical methods for inverting the measurement rely on delay-and-sum algorithms that rapidly produce reconstructions but with significant artifacts. Recently, model-based reconstruction approaches using a linear forward model have been shown to significantly improve image quality compared to the conventional approach. Nonetheless, even these techniques result in artifacts for complex objects because of the inherent non-linearity of the ultrasound forward model.In this paper, we propose a non-iterative model-based reconstruction method for inverting measurements that are based on non-linear forward models for ultrasound imaging. Our approach involves obtaining an approximate estimate of the reconstruction using a simple linear back-projection and training a deep neural network to refine this to the actual reconstruction. As a result, we apply our method to simulated and experimental ultrasound data to demonstrate dramatic improvements in image quality compared to the delay-and-sum approach and the linear model-based reconstruction approach.

97 MATHEMATICS AND COMPUTING↗

SHARP-Net: Platform for Self-Healing and Attack Resilient PMU Networks

Synchrophasor technology plays a pivotal role in developing the next generation of wide-area monitoring, protection, and control in the smart grid environment. As technology and communications infrastructures evolve, however, so do the attack surfaces in the synchrophasor network that can be exploited by advanced persistent threat (APT) actors to affect power system stability and reliability. In this paper, we propose a novel platform for developing a self-healing and attack-resilient PMU network (SHARP-Net) by instituting a state-of-the-art intrusion detection system (IDS) with an intrusion mitigation system (IMS) and an alert management system (AMS). In particular, the proposed platform detects anomalies during cyberattacks on phasor data concentrators (PDCs) based on the rules defined in the IDS, then the generated alerts are published to the IMS through the AMS. The proposed IMS proceeds to take automated corrective responses to mitigate cyberattacks by reconfiguring the synchrophasor network to isolate the compromised PDCs, and it orchestrates new PDCs to prevent the future propagation of attacks. Further, the IMS restores the system's observability by reconnecting the new PDCs to make the grid attack-resilient. In this work, the SHARP-Net platform is developed by using Python-based libraries, minimega's software-defined network, and virtual machine orchestration. We implement and validate the proposed SHARP-Net architecture by testing a PMU network in the smart grid environment. SHARP-Net showed promising performance in detecting cyberattacks and mitigating them through the network reconfiguration.

computer architecture↗

Bias Correction and Statistical Downscaling of Future Solar Irradiance Projections Using the NSRDB

Assessing renewable energy resources under future climate scenarios has been highlighted to understand potential impacts of future climate change in renewable generation on the power sector. Climate model projection has been recognized by the renewable energy community as a useful data set to analyze the impacts of future climate change on renewable resources. However, future climate projections generated from general circulation models (GCMs) contain inherent biases that need to be corrected for accurate analysis of future projections of climate variables. In addition, the coarse spatiotemporal resolution of GCMs needs to be improved for regional climate studies. In this work, we develop statistical methods to downscale future projections of global horizontal irradiance (GHI) in a computationally efficient way. Our approach builds statistical downscaling models that correct bias of climate projection of GHI and downscale the future GHI projection from daily-scale to hourly-scale. The National Solar Radiation Database (NSRDB) is used to calibrate the statistical models and validate the downscaled GHI projections across the contiguous United State (CONUS). Preliminary results show that the statistical approach efficiently downscales climate projections of GHI with a nBIAS of 3%, nMAE of 34 % and nRMSE of 46% calculated against NSRDB for CONUS. This study describes the implemented methodology and initial results as well as future research to create high-resolution climate data sets for solar energy applications.

analytical models↗

Characterization, Analysis, and Implementation of Integrated Bandstop Structures on Ultra-Wideband Archimedean Spiral Antenna

A subsection of the radiating spiral arm, placed in parallel, induces a bandstop response at a notch frequency proportional to the resonant length of the strip. Detailed parametric study on the effect of variation of design parameters for an Archimedean spiral antenna and the resonant parallel strip (RPS) is presented. Empirical analysis on phase velocity on the radiating spiral arms allows characterization of RPS in terms of its resonant length. Identified systematic relation between design parameters and filter response is applied to design an antenna for the 3.1-10.5 GHz operating band with the notch response over the IEEE 802.11a band, 5.15-5.95 GHz. Successful implementation is then demonstrated through performance comparison between simulated and experimental results.

42 ENGINEERING↗

Monaural Music Source Separation using Convolutional Sparse Coding

We present a comprehensive performance study of a new time-domain approach for estimating the components of an observed monaural audio mixture. Unlike existing time-frequency approaches that use the product of a set of spectral templates and their corresponding activation patterns to approximate the spectrogram of the mixture, the proposed approach uses the sum of a set of convolutions of estimated activations with prelearned dictionary filters to approximate the audio mixture directly in the time domain. The approximation problem can be solved by an efficient convolutional sparse coding algorithm. The effectiveness of this approach for source separation of musical audio has been demonstrated in our prior work, but under rather restricted and controlled conditions, requiring the musical score of the mixture being informed a priori and little mismatch between the dictionary filters and the source signals. In this paper, we report an evaluation that considers wider, and more practical, experimental settings. This includes the use of an audio-based multipitch estimation algorithm to replace the musical score, and an external dataset of audio single notes to construct the dictionary filters. Here, our result shows that the proposed approach remains effective with a larger dictionary, and compares favorably with the state-of-the-art nonnegative matrix factorization approach. However, in the absence of the score and in the case of a small dictionary, our approach may not be better.

97 MATHEMATICS AND COMPUTING↗