Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Order parameters”

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 811 records · Page 45

Analysis of Air Traffic Track Data with the AutoBayes Synthesis System

The Next Generation Air Traffic System (NGATS) is aiming to provide substantial computer support for the air traffic controllers. Algorithms for the accurate prediction of aircraft movements are of central importance for such software systems but trajectory prediction has to work reliably in the presence of unknown parameters and uncertainties. We are using the AutoBayes program synthesis system to generate customized data analysis algorithms that process large sets of aircraft radar track data in order to estimate parameters and uncertainties. In this paper, we present, how the tasks of finding structure in track data, estimation of important parameters in climb trajectories, and the detection of continuous descent approaches can be accomplished with compact task-specific AutoBayes specifications. We present an overview of the AutoBayes architecture and describe, how its schema-based approach generates customized analysis algorithms, documented C/C++ code, and detailed mathematical derivations. Results of experiments with actual air traffic control data are discussed.

Schumann, Johann Martin Philip↗

Multi-Fidelity Scheme for Accelerating 4D Finite Element Analysis for Mircoreactors

Next-generation microreactors are currently being designed to be operated terrestrial and extraterrestrial for remote surface power production. These systems will provide an alternative source of carbon-free energy that is versatile and can be utilized for various applications. These applications of microreactors have prompted the usage of high-fidelity unstructured finite element (FE) based approaches to provide time-dependent solutions for multiple physics fields. Computing these high-fidelity (4-D) solutions for multiple design iterations, physics fields, and transient events requires an immense computational cost. These high-fidelity solutions are computational expensive and the cost can be reduced through the implementation of less accurate low-fidelity solutions. Unlike the current fleet of commercial nuclear reactors, these next-gen systems present challenges due to the material and physical limitations required. Due to these constraints, the brute force technique of parameterizing important system characteristics determines whether a design meets the project objective. A system such as a nuclear reactor could have thousands of design parameters that affect system performance. In order to analyze the entire parameter space of such a complex system would require millions of CPU hours and countless design iterations. This costly approach is not practical due to regulatory and budget limitations. In this proposal, an approach that utilizes hybrid high-fidelity and low-fidelity FE models to reduce the computational cost of evaluating these applications in 4-D will be presented. The accuracy of the high-fidelity model and the computational efficiency of the low-fidelity model are taken advantage of to produce a solution that closely resembles the full order high-fidelity solution. By utilizing an unconverged coarse FE mesh, operation limits such as temperature, structural loading, etc., can be evaluated in an accelerated fashion and then can be spatially interpolated onto a finer mesh. The resulting error arising from the coarse mesh can be mitigated with a discrepancy function that actively quantifies and corrects the error in the coarse solution. This discrepancy function can be periodically updated with high-fidelity calculations across the temporal domain thus requiring less iterations on the fine mesh. As a result, the computational cost can be reduced for the evaluation and design iteration of microreactors. The proposal for this research contains three sections: Section 2 provides a literature review, Section 3 outlines the methodology for the proposed multi-fidelity scheme, and Section 4 displays preliminary results of the proposed research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Genetic Algorithm-Guided, Adaptive Model Order Reduction of Flexible Aircrafts

This paper presents a methodology for automated model order reduction (MOR) of flexible aircrafts to construct linear parameter-varying (LPV) reduced order models (ROM) for aeroservoelasticity (ASE) analysis and control synthesis in broad flight parameter space. The novelty includes utilization of genetic algorithms (GAs) to automatically determine the states for reduction while minimizing the trial-and-error process and heuristics requirement to perform MOR; balanced truncation for unstable systems to achieve locally optimal realization of the full model; congruence transformation for "weak" fulfillment of state consistency across the entire flight parameter space; and ROM interpolation based on adaptive grid refinement to generate a globally functional LPV ASE ROM. The methodology is applied to the X-56A MUTT model currently being tested at NASA/AFRC for flutter suppression and gust load alleviation. Our studies indicate that X-56A ROM with less than one-seventh the number of states relative to the original model is able to accurately predict system response among all input-output channels for pitch, roll, and ASE control at various flight conditions. The GA-guided approach exceeds manual and empirical state selection in terms of efficiency and accuracy. The adaptive refinement allows selective addition of the grid points in the parameter space where flight dynamics varies dramatically to enhance interpolation accuracy without over-burdening controller synthesis and onboard memory efforts downstream. The present MOR framework can be used by control engineers for robust ASE controller synthesis and novel vehicle design.

Numerical Analysi↗

Relativistic Hartree–Fock chiral Lagrangians with confinement, nucleon finite size and short-range effects

