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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Root-driven weathering impacts on mineral-organic associations in deep soils over pedogenic time scales

Plant roots are primary weathering agents in deep soils, yet the impact of resulting mineral transformations on the vast deep soil carbon (C) reservoir are largely unknown. Root-driven weathering of primary minerals causes the formation of reactive secondary minerals, which may protect mineral-organic associations (MOAs) for centuries or millennia. Conversely, root-driven weathering may also transform secondary minerals, potentially enhancing the bioavailability of C previously protected in MOAs. Here we examined the impact of root-driven weathering on MOAs and their capacity to store C over pedogenic time scales. To accomplish this, we compared soil that experienced root-driven weathering, resulting in the formation of discrete rhizosphere zones in deep soil horizons (100-160 cm) of the Santa Cruz Marine Terrace chronosequence (65ka-226ka), with adjacent soil that experienced no root growth. Using a combination of radiocarbon, mass spectrometry, Mössbauer spectroscopy, and X-ray spectromicroscopy approaches, we characterized transformations of mineral-organic associations in relation to changes in C content, turnover and chemistry across four soils ranging in age from 65 to 226 ka. We found that the onset of root-driven weathering (65-90ka) increased the amount of C associated with poorly crystalline iron (Fe) and aluminum (Al) phases, particularly highly-disordered np-goethite. This increase coincided with greater C concentrations, longer C residence times, and a greater abundance of microbially-derived C. Continued root-driven weathering (137-226ka) did not significantly change the amount of C associated with crystalline Fe and Al phases, but resulted in a decline in the amount of C associated with poorly crystalline metal phases. This decline coincided with a decrease in C concentrations and potential turnover rates, and a shift toward plant-derived C. In contrast, soil not affected by root-driven weathering showed comparatively low amounts of C bound in poorly crystalline Fe and Al phases regardless of soil age and, correspondingly, lower C concentrations and estimated residence times. Furthermore, our results demonstrate that root driven-formation and disruption of poorly crystalline Fe and Al phases is a direct control on both C accrual and loss in deep soil. This finding suggests that root impacts on soil C storage are dependent on soil weathering stage, a consideration that is critical for future predictions of the vulnerability of deep soil C to global change.

58 GEOSCIENCES↗
DOE OSTI2025

Data from: 'Abiotic influences on continuous conifer forest structure across a subalpine watershed'

This package archives the core data used for analysis and inference in 'Abiotic influences on continuous conifer forest structure across a subalpine watershed' (Worsham et al., 2025). All data were collected in the East River, Washington Gulch, Slate River, and Coal Creek watersheds of Colorado. In the paper, we quantified the relative influence of climate, topographic, edaphic, and geologic factors on conifer stand structure and composition, and their functional relationships, at the watershed scale. We used waveform LiDAR data to derive spatially continuous stand structure metrics. We fused these with a species-level classification map to estimate tree species abundance. We applied generalized additive and generalized boosted models to evaluate the covariability of structural and compositional metrics with abiotic variables. The package contains the essential products required for reproducing our analysis and the tables and figures reported in the publication. The products comprise four classes: (1) geospatial data, (2) tabular data used for inferential analysis, (3) tabular data describing analytical results and performance statistics, and (4) a data user guide. (1) includes discretized waveform LiDAR data, locations and attributes of individual tree crowns, sampling locations and domain boundaries, a canopy height model, and raster files of estimated forest structural and compositional metrics at 100 m grid scale. (2) includes all response and explanatory variable values applied in inferential models. Response variables include conifer forest stand density, basal area, 95th percentile height, quadratic mean diameter, and others. Explanatory variables include climatic water deficit, actual evapotranspiration, elevation, heat load, soil available water content, and others. (3) includes results of training and testing several individual tree detection (ITD) algorithms, as well as inferential modeling results. (4) is a PDF user guide for this data package, including detailed descriptions and data dictionaries for all files. The data package root contains 17 assets: 8 compressed tape archive (.tar.gz) files, 5 comma-separated values (.csv) files, 3 Geographic Tagged Image File Format (GeoTIFF) (.tif) files, and 1 Portable Document Format (.pdf) file. The compressed .tar.gz archives contain ESRI shapefiles (.shp) .tif, compressed LASer (.laz), and .csv files. The archives must first be decompressed using the widely distributed command-line software utility TAR. All other files, including constituent files within the .tar.gz archives, can be opened in the open-source R statistical computing environment. Alternatively, .csv files may also be read in any simple text editor software or Microsoft Excel. Geospatial files including .shp and .tif files can also be opened in GIS software, such as QGIS (open-source) or ESRI ArcGIS (proprietary). The .pdf Data User Guide can be read with Adobe Acrobat Reader or other compatible readers.

