Search NASA⌕ Search

SEARCH · Search NASA

Results for “Continuity equation”

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 109 records · Page 6

Exploring the unknown Lambda-neutron interaction

No published \Lambda Λ n scattering data exist. A relativistic heavy-ion experiment has suggested that a \Lambda Λ nn bound state was seen. However, several theoretical analyses have cast serious doubt on the bound-state assertion. Nevertheless, there could exist a three-body \Lambda Λ nn resonance. Such a resonance could be used to constrain the \Lambda Λ n interaction. We discuss \Lambda Λ nn calculations using nn and \Lambda Λ n pairwise interactions of rank-one, separable form that fit effective range parameters of the nn system and those hypothesized for the as yet unobserved \Lambda Λ n system based upon four different \Lambda Λ N potentials. The use of rank-one separable potentials allows one to analytically continue the \Lambda Λ nn Faddeev equations onto the second complex energy plane in search of resonance poles, by examining the eigenvalue spectrum of the kernel of the Faddeev equations. Although each of the potential models predicts a \Lambda Λ nn sub-threshold resonance pole, scaling of the \Lambda Λ n interaction by as little as \sim ∼ 5% does produce a physical resonance. This suggests that one may use photo-(electro-)production of the \Lambda Λ nn system from tritium as a tool to examine the strength of the \Lambda Λ n interaction.

Gibson, Benjamin↗

Neural lumped parameter differential equations with application in friction-stir processing

Lumped parameter methods aim to simplify the evolution of spatially-extended or continuous physical systems to that of a “lumped” element representative of the physical scales of the modeled system. For systems where the definition of a lumped element or its associated physics may be unknown, modeling tasks may be restricted to full-fidelity physics simulations. Here, in this work, we consider data-driven modeling tasks with limited point-wise measurements of otherwise continuous systems. We build upon the notion of the Universal Differential Equation (UDE) to construct data-driven models for reducing dynamics to that of a lumped parameter and inferring its properties. The flexibility of UDEs allow for composing various known physical priors suitable for application-specific modeling tasks, including lumped parameter methods. The motivating example for this work is the plunge and dwell stages for friction-stir welding; specifically, (i) mapping power input into the tool to a point-measurement of temperature and (ii) using this learned mapping for process control.

97 MATHEMATICS AND COMPUTING↗

Optimization-based, property-preserving finite element methods for scalar advection equations and their connection to Algebraic Flux Correction

In this paper, we continue our efforts to exploit optimization and control ideas as a common foundation for the development of property-preserving numerical methods. Here we focus on a class of scalar advection equations whose solutions have fixed mass in a given Eulerian region and constant bounds in any Lagrangian volume. Our approach separates discretization of the equations from the preservation of their solution properties by treating the latter as optimization constraints. This relieves the discretization process from having to comply with additional restrictions and makes stability and accuracy the sole considerations in its design. A property-preserving solution is then sought as a state that minimizes the distance to an optimally accurate but not property-preserving target solution computed by the scheme, subject to constraints enforcing discrete proxies of the desired properties. Furthermore, we consider two such formulations in which the optimization variables are given by the nodal solution values and suitably defined nodal fluxes, respectively. A key result of the paper reveals that a standard Algebraic Flux Correction (AFC) scheme is a modified version of the second formulation obtained by shrinking its feasible set to a hypercube. In conclusion, we present numerical studies illustrating the optimization-based formulations and comparing them with AFC

97 MATHEMATICS AND COMPUTING↗

Generalized Grad–Shafranov equation for non-axisymmetric MHD equilibria

The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad–Shafranov equation, which is a single elliptic partial differential equation in two space dimensions. By revealing a hidden symmetry, here we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad–Shafranov equation. In contrast to solutions of the classical Grad–Shafranov equation, solutions of the generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad–Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

What is (quantitative) system dynamics modeling? Defining characteristics and the opportunities they create

A clear definition of system dynamics modeling can provide shared understanding and clarify the impact of the field. We introduce a set of characteristics that define quantitative system dynamics, selected to capture core philosophy, describe theoretical and practical principles, and apply to historical work but be flexible enough to remain relevant as the field progresses. The defining characteristics are: (1) models are based on causal feedback structure, (2) accumulations and delays are foundational, (3) models are equation-based, (4) concept of time is continuous, and (5) analysis focuses on feedback dynamics. We discuss the implications of these principles and use them to identify research opportunities in which the system dynamics field can advance. These research opportunities include causality, disaggregation, data science and AI, and contributing to scientific advancement. Progress in these areas has the potential to improve both the science and practice of system dynamics.

97 MATHEMATICS AND COMPUTING↗

MOOSE Navier–Stokes module

The MOOSE Navier–Stokes module solves mass, momentum, energy, and passive scalar conservation equations in the context of fluid flow. The module supports solution of these equations in both free flow and porous medium contexts and for a range of fluid compressibility. The conservation equations can be discretized in space using continuous Galerkin finite elements or with cell centered finite volumes.

