Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Continuous derivative”

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 55 records · Page 3

Electromagnetic momentum and the energy–momentum tensor in a linear medium with magnetic and dielectric properties

In a continuum setting, the energy–momentum tensor embodies the relations between conservation of energy, conservation of linear momentum, and conservation of angular momentum. The well-defined total energy and the well-defined total momentum in a thermodynamically closed system with complete equations of motion are used to construct the total energy–momentum tensor for a stationary simple linear material with both magnetic and dielectric properties illuminated by a quasimonochromatic pulse of light through a gradient-index antireflection coating. The perplexing issues surrounding the Abraham and Minkowski momentums are bypassed by working entirely with conservation principles, the total energy, and the total momentum. We derive electromagnetic continuity equations and equations of motion for the macroscopic fields based on the material four-divergence of the traceless, symmetric total energy–momentum tensor. We identify contradictions between the macroscopic Maxwell equations and the continuum form of the conservation principles. We resolve the contradictions, which are the actual fundamental issues underlying the Abraham–Minkowski controversy, by constructing a unified version of continuum electrodynamics that is based on establishing consistency between the three-dimensional Maxwell equations for macroscopic fields, the electromagnetic continuity equations, the four-divergence of the total energy–momentum tensor, and a four-dimensional tensor formulation of electrodynamics for macroscopic fields in a simple linear medium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stabilization of Softwood-Derived Pyrolysis Oils for Continuous Bio-oil Hydroprocessing

The use of fast pyrolysis as a potential renewable liquid transportation fuel alternative to crude oil depends on successful catalytic upgrading to produce a refinery-ready product with oxygen content and qualities (i.e. specific functional group or compound content) that is compatible with the product’s proposed insertion point. Catalytic upgrading of bio-oil requires high temperature and pressure, while similar to crude oil hydrotreating, is not as straightforward for the thermally unstable pyrolysis oil. For years, a two-temperature zone, downflow trickle bed reactor was the state-of-the art for continuous operation. However, pressure excursion due to plug formation still occurred, typically at the high temperature transition zone, leading to a process shutdown within 140 h. Recently, a bio-oil pre-treatment process, together with a robust commercial catalyst, was found to be enabling the continuous operation of the two-zone hydroprocessing system. Here, we report the results on pre-treating bio-oil at 413 K and 8.4 MPa of flowing H2 (500 L H2/L bio-oil, 0.5 L bio-oil/L catalyst bed) and the attempts to characterize this oil product to understand the chemistry which enabled the long-term processing of bio-oil.

carbonyls↗

Quantized vortices around wavefront nodes, 2

Quantized vortices can occur around nodal points in wavefunctions. The derivation depends only on the wavefunction being single valued, continuous, and having continuous first derivatives. Since the derivation does not depend upon the dynamical equations, the quantized vortices are expected to occur for many types of waves such as electromagnetic and acoustic. Such vortices have appeared in the calculations of the H + H2 molecular collisions and play a role in the chemical kinetics. In a companion paper, it is shown that quantized vortices occur when optical waves are internally reflected from the face of a prism or particle beams are reflected from potential energy barriers.

Hirschfelder, J. O.↗

Local distributed estimation

Based on partial differential equations of motion the closed form solution for the optimal estimation of a spatially continuous state vector is derived, using a continuously distributed sensor. Local control is shown to be the feedback that minimizes a quadratic performance index of sensor and process disturbances. A detailed example of the control of a string in tension is presented.

Schaechter, D. B.↗

Enhanced Elliptic Grid Generation

An enhanced method of elliptic grid generation has been invented. Whereas prior methods require user input of certain grid parameters, this method provides for these parameters to be determined automatically. "Elliptic grid generation" signifies generation of generalized curvilinear coordinate grids through solution of elliptic partial differential equations (PDEs). Usually, such grids are fitted to bounding bodies and used in numerical solution of other PDEs like those of fluid flow, heat flow, and electromagnetics. Such a grid is smooth and has continuous first and second derivatives (and possibly also continuous higher-order derivatives), grid lines are appropriately stretched or clustered, and grid lines are orthogonal or nearly so over most of the grid domain. The source terms in the grid-generating PDEs (hereafter called "defining" PDEs) make it possible for the grid to satisfy requirements for clustering and orthogonality properties in the vicinity of specific surfaces in three dimensions or in the vicinity of specific lines in two dimensions. The grid parameters in question are decay parameters that appear in the source terms of the inhomogeneous defining PDEs. The decay parameters are characteristic lengths in exponential- decay factors that express how the influences of the boundaries decrease with distance from the boundaries. These terms govern the rates at which distance between adjacent grid lines change with distance from nearby boundaries. Heretofore, users have arbitrarily specified decay parameters. However, the characteristic lengths are coupled with the strengths of the source terms, such that arbitrary specification could lead to conflicts among parameter values. Moreover, the manual insertion of decay parameters is cumbersome for static grids and infeasible for dynamically changing grids. In the present method, manual insertion and user specification of decay parameters are neither required nor allowed. Instead, the decay parameters are determined automatically as part of the solution of the defining PDEs. Depending on the shape of the boundary segments and the physical nature of the problem to be solved on the grid, the solution of the defining PDEs may provide for rates of decay to vary along and among the boundary segments and may lend itself to interpretation in terms of one or more physical quantities associated with the problem.

