Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Regularization”

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 271 records · Page 15

From Latent Dynamics to Meaningful Representations

While representation learning has been central to the rise of machine learning and artificial intelligence, a key problem remains in making the learnt representations meaningful. For this the typical approach is to regularize the learned representation through prior probability distributions. However such priors are usually unavailable or are ad hoc. To deal with this, recent efforts have shifted towards leveraging the insights from physical principles to guide the learning process. In this spirit, we propose a purely dynamics-constrained representation learning framework. Instead of relying on predefined probabilities, we restrict the latent representation to follow overdamped Langevin dynamics with a learnable transition density — a prior driven by statistical mechanics. We show this is a more natural constraint for representation learning in stochastic dynamical systems, with the crucial ability to uniquely identify the ground truth representation. We validate our framework for different systems including a real-world fluorescent DNA movie dataset. Here, we show that our algorithm can uniquely identify orthogonal, isometric and meaningful latent representations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spacecraft Availability Enhancement by In-flight Testing of Spare Parts

This paper describes an approach to improve availability by testing redundant parts at pre-determined intervals. The purpose of the testing is to detect nonfunctional back-up equipment and develop work-around measures or replacement spacecraft before a failure of the primary equipment reduces availability. The work reported here is an outgrowth of the NASA Space Network's program for the maintenance and replenishment of the Tracking and Data Relay Satellite System. The approach is based on a standby factor and a cyclic stress factor. The standby factor accounts for the effects of adverse storage conditions encountered as part of the in-flight environment. The stress factor accounts for the effects of physical or thermal cycling due to the application of force or power that characterizes the operation or use of a component. By the quantitative consideration of standby and cycling risks, a regular testing interval can be calculated - an interval for testing that furnishes information for availability planning but does not subject the spacecraft to undue risk. This paper includes quantitative solutions for appropriate testing intervals for equipment configurations and failure rates that are representative of a Tracking and Relay Data Satellite. The effect of a single test on the availability of certain equipment is also illustrated.

Cutler, Robert M.↗

Persistent Assets in Zero-G and on Planetary Surfaces: Enabled by Modular Technology and Robotic Operations

Space operations are on the cusp of a revolutionary new operational paradigm that leverages modular systems and recurring robotic visits to “Persistent Assets” enabling asset maintenance, repair, and enhancement. A “Persistent Asset” is defined here as any zero-g or planetary surface system that benefits from in-space assembly (ISA) or multiple visits for servicing, repairs, and upgrades. This term is an extension of the term “Persistent Platform” used by Ms. Pam Melroy at the Defense Advanced Research Projects Agency to describe the vision of the Agency for a geosynchronous Earth orbiting platform. In this paper, the term “Persistent Asset” is introduced to encompass not only zero-g systems; such as telecommunication platforms, Earth observing science platform, Department of Defense platforms, and scientific telescope systems, but also planetary surface systems that support missions such as human outposts, science stations, and in-situ resource utilization systems. In contrast to the current state of the art, where space systems are typically launched as a single unit and operated without any further physical intervention after launch; future systems will be maintained, enhanced and reconfigured in-situ as new technology becomes available or mission needs change. Visits to the persistent asset can be regularly scheduled or dictated by funding constraints enabling a pay-as-you-go approach, which is largely independent of time-constraints and able to exploit launches of opportunity. In this paper, historical in-space assembly activities which relied heavily on astronaut extra-vehicular activity will be reviewed as well as early robotic assembly activities. These approaches will be contrasted with emerging modular approaches supporting realization of a new Persistent Asset operational paradigm. The paper will define attributes of the Persistent Asset paradigm and illustrate advantages by applying the paradigm to two relevant applications: 1) a large space telescope backing structure and 2) the backbone structure for a solar electric transport vehicle. Finally, recently developed unique Persistent Asset elements (modules, and interface approaches) will be described.

Doggett, William R.↗

Persistent Assets in Zero-G and on Planetary Surfaces: Enabled by Modular Technology and Robotic Operations

