Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First law”

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

High-Order Space-Time Methods for Conservation Laws

Current high-order methods such as discontinuous Galerkin and/or flux reconstruction can provide effective discretization for the spatial derivatives. Together with a time discretization, such methods result in either too small a time step size in the case of an explicit scheme or a very large system in the case of an implicit one. To tackle these problems, two new high-order space-time schemes for conservation laws are introduced: the first is explicit and the second, implicit. The explicit method here, also called the moment scheme, achieves a Courant-Friedrichs-Lewy (CFL) condition of 1 for the case of one-spatial dimension regardless of the degree of the polynomial approximation. (For standard explicit methods, if the spatial approximation is of degree p, then the time step sizes are typically proportional to 1/p(exp 2)). Fourier analyses for the one and two-dimensional cases are carried out. The property of super accuracy (or super convergence) is discussed. The implicit method is a simplified but optimal version of the discontinuous Galerkin scheme applied to time. It reduces to a collocation implicit Runge-Kutta (RK) method for ordinary differential equations (ODE) called Radau IIA. The explicit and implicit schemes are closely related since they employ the same intermediate time levels, and the former can serve as a key building block in an iterative procedure for the latter. A limiting technique for the piecewise linear scheme is also discussed. The technique can suppress oscillations near a discontinuity while preserving accuracy near extrema. Preliminary numerical results are shown

Huynh, H. T.↗

Adaptive Control Law for PID

This paper is the final report for my Spring 2019 internship at the Kennedy Space Center in Cape Canaveral, Florida. The official title of my internship is 'Launch Vehicle Control Study internship,' and I spent the spring working in GMRO- Granular Mechanics and Regolith Operations. There are two components to my project involving developing an adaptive control law for a PID controller. The first component involves performing a deep study and analysis of adaptive control laws, with the goal of developing stability proofs about the control law and the gain and phase margins. The second component involves analyzing the data received after implementing the adaptive control law. Unfortunately, as of the time of writing this paper, the data is not considered ‘clean' enough to analyze. Therefore, in this paper I will give an overview of adaptive control law stability proofs and will write about the data analysis separately. Section II includes several general definitions, and the later sections have additional definitions at the end of each section.A PID (proportional integral derivative) controller is a common control algorithm which isimplemented in NASA’s rocket launch system. The idea behind the PID controller is to calculatethe error (current position – desired position) and drive it to zero by using proportional, integral,and derivative influences on the controller. For example, when a sailor steers a ship which isheading towards location 𝑥 but would like to turn the ship to reach location 𝑥∗, he/she wouldinitially turn the wheel significantly, and as the ship proceeds towards 𝑥∗ the sailor would slowlyshift the wheel back to its original position, and thereby drive the error, 𝑥∗ − 𝑥, to zero.

Eisenberg, Yael↗

A Control Model: Interpretation of Fitts' Law

The analytical results for several models are given: a first order model where it is assumed that the hand velocity can be directly controlled, and a second order model where it is assumed that the hand acceleration can be directly controlled. Two different types of control-laws are investigated. One is linear function of the hand error and error rate; the other is the time-optimal control law. Results show that the first and second order models with the linear control-law produce a movement time (MT) function with the exact form of the Fitts' Law. The control-law interpretation implies that the effect of target width on MT must be a result of the vertical motion which elevates the hand from the starting point and drops it on the target at the target edge. The time optimal control law did not produce a movement-time formula simular to Fitt's Law.

Connelly, E. M.↗

Neural Scaling Laws for Jet Generation

Recently observed empirical scaling laws describe the performance of foundation-type models as three independent key quantities -- dataset size, compute, and model parameters -- are modified. Extracting these scaling laws informs the training of large complex models for which the tuning of hyperparameters in traditional ways is not feasible. This work for the first time explores if scaling laws can also be observed for the task of particle jet generation -- both relevant as a pre-training objective for foundation models and as in-situ simulation by itself. We indeed replicate the key logarithmic scaling law behavior for model-size scaling. Beyond studying the next token prediction validation loss of the generative model, we also study the sliced Wasserstein distance of five physical quantities that are not immediately available to the model during training. Our study shows that this quantity is monotonically related to the next token prediction validation loss, meaning that this loss is indeed a good proxy for the physics performance. For the scaling with dataset size and compute, we observe substantially weaker scaling behavior of both the loss and the sliced Wasserstein distance. We analyze this behavior by introducing the concept of a learnable window, and argue that autoregressive next token prediction on jet constituents exhibits comparatively rapid saturation relative to language-model studies. We discuss possible origins of this behavior, including the stochastic nature of QCD radiation and differences between generative and supervised learning tasks in collider physics.

