Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic equivalent”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impedance Scan of Inverter-Based Resources and Diesel Generator for Stability Analysis

Impedance-based methods are widely used for power system stability analysis with inverter-based resources (IBRs), e.g., assessing dynamic interactions between the power grid and an IBR, control interactions between multiple IBRs, and the sub-synchronous oscillation and damping phenomenon. Since it is difficult to get a numerical model 100% matching with the hardware IBR, using the hardware inverter directly to obtain its output impedance has become a prominent approach nowadays. Therefore, this article presents the impedance scan using hardware IBRs, and also a hardware diesel generator as it still stays with the grid before the grid completely goes to renewable. The devices under test (DuTs) for the impedance scan includes two 3-..phi.., 480 V, 60 Hz commercial grid-forming IBRs (one of 250 kVA and another of 125 kVA rating) in series with ..delta..-Y transformers, one 3-..phi.., 480 V, 60 Hz commercial grid-following IBR (of 125 kVA rating), and a 3-..phi.., 480 V, 60 Hz commercial diesel generator (of 187.5 kVA rating). Using voltage signals perturbed with sub-, inter-, and higher harmonic components, and measuring the current response, the positive-sequence impedances are computed via an offline-based post-analysis. Moreover, best-fit transfer functions are estimated that closely resemble the measured data points of the positive-sequence impedances. Based on the observations from various outcomes of the hardware experiments, this article also provides some fundamental insights on the equivalent positive-sequence impedance of a combination of multiple hardware components by comparing the estimated and the empirically computed impedances. A comparative insight on the damping capability of the DuTs using the positive-sequence impedances of the hardware is also discussed.

current measurement↗
DOE OSTI2023

IDAES-PSE 2.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.0.0 Release Highlights Removal of deprecated features from IDAES v1 Update to Pyomo v6.5 – this required a number of updates to support the new NL solver writer and to address some changes in Pyomo Creation of new testing suite for backward compatibility, model robustness and verification More general implementation of the Helmholtz EoS. This brings some new features like standard property diagrams, choice of mass or mole basis, and new state variable options Standardizing names in Heat Exchanger models (breaking change from v2.0.0a2): Control Volumes named hot_side and cold_side Ports names hot_side_inlet, hot_side_outlet, cold_side_inlet and cold_side_outlet Config Blocks names hot_side_config and cold_side_config Config arguments for user provided names for each side: hot_side_name and cold_side_name. Updating Keras surrogate tool to use v1.1 of OMLT New prototype API for model initialization (idaes.core.initialization) The new API uses "Model Initializer" objects instead of class methods, allowing for the definition of multiple initialization routines for a single model A number of common, model agnostic initialization routines have also been defined, including initialization from data, block-decomposition and a general hierarchical approach equivalent to the existing method for common unit models New metadata for thermophysical properties – valid_range This can be used to record the range of values over which a property value can be trusted, such as the range of experimental data used to regress parameters A number of new utility functions have been added to check for properties with values outside the valid range and to set bounds based on this metadata Updated construction of balance expressions in Control Volumes to remove unneeded terms In the past, unneeded terms were added as a constant 0 term, however they will now be dropped entirely from the expression This was necessary due to more strict unit checking in the new Pyomo solver writer which no longer ignores 0 terms Updates to metadata for thermophysical properties to better define known properties and units of measurement This results in more strict enforcement of standard naming for thermophysical and reaction properties Users can still define custom properties, but these must be done explicitly using the define_custom_properties() method instead of being implicitly created by add_property() Updated convergence tester utility tool to support definition of benchmark files (JSON format) and comparison of performance to benchmarks Set default iteration limit for IPOPT in IDAES config to 200 iterations Update scaling of example models to work with new Pyomo NL solver writer Improve testing of extensions and examples infrastructure to avoid need for downloading files Updated distillation column to centralize common functionality and remove a number of Pyomo warnings

IDAES↗
NASA NTRS2002

Length Scale and Gravity Effects on Microgravity Boiling Heat Transfer