Space operations are on the cusp of a revolutionary new operational paradigm that leverages modular systems and recurring robotic visits to “Persistent Assets” enabling asset maintenance, repair, and enhancement. A “Persistent Asset” is defined here as any zero-g or planetary surface system that benefits from in-space assembly (ISA) or multiple visits for servicing, repairs, and upgrades. This term is an extension of the term “Persistent Platform” used by Ms. Pam Melroy at the Defense Advanced Research Projects Agency to describe the vision of the Agency for a geosynchronous Earth orbiting platform. In this paper, the term “Persistent Asset” is introduced to encompass not only zero-g systems; such as telecommunication platforms, Earth observing science platform, Department of Defense platforms, and scientific telescope systems, but also planetary surface systems that support missions such as human outposts, science stations, and in-situ resource utilization systems. In contrast to the current state of the art, where space systems are typically launched as a single unit and operated without any further physical intervention after launch; future systems will be maintained, enhanced and reconfigured in-situ as new technology becomes available or mission needs change. Visits to the persistent asset can be regularly scheduled or dictated by funding constraints enabling a pay-as-you-go approach, which is largely independent of time-constraints and able to exploit launches of opportunity. In this paper, historical in-space assembly activities which relied heavily on astronaut extra-vehicular activity will be reviewed as well as early robotic assembly activities. These approaches will be contrasted with emerging modular approaches supporting realization of a new Persistent Asset operational paradigm. The paper will define attributes of the Persistent Asset paradigm and illustrate advantages by applying the paradigm to two relevant applications: 1) a large space telescope backing structure and 2) the backbone structure for a solar electric transport vehicle. Finally, recently developed unique Persistent Asset elements (modules, and interface approaches) will be described.

William R Doggett↗

Construction of two-dimensional topological field theories with non-invertible symmetries

We construct the defining data of two-dimensional topological field theories (TFTs) enriched by non-invertible symmetries/topological defect lines. Simple formulae for the three-point functions and the lasso two-point functions are derived, and crossing symmetry is proven. The key ingredients are open-to-closed maps and a boundary crossing relation, by which we show that a diagonal basis exists in the defect Hilbert spaces. We then introduce regular TFTs, provide their explicit constructions for the Fibonacci, Ising and Haagerup $\mathscr{H}$ 3 fusion categories, and match our formulae with previous bootstrap results. We end by explaining how non-regular TFTs are obtained from regular TFTs via generalized gauging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy and momentum conservation in the Euler–Poincaré formulation of local Vlasov–Maxwell-type systems

The action principle by Low for the classic Vlasov–Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler–Poincaré formulation of Vlasov–Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while Cendra et al, it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler–Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov–Darwin system Sugama et al. The present exposition discusses a generic class of local Vlasov–Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler–Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. Here, we also discuss what happens if no symmetries exist. Finally, two explicit examples will be given—the classic Vlasov–Maxwell and the drift-kinetic Vlasov–Maxwell—and the results expressed in the language of regular vector calculus for familiarity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MOOSE Framework Meshing Enhancements to Support Reactor Analysis

MOOSE-based physics codes require an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The use of difficult-touse, external licensed software is often required to generate high quality meshes for reactor geometries. High-fidelity geometry modeling also requires elaborate tracking of groups of elements for material property assignment and output reporting which can be considerably complex for the user to identify and maintain. Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, several meshingrelated enhancements have been developed for the MOOSE framework to address user challenges in creating finite element meshes for advanced reactor geometries. MOOSE mesh generators have been developed to mesh hexagonal geometries (pins, ducted assemblies, and cores) commonly found in liquid-metal cooled fast reactor concepts. The mesh generator used for hexagonal pin cells is generic for regular polygons and therefore may also be used for Cartesian pin cells. Hexagonal pin cells can be stitched into ducted assemblies, and assemblies can be stitched together into a core. The user may specify region ids, region names, and other preferences on the mesh. This control is useful for later material mapping in the MOOSE-based physics codes input. A capability was also developed for meshing rotating control drums including determination of material volume fractions in each mesh element as a function of time. Control drum meshes may be stitched to other hexagonal assemblies to create a core configuration. Additional mesh generators were developed that wrap around the hexagonal meshing capabilities and utilize “extra element integer” ID values on each element. In regular Cartesian or hexagonal assemblies or cores, the bookkeeping of element groups for both material assignment and output reporting can now be automated through assignment of pin, assembly, core, axial and depletion id values stored as extra element integers. The extra element tags on the mesh greatly speed the reactor analyst’s efforts to map materials to meshes, track depletion zones, and parse output such as axial pin power distributions. At the highest level, pin, assembly, and core mesh generators (with this reactor terminology) have also been developed to easily generate regular Cartesian and hexagonal cores, including axial extrusion. These reactor geometry builders call upon the previously mentioned capabilities to produce analysis-ready 3D meshes including material assignments. Open source mesh triangulation capabilities were also investigated for integration into the MOOSE framework to address the need for meshing the core periphery region which extends from the irregular outer assembly border to a cylindrical boundary. Options are limited due to licensing constraints, and the recommendation is pursue building a native MOOSE Delaunay triangulator routine with full functionality. Finally, a series of verification problems were performed with NEAMS physics tools. All developed capabilities will be available in the new open-source “Reactor” module of the MOOSE framework, which is accessible to any MOOSE-based NEAMS physics tool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Pion valence quark distribution at physical pion mass of N f = 2 + 1 + 1 lattice QCD