Amram, Oz [Fermilab]↗

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗

Crater size estimates for large-body terrestrial impact

Calculating the effects of impacts leading to global catastrophes requires knowledge of the impact process at very large size scales. This information cannot be obtained directly but must be inferred from subscale physical simulations, numerical simulations, and scaling laws. Schmidt and Holsapple presented scaling laws based upon laboratory-scale impact experiments performed on a centrifuge (Schmidt, 1980 and Schmidt and Holsapple, 1980). These experiments were used to develop scaling laws which were among the first to include gravity dependence associated with increasing event size. At that time using the results of experiments in dry sand and in water to provide bounds on crater size, they recognized that more precise bounds on large-body impact crater formation could be obtained with additional centrifuge experiments conducted in other geological media. In that previous work, simple power-law formulae were developed to relate final crater diameter to impactor size and velocity. In addition, Schmidt (1980) and Holsapple and Schmidt (1982) recognized that the energy scaling exponent is not a universal constant but depends upon the target media. Recently, Holsapple and Schmidt (1987) includes results for non-porous materials and provides a basis for estimating crater formation kinematics and final crater size. A revised set of scaling relationships for all crater parameters of interest are presented. These include results for various target media and include the kinematics of formation. Particular attention is given to possible limits brought about by very large impactors.

Schmidt, Robert M.↗

Spectral Single and Double Power-law Formation by Sequential Particle Acceleration in Flare Magnetic Islands

Spectral single and double power laws are common in high-energy phenomena, such as solar flares and solar energetic particles, including ground level enhancement events. It is not clear what determines the energy breaks and spectral indexes of these power laws. Here, we will describe a first-principles model of pitch-angle and energy distribution function evolution, which produces power laws and provides physical interpretation for such spectral features (Guidoni et al. 2022, ApJ). In this model, a prescribed fraction of particles sequentially “hops” between shrinking magnetic islands (accelerators) formed by flare reconnection. Each accelerator increases particles’ energies by a modest amount, but particles must visit only a few accelerators to increase their energies by orders of magnitude. Data from global magnetohydrodynamic simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. We will also describe the fully analytic method for forming and interpreting power laws, which requires only a few constrained physical parameters of the acceleration region and is independent of the acceleration model, as well as preliminary results extending the analytical model to the formation of double power laws.

S E Guidoni↗

Analytical guidance law development for aerocapture at Mars

During the first part of this reporting period research has concentrated on performing a detailed evaluation, to zero order, of the guidance algorithm developed in the first period taking the numerical approach developed in the third period. A zero order matched asymptotic expansion (MAE) solution that closely satisfies a set of 6 implicit equations in 6 unknowns to an accuracy of 10(exp -10), was evaluated. Guidance law implementation entails treating the current state as a new initial state and repetitively solving the MAE problem to obtain the feedback controls. A zero order guided solution was evaluated and compared with optimal solution that was obtained by numerical methods. Numerical experience shows that the zero order guided solution is close to optimal solution, and that the zero order MAE outer solution plays a critical role in accounting for the variations in Loh's term near the exit phase of the maneuver. However, the deficiency that remains in several of the critical variables indicates the need for a first order correction. During the second part of this period, methods for computing a first order correction were explored.

Calise, A. J.↗

Safe Tracking Control of an Uncertain Euler-Lagrange System with Full-State Constraints using Barrier Functions

This paper presents a novel, safe tracking control design method that learns the parameters of an uncertain Euler-Lagrange (EL) system online using adaptive learning laws. A barrier function (BF) is first used to transform the full-state constrained EL-dynamics into an equivalent unconstrained dynamics. An adaptive tracking controller is then developed along with the parameter update law in the transformed state space such that the states remain bounded for all time within a prescribed bound. A stability analysis is developed that considers the EL-dynamics’ uncertainty, yielding a semi-globally uniformly ultimately bounded (SGUUB) tracking error and the parameter estimation error. The controller design is validated in simulations using a two-link planar manipulator. The results show the proposed method’s ability to track the reference trajectory while remaining inside each of the predefined state bounds.

Robots↗

Facile Access to Organostibines via Selective Organic Superbase Catalyzed Antimony‐Carbon Protonolysis