97 MATHEMATICS AND COMPUTING↗

Institutional Computing: Annual Progress Report, Precision early universe simulations to constrain nuclear reactions and beyond standard model physics

The current one-year allocation (w20 qburst) is the first year in a two year allocation of the same name. The project has supported two publications in peer-review journals, another manuscript in preparation for publication, a white paper for the Astro2020 Decadal Survey, and a number of talks (see below). Code development continues on implementing a neutrino quantum kinetic equation solver into BURST, now called QBURST. While implementing coherent and forward-scattering processes into the dynamical neutrino evolution equations, multiple numerical issues have appeared. We have remedied a number of these issues through extensive debugging and testing, which include streamlining interpolation and implementing new ODE integration techniques. Further testing is required to solve the problem at a precision of better than one part in 10 6 . We expect that we will understand the numerical difficulties and implement solutions by end of year 2021.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗

Vibration-Based Sensor Design: A Grey-Box Approach

Knowledge of the internal structure of an object or device under investigation proceeds from the basic idea of constructing its dynamic behavioral relations governed by a set of differential/algebraic equations that characterize its response. These equations can be partial differential equations leading to finite element or finite difference relations requiring a complex numerical solution on a super computer or ordinary differential equations requiring sophisticated numerical integration techniques to obtain the desired solution. Discrete dynamic systems evolving from digitized data acquisition are typically captured by sampled-data (continuous-to-discrete) representations characterized by a set of difference equations specifying the underlying system dynamics. In any case, with a mathematical description in hand, Grey-Box modeling techniques have evolved, concerned with the estimation of model parameters embedded in a prescribed set of equations (the system) governing its behavior, while capturing the underlying physical phenomenology of the problem at hand.

97 MATHEMATICS AND COMPUTING↗

Bifurcation Analysis Reveals Solution Structures of Phase Field Models

The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems. Here, in this study, we develop a new analytical approach based on the bifurcation analysis to explore the mathematical solution structure of phase field models. Revealing such solution structures not only is of great mathematical interest but also may provide guidance to experimentally or computationally uncover new morphological evolution phenomena in materials undergoing electronic and structural phase transitions. To elucidate the idea, we apply this analytical approach to three representative phase field equations: the Allen-Cahn equation, the Cahn-Hilliard equation, and the Allen-Cahn-Ohta-Kawasaki system. The solution structures of these three phase field equations are also verified numerically by the homotopy continuation method.

97 MATHEMATICS AND COMPUTING↗

PyOECP: A flexible open-source software library for estimating and modeling the complex permittivity based on the open-ended coaxial probe (OECP) technique

Here, we present PyOECP, a Python-based flexible open-source software for estimating and modeling the complex permittivity obtained from the open-ended coaxial probe (OECP) technique. The transformation of the measured reflection coefficient to complex permittivity is performed based on three different methods. The software library contains the dielectric spectra of common reference liquids, which can be used to transform the reflection coefficient into the dielectric spectra. Several Python routines that are commonly employed (e.g., SciPy and NumPy) in the field of science and engineering are required only so that the users can alter the software structure depending on their needs. The modeling algorithm exploits the Markov Chain Monte Carlo method for the data regression. The discrete relaxation models can be built by a proper combination of well-known relaxation models. In addition to these models, electrode polarization, a typical measurement artifact for interpreting dielectric spectra, can be incorporated into the modeling algorithm. A continuous relaxation model, which solves the Fredholm integral equation of the first kind (a mathematically ill-posed problem), is also included. This open-source software enables users to freely adjust the physical parameters to obtain physical insight into their materials under test and will be consistently updated for more accurate measurement and interpretation of dielectric spectra in an automated manner. This work describes the theoretical and mathematical background of the software, lays out the workflow, and validates the software functionality based on both synthetic and empirical data included in the software.

97 MATHEMATICS AND COMPUTING↗

A Continental-Scale Estimate of Soil Organic Carbon Change at NEON Sites and Their Environmental and Edaphic Controls

Current carbon cycle models focus on the effects of climate and land-use change on primary productivity and microbial-mineral dependent carbon turnover in the topsoil, while less attention has been paid to vertical soil processes and soil-dependent response to land-use change along the profile. In this study, a spatial-temporal analysis was used to estimate soil organic carbon (SOC) change in topsoil/A horizon and subsoil/B horizon at National Ecological Observatory Network (NEON) sites, USA over 30 years. To separate the effects of land-use, environmental, and edaphic factors on SOC change, space-for-time substitution was used in combination with the Continuous Change Detection and Classification algorithm and Structural Equation Modeling. Results showed that (a) under natural vegetation, Spodosols and Inceptisols found in the eastern NEON sites had substantial topsoil SOC accumulation (+0.4 to +1.2 Mg C ha –1 year –1 ), while Inceptisols and Andisols in the west had a comparable magnitude of topsoil SOC loss (–0.5 to –1.8 Mg C ha –1 year –1 ); (b) Mollisols and Alfisols in the Central Plains sites were susceptible to significant SOC loss under farming and grazing; (c) Runoff/erosion and leaching potential, vertical translocation, and mineral sorption were the most important factors controlling SOC variation across the NEON sites. Our work could be used to parameterize ecosystem models simulating SOC change.