2018 NEON and 2025 CHESS Campaigns↗
NASA NTRS2018

Modified Geometric Truncation of the Scattering Phase Function

Phase function of light scattering on large atmospheric particles has very strong peak in forward direction constituting a challenge for accurate numerical calculations of radiance required in remote sensing problems. Scaling transformation replaces original phase function with a sum of the delta function and a new regular smooth phase function. Geometric truncation is one of the ways to construct such a smooth function. The replacement phase function coincides with the original one outside the forward cone and preserves the asymmetry parameter. It has discontinuity at the cone. Another simple functional form of the replacement phase function within the cone is suggested. It enables continuity and allows for a number of modifications. Three of them are considered in this study: preserving asymmetry parameter, providing continuity of the 1st derivative of the phase function, and preserving mean scattering angle. Yet another problem addressed in this study is objective selection of the width of the forward cone. That angle affects truncation fraction and values of the phase function within the cone. A heuristic approach providing unambiguous criterion of selection of the truncation angle is proposed. The approach has easy numerical implementation. Suggested modifications were tested on cloud phase function using discrete ordinates and Monte Carlo methods. It was shown that the modifications provide better accuracy of the radiance computation compare to the original geometric truncation with discrete ordinates while continuous derivative approach provides significant gain in computer time with Monte Carlo simulations.

Alexander Radkevich↗
DOE OSTI2018

Nilpotent symmetries of a 4D model of Hodge theory: Augmented (anti-)chiral superfield formalism

We derive the continuous nilpotent symmetries of the four (3+1)-dimensional (4D) model of the Hodge theory (i.e. 4D Abelian 2-form gauge theory) by exploiting the beauty and strength of the symmetry invariant restrictions on the (anti-)chiral superfields. The above off-shell nilpotent symmetries are the Becchi–Rouet–Stora–Tyutin (BRST), anti-BRST and (anti-)co-BRST transformations which turn up beautifully due to the (anti-)BRST and (anti-)co-BRST invariant restrictions on the (anti-)chiral superfields that are defined on the (4, 1)-dimensional (anti-)chiral super-submanifolds of the general (4, 2)-dimensional supermanifold on which our ordinary 4D theory is generalized. The latter supermanifold is characterized by the superspace coordinates Z{sup M}=(x{sup μ}, θ, θ-bar ) where x{sup μ} (μ=0,1,2,3) are the bosonic coordinates and a pair of Grassmannian variables θ and θ-bar are fermionic in nature as they obey the standard relationships: θ{sup 2}= θ-bar {sup 2}=0, θ θ-bar + θ-bar θ=0). The derivation of the proper (anti-)co-BRST symmetries and proof of the absolute anticommutativity property of the conserved (anti-)BRST and (anti-) co-BRST charges are novel results of our present investigation (where only the (anti-)chiral superfields and their super-expansions have been taken into account).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Lithium and boron isotopic compositions of olivine in chondrules from carbonaceous and ordinary chondrite meteorites: Implications for the origin of solar 11 B/ 10 B ratio

We report the origin of solar 11 B/ 10 B ratio ~4 remains an open question. It has been thought that a significant portion of boron in the Solar System was derived from continuous spallation nucleosynthesis during interactions between Galactic Cosmic Rays and C—N—O nuclei in the interstellar medium. However, because GCR-produced boron is characterized by 11 B/ 10 B ~2.5, an endmember with 11 B/ 10 B > 4 is required to account for the solar 11 B/ 10 B ratio. Two leading hypotheses for the sources of 11 B-rich components include low energy spallation in the Sun’s parental molecular cloud and the neutrino process during the supernova explosions. In this study, lithium and boron elemental and isotopic compositions of seven porphyritic olivine chondrules and one isolated olivine crystal from four meteorites (Allende, Yamato 81020, Asuka 12236, and QUE 97008) were determined in-situ to help constrain which of the two nucleosynthetic mechanisms has most likely supplied the forming Solar System with extra 11 B. Apparent isotopic variations in Li/Si, B/Si, δ 7 Li and δ 11 B were found in chondrule olivine crystals, but only three chondrules exhibit statistically resolved δ 11 B heterogeneities. Using these three chondrules as a basis for discussion, we evaluated the processes that could potentially cause elemental and isotopic variations of Li and B. We argue that the data can be best understood in the context of condensation of lithium-boron-rich material onto chondrule precursors in (an) initially heterogeneous gaseous reservoir(s), which could be realized if the Solar System derived 11 B-rich components from a supernova or supernovae.