Boiling is a complex phenomenon where hydrodynamics, heat transfer, mass transfer, and interfacial phenomena are tightly interwoven. An understanding of boiling and critical heat flux in microgravity environments is of importance to space based hardware and processes such as heat exchange, cryogenic fuel storage and transportation, electronic cooling, and material processing due to the large amounts of heat that can be removed with relatively little increase in temperature. Although research in this area has been performed in the past four decades, the mechanisms by which heat is removed from surfaces in microgravity are still unclear. In earth gravity, buoyancy is an important parameter that affects boiling heat transfer through the rate at which bubbles are removed from the surface. A simple model describing the bubble departure size based on a quasistatic force balance between buoyancy and surface tension is given by the Fritz [I] relation: Bo(exp 1/2) = 0.0208 theta where Bo is the ratio between buoyancy and surface tension forces. For small, rapidly growing bubbles, inertia associated with the induced liquid motion can also cause bubble departure. In microgravity, the magnitude of effects related to natural convection and buoyancy are small and physical mechanisms normally masked by natural convection in earth gravity such as Marangoni convection can substantially influence the boiling and vapor bubble dynamics. CHF (critical heat transfer) is also substantially affected by microgravity. In 1 g environments, Bo has been used as a correlating parameter for CHF. Zuber's CHF model for an infinite horizontal surface assumes that vapor columns formed by the merger of bubbles become unstable due to a Helmholtz instability blocking the supply of liquid to the surface. The jets are spaced lambda(sub D) apart, where lambda(sub D) = 2pi square root of 3[(sigma)/(g(rho(sub l) - rho(sub v)](exp 1/2) = 2pi square root of 3 L Bo(exp -1/2) = square root of 3 lambda(sub c) and is the wavelength that amplifies most rapidly. The critical wavelength, lambda(sub c), is the wavelength below which a vapor layer underneath a liquid layer is stable. For heaters with Bo smaller than about 3 (heaters smaller than lambda(sub D)), the above model is not applicable, and surface tension effects dominate. Bubble coalescence is thought to be the mechanism for CHF under these conditions. Small Bo can result by decreasing the size of a heater in earth gravity, or by operating a large heater in a lower gravity environment. In the microgravity of space, even large heaters can have low Bo, and models based on Helmholtz instability should not be applicable. The macrolayer model of Haramura and Katto is dimensionally equivalent to Zuber's model and has the same dependence on gravity, so it should not be applicable as well. The goal of this work is to determine how boiling heat transfer mechanisms in a low-g environment are altered from those at higher gravity levels. Boiling data using a microheater array was obtained under gravity environments ranging from 1.8 g to 0.02 g with heater sizes ranging from 2.7 mm to 1 mm. The boiling behavior for 2.7 mm at 0.02 g looked quite similar to boiling on the 1 mm heater at 1 g-the formation of a large primary bubble surrounded by smaller satellite bubbles was observed under both conditions. The similarity suggests that for heaters smaller than some fraction of I(sub c), coalescence and surface tension dominate boiling heat transfer. It also suggests that microgravity boiling can be studied by studying boiling on very small heaters.

Kim, Jungho↗
DOE OSTI2021

Influence of Surface and Structural Variations in Donor–Acceptor–Donor Sensitizers on Photoelectrocatalytic Water Splitting

Conjugated organic chromophores composed of linked donor (D) and acceptor (A) moieties have attracted considerable attention for photoelectrochemical applications. In this work, we compare the optoelectronic properties and photoelectrochemical performance of two D–A–D structural isomers with thiophene-X-carboxylic acid (X denotes 3 and 2 positions) derivatives and 2,1,3-benzothiadiazole as the D and A moieties, respectively. 5,5′-(Benzo[c][1,2,5]thiadiazole-4,7-diyl)bis(thiophene-3-carboxylic acid), BTD1, and 5,5′-(benzo[c][1,2,5]thiadiazole-4,7-diyl)bis(thiophene-2-carboxylic acid), BTD2, were employed in the study to understand how structural isomers affect surface attachments within chromophore–catalyst assemblies and their influence on charge-transfer dynamics. Crystal structures revealed that varying the position of the −COOH anchoring group causes the molecules to either contort out of a plane (BTD1) or adopt a near-perfect planar conformation (BTD2). BTD1 and BTD2 were co-loaded with either a water oxidation catalyst, [Ru(2,6-bis(1-methylbenzimidazol-2-yl)pyridine)-(4,4′-((HO) 2 OPCH 2 )2-2,2′-bipyridine)(OH 2 )] 2 , RuCt 2+ , or proton reduction catalyst [Ni(P 2 Ph N 2 C 6 H 4 CH 2 PO 3 H 2 ) 2 ] 2+ , NiCt 2+ , on oxide electrodes to facilitate photodriven water splitting reactions. Emission quenching measurements indicate that both BTD1 and BTD2 inject electrons into n-type SnO 2 |TiO 2 electrodes and holes into p-type NiO semiconductors from their respective excited states at high efficiencies >60%. Photocurrent densities of chromophore–catalyst assemblies obtained using linear sweep voltammetry (LSV) show that BTD2-sensitized photoanodes generate significantly more photocurrent than BTD1-sensitized electrodes; however, both exhibit similar performances at the photocathode. Photoelectrocatyltic measurements demonstrate that both BTD1 and BTD2 performed similarly, generating Faradaic efficiencies of 39 and 38% at the anode or 61 and 79% at the cathode. Transient absorption measurements suggest that the differences between the LSV and photoelectrocatalytic measurements result from the differences in quantum yields of the photogenerated redox equivalents, which is also a reflection of the varying metal oxide surface conformation. Our findings suggest that BTD2 should be investigated further in photocathodic studies since it has the structural advantage of being incorporated into diverse types of chromophore–catalyst assemblies.