54 ENVIRONMENTAL SCIENCES↗

PyOECP: A flexible open-source software for the measurement and modeling of the dielectric spectrum from the reflection coefficient

We present PyOECP, a Python-based flexible open-source software for estimating and modeling the complex permittivity obtained from the open-ended coaxial probe (OECP) technique. The transformation of the measured reflection coefficient to complex permittivity is performed based on two different models, including the capacitance model and the antenna model. The software library contains dielectric spectra of common reference liquids, which can be used to transform the reflection coefficient into the dielectric spectra. Several Python routines that are commonly employed (e.g., SciPy and NumPy) are only required so that the users can alter the software structure depending on their needs. The modeling algorithm exploits the Markov chain Monte Carlo method for the data regression. The discrete relaxation models can be built by a proper combination of well-known relaxation models. In addition to these models, the electrode polarization, which is a common measurement artifact for interpreting the dielectric spectra, can be incorporated in the modeling algorithm. A continuous relaxation model, which solves the Fredholm integral equation of the first kind (a mathematically ill-posed problem) is also included.

Yoon, Tae-Jun↗

REopt Lite User Manual

REopt Lite evaluates the economic viability of grid-connected solar photovoltaics, wind, combined heat and power, and electric and thermal storage at commercial and small industrial sites. It allows building owners to identify the system sizes and dispatch strategies that minimize the site’s life cycle cost of energy. REopt Lite also estimates the amount of time on-site generation and storage can sustain the site's critical load during a grid outage and allows the user the choice of optimizing for energy resilience. It is primarily used to inform project development decisions and to support research on the factors that drive project feasibility for market development and policy analysis. It is available through a web interface, application programing interface, and open-source code. This user manual provides an overview of the model, including its capabilities and typical applications; inputs and outputs; economic calculations; technology descriptions; and model parameters, variables, and equations. The model is highly flexible and is continually evolving to meet the needs of each analysis. Therefore, this report is not an exhaustive description of all capabilities, but rather a summary of the core components of the model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Exact quantization and analytic continuation

In this paper we give a streamlined derivation of the exact quantization condition (EQC) on the quantum periods of the Schrödinger problem in one dimension with a general polynomial potential, based on Wronskian relations. We further generalize the EQC to potentials with a regular singularity, describing spherical symmetric quantum mechanical systems in a given angular momentum sector. We show that the thermodynamic Bethe ansatz (TBA) equations that govern the quantum periods undergo nontrivial monodromies as the angular momentum is analytically continued between integer values in the complex plane. The TBA equations together with the EQC are checked numerically against Hamiltonian truncation at real angular momenta and couplings, and are used to explore the analytic continuation of the spectrum on the complex angular momentum plane in examples.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Synthesis of Acid Hydrazides from Carboxylic Acids in Continuous Flow

A concise methodology for the synthesis of acid hydrazides from carboxylic acids via a continuous flow process in yields ranging from 65 to 91% is described. Using short residence times of between 13 and 25 mins, the conditions proved amenable to a variety of mono- and diacids, including aliphatic, aromatic, and heteroaromatic functionalities. Furthermore, to demonstrate scalability, a large-scale (200 g) synthesis of azelaic dihydrazide was conducted over a continuous run of 9 h with an 86% overall yield, equating to a 22 g/h output.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Second-Order Invariant Domain Preserving ALE Approximation of Euler Equations

Abstract An invariant domain preserving arbitrary Lagrangian-Eulerian method for solving non-linear hyperbolic systems is developed. The numerical scheme is explicit in time and the approximation in space is done with continuous finite elements. The method is made invariant domain preserving for the Euler equations using convex limiting and is tested on various benchmarks.

Guermond, Jean-Luc↗

On the Cauchy problem for Boltzmann equation modeling a polyatomic gas

In the present article, we consider the Boltzmann equation that models a polyatomic gas by introducing one additional continuous variable, referred to as microscopic internal energy. We establish existence and uniqueness theory in the space homogeneous setting for the full non-linear case, under an extended Grad-type assumption on transition probability rates, which comprises hard potentials for both the relative speed and internal energy with the rate in the interval 0,2, multiplied by an integrable angular part and integrable partition functions. The Cauchy problem is resolved by means of an abstract ordinary differential equation (ODE) theory in Banach spaces for the initial data with finite and strictly positive gas mass and energy, finite momentum, and additionally finite K* polynomial moment, with K* depending on the rate of the transition probability and the structure of a polyatomic molecule or its internal degrees of freedom. Moreover, we prove that polynomially and exponentially weighted Banach space norms associated with the solution are both generated and propagated uniformly in time.

Physics↗