58 GEOSCIENCES↗
DOE OSTI2024

Stress Field Dynamics and Fault Slip Potential in the Paradox Basin

Abstract The Paradox Basin, straddling Utah, Colorado, Arizona, and New Mexico is characterized by an intricate amalgamation of evaporites and clastic layers and is dominated by prominent salt walls and related subsurface structures. Our research offers a new examination of the stress distribution across the basin, deriving from continuous and discrete stress measurements conducted in boreholes in the region and focal mechanism analysis, emphasizing variations over salt structures. Integrating Coulomb failure criteria with probabilistic methods, we assess potential fault movements resulting from fluid pressure alterations. Our approach provides a comprehensive understanding of the Paradox Basin's state of stress, showing a continuous change of the maximum horizontal stress orientation from N‐S at the Wasatch Fault Zone to WNW‐ESE in the northern part of the Paradox Basin and to WSW‐ENE in the southern part of the basin. Further East, into the Colorado Plateau and the Uncompahgre Uplift, the S H max orientation becomes E‐W. Decoding stress orientation dynamics has enabled critical insights into fault slip potential, especially in the basin's northern region. The salt wall faults are less likely to slip, and the Paradox Formation's evaporite and clastic rock sequence can serve as a potential low seismic risk target for carbon storage and hydrocarbon extraction.

Geochemistry & Geophysics↗
DOE OSTI2015

Vertical deformation at western part of Sumatra

This research tries to make advancement in GPS signal processing to estimate the interseismic vertical deformation field at western part of Sumatra Island. The data derived by Continuous Global Positioning System (CGPS) from Badan Informasi Geospasial (BIG) between 2010 and 2012. GPS Analyze at Massachusetts Institute of Technology (GAMIT) software and Global Kalman Filter (GLOBK) software are used to process the GPS signal to estimate the vertical velocities of the CGPS station. In order to minimize noise due to atmospheric delay, Vienna Mapping Function 1 (VMF1) is used as atmospheric parameter model and include daily IONEX file provided by the Center for Orbit Determination in Europe (CODE) as well. It improves GAMIT daily position accuracy up to 0.8 mm. In a second step of processing, the GLOBK is used in order to estimate site positions and velocities in the ITRF08 reference frame. The result shows that the uncertainties of estimated displacement velocity at all CGPS stations are smaller than 1.5 mm/yr. The subsided deformation patterns are seen at the northern and southern part of west Sumatra. The vertical deformation at northern part of west Sumatra indicates postseismic phase associated with the 2010 and 2012 Northern Sumatra earthquakes and also the long-term postseismic associated with the 2004 and 2005 Northern Sumatra earthquakes. The uplifted deformation patterns are seen from Bukit Tinggi to Seblat which indicate a long-term interseismic phase after the 2007 Bengkulu earthquake and 2010 Mentawai earthquake. GANO station shows a subsidence at rate 12.25 mm/yr, indicating the overriding Indo-Australia Plate which is dragged down by the subducting Southeast Asian Plate.

58 GEOSCIENCES↗
DOE OSTI2018

A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons

A conservative scheme has been formulated and confirmed for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from the drift kinetic equation is used to time advance the electron density perturbation by using the perturbed mechanical flow calculated from the parallel vector potential, and the parallel vector potential is solved by using the perturbed canonical flow from the perturbed distribution function. In gyrokinetic particle simulations using this new scheme, the shear Alfvén wave dispersion relation in the shearless slab and continuum damping in the sheared cylinder have been recovered. The new scheme overcomes the stringent requirement in the conventional perturbative simulation method that perpendicular grid size needs to be as small as electron collisionless skin depth even for the long wavelength Alfvén waves. The new scheme also circumvents the problem in the typical approach that an unphysically large parallel electric field arises due to the inconsistency between electrostatic potential calculated from the perturbed density and vector potential calculated from the perturbed canonical flow. Lastly, the gyrokinetic particle simulations of the Alfvén waves in sheared cylinder have superior numerical properties compared with the fluid simulations, which suffer from numerical difficulties associated with singular mode structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

THE GALEX ARECIBO SDSS SURVEY. VII. THE BIVARIATE NEUTRAL HYDROGEN-STELLAR MASS FUNCTION FOR MASSIVE GALAXIES