A relativistic Hartree–Fock Lagrangian including a chiral potential and nucleon polarisation is investigated in hopes of providing a better description of dense nuclear matter. Here, we fully consider the contribution of the exchange Fock term to the energy and the self-energies, and in addition we investigate the nucleon’s compositeness and finite size effects (confinement and form factors) and short range correlations modeled by a Jastrow ansatz. These effects are added step by step, such that their impact on the dense matter properties can be analysed in details. The parameters of the model are adjusted to reproduce fundamental properties related to the QCD theory at low energy, such as the chiral symmetry breaking, nucleon’s quark substructure and Lattice-QCD predictions, as well as two empirical properties at saturation: the binding energy and the density. All other empirical parameters, e.g., symmetry energy and its slope, incompressibility modulus, effective mass, as well as spin–isospin Landau–Midgal parameter are predictions of the models and can be used to evaluate the gain of the different approximation schemes in describing nuclear properties. Bayesian statistics is employed in order to propagate parameter uncertainties into predictions for the nuclear matter properties. We show that the splitting of the effective Landau mass is largely influenced by the value of the ρ T coupling, and we show that the fit to the symmetry energy, which induces an increase of the coupling constant g ρ by about 20–25% compared to the case where it is fixed by the quark model, provides a very good EoS compatible with the present nuclear physics knowledge.

74 ATOMIC AND MOLECULAR PHYSICS↗

Analytic expressions for perturbations and partial derivatives of range and range rate of a spacecraft with respect to the coefficient of the second harmonic

Closed-form analytic expressions for the time variations of instantaneous orbital parameters and of the topocentric range and range rate of a spacecraft moving in the gravitational field of an oblate large body are derived using a first-order variation of parameters technique. In addition, the closed-form analytic expressions for the partial derivatives of the topocentric range and range rate are obtained, with respect to the coefficient of the second harmonic of the potential of the central body (J sub 2). The results are applied to the motion of a point-mass spacecraft moving in the orbit around the equatorially elliptic, oblate sun, with J sub 2 approximately equal to .000027.

Georgevic, R. M.↗

Effectively Using Remote I/O For Work Composition in Distributed Workflows

Distributed scientific workflows are becoming more important with the interest in incorporating AI into their loops. A critical programming and performance question is how to compose workflow tasks when data is produced on one system but must be consumed on another. Since the dominant technique is composition with remote I/O, this paper explores its performance expectations. We describe BigFlowSim, a workflow I/O simulator that captures key implementation choices for remote I/O, including intensity, reuse, locality, access pattern, and data movement.With BigFlowSim, we generate a synthetic benchmark. We quantify the effects of each parameter with a performance sensitivity study. We explain trends in terms of data movement reduction and show that, under certain conditions, it is possible to establish a total order among most parameters. We apply these insights to a high energy physics workflow, Belle II Monte Carlo and simulate several I/O optimizations. Speedups range from 5% to 2×, without changing compute time.

Friese, Ryan D.↗

Constructing chiral effective field theory around unnatural leading-order interactions

A momentum-dependent formulation based on a stationary spin-0 and isospin-1 dibaryon field is proposed to improve convergence of chiral effective field theory in the 1 S 0 channel of NN scattering. Although the two-parameter leading-order interaction appears to be unnatural, it nevertheless has the necessary features of an effective field theory. A rapid order-by-order convergence is found in 1 S 0 . As an application beyond the two-body level, the triton binding energy is studied and compared to standard chiral effective field theory with partly perturbative pions. Here, the consistency of the chiral Lagrangian for the new formulation is examined by working out the pionic radiative corrections, and consequences of nontrivial chiral-connection terms are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Characterization data of an (AlFeNiTiVZr) 1-x Cr x multi-principal element alloy continuous composition spread library

The data provided in this article is related to the research article entitled “Phase stabilization and oxidation of a continuous composition spread multi-principal element (AlFeNiTiVZr) 1-x Cr x alloy”. This data article describes the high-throughput synthesis and characterization processes of an (AlFeNiTiVZr) 1-x Cr x alloy system. Continuous composition spread (CCS) thin-film libraries were synthesized by co-depositing an AlFeNiTiVZr metal alloy target and Cr target via magnetron sputtering. Post-processing was performed on the sample libraries with a vacuum anneal at 873 K and an air anneal at 873 K. Compositional data was determined via WDS in order to verify parameters provided by an in-house sputter model. Crystallographic data was captured via synchrotron diffraction and diffractograms were compared as a function of the change in Cr concentration. These measurements were taken in order to observe phase behavior after oxidation throughout the composition library. Furthermore, vibrational spectrographic data is provided of the oxidized library to show surface speciation along the composition gradient of the alloy system. The structural and oxidative behavior of the (AlFeNiTiVZr) 1-x Cr x alloy can be analysed using the data provided in this article. Additionally, this characterization dataset can be utilized in machine learning algorithms for determining important features and parameters for future hypothesis generation of functional multi-principal element alloys (MPEAs).