Kaul, Upender K.↗

A prefiltering version of the Kalman filter with new numerical integration formulas for Riccati equations

A prefiltering version of the Kalman filter is derived for both discrete and continuous measurements. The derivation consists of determining a single discrete measurement that is equivalent to either a time segment of continuous measurements or a set of discrete measurements. This prefiltering version of the Kalman filter easily handles numerical problems associated with rapid transients and ill-conditioned Riccati matrices. Therefore, the derived technique for extrapolating the Riccati matrix from one time to the next constitutes a new set of integration formulas which alleviate ill-conditioning problems associated with continuous Riccati equations. Furthermore, since a time segment of continuous measurements is converted into a single discrete measurement, Potter's square root formulas can be used to update the state estimate and its error covariance matrix. Therefore, if having the state estimate and its error covariance matrix at discrete times is acceptable, the prefilter extends square root filtering with all its advantages, to continuous measurement problems.

Womble, M. E.↗

Cascading Failures as Continuous Phase-Space Transitions

In network systems, a local perturbation can amplify as it propagates, potentially leading to a large-scale cascading failure. We derive a continuous model to advance our understanding of cascading failures in power-grid networks. The model accounts for both the failure of transmission lines and the desynchronization of power generators and incorporates the transient dynamics between successive steps of the cascade. In this framework, we show that a cascade event is a phase-space transition from an equilibrium state with high energy to an equilibrium state with lower energy, which can be suitably described in a closed form using a global Hamiltonian-like function. From this function, we show that a perturbed system cannot always reach the equilibrium state predicted by quasi-steady-state cascade models, which would correspond to a reduced number of failures, and may instead undergo a larger cascade. We also show that, in the presence of two or more perturbations, the outcome depends strongly on the order and timing of the individual perturbations. These results offer new insights into the current understanding of cascading dynamics, with potential implications for control interventions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Pressure updating methods for the steady-state fluid equations

We consider the steady state equations for a compressible fluid. Since we wish to solve for a range of speeds we must consider the equations in conservation form. For transonic speeds these equations are of mixed type. Hence, the usual approach is to add time derivatives to the steady state equations and then march these equations in time. One then adds a time derivative of the density to the continuity equation, a derivative of the momentum to the momentum equation and a derivative of the total energy to the energy equation. This choice is dictated by the time consistent equations. However, since we are only interested in the steady state this is not necessary. Thus we shall consider the possibility of adding a time derivative of the pressure to the continuity equation and similar modifications for the energy equation. This can then be generalized to adding combinations of time derivatives to each equation since these vanish in the steady state. When using acceleration techniques such as residual smoothing, multigrid, etc. these are applied to the pressure rather than the density. Hence, the code duplicates the behavior of the incompressible equations for low speeds.

Fiterman, A.↗

Geodetic Strain Analysis Tool

A geodetic software analysis tool enables the user to analyze 2D crustal strain from geodetic ground motion, and create models of crustal deformation using a graphical interface. Users can use any geodetic measurements of ground motion and derive the 2D crustal strain interactively. This software also provides a forward-modeling tool that calculates a geodetic velocity and strain field for a given fault model, and lets the user compare the modeled strain field with the strain field obtained from the user s data. Users may change parameters on-the-fly and obtain a real-time recalculation of the resulting strain field. Four data products are computed: maximum shear, dilatation, shear angle, and principal components. The current view and data dependencies are processed first. The remaining data products and views are then computed in a round-robin fashion to anticipate view changes. When an analysis or display parameter is changed, the affected data products and views are invalidated and progressively re-displayed as available. This software is designed to facilitate the derivation of the strain fields from the GPS and strain meter data that sample it to facilitate the understanding of the strengths and weaknesses of the strain field derivation from continuous GPS (CGPS) and other geodetic data from a variety of tectonic settings, to converge on the "best practices" strain derivation strategy for the Solid Earth Science ESDR System (SESES) project given the CGPS station distribution in the western U.S., and to provide SESES users with a scientific and educational tool to explore the strain field on their own with user-defined parameters.

Kedar, Sharon↗

The use of precession modulation for nutation control in spin-stabilized spacecraft