We present the bivariate neutral atomic hydrogen (H I)-stellar mass function (HISMF) φ(M{sub H{sub I}}, M{sub *}) for massive (log M{sub *}/M{sub ☉} \gt 10) galaxies derived from a sample of 480 local (0.025 < z < 0.050) galaxies observed in H I at Arecibo as part of the GALEX Arecibo SDSS Survey. We fit six different models to the HISMF and find that a Schechter function that extends down to a 1% H I gas fraction, with an additional fractional contribution below that limit, is the best parameterization of the HISMF. We calculate Ω{sub H{sub I,{sub M{sub *>10{sup 1}{sup 0}}}}} and find that massive galaxies contribute 41% of the H I density in the local universe. In addition to the binned HISMF, we derive a continuous bivariate fit, which reveals that the Schechter parameters only vary weakly with stellar mass: M{sub H{sub I}{sup *}}, the characteristic H I mass, scales as M{sub *}{sup 0.39}; α, the slope of the HISMF at moderate H I masses, scales as M{sub *}{sup 0.07}; and f, the fraction of galaxies with H I gas fraction greater than 1%, scales as M{sub *}{sup -0.24}. The variation of f with stellar mass should be a strong constraint for numerical simulations. To understand the physical mechanisms that produce the shape of the HISMF, we redefine the parameters of the Schechter function as explicit functions of stellar mass and star formation rate (SFR) to produce a trivariate fit. This analysis reveals strong trends with SFR. While M{sub H{sub I}{sup *}} varies weakly with stellar mass and SFR (M{sub H{sub I}{sup *}} ∝ M{sub *}{sup 0.22}, M{sub H{sub I}{sup *}} ∝ SFR{sup –0.03}), α is a stronger function of both stellar mass and especially SFR (α ∝ M{sub *}{sup 0.47}, α ∝ SFR{sup 0.95}). The HISMF is a crucial tool that can be used to constrain cosmological galaxy simulations, test observational predictions of the H I content of populations of galaxies, and identify galaxies whose properties deviate from average trends.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Crustal Strain Rates in the Western United States and Their Relationship with Earthquake Rates

Abstract We present a suite of strain rate models for the western United States based on geologic and geodetic data. The geologic data consist of Quaternary fault-slip rates and the geodetic data consists of a new compilation of Global Positioning System (GPS) velocities derived from continuous, semicontinuous, and campaign measurements. We remove postseismic deformation from the GPS time series in order for our geodetic strain rate model to best capture the interseismic strain accumulation rate. We present models based on either geologic or geodetic data, but also create a hybrid model. Although there are some differences between the models, the large-scale features are the same, with the noticeable exception for the Pacific Northwest where interseismic strain is naturally more distributed than the long-term strain release. We also present a map of earthquake rate densities based on mainshocks, and the result has similar spatial features similar to the strain rate models (at least in the southwestern United States). We perform a general correlation analysis between strain rate and seismicity rate (south of Cascadia) and find a change in linearity between seismicity and strain rates from slow to faster deforming areas with seismicity rates relatively lower for the latter. The extent of that change depends a bit on assumptions made on the declustering and completeness of the catalog, but the finding of a change in slope is robust across the different strain rate models. Linearity for all areas is only expected when Gutenberg–Richter parameters and parameters involved in the conversion from strain to moment rate are uniform across the study area. We discuss these qualifications, but find no single satisfactory explanation for our observation. Moreover, when considering a rather short time and space, theoretical considerations of sampling from a power-law distribution actually predict there to be a power law instead of a linear relationship, generally consistent with our observation.

Geochemistry & Geophysics↗
DOE OSTI2022

Kink solutions with power law tails

We present a comprehensive review about the various facets of kink solutions with a power law tail, which have received considerable attention during the last few years. This area of research is in its early stages; although several aspects have become clear by now, there are a number of issues which have only been partially understood or not understood at all. We first discuss the aspects which are reasonably well known and then address in some detail the issues which are only partially or not understood at all. We present a wide class of higher (than sixth) order field theory models admitting implicit kink as well as mirror kink solutions where the two tails facing each other have a power law or a power-tower type fall off, whereas the other two ends not facing each other could have either an exponential or a power law tail. The models admitting implicit kink solutions where the two ends facing each other have an exponential tail while the other two ends have a power law tail are also discussed. Moreover, we present several field theory models which admit explicit kink solutions with a power law fall off; we note that in all these polynomial models while the potential V ( ϕ ) is continuous, its derivative is discontinuous. We also discuss one of the most important and only partially understood issues of the kink–kink and the kink–antikink forces in case the tails facing each other have a power law fall off. Finally, we briefly discuss the kink–antikink collisions at finite velocity and present some open questions.

power-tower tail↗
NASA NTRS1974

A smoothing algorithm using cubic spline functions