Abstract We present a state-of-the-art calculation of the unpolarized pion valence-quark distribution in the framework of large-momentum effective theory (LaMET) with improved handling of systematic errors as well as two-loop perturbative matching. We use lattice ensembles generated by the MILC collaboration at lattice spacinga≈ 0.09 fm, lattice volume 64 3 × 96,N f = 2 + 1 + 1 flavors of highly-improved staggered quarks and a physical pion mass. The LaMET matrix elements are calculated with pions boosted to momentumP z ≈ 1.72 GeV with high-statistics ofO(10 6 ) measurements. We study the pion PDF in both hybrid-ratio and hybrid-regularization-independent momentum subtraction (hybrid-RI/MOM) schemes and also compare the systematic errors with and without the addition of leading-renormalon resummation (LRR) and renormalization-group resummation (RGR) in both the renormalization and lightcone matching. The final lightcone PDF results are presented in the modified minimal-subtraction scheme at renormalization scaleμ= 2.0 GeV. We show that thex-dependent PDFs are compatible between the hybrid-ratio and hybrid-RI/MOM renormalization with the same improvements. We also show that systematics are greatly reduced by the simultaneous inclusion of RGR and LRR and that these methods are necessary if improved precision is to be reached with higher-order terms in renormalization and matching.

Physics↗

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↗

Optical studies of shock generated transient supersonic base flows

A shock tube employing interferometric and schlieren techniques is used to study transient base flow phenomena following shock wave passage over two plane bluff bodies: a hemicircular cylinder and a cylinder with the Galileo Jovian probe profile. An attempt is made to understand the physics of transition from transient to steady state flow, and to provide code verification for a study employing the Illiac IV computer. Transient base flow interactions include a series of shock diffraction, regular, and Mach reflections, coupled with boundary layer development, separation, and recompression. Vorticity generation and transport underlie these features. The quantitative verification of the computer code includes comparisons of transient pressure and density fields, near wake geometries, and bow shock standoff distances.

Liang, P.-Y.↗

Gigatraj: An Atmospheric Trajectory Model

Atmospheric trajectory models have a long history of success in tracking air motions in the lower stratosphere and upper troposphere over periods of up to a few days. Parcels have been traced backwards from observations to identify whatever phenomena (strong convection, volcanic eruptions, rocket launches, etc.) put their signature on them. Parcels have also been initialized at a known event and traced forward to examine their subsequent physical and chemical evolution. We describe a new trajectory model, "gigatraj," that aims to increase exibility by (a) making it straightforward to use new meteorological data sources, including those not based on regular latitude-longitude grids; (b) enabling a run-time choice of vertical coordinate system for kinematic and/or quasi-isentropic calculations; (c) allowing for the output of arbitrary meteorological products, selectable by the user and interpolated to the parcels' locations and times. The model can be run in a serial or parallel processing environment, so that large numbers of parcels can be traced in a reasonable time. Information is presented on model accuracy and performance. The former is demonstrated by runs using both test-pattern winds (comparing expected paths with actual output) and real-world winds (comparing forward and backward runs to characterize how well parcels retrace their paths). Sample cases are also shown, including a reverse domain lling (RDF) calculation illustrating a tropopause fold event. Model output can be displayed using the new Visualization And Lagrangian dynamics Immersive eXtended Reality (VALIXR) system, and an example will be shown. In addition, we describe work to incorporate a version of gigatraj into the Goddard Earth Observing System (GEOS) of NASA's Global Modeling and Assimilation O ce (GMAO) at Goddard Space Flight Center. This enables trajectory calculations to be performed within the running GEOS model at the latter's native time resolution, instead of using the winds from every few hours. It also provides access to all of GEOS's internal variables as they are calculated. This module may be useful, for example, for tracking rapid chemical changes in a Lagrangian framework.

dynamics↗

Multiscale Neural Networks for Approximating Green’s Functions

Neural networks (NNs) have been widely used to solve partial differential equations (PDEs) in the applications of physics, biology, and engineering. One effective approach for solving PDEs with a fixed differential operator is learning Green’s functions. However, Green’s functions are notoriously difficult to learn due to their poor regularity, which typically requires larger NNs and longer training times. In this work, we address these challenges by leveraging multiscale NNs to learn Green’s functions. Through theoretical analysis using multiscale Barron space methods and experimental validation, we show that the multiscale approach significantly reduces the necessary NN size and accelerates training.

97 MATHEMATICS AND COMPUTING↗

High resolution spectrograph for the Space Telescope