Chromophores↗
NASA NTRS2002

Polar Wander on Triton and Pluto Due to Volatile Migration

Polar wander may occur on Triton and Pluto because of volatile migration. Triton, with its low obliquity, can theoretically sublimate volatiles (mostly nitrogen) at the rate of approximately 10(exp 14) kilograms per year from the equatorial regions and deposit them at the poles. Assuming Triton to be rigid on the sublimation timescale, after approximately 10(exp 5) years the polar caps would become large enough to cancel the rotational flattening, with a total mass equivalent to a global layer approximately 120-250 m in depth. At this point the pole wanders about the tidal bulge axis, which is the line joining Triton and Neptune. Rotation about the bulge axis might be expected to disturb the leading side/trailing side cratering statistics. Because no such disturbance is observed, it may be that Triton's mantle viscosity is too high but its surface volatile inventory is too low to permit wander. On the other hand, its mantle viscosity might be low, so that any uncompensated cap load might be expected to wander toward the tidal bulge axis. In this case, the axis of wander passes through the equator from the leading side to the trailing side; rotation about this wander axis would not disturb the cratering statistics. Low-viscosity polar wander may explain the bright southern hemisphere: this is the pole which is wandering toward the equator. In any case the permanent polar caps may be geologically very young. Polar wander may possibly take place on Pluto, due to its obliquity oscillations and perihelion-pole geometry. However, Pluto is probably not experiencing any wander at present. The Sun has been shining strongly on the poles over the last half of the obliquity cycle, so that volatiles should migrate to the equator, stabilizing the planet against wander. Spacecraft missions to Triton and Pluto which measure the dynamical flattening could give information about the accumulation of volatiles at the poles. Such information is best obtained by measuring gravity and topography from orbiters, as was done for Mars with the highly successful Mars Global Surveyor.

Rubincam, David Parry↗
NASA NTRS1993

Stable isotopic biogeochemistry of carbon and nitrogen in a perennially ice-covered Antarctic lake

Lake Hoare (77 degrees 38' S, 162 degrees 53' E) is an amictic, oligotrophic, 34-m-deep, closed-basin lake in Taylor Valley, Antarctica. Its perennial ice cover minimizes wind-generated currents and reduces light penetration, as well as restricts sediment deposition into the lake and the exchange of atmospheric gases between the water column and the atmosphere. The biological community of Lake Hoare consists solely of microorganisms -- both planktonic populations and benthic microbial mats. Lake Hoare is one of several perennially ice-covered lakes in the McMurdo Dry Valleys that represent the end-member conditions of cold desert and saline lakes. The dry valley lakes provide a unique opportunity to examine lacustrine processes that operate at all latitudes, but under an extreme set of environmental conditions. The dry valley lakes may also offer a valuable record of catchment and global changes in the past and present. Furthermore, these lakes are modern-day equivalents of periglacial lakes that are likely to have been common during periods of glacial maxima at temperate latitudes. We have analyzed the dissolved inorganic carbon (DIC) of Lake Hoare for delta 13C and the organic matter of the sediments and sediment-trap material for delta 13C and delta 15N. The delta 13C of the DIC indicates that 12C is differentially removed in the shallow, oxic portions of the lake via photosynthesis. In the anoxic portions of the lake (27-34 m) a net addition of 12C to the DIC pool occurs via organic matter decomposition. The dissolution of CaCO3 at depth also contributes to the DIC pool. Except near the Canada Glacier where a substantial amount of allochthonous organic matter enters the lake, the organic carbon being deposited on the lake bottom at different sites is isotopically similar, suggesting an autochthonous source for the organic carbon. Preliminary inorganic carbon flux calculations suggest that a high percentage of the organic carbon fixed in the water column is remineralized as it falls through the water column. At nearby Lake Fryxell, the substantial (relative to Lake Hoare) glacial meltstream input overprints Fryxell's shallow-water biological delta 13C signal with delta 13C-depleted DIC. In contrast, Lake Hoare is not significantly affected by surface-water input and mixing, and therefore the delta 13C patterns observed arise primarily from biological dynamics within the lake. Organic matter in Lake Hoare is depleted in 15N, which we suggest is partially the result of the addition of relatively light inorganic nitrogen into the lake system from terrestrial sources.