Two algorithms are presented for smoothing arbitrary sets of data. They are the explicit variable algorithm and the parametric variable algorithm. The former would be used where large gradients are not encountered because of the smaller amount of calculation required. The latter would be used if the data being smoothed were double valued or experienced large gradients. Both algorithms use a least-squares technique to obtain a cubic spline fit to the data. The advantage of the spline fit is that the first and second derivatives are continuous. This method is best used in an interactive graphics environment so that the junction values for the spline curve can be manipulated to improve the fit.

Smith, R. E., Jr.↗
NASA NTRS1974

Quantum mechanical streamlines. I - Square potential barrier

Exact numerical calculations are made for scattering of quantum mechanical particles hitting a square two-dimensional potential barrier (an exact analog of the Goos-Haenchen optical experiments). Quantum mechanical streamlines are plotted and found to be smooth and continuous, to have continuous first derivatives even through the classical forbidden region, and to form quantized vortices around each of the nodal points. A comparison is made between the present numerical calculations and the stationary wave approximation, and good agreement is found between both the Goos-Haenchen shifts and the reflection coefficients. The time-independent Schroedinger equation for real wavefunctions is reduced to solving a nonlinear first-order partial differential equation, leading to a generalization of the Prager-Hirschfelder perturbation scheme. Implications of the hydrodynamical formulation of quantum mechanics are discussed, and cases are cited where quantum and classical mechanical motions are identical.

Hirschfelder, J. O.↗
NASA NTRS1976

Structural optimization via a design space hierarchy

Mathematical programming techniques provide a general approach to automated structural design. An iterative method is proposed in which design is treated as a hierarchy of subproblems, one being locally constrained and the other being locally unconstrained. It is assumed that the design space is locally convex in the case of good initial designs and that the objective and constraint functions are continuous, with continuous first derivatives. A general design algorithm is outlined for finding a move direction which will decrease the value of the objective function while maintaining a feasible design. The case of one-dimensional search in a two-variable design space is discussed. Possible applications are discussed. A major feature of the proposed algorithm is its application to problems which are inherently ill-conditioned, such as design of structures for optimum geometry.

Vanderplaats, G. N.↗
NASA NTRS1979

Simple spline-function equations for fracture mechanics calculations

The paper presents simple spline-function equations for fracture mechanics calculations. A spline function is a sequence of piecewise polynomials of degree n greater than 1 whose coefficients are such that the function and its first n-1 derivatives are continuous. Second-degree spline equations are presented for the compact, three point bend, and crack-line wedge-loaded specimens. Some expressions can be used directly, so that for a cyclic crack propagation test using a compact specimen, the equation given allows the cracklength to be calculated from the slope of the load-displacement curve. For an R-curve test, equations allow the crack length and stress intensity factor to be calculated from the displacement and the displacement ratio.

Orange, T. W.↗
NASA NTRS1981

Analysis of capacitive heat exchangers, part 1

Governing equations are developed for a simple capacitive heat exchanger. This type of heat exchanger consists of hot spherical particles falling through an ascending cold gas stream. The assumptions made in deriving the continuity, momentum and energy equations are clearly stated. The analysis yields a system of first order, ordinary, nonlinear equations which form a complex boundary value problem. The method of solution is presented together with a comparison between the performance of capacitive heat exchangers and conventional counter flow ones.

Schonfeld, D.↗
NASA NTRS1982

Global differential geometry: An introduction for control engineers

The basic concepts and terminology of modern global differential geometry are discussed as an introduction to the Lie theory of differential equations and to the role of Grassmannians in control systems analysis. To reach these topics, the fundamental notions of manifolds, tangent spaces, vector fields, and Lie algebras are discussed and exemplified. An appendix reviews such concepts needed for vector calculus as open and closed sets, compactness, continuity, and derivative. Although the content is mathematical, this is not a mathematical treatise but rather a text for engineers to understand geometric and nonlinear control.

Doolin, B. F.↗
NASA NTRS1989

Numerical simulation of thermocapillary bubble migration under microgravity for large Reynolds and Marangoni numbers

A numerical procedure in which the Navier-Stokes equations are discretized using tightly coupled discretizations of pressure derivatives and continuity equations is used here to extend the range of known terminal velocities of gaseous bubbles in liquids well beyond that in previous investigations. Computations performed for Reynolds numbers up to 2000 and Marangoni numbers up to 1000 show only a modest variation of the scaled bubble velocity between 0.16 and 0.5. The bubble velocity is influenced more by the Marangoni number than by the Reynolds number.

Balasubramaniam, R.↗