The high resolution spectrograph (HRS) for ultraviolet astronomy with the Space Telescope will provide a spectral resolution of approximately 120,000 over a nominal wavelength range of 110-320 nm, together with a spatial resolution of about 0.25 arc seconds. The two detectors will consist of 512-element Digicons with cesium telluride and cesium iodide photocathodes, respectively. Photoelectrons in transit between the photocathodes and the diodes within the Digicons can be deflected in two axes with 12-bit resolution. This feature facilitates a design that emphasizes reliability since (once a hermetic seal is opened in orbit), only two moving parts, a grating carrousel and a shutter, are required for regular operation of the HRS. The instrument will be controlled by a computer in the spacecraft. The scientific objectives of the HRS investigation relate to interstellar matter in our own and nearby galaxies, physical processes of stellar mass loss and mass transfer, chemical abundances, bright quasars and Seyfert galaxy nuclei, and solar system phenomena.

Brandt, J. C.↗

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗

Effects of regular exercise training on skeletal muscle contractile function

Skeletal muscle function is critical to movement and one's ability to perform daily tasks, such as eating and walking. One objective of this article is to review the contractile properties of fast and slow skeletal muscle and single fibers, with particular emphasis on the cellular events that control or rate limit the important mechanical properties. Another important goal of this article is to present the current understanding of how the contractile properties of limb skeletal muscle adapt to programs of regular exercise.

NASA Discipline Musculoskeletal↗

In-Situ Measurement of Hall Thruster Erosion Using a Fiber Optic Regression Probe

One potential life-limiting mechanism in a Hall thruster is the erosion of the ceramic material comprising the discharge channel. This is especially true for missions that require long thrusting periods and can be problematic for lifetime qualification, especially when attempting to qualify a thruster by analysis rather than a test lasting the full duration of the mission. In addition to lifetime, several analytical and numerical models include electrode erosion as a mechanism contributing to enhanced transport properties. However, there is still a great deal of dispute over the importance of erosion to transport in Hall thrusters. The capability to perform an in-situ measurement of discharge channel erosion is useful in addressing both the lifetime and transport concerns. An in-situ measurement would allow for real-time data regarding the erosion rates at different operating points, providing a quick method for empirically anchoring any analysis geared towards lifetime qualification. Erosion rate data over a thruster s operating envelope would also be useful in the modeling of the detailed physics inside the discharge chamber. There are many different sensors and techniques that have been employed to quantify discharge channel erosion in Hall thrusters. Snapshots of the wear pattern can be obtained at regular shutdown intervals using laser profilometry. Many non-intrusive techniques of varying complexity and sensitivity have been employed to detect the time-varying presence of erosion products in the thruster plume. These include the use quartz crystal microbalances, emission spectroscopy, laser induced flourescence, and cavity ring-down spectroscopy. While these techniques can provide a very accurate picture of the level of eroded material in the thruster plume, it is more difficult to use them to determine the location from which the material was eroded. Furthermore, none of the methods cited provide a true in-situ measure of erosion at the channel surface while the thruster is in operation (i.e. none yield a continuous channel erosion measurement). A recent fundamental sensor development effort has led to a novel regression, erosion, and ablation sensor technology (REAST). The REAST sensor allows for measurement of real-time surface erosion rates at a discrete surface location. The sensor was tested using a linear Hall thruster geometry (see Fig. 1), which served as a means of producing plasma erosion of a ceramic discharge chamber. The mass flow rate, discharge voltage, and applied magnetic field strength could be varied, allowing for erosion measurements over a broad thruster operating envelope. Results are presented demonstrating the ability of the REAST sensor to capture not only the insulator erosion rates but also changes in these rates as a function of the discharge parameters.

Polzin, Kurt↗

Partial muon capture rates in A = 3 and A = 6 nuclei with chiral effective field theory

Searches for neutrinoless-double beta decay rates are crucial in addressing questions within fundamental symmetries and neutrino physics. The rates of these decays depend not only on unknown parameters associated with neutrinos, but also on nuclear properties. In order to reliably extract information about the neutrino, one needs an accurate treatment of the complex many-body dynamics of the nucleus. Neutrinoless-double beta decays take place at momentum transfers on the order of 100 MeV/$c$ and require both nuclear electroweak vector and axial current matrix elements. Muon capture, a process in the same momentum transfer regime, has readily available experimental data to validate these currents. In this work, we present results of {\it ab initio} calculations of partial muon capture rates for $^3$He and $^6$Li nuclei using variational and Green's Function Monte Carlo computational methods. We estimate the impact of the three-nucleon interactions, the cutoffs used to regularize two-nucleon ($2N$) interactions, and the energy range of $2N$ scattering data used to fit these interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