NASA Program Exobiology↗
NASA NTRS2006

Semiautonomous Avionics-and-Sensors System for a UAV

Unmanned Aerial Vehicles (UAVs) autonomous or remotely controlled pilotless aircraft have been recently thrust into the spotlight for military applications, for homeland security, and as test beds for research. In addition to these functions, there are many space applications in which lightweight, inexpensive, small UAVS can be used e.g., to determine the chemical composition and other qualities of the atmospheres of remote planets. Moreover, on Earth, such UAVs can be used to obtain information about weather in various regions; in particular, they can be used to analyze wide-band acoustic signals to aid in determining the complex dynamics of movement of hurricanes. The Advanced Sensors and Electronics group at Langley Research Center has developed an inexpensive, small, integrated avionics-and-sensors system to be installed in a UAV that serves two purposes. The first purpose is to provide flight data to an AI (Artificial Intelligence) controller as part of an autonomous flight-control system. The second purpose is to store data from a subsystem of distributed MEMS (microelectromechanical systems) sensors. Examples of these MEMS sensors include humidity, temperature, and acoustic sensors, plus chemical sensors for detecting various vapors and other gases in the environment. The critical sensors used for flight control are a differential- pressure sensor that is part of an apparatus for determining airspeed, an absolute-pressure sensor for determining altitude, three orthogonal accelerometers for determining tilt and acceleration, and three orthogonal angular-rate detectors (gyroscopes). By using these eight sensors, it is possible to determine the orientation, height, speed, and rates of roll, pitch, and yaw of the UAV. This avionics-and-sensors system is shown in the figure. During the last few years, there has been rapid growth and advancement in the technological disciplines of MEMS, of onboard artificial-intelligence systems, and of smaller, faster, and smarter wireless telemetry systems. The major attraction of MEMS lies in orders-of-magnitude reductions of power requirements relative to traditional electronic components that perform equivalent functions. In addition, the compactness of MEMS, relative to functionally equivalent traditional electronics systems, makes MEMS attractive for UAV applications. Recent advances in MEMS have made it possible to produce pressure, acceleration, humidity, and temperature sensors having masses in subgram range and possessing sensitivities and accuracies comparable to those of larger devices.

Shams, Qamar↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Characterizing Uranus with an Ice giant Planetary Origins Probe (Ice-POP)

