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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Large-Eddy Simulation of Laminar-Turbulent Transition in a Swept-Wing Boundary Layer

A continuous adjoint approach for obtaining sensitivity derivatives on unstructured grids is developed and analyzed. The derivation of the costate equations is presented, and a second-order accurate discretization method is described. The relationship between the continuous formulation and a discrete formulation is explored for inviscid, as well as for viscous flow. Several limitations in a strict adherence to the continuous approach are uncovered, and an approach that circumvents these difficulties is presented. The issue of grid sensitivities, which do not arise naturally in the continuous formulation, is investigated and is observed to be of importance when dealing with geometric singularities. A method is described for modifying inviscid and viscous meshes during the design cycle to accommodate changes in the surface shape. The accuracy of the sensitivity derivatives is established by comparing with finite-difference gradients and several design examples are presented.

Huai, X.↗
NASA NTRS1993

Theoretical calculations on the electron absorption spectra of selected Polycyclic Aromatic Hydrocarbons (PAH) and derivatives

As a theoretical component of the joint effort with the laboratory of Dr. Lou Allamandola to search for potential candidates for interstellar organic carbon compound that are responsible for the visible diffuse interstellar absorption bands (DIB's), quantum mechanical calculations were performed on the electron absorption spectra of selected polycyclic aromatic hydrocarbons (PAH) and derivatives. In the completed project, 15 different species of naphthalene, its hydrogen abstraction and addition derivatives, and corresponding cations and anions were studied. Using semiempirical quantum mechanical method INDO/S, the ground electronic state of each species was evaluated with restricted Hartree-Fock scheme and limited configuration interaction. The lowest energy spin state for each species was used for electron absorption calculations. Results indicate that these calculations are accurate enough to reproduce the spectra of naphthalene cation and anion observed in neon matrix. The spectral pattern of the hydrogen abstraction and addition derivatives predicted based on these results indicate that the electron configuration of the pi orbitals of these species is the dominant determinant. A combined list of 19 absorptions calculated from 4500 A to 10,400 A were compiled and suggested as potential candidates that are relevant for the DIB's absorptions. Continued studies on pyrene and derivatives revealed the ground state symmetries and multiplicities of its neutral, anionic, and cationic species. Spectral calculations show that the cation (B(sub 3g)-2) and the anion (A(sub u)-2) are more likely to have low energy absorptions in the regions between 10 kK and 20 kK, similar to naphthalene. These absorptions, together with those to be determined from the hydrogen abstraction and addition derivatives of pyrene, can be used to provide additional candidates and suggest experimental work in the search for interstellar compounds that are responsible for DIB's.

Du, Ping↗
NASA NTRS1985

Noise correction on LANDSAT images using a spline-like algorithm

Many applications using LANDSAT images face a dilemma: the user needs a certain scene (for example, a flooded region), but that particular image may present interference or noise in form of horizontal stripes. During automatic analysis, this interference or noise may cause false readings of the region of interest. In order to minimize this interference or noise, many solutions are used, for instane, that of using the average (simple or weighted) values of the neighboring vertical points. In the case of high interference (more than one adjacent line lost) the method of averages may not suit the desired purpose. The solution proposed is to use a spline-like algorithm (weighted splines). This type of interpolation is simple to be computer implemented, fast, uses only four points in each interval, and eliminates the necessity of solving a linear equation system. In the normal mode of operation, the first and second derivatives of the solution function are continuous and determined by data points, as in cubic splines. It is possible, however, to impose the values of the first derivatives, in order to account for shapr boundaries, without increasing the computational effort. Some examples using the proposed method are also shown.

Vijaykumar, N. L.↗
DOE OSTI2026

Data for Clumping Index Estimation With 30°-tilted Cameras in Row Crops: Evaluation of Methods and Segment Size Effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

Modeling↗
DOE OSTI2025

Moment method and continued fraction expansion in Floquet operator Krylov space

Recursion methods such as Krylov techniques map complex dynamics to an effective noninteracting problem in one dimension. For example, the operator Krylov space for Floquet dynamics can be mapped to the dynamics of an edge operator of the one-dimensional Floquet inhomogeneous transverse field Ising model (ITFIM), where the latter, after a Jordan-Wigner transformation, is a Floquet model of noninteracting Majorana fermions and the couplings correspond to Krylov angles. We present an application of this showing that a moment method exists where given an autocorrelation function, one can construct the corresponding Krylov angles and from that the corresponding Floquet ITFIM. Consequently, when no solutions for the Krylov angles are obtained, it indicates that the autocorrelation is not generated by unitary dynamics. We highlight this by studying certain special cases: stable m-period dynamics derived using the method of continued fractions, exponentially decaying, and power-law decaying stroboscopic dynamics. Remarkably, our examples of stable m-period dynamics correspond to m-period edge modes for the Floquet ITFIM where, deep in the chain, the couplings correspond to a critical phase. Furthermore our results pave the way to engineer Floquet systems with desired properties of edge modes and also provide examples of persistent edge modes in gapless Floquet systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1990

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Park, K. C.↗
NASA NTRS2002

Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives

In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.

Yan, Jue↗
DOE OSTI2015

Dynamic density functional theory with hydrodynamic interactions: Theoretical development and application in the study of phase separation in gas-liquid systems

Building on recent developments in dynamic density functional theory, we have developed a version of the theory that includes hydrodynamic interactions. This is achieved by combining the continuity and momentum equations eliminating velocity fields, so the resulting model equation contains only terms related to the fluid density and its time and spatial derivatives. The new model satisfies simultaneously continuity and momentum equations under the assumptions of constant dynamic or kinematic viscosity and small velocities and/or density gradients. We present applications of the theory to spinodal decomposition of subcritical temperatures for one-dimensional and three-dimensional density perturbations for both a van der Waals fluid and for a lattice gas model in mean field theory. In the latter case, the theory provides a hydrodynamic extension to the recently studied dynamic mean field theory. We find that the theory correctly describes the transition from diffusive phase separation at short times to hydrodynamic behaviour at long times.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2023

Assessment of the Impact of Discontinuity in Satellite Instruments and Retrievals on Global PM2.5

Exposure to fine particulate matter (PM2.5) is the leading environmental risk factor for mortality globally. Satellite-derived estimates of surface PM2.5 developed from a combination of satellites, simulations, and ground monitor data are relied upon for health impact studies. The ability to develop satellite-derived PM2.5 estimates requires the continued availability of aerosol optical depth (AOD) sources. This work examines the impact of the addition or loss of satellite AOD data sources on global PM2.5 estimation and the impact of continuing the longterm record with AOD from the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi-National Polar orbiting Partnership (S-NPP) satellite after the loss of the MODIS (MODerate resolution Imaging Spectroradiometer) and MISR (Multi-angle Imaging Spectroradiometer) instruments on board the Terra and Aqua satellites. We find that the addition of VIIRS S-NPP AOD products to geophysical PM2.5 estimates from satellites and simulations causes regional differences that correspond to differences in the VIIRS and MODIS Deep Blue AOD algorithms and sampling. Changes in long-term trends and timeseries due to the addition or loss of AOD data sources are generally within their uncertainties. Statistical fusion with ground monitor data partially corrects for changes due to sampling differences when introducing the VIIRS AOD products, but uncertainty remains over desert regions where ground monitor coverage is sparse. This work provides promise for the sustained development of global satellite-derived PM2.5 estimates, despite discontinuities in instruments and retrieval methods.

Remote sensing↗
DOE OSTI2023

Necessary and sufficient conditions for quadratic stabilizability of switched systems on non-uniform time domains

In this paper, we consider the quadratic stabilizability via state feedback for a particular class of switched systems that evolve on a non-uniform time domain by introducing time scales theory. The system considered switches between a continuous-time subsystem with variable lengths and a discrete-time subsystem with variable discrete step sizes. Necessary and sufficient conditions are derived to guarantee the quadratic stability of this class of switched systems via a switching state feedback law based on the existence of a common positive definite matrix satisfying the quadratic stabilizability condition by considering that the two subsystems are unstable. By state feedback, we mean that the switching among subsystems depends on the system states. Current results for this kind of state switching feedback control are derived only for switched systems evolving on a continuous time domain or a discrete time domain with fixed step’s size. These results are not applicable for the particular class of switched systems where there is a mixing between the continuous and discrete dynamics. This motivates the derivation of a new and more general state feedback control law for switched systems in this work. Here, a numerical example illustrating the results is presented.

42 ENGINEERING↗
DOE OSTI2023

Electronic and Thermal Properties of the Cation Substitution-Derived Quaternary Chalcogenide CuInSnSe 4

Quaternary chalcogenides continue to be of interest for a variety of technological applications, with physical properties stemming from their structural complexity and stoichiometric variation. In certain structure types, partial vacancies on specific lattice positions present an opportunity to investigate electrical and thermal properties in light of these lattice defects. In this work, we investigated the structural, thermal, and electronic properties of CuInSnSe 4 , a material that belongs to a relatively unexplored class of quaternary chalcogenides with a defect adamantine crystal structure. First-principles calculations together with experimental measurements revealed a chalcopyrite-like structure with inherent vacancies and characteristic s–p and p–d orbital hybridizations in the electronic structure of the material. Cation disorder and lattice anharmonicity result in very low thermal conductivity with values significantly lower than those for related compositions. In conclusion, this work reveals the fundamental physical properties of a previously uninvestigated quaternary chalcogenide and may aid investigations of similar as well as other quaternary chalcogenide compositions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

TORUS Annual Continuation and Progress Report

The TORUS collaboration derives its name from the research it focuses on, namely the Theory of Reactions for Unstable iSotopes. It is a Topical Collaboration in Nuclear Theory, and funded by the Nuclear Theory Division of the Office of Nuclear Physics in the Office of Science of the Department of Energy. The funding supports one postdoctoral researcher for the years 1 through 4. The collaboration brings together as Principal Investigators a large fraction of the nuclear reaction theorists currently active within the USA. The mission of the TORUS Topical Collaboration is to develop new methods that will advance nuclear reaction theory for unstable isotopes by using three-body techniques to improve directreaction calculations. This multi-institution collaborative effort is directly relevant to three areas of interest: the properties of nuclei far from stability; microscopic studies of nuclear input parameters for astrophysics, and microscopic nuclear reaction theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2021

Improving the CERES SYN Cloud and Flux Products by Identifying GOES-17 Scan Anomalies Using a Convolutional Neural Network

The NASA Clouds and the Earth’s Radiant Energy System (CERES) project relies on top-of-atmosphere (TOA) broadband fluxes derived from geostationary (GEO) satellite imagery to account for the diurnal flux variations between the CERES observation intervals, and thereby produce a synoptic gridded (SYN1deg) product based on continuous temporal observations. Consistent broadband flux derivation depends on accurate radiative property measurements and cloud retrievals, which largely determine the radiance-to-flux conversion process. Therefore, it is important to ensure a high quality of cloud property input in order to maintain a reliable broadband flux record. In Edition 4 of the CERES SYN1deg product, a robust automated image anomaly detection algorithm based on inter-line and inter-pixel differences, spatial variance, and 2-D Fourier analysis has been successful in identifying imagery with linear artifacts, but the line-by-line inspection and cleaning process must still be performed by a human. Therefore, further automation of this quality assurance process is warranted, especially considering the excessive amount of additional cleaning necessitated by the GOES-17 Advance Baseline Imager (ABI) cooling system anomaly. As such, this article highlights advancement of the CERES GEO image artifact cleaning approach based on a convolutional neural network (CNN) for classification of bad scanlines. Once trained, the CNN approach is a computationally inexpensive means to ensure greater consistency in cloud retrievals, and therefore broadband flux derivation, based on GOES-17 measurements.

Benjamin Scarino↗
NASA NTRS1995

Unified approach for incompressible flows

A unified approach for solving incompressible flows has been investigated in this study. The numerical CTVD (Centered Total Variation Diminishing) scheme used in this study was successfully developed by Sanders and Li for compressible flows, especially for the high speed. The CTVD scheme possesses better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that the CTVD scheme can equally well apply to solve incompressible flows. Because of the mathematical difference between the governing equations for incompressible and compressible flows, the scheme can not directly apply to the incompressible flows. However, if one can modify the continuity equation for incompressible flows by introducing pseudo-compressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of the algorithm to incompressible flows thus becomes feasible. In this study, the governing equations for incompressible flows comprise continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the physical and numerical boundary conditions are properly implemented by the characteristic boundary conditions. Accordingly, a CFD code has been developed for this research and is currently under testing. Flow past a circular cylinder was chosen for numerical experiments to determine the accuracy and efficiency of the code. The code has shown some promising results.

Chang, Tyne-Hsien↗
DOE OSTI2013

Preparing the Way for New Policy Regarding Adaptation of US Electricity Infrastructure to Climate Change

The following pages represent the status of policy regarding adaptation of the electric grid to climate change and proposed directions for new policy development. While strides are being made to understand the current climate and to predict hazards it may present to human systems, both the science and the policy remain at present in an analytical state. The policy proposed in this document involves first continued computational modeling of outcomes which will produce a portfolio of options to be considered in light of specific region-related risks. It is proposed that the modeling continue not only until reasonable policy at various levels of jurisdiction can be derived from its outcome but also on a continuing basis so that as improvements in the understanding of the state and trajectory of climate science along with advancements in technology arise, they can be incorporated into an appropriate and evolving policy.

Allen, Melissa R.↗
DOE OSTI2018

Submarine exposure guideline recommendations for carbon dioxide based on the prenatal developmental effects of exposure in rats

Background To protect crewmember health, the U.S. Navy sets exposure limits for more than 200 components of submarine atmospheres. The addition of females to nuclear submarines required a reevaluation of these exposure limits, originally established for all‐male crews. In the case of carbon dioxide (CO 2 ), the only available data suitable for deriving an exposure limit were from a 2010 study sponsored by the British Royal Navy that reported a debatable interpretation casting doubt on whether current U.S. Navy exposure limits served to protect fetal developmental health. Methods About 120 time‐mated female Sprague–Dawley rats (Crl: CD[SD]) were exposed to CO 2 at levels of 1.5%, 2.0%, 2.5%, and 3.0% from gestation days 6 to 20. Dams were euthanized and fetuses were examined. Results Findings with implications for exposure limits for CO 2 during pregnancy were an increased mean litter proportion of early resorptions and a lower mean litter proportion of viable fetuses in the 3.0% CO 2 group. Conclusion The results yield a No Observed Adverse Effect Level (NOAEL) of 2.5% and a Lowest Observed Adverse Effect Level (LOAEL) of 3.0%. The results reasonably allow a point of departure of 2.5% CO 2 for deriving an exposure recommendation. An interspecies uncertainty factor was applied to derive a recommended 90‐day continuous exposure limit (CEL) of 0.8% for CO 2 . As reproductive endpoints that are developmental in nature must be assumed to result from a single exposure at a critical point during gestation, it is further recommended that the 24‐hr emergency exposure limit (EEL) also be 0.8%.

Howard, William R.↗
DOE OSTI2024

Ladder symmetries and Love numbers of Reissner-Nordström black holes

It is well known that asymptotically flat black holes in general relativity have vanishing tidal Love numbers. In the case of Schwarzschild and Kerr black holes, this property has been shown to be a consequence of a hidden structure of ladder symmetries for the perturbations. In this work, we extend the ladder symmetries to non-rotating charged black holes in general relativity. As opposed to previous works in this context, we adopt a more general definition of Love numbers, including quadratic operators that mix gravitational and electromagnetic perturbations in the point-particle effective field theory. We show that the calculation of a subset of those couplings in full general relativity is affected by an ambiguity in the split between source and response, which we resolve through an analytic continuation. As a result, we derive a novel master equation that unifies scalar, electromagnetic and gravitational perturbations around Reissner-Nordström black holes. The equation is hypergeometric and can be obtained from previous formulations via nontrivial field redefinitions, which allow to systematically remove some of the singularities and make the presence of the ladder symmetries more manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