The selective formation of antimony-carbon bonds via organic superbase catalysis under metal- and salt-free conditions is reported. Here, this novel approach utilizes electron-deficient stibine, Sb(C 6 F 5 ) 3 , to give upon base-catalyzed reactions with weakly acidic aromatic and heteroaromatic hydrocarbons access to a range of new aromatic and heteroaromatic stibines, respectively, with loss of C 6 HF 5 . Also, the significantly less electron-deficient stibines, Ph 2 SbC 6 F 5 and PhSb(C 6 F 5 ) 2 smoothly underwent base-catalyzed exchange reactions with a range of terminal alkynes to generate the stibines of formulae PhSb(C≡CPh) 2 , and Ph 2 SbC≡CR [R=C 6 H 5 , C 6 H 4 -NO 2 , COOEt, CH 2 Cl, CH 2 NEt 2 , CH 2 OSiMe 3 , Sb(C 6 H 5 ) 2 ], respectively. These formal substitution reactions proceed with high selectivity as only the C 6 F 5 groups serve as a leaving group to be liberated as C 6 HF 5 upon formal proton transfer from the alkyne. Kinetic studies of the base-catalyzed reaction of Ph 2 SbC 6 F 5 with phenyl acetylene to form Ph 2 SbC≡CPh and C 6 HF 5 suggested the empirical rate law to exhibit a first-order dependence with respect to the base catalyst, alkyne and stibine. DFT calculations support a pathway proceeding via a concerted σ-bond metathesis transition state, where the base catalyst activates the Sb-C 6 F 5 bond sequence through secondary bond interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock Sensitivity of energetic materials

Viscoplastic deformation is examined as the principal source of hot energy. Some shock sensitivity data on a proposed model is explained. A hollow sphere model is used to approximate complex porous matrix of energetic materials. Two pieces of shock sensitivity data are qualitatively compared with results of the proposed model. The first is the p2 tau law. The second is the desensitization of energetic materials by a ramp wave applied stress. An approach to improve the model based on experimental observations is outlined.

Kim, K.↗

A high temperature thermal diffusivity determination procedure for solids and liquids

A new method for measuring the thermal diffusivity of materials at high temperatures is presented. The method is applicable to solids on Earth, and to liquids in the reduced gravity environment of space. It is especially suited to levitated liquid metals at elevated temperatures where thermal diffusivity data is not available. The method is applied in two parts, such that lumped analysis is valid in the first part, and Fourier's law of conduction is valid in the second. In both parts, the spherical specimen is assumed to have been heated to a desired temperature and cooled. An inverse conduction problem is then formulated and solved using Laplace transformation techniques. Using this solution, sample sizes, and experimentally obtained surface temperature history, the thermal diffusivity is determined by minimizing a function that satisifies the heat balance at the surface. Minimization is performed using a modified quasilinearization algorithm. Accuracy is very sensitive to error in the temperature data and increases with better curve fits to the temperature data. An error analysis is also performed, and the effect of errors in the various parameters on the evaluated thermal diffusivity is determined. An experimental study for solids on Earth is suggested, before development for implementation in space.

Bayazitoglu, Yildiz↗

Space-Plane Spreadsheet Program

Basic Hypersonic Data and Equations (HYPERDATA) spreadsheet computer program provides data gained from three analyses of performance of space plane. Equations used to perform analyses derived from Newton's second law of physics, derivation included. First analysis is parametric study of some basic factors affecting ability of space plane to reach orbit. Second includes calculation of thickness of spherical fuel tank. Third produces ratio between volume of fuel and total mass for each of various aircraft. HYPERDATA intended for use on Macintosh(R) series computers running Microsoft Excel 3.0.

Mackall, Dale↗

The formation of ore mineral deposits on the Moon: A feasibility study