We now know from studies of planetary transits and microlensing that Neptune-mass planets are ubitquitous and may be the most common class of planets in the Galaxy. As such it is crucial that we understand the formation and evolution of the ice giant planets in our own solar system so that we can better understand planet formation throughout the galaxy. An entry probe mission to Uranus would help accomplish this goal. In fact the Planetary Decadal Survey recommended a Uranus orbiter with entry probe but did not explore in detail the specifications for the entry probe. NASA Ames is currently studying thermal protection system requirements for such a mission and this has led to questions regarding the minimum interesting science payload of such an entry probe. The single most important in-situ measurement for an ice giant entry probe is a measurement of atmospheric composition. For Uranus this would specifically include the methane and noble gas abundances. An in situ measurement of the methane abundance, from below the methane cloud, would constrain the atmospheric carbon abundance, which is believed to be roughly 30 to 50 times solar. There are hints from the transiting planets that extrasolar ice giants show comparable or even greater enhancements of heavy elements compared to their primary stars. However the origin of this carbon enhancement is controversial. Is Uranus a "failed core" of a larger gas giant or was the atmosphere enhanced by accretion of icy planetesimals' Constraining atmospheric abundances of C and perhaps S or even N from below 5 bars would provide badly needed data to address such issues. A measurement of the N abundance would provide clues on the origin of the planetesimals that formed Uranus. Low N-abundance indicates planetesimals from 'warmer' regions where N was mainly in form of NH3, whereas a strong enrichment could indicate planetesimals / cometary material from the colder outer regions of the nebula. Furthermore CO and HCN have been detected in Neptune but not in Uranus. A measurement of the abundance of either would constrain the source mechanisms for these molecules (exogenic or internal). A major surprise from the Galileo Entry Probe was that the heavier noble gases Ar, Kr, and Xe are enhanced in Jupiter's atmosphere at a level comparable to what was seen for the chemically active volatiles N, C, and S. It had been generally expected that Ar, Kr, and Xe would be present in solar abundances, as all were expected to accrete with hydrogen during the gravitational capture of nebular gases. Enhanced abundances of Ar, Kr, and Xe is equivalent to saying that these noble gases have been separated from hydrogen. There are several mechanisms that could accomplish this but these hypotheses require further testing. Measurement of noble gas abundances in an ice giant would constrain the planetary formation and nebular mechanisms responsible for this enhancement. Standard three-layer models of Uranus find that the outer, predominantly H/He layer of Uranus does not reach pressures high enough (approximately 1 Mbar) for H2 to transition to liquid metallic hydrogen. However, valid models can also be constructed with a smaller intermediate water-rich layer, with hydrogen then reaching the metallic hydrogen phase. If this occurs, He should phase separate from the hydrogen and ``rain out," taking along a substantial abundance of Ne, as suggested for Jupiter (and likely also for Saturn). Hence He and Ne depletions could be probes of the planet's structure in the much deeper interior. A determination of Uranus' atmospheric abundances, particularly of the noble gasses, is thus critical to understanding the formation of Uranus, and giant planets in general. These measurements can only be performed with an entry probe. The second key measurement would be a temperature-pressure sounding to provide ground truth for remote measurements of atmospheric temperature and composition and to constrain the internal heat flow. This would also establish that the methane abundance measurements have indeed been made below any possible methane cloud. Finally an ultra stable oscillator would measure wind speeds and constrain atmospheric dynamics. In our presentation we will discuss the importance of all of these measurements and argue that an entry probe is a crucial component of any ice giant mission.

Marley, Mark S.↗
DOE OSTI2021

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

Dynamic interactions between hypersonic vehicle aerodynamics and propulsion system performance

Described here is the development of a flexible simulation model for scramjet hypersonic propulsion systems. The primary goal is determination of sensitivity of the thrust vector and other system parameters to angle of attack changes of the vehicle. Such information is crucial in design and analysis of control system performance for hypersonic vehicles. The code is also intended to be a key element in carrying out dynamic interaction studies involving the influence of vehicle vibrations on propulsion system/control system coupling and flight stability. Simple models are employed to represent the various processes comprising the propulsion system. A method of characteristics (MOC) approach is used to solve the forebody and external nozzle flow fields. This results in a very fast computational algorithm capable of carrying out the vast number of simulation computations needed in guidance, stability, and control studies. The three-dimensional fore- and aft body (nozzle) geometry is characterized by the centerline profiles as represented by a series of coordinate points and body cross-section curvature. The engine module geometry is represented by an adjustable vertical grid to accommodate variations of the field parameters throughout the inlet and combustor. The scramjet inlet is modeled as a two-dimensional supersonic flow containing adjustable sidewall wedges and multiple fuel injection struts. The inlet geometry including the sidewall wedge angles, the number of injection struts, their sweepback relative to the vehicle reference line, and strut cross-section are user selectable. Combustion is currently represented by a Rayleigh line calculation including corrections for variable gas properties; improved models are being developed for this important element of the propulsion flow field. The program generates (1) variation of thrust magnitude and direction with angle of attack, (2) pitching moment and line of action of the thrust vector, (3) pressure and temperature distributions throughout the system, and (4) performance parameters such as thrust coefficient, specific impulse, mass flow rates, and equivalence ratio. Preliminary results are in good agreement with available performance data for systems resembling the NASP vehicle configuration.

Flandro, G. A.↗
DOE OSTI2024

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2014

Interannual Variation in Phytoplankton Class-Specific Primary Production at a Global Scale