36 MATERIALS SCIENCE↗

Adaptive tracking for complex systems using reduced-order models

Reduced-order models are considered in the context of parameter adaptive controllers for tracking workspace trajectories. A dual-arm manipulation task is used to illustrate the methodology and provide simulation results. A parameter adaptive controller is designed to track the desired position trajectory of a payload using a four-parameter model instead of a full-order, nine-parameter model. Several simulations with different payload-to-arm mass ratios are used to illustrate the capabilities of the reduced-order model in tracking the desired trajectory.

Carignan, Craig R.↗

Adaptive tracking for complex systems using reduced-order models

Reduced-order models are considered in the context of parameter adaptive controllers for tracking workspace trajectories. A dual-arm manipulation task is used to illustrate the methodology and provide simulation results. A parameter adaptive controller is designed to track a payload trajectory using a four-parameter model instead of the full-order, nine-parameter model. Several simulations with different payload-to-arm mass ratios are used to illustrate the capabilities of the reduced-order model in tracking the desired trajectory.

Carnigan, Craig R.↗

Learning and meta-learning of stochastic advection–diffusion–reaction systems from sparse measurements

Physics-informed neural networks (PINNs) were recently proposed in [18] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution, while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a non-linear advection–diffusion–reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimise the hyper-parameters of sPINN by using the Bayesian optimisation method (meta-learning) and compare the results with the empirically selected hyper-parameters of sPINN. In particular, for the first part in solving the inverse deterministic ADR, we assume that we only have a few high-fidelity measurements, whereas the rest of the data is of lower fidelity. Hence, the PINN is trained using a composite multi-fidelity network, first introduced in [12], that learns the correlations between the multi-fidelity data and predicts the unknown values of diffusivity, transport velocity and two reaction constants as well as the concentration field. For the stochastic ADR, we employ a Karhunen–Loève (KL) expansion to represent the stochastic diffusivity, and arbitrary polynomial chaos (aPC) to represent the stochastic solution. Correspondingly, we design multiple NNs to represent the mean of the solution and learn each aPC mode separately, whereas we employ a separate NN to represent the mean of diffusivity and another NN to learn all modes of the KL expansion. For the inverse problem, in addition to stochastic diffusivity and concentration fields, we also aim to obtain the (unknown) deterministic values of transport velocity and reaction constants. The available data correspond to 7spatial points for the diffusivity and 20 space–time points for the solution, both sampled 2000 times. We obtain good accuracy for the deterministic parameters of the order of 1–2% and excellent accuracy for the mean and variance of the stochastic fields, better than three digits of accuracy. In the second part, we consider the previous stochastic inverse problem, and we use Bayesian optimisation to find five hyper-parameters of sPINN, namely the width, depth and learning rate of two NNs for learning the modes. Here, we obtain much deeper and wider optimal NNs compared to the manual tuning, leading to even better accuracy, i.e., errors less than 1% for the deterministic values, and about an order of magnitude less for the stochastic fields.

97 MATHEMATICS AND COMPUTING↗

Thermalization and annihilation of dark matter in neutron stars

The capture of dark matter, and its subsequent annihilation, can heat old, isolated neutron stars. In order for kinetic heating to be achieved, the captured dark matter must undergo sufficient scattering to deposit its kinetic energy in the star. We find that this energy deposit typically occurs quickly, for most of the relevant parameter space. In order for appreciable annihilation heating to also be achieved, the dark matter must reach a state of capture-annihilation equilibrium in the star. We show that this can be fulfilled for all types of dark matter - baryon interactions. This includes cases where the scattering or annihilation cross sections are momentum or velocity suppressed in the non-relativistic limit. Importantly, we find that capture-annihilation equilibrium, and hence maximal annihilation heating, can be achieved without complete thermalization of the captured dark matter. For scattering cross sections that saturate the capture rate, we find that capture-annihilation equilibrium is typically reached on a timescale of less than 1 year for vector interactions and 10 4 years for scalar interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Novel use of Direct Simulation Monte-Carlo to Model Dynamics of COVID-19 Pandemic Spread