The relations which determine the nutation effects induced in a spinning spacecraft by periodic precession thrust pulses are derived analytically. By utilizing the idea that nutation need only be observed just before each precession thrust pulse, a difficult continuous-time derivation is replaced by a simple discrete-time derivation using z-transforms. The analytic results obtained are used to develop two types of modulated precession control laws which use the precession maneuver to concurrently control nutation. Results are illustrated by digital simulation of an actual spacecraft configuration.

Taylor, J. M.↗

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

International Research Results and Accomplishments From the International Space Station - A New Compilation

In 2016, the International Space Station (ISS) partnership published the first-ever compilation of international ISS research publications resulting from research performed on the ISS through 2011 (Expeditions 0 through 30). International Space Station Research Accomplishments: An Analysis of Results. From 2000-2011 is a collection of over 1,200 journal publications that describe ISS research in the areas of biology and biotechnology; Earth and space science; educational activities and outreach; human research; physical sciences; technology development and demonstration; and, results from ISS operations. This paper will summarize the ISS results publications obtained through 2011 on behalf of the ISS Program Science Forum that is made up of senior science representatives across the international partnership. NASA's ISS Program Science office maintains an online experiment database (www.nasa.gov/iss- science) that tracks and communicates ISS research activities across the entire ISS partnership, and it is continuously updated by cooperation and linking with the results tracking activities of each partner. It captures ISS experiment summaries and results and includes citations to the journals, conference proceedings, and patents as they become available. This content is obtained through extensive and regular journal and patent database searches, and input provided by the ISS international partners ISS scientists themselves. The International Space Station Research Accomplishments: An Analysis of Results From 2000-2011 is a testament to the research that was underway even as the ISS laboratory was being built. It rejects the scientific knowledge gained from ISS research, and how it impact the fields of science in both space and traditional science disciplines on Earth. Now, during a time when utilization is at its busiest, and with extension of the ISS through at least 2024, the ISS partners work together to track the accomplishments and the new knowledge gained in a way that will impact humanity like no laboratory on Earth. Examples of the highest pro le publications to date from each discipline will also be presented. As ISS research activities and operations continue, scientific data derived from earlier experiments will continuously be re-examined, refined, and assembled with new data and findings, including data from other fields never considered. New results will be produced, allowing breakthroughs in new areas of research and innovative solutions to problems on Earth. The ISS Program Science Forum will continue to capture and report on these results in the form of journal publications, conference proceedings, and patents. We anticipate that successful ISS research will continue to contribute to the science literature in a way that helps to formulate new hypotheses and conclusions that will enable science advancements across a wide range of scientific disciplines both in space and on Earth.

Ruttley, Tara↗

Equations of state, phase relations, and oxygen fugacity of the Ru-RuO 2 buffer at high pressures and temperatures

Experimental studies and measurements of inclusions in diamonds show that ferric iron components are increasingly stabilized with depth in the mantle. To determine the thermodynamic stability of such components, their concentration needs to be measured at known oxygen fugacities. The metal-oxide pair Ru and RuO 2 are ideal as an internal oxygen fugacity buffer in high-pressure experiments. Both phases remain solid to high temperatures and react minimally with silicates, only exchanging oxygen. To calculate oxygen fugacities at high pressure and temperature, however, requires information on the phase relations and equation of state properties of the solid phases. We have made in situ synchrotron X-ray diffraction measurements in a multi-anvil press on mixtures of Ru and RuO 2 to 19.4 GPa and 1473 K with which we have determined phase relations of the RuO 2 phases and derived thermal equations of state (EoS) parameters for both Ru and RuO 2 . Rutile-structured RuO 2 was found to undergo two phase transformations, first at ~7 GPa to an orthorhombic structure and then above 12 GPa to a cubic structure. The phase boundary of the cubic phase was constrained for the first time at high pressure and temperature. We have derived a continuous Gibbs free energy expression for the tetragonal and orthorhombic phases of RuO 2 by fitting the second-order phase transition boundary and P-V-T data for both phases, using a model based on Landau theory. The transition between the orthorhombic and cubic phases was then used along with EoS terms derived for both phases to determine a Gibbs free energy expression for the cubic phase. We have used these data to calculate the oxygen fugacity of the Ru + O 2 = RuO 2 equilibrium, which we have parameterized as a single polynomial across the stability fields of all three phases of RuO 2 . The expression is log 10 f O 2 (Ru – RuO 2 ) = (7.782 – 0.00996P + 0.001932P 2 – 3.76 × 10 –5 P 3 ) + (–13 763 + 592P – 3.955P 2 )/T + (–1.05 × 10 6 – 4622P)/T 2 , which should be valid from room pressure up to 25 GPa and 773–2500 K, with an estimated uncertainty of 0.2 log units. Our calculated fO 2 is shown to be up to 1 log unit lower than estimates that use previous expressions or ignore EoS terms.

Armstrong, Katherine↗