Most of the ore deposits on Earth are the direct result of formation by hydrothermal solutions. Analogous mineral concentrations do not occur on the Moon, however, because of the absence of water. Stratified ore deposits form in layered instrusives on Earth due to fractional crystallization of magma and crystal settling of high-density minerals, particularly chromium in the mineral chromite. We have evaluated the possibility of such mineral deposition on the Moon, based upon considerations of 'particle settling velocities' in lunar vs. terrestrial magmas. A first approximation of Stoke's Law would seem to indicate that the lower lunar gravity (1/6 terrestrial) would result in slower crystal settling on the Moon. However, the viscosity of the silicate melt is the most important factor affecting the settling velocity. The viscosities of typical lunar basaltic melts are 10-100 times less than their terrestrial analogs. These lower viscosities result from two factors: (1) lunar basaltic melts are typically higher in FeO and lower in Al2O3, Na2O, and K2O than terrestrial melts; and (2) lunar igneous melts and phase equilibria tend to be 100-150 C higher than terrestrial, largely because of the general paucity of water and other volatile phases on the Moon. Therefore, particle settling velocities on the Moon are 5-10 times greater than those on Earth. It is highly probable that stratiform ore deposits similar to those on Earth exist on the Moon. The most likely ore minerals involved are chromite, ilmenite, and native FeNi metal. In addition, the greater settling velocities of periodotite in lunar magmas indicate that the buoyancy effects of the melt are less than on Earth. Consequently, the possibility is considerably less than on Earth of deep-seated volcanism transporting upper mantle/lower crustal xenoliths to the surface of the Moon, such as occurs in kimberlites on Earth.

Taylor, Lawrence A.↗

An Astronomical Test of CCD Photometric Precision

This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques. we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approximations of the exact error representation formula are then proposed and numerically evaluated for Ringleb flow, an exact solution of the 2-D Euler equations.

Koch, David↗

A Posteriori Error Estimation for Discontinuous Galerkin Approximations of Hyperbolic Systems

This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques, we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approximations of the exact error representation formula are then proposed and numerically evaluated for Ringleb flow, an exact solution of the 2-D Euler equations.

Larson, Mats G.↗

Predicting Rocket or Jet Noise in Real Time

A semi-empirical theoretical model and a C++ computer program that implements the model have been developed for use in predicting the noise generated by a rocket or jet engine. The computer program, entitled the Realtime Rocket and Jet Engine Noise Analysis and Prediction Software, is one of two main subsystems of the Acoustic Prediction/Measurement Tool, which comprises software, acoustic instrumentation, and electronic hardware combined to afford integrated capabilities for real-time prediction and measurement of noise emitted by rocket and jet engines. [The other main subsystem, consisting largely of acoustic instrumentation and electronic hardware, is described in Wireless Acoustic Measurement System, which appears elsewhere in this section.] The theoretical model was derived from the fundamental laws of fluid mechanics, as first was done by M. J. Lighthill in his now famous theory of aerodynamically generated sound. The far-field approximation of the Lighthill theory is incorporated into this model. Many other contributions from various researchers have also been introduced into the model. The model accounts for two noise components: shear noise and self noise. The final result of the model is expressed in terms of a volume integral of the acoustic intensities attributable to these two components, subject to various directivity coefficients. The computer program was written to solve the volume integral. The inputs required by the program are two data files from a computational fluid dynamics (CFD) simulation of the flow of interest: the computational-grid file and the solution file. The CFD solution should be one that has been obtained for conditions that closely approximate those of an experimental test that is yet to be performed. In the current state of development of the model and software, it is recommended that the observation points lie along a radius at an angle >60 from the jet axis. The software provides, and is driven via, a graphical user interface, which facilitates its use. Optionally, the program accepts additional input in the form of data on the measured sound pressure level as a function of frequency at a given far-field location, preferably at an angle of 90 from the jet axis. The user is prompted to use default empirical constants or to choose constants based the measurement data. The user can view the results and compare them with other computational or experimental data. Once satisfied with the results, the user can save a graph of the results in a file that can be imported into documents.

Frendi, Kader↗

Low Activity Waste Glass Optimization with Property Models from Machine Learning, Part 2: Experimental Validation and Active Learning

The United States Department of Energy is responsible for managing legacy nuclear waste stored in underground tanks at the Hanford Site. To treat the waste, it is planned as the current baseline to separately vitrify low-activity waste (LAW) and high-level waste fractions. Previously, machine learning (ML) based glass property models (e.g., chemical durability, viscosity, electrical conductivity and SO3 solubility) were developed with prediction uncertainties. A waste glass optimization approach was then established to enable the capability of using these ML models in LAW glass formulation. In this study, the previous ML models were first experimentally validated, and the results were incorporated back into the database to update the ML models. The updated models and formulations showed increased waste loading while reducing the failure rate, demonstrating improved predictive accuracy, reduced uncertainties, and the effectiveness of active learning in guiding high-dimensional, nonlinear LAW glass design. This represents the first experimental validation of ML based LAW glass formulation, with practical benefits such as higher waste loading, shorter mission duration, and lower operational risk.

Lu, Xiaonan (ORCID:0000000179708148)↗