In this report, we evaluate a novel method for modeling the spread of COVID-19 pandemic. In this new approach we leverage methods and algorithms developed for fully-kinetic plasma physics simulations using Particle-In-Cell (PIC) Direct Simulation Monte-Carlo (DSMC) models. This approach then leverages Sandia-unique simulation capabilities, and High-Performance Computer (HPC) resources and expertise in particle-particle interactions using stochastic processes. Our hypothesis is that this approach would provide a more efficient platform with assumptions based on physical data that would then enable the user to assess the impact of mitigation strategies and forecast different phases of infection. This work addresses key scientific questions related to the assumptions this new approach must make to model the interactions of people using algorithms typically used for modeling particle interactions in physics codes (kinetic plasma, gas dynamics). The model developed uses rational/physical inputs while also providing critical insight; the results could serve as inputs to, or alternatives for, existing models. The model work presented was developed over a four-week time frame, thus far showing promising results and many ways in which this model/approach could be improved. This work is aimed at providing a proof-of-concept for this new pandemic modeling approach, which could have an immediate impact on the COVID-19 pandemic modeling, while laying a basis to model future pandemic scenarios in a manner that is timely and efficient. Additionally, this new approach provides new visualization tools to help epidemiologists comprehend and articulate the spread of this and other pandemics as well as a more general tool to determine key parameters needed in order to better predict pandemic modeling in the future. In the report we describe our model for pandemic modeling, apply this model to COVID-19 data for New York City (NYC), assess model sensitivities to different inputs and parameters and , finally, propagate the model forward under different conditions to assess the effects of mitigation and associated timing. Finally, our approach will help understand the role of asymptomatic cases, and could be extended to elucidate the role of recovered individuals in the second round of the infection, which is currently being ignored.

59 BASIC BIOLOGICAL SCIENCES↗

Competing orders at higher-order Van Hove points

Van Hove points are special points in the energy dispersion, where the density of states exhibits analytic singularities. When a Van Hove point is close to the Fermi level, tendencies towards density wave orders, Pomeranchuk orders, and superconductivity can all be enhanced, often in more than one channel, leading to a competition between different orders and unconventional ground states. In this study we consider the effects from higher-order Van Hove points, around which the dispersion is flatter than near a conventional Van Hove point, and the density of states has a power-law divergence. We argue that such points are present in intercalated graphene and other materials. We use an effective low-energy model for electrons near higher-order Van Hove points and analyze the competition between different ordering tendencies using an unbiased renormalization-group approach. For purely repulsive interactions, we find that two key competitors are ferromagnetism and chiral superconductivity. For a small attractive interaction, we find an unconventional spin Pomeranchuk order, in wich the spin oder parameter winds around the Fermi surface. The supermetal state, predicted for a single higher-order Van Hove point, is an unstable fixed point in our case.

36 MATERIALS SCIENCE↗

Cosmological forecasts for future galaxy surveys with the linear point standard ruler: Toward consistent BAO analyses far from a fiducial cosmology

The linear point (LP) standard ruler was identified as the basis of a purely geometric method for exploiting the baryon acoustic oscillations (BAO). The LP exploits the BAO feature imprinted in the galaxy two-point correlation function to measure cosmological distances independent of any specific cosmological model. We forecast the expected precision of future and ongoing spectroscopic galaxy surveys to measure distances leveraging the linear point. We investigate the cosmological implications of our forecasted results. We focus in particular on a relevant working example: the detection of the late-time cosmic acceleration independent of other cosmological probes. Further, our findings show that, even within the Λ⁢CDM standard cosmological paradigm, estimated distances need to be reliable over a very wide parameter range in order to realize their maximum utility. This is particularly relevant if we aim to properly characterize cosmological tensions. The LP is a promising candidate approach to achieve this reliability. In contrast, widely employed procedures in BAO analysis estimate distances keeping fixed cosmological parameters to fiducial values close to cosmic-microwave-background constraints in flat-Λ⁢CDM. It is unclear whether they are purely geometric methods. Moreover, they rely on untested extrapolations to explore the parameter space away from those fiducial flat-Λ⁢CDM values. We recommend that all BAO methodologies be validated across the full range of models and parameters over which their results are quoted, first by means of linear predictions and then N -body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Periodic solutions of second-order nonlinear difference equations containing a small parameter. II - Equivalent linearization

The classical method of equivalent linearization is extended to a particular class of nonlinear difference equations. It is shown that the method can be used to obtain an approximation of the periodic solutions of these equations. In particular, the parameters of the limit cycle and the limit points can be determined. Three examples illustrating the method are presented.

Mickens, R. E.↗