We used the NASA Ocean Biogeochemical Model (NOBM) combined with remote sensing data via assimilation to evaluate the contribution of 4 phytoplankton groups to the total primary production. First we assessed the contribution of each phytoplankton groups to the total primary production at a global scale for the period 1998-2011. Globally, diatoms were the group that contributed the most to the total phytoplankton production (50, the equivalent of 20 PgC y-1. Coccolithophores and chlorophytes each contributed to 20 (7 PgC y-1 of the total primary production and cyanobacteria represented about 10 (4 PgC y(sub-1) of the total primary production. Primary production by diatoms was highest in high latitude (45) and in major upwelling systems (Equatorial Pacific and Benguela system). We then assessed interannual variability of this group-specific primary production over the period 1998-2011. Globally the annual relative contribution of each phytoplankton groups to the total primary production varied by maximum 4 (1-2 PgC y-1. We assessed the effects of climate variability on the class-specific primary production using global (i.e. Multivariate El Nio Index, MEI) and regional climate indices (e.g. Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO) and North Atlantic Oscillation (NAO)). Most interannual variability occurred in the Equatorial Pacific and was associated with climate variability as indicated by significant correlation (p 0.05) between the MEI and the class-specific primary production from all groups except coccolithophores. In the Atlantic, climate variability as indicated by NAO was significantly correlated to the primary production of 2 out of the 4 groups in the North Central Atlantic (diatomscyanobacteria) and in the North Atlantic (chlorophytes and coccolithophores). We found that climate variability as indicated by SAM had only a limited effect on the class-specific primary production in the Southern Ocean. These results provide a modeling and data assimilation perspective to phytoplankton partitioning of primary production and contribute to our understanding of the dynamics of the carbon cycle in the oceans at a global scale.

Chl-a↗
NASA NTRS2000

Gravity Effects in Diffusive Coarsening of Bubble Lattices: von Neumann's Law

von Neumann modelled the evolution of two-dimensional soap froths as a purely diffusive phenomenon; the area growth of a given cell was found to depend only on the geometry of the bubble lattice. In the model, hexagons are stable, pentagons shrink and heptagons grow. The simplest equivalent to the area growth law is / approximately t(sub beta). The result depends on assuming (1) an incompressible gas; (2) bubble walls which meet at 120 deg and (3) constant wall thickness and curvature. Each assumption is borne out in experiments except the last one: bubble wall thickness between connecting cells varies in unit gravity because of gravity drainage. The bottom part of the soap membrane is thickened, the top part is thinned, such that gas diffusion across the membrane shows a complex dependence on gravity. As a result, experimental tests of von Neumann's law have been influenced by effects of gravity; fluid behavior along cell borders can give non-uniform wall thicknesses and thus alter the effective area and gas diffusion rates between adjacent bubbles. For area plotted as a function of time, Glazier (J.A. Glazier, S.P. Gross, and I. Stavans, Phys. Rev. A. 36, 306 (1987); J. Stavans, J.A, Glazier, Phys. Rev. Lett. 62, 1318 (1989).) suggest that in some cases their failure to observe von Neumann's predicted growth exponent ((sup beta)theor(sup =1; beta)exp(sup =0.70 + 0.10)) may have been the result of such "fluid drainage onto the lower glass plate". Additional experiments which varied plate spacing gave different beta exponents in a fashion consistent with this suggestion. During preliminary long duration experiments (approximately 100 h) aboard Spacelab-J, a low-gravity test of froth coarsening has examined (1) power law scaling of von Neumann's law (beta values) in the appropriate diffusive limits; (2) new bubble lattice dynamics such as greater fluid wetting behavior on froth membranes in low gravity; and (3) explicit relations for the gravity dependence of the second moment (or disorder parameter) governing the geometric spread in cell-sidedness around the mean of perfect hexagonal filling. By reducing the gravity-induced distortion in lattice wall thickness, the diffusion-limited regime of bubble coarsening becomes available for performing critical tests of network dynamics.

Noever, David A.↗
NASA NTRS2013

Interannual Variation in Phytoplankton Primary Production at a Global Scale

We used the NASA Ocean Biogeochemical Model (NOBM) combined with remote sensing data via assimilation to evaluate the contribution of four phytoplankton groups to the total primary production. First, we assessed the contribution of each phytoplankton groups to the total primary production at a global scale for the period 1998-2011. Globally, diatoms contributed the most to the total phytoplankton production ((is)approximately 50%, the equivalent of 20 PgC·y1). Coccolithophores and chlorophytes each contributed approximately 20% ((is) approximately 7 PgC·y1) of the total primary production and cyanobacteria represented about 10% ((is) approximately 4 PgC·y1) of the total primary production. Primary production by diatoms was highest in the high latitudes ((is) greater than 40 deg) and in major upwelling systems (Equatorial Pacific and Benguela system). We then assessed interannual variability of this group-specific primary production over the period 1998-2011. Globally the annual relative contribution of each phytoplankton groups to the total primary production varied by maximum 4% (1-2 PgC·y1). We assessed the effects of climate variability on group-specific primary production using global (i.e., Multivariate El Niño Index, MEI) and "regional" climate indices (e.g., Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO) and North Atlantic Oscillation (NAO)). Most interannual variability occurred in the Equatorial Pacific and was associated with climate variability as indicated by significant correlation (p (is) less than 0.05) between the MEI and the group-specific primary production from all groups except coccolithophores. In the Atlantic, climate variability as indicated by NAO was significantly correlated to the primary production of 2 out of the 4 groups in the North Central Atlantic (diatoms/cyanobacteria) and in the North Atlantic (chlorophytes and coccolithophores). We found that climate variability as indicated by SAM had only a limited effect on group-specific primary production in the Southern Ocean. These results provide a modeling and data assimilation perspective to phytoplankton partitioning of primary production and contribute to our understanding of the dynamics of the carbon cycle in the oceans at a global scale.

Coccolithophores↗
DOE OSTI2020

Canonical Transformation of Potential Model Hamiltonian Mechanics to Geometrical Form I

Using the methods of symplectic geometry, we establish the existence of a canonical transformation from potential model Hamiltonians of standard form in a Euclidean space to an equivalent geometrical form on a manifold, where the corresponding motions are along geodesic curves. The advantage of this representation is that it admits the computation of geodesic deviation as a test for local stability, shown in recent previous studies to be a very effective criterion for the stability of the orbits generated by the potential model Hamiltonian. We describe here an algorithm for finding the generating function for the canonical transformation and describe some of the properties of this mapping under local diffeomorphisms. We give a convergence proof for this algorithm for the one-dimensional case, and provide a precise geometric formulation of geodesic deviation which relates the stability of the motion in the geometric form to that of the Hamiltonian standard form. We apply our methods to a simple one-dimensional harmonic oscillator and conclude with a discussion of the relation of bounded domains in the two representations for which Morse theory would be applicable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2023

A Multi-Pixel Split-Window Approach to Sea Surface Temperature Retrieval from Thermal Imagers with Relatively High Radiometric Noise: Preliminary Studies

In the following decade(s), a set of satellite missions carrying thermal infrared (TIR) imagers with a relatively high noise equivalent differential temperature (NEdT) are expected, e.g., the high resolution TIR imagers flying on the future Thermal infraRed Imaging Satellite for High-resolution Natural resource Assessment (TRISHNA), Land Surface Temperature Monitoring (LSTM) and NASA-JPL/ASI Surface Biology and Geology Thermal (SBG) missions or the secondary payload on board the ESA Earth Explorer 10 Harmony. The instruments on board these missions are expected to be characterized by an NEdT of ⪆0.1 K. In order to reduce the impact of radiometric noise on the retrieved sea surface temperature (SST), this study investigates the possibility of applying a multi-pixel atmospheric correction based on the hypotheses that (i) the spatial variability scales of radiatively active atmospheric variables are, on average, larger than those of the SST and (ii) the effect of atmosphere is accounted for via the split window (SW) difference. Based on 32 Sentinel 3 SLSTR case studies selected in oceanic regions where SST features are mainly driven by meso to sub-mesoscale turbulence (e.g., corresponding to major western boundary currents), this study documents that the local spatial variability of the SW difference term on the scale of ≃3 × 3 km2 is comparable with the noise associated with the SW difference. Similarly, the power spectra of the SW term are shown to have, for small scales, the behavior of white noise spectra. On this basis, we suggest to average the SW term and to use it for the atmospheric correction procedure to reduce the impact of radiometric noise. In principle, this methodology can be applied on proper scales that can be dynamically defined for each pixel. The applicability of our findings to high-resolution TIR missions is discussed and an example of an application to ECOSTRESS data is reported.

sea surface temperature↗
NASA NTRS2008

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