Search NASA⌕ Search

SEARCH · Search NASA

Results for “STEP FUNCTION”

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 325 records · Page 18

Operator dynamics in Floquet many-body systems

We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading. Published by the American Physical Society 2025

Yoshimura, Takato (ORCID:0000000309159846)↗

Estimating Irrigation Water Requirements using MODIS Vegetation Indices and Inverse Biophysical Modeling

An inverse process approach using satellite-driven (MODIS) biophysical modeling was used to quantitatively assess water resource demand in semi-arid and arid agricultural lands by comparing the carbon and water flux modeled under both equilibrium (in balance with prevailing climate) and non-equilibrium (irrigated) conditions. Since satellite observations of irrigated areas show higher leaf area indices (LAI) than is supportable by local precipitation, we postulate that the degree to which irrigated lands vary from equilibrium conditions is related to the amount of irrigation water used. For an observation year we used MODIS vegetation indices, local climate data, and the SiB2 photosynthesis-conductance model to examine the relationship between climate and the water stress function for a given grid-cell and observed leaf area. To estimate the minimum amount of supplemental water required for an observed cell, we added enough precipitation to the prevailing climatology at each time step to minimize the water stress function and bring the soil to field capacity. The experiment was conducted on irrigated lands on the U.S. Mexico border and Central Asia and compared to estimates of irrigation water used.

Imhoff, Marc L.↗

Tunneling calculations for GaAs-Al(x)Ga(1-x)As graded band-gap sawtooth superlattices

The transmission resonance spectra and tunneling current-voltage characteristics for direct conduction band electrons in sawtooth GaAs-Al(x)Ga(1-x)As superlattices are computed. Only direct-gap interfaces are considered. It is found that sawtooth superlattices exhibit resonant tunneling similar to that in step superlattices, manifested by correlation of peaks and regions of negative differential resistance in the current-voltage curves with transmission resonances. The Stark shift of the resonances of step-barrier superlattices is a linear function of the field, whereas in sawtooth superlattices under strong fields the shift is not a simple function of the field. This follows from the different ways in which the two structures deform under uniform electric fields: the sawtooth deforms into a staircase, at which field strength all barriers to tunneling are eradicated. The step-barrier superlattice always presents some barrier to tunneling, no matter how high the electric field strength.

Forrest, Kathrine↗

Bidirectional Suzuki Catalyst Transfer Polymerization of Poly( p -phenylene)

Suzuki catalyst transfer polymerization (SCTP) has emerged as an effective method for accessing length-controlled π-conjugated poly(p-phenylene). Regio-controlled functional group sequencing along the backbone and the chain ends of synthetic polymers remains a challenge for the successful integration of organic semiconductors in sensors and electronic devices. Here, we report a bidirectional SCTP system based on dinuclear palladium(II) initiators. Functional groups transferred during the initiation and termination steps unlock independently addressable synthetic handles at the polymer core and respective chain ends. These functional groups open opportunities for late-stage regio-controlled and chemoselective derivatizations. Control over key polymer parameters, including molecular weight (M n ), dispersity ( Đ = M w /M n ), degree of polymerization ( DP ), termination efficiency, and functional group interconversion, is corroborated by size exclusion chromatography (SEC) and NMR spectroscopy. The modular design of SCTP initiators, in combination with commercially available terminating groups, represents a highly flexible toolbox for late-stage polymer conjugation, e.g., chemoselective anchoring groups for integration with functional electronics or bioorthogonal conjugation for molecular sensing.

Aldehydes↗

On the numerical solution of the dynamically loaded hydrodynamic lubrication of the point contact problem

The transient analysis of hydrodynamic lubrication of a point-contact is presented. A body-fitted coordinate system is introduced to transform the physical domain to a rectangular computational domain, enabling the use of the Newton-Raphson method for determining pressures and locating the cavitation boundary, where the Reynolds boundary condition is specified. In order to obtain the transient solution, an explicit Euler method is used to effect a time march. The transient dynamic load is a sinusoidal function of time with frequency, fractional loading, and mean load as parameters. Results include the variation of the minimum film thickness and phase-lag with time as functions of excitation frequency. The results are compared with the analytic solution to the transient step bearing problem with the same dynamic loading function. The similarities of the results suggest an approximate model of the point contact minimum film thickness solution.

Lim, Sang G.↗

Use of Functionalized Carbon Nanotubes for Covalent Attachment of Nanotubes to Silicon

The purpose of the invention is to covalently attach functionalized carbon nanotubes to silicon. This step allows for the introduction of carbon nanotubes onto all manner of silicon surfaces, and thereby introduction of carbon nano - tubes covalently into silicon-based devices, onto silicon particles, and onto silicon surfaces. Single-walled carbon nanotubes (SWNTs) dispersed as individuals in surfactant were functionalized. The nano - tube was first treated with 4-t-butylbenzenediazonium tetrafluoroborate to give increased solubility to the carbon nanotube; the second group attached to the sidewall of the nanotube has a silyl-protected terminal alkyne that is de-protected in situ. This gives a soluble carbon nanotube that has functional groups appended to the sidewall that can be attached covalently to silicon. This reaction was monitored by UV/vis/NJR to assure direct covalent functionalization.

Tour, James M.↗

Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations

A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.

Changwoon Jang↗

On the relation between upwind-differencing schemes of Godunov, Engquist-Osher and Roe

The upwind differencing first order schemes of Godunov, Engquist-Osher and Roe are discussed on the basis of the inviscid Burgers equations. The differences between the schemes are interpreted as differences between the approximate Riemann solutions on which their numerical flux functions are based. Special attention is given to the proper formulation of these schemes when a source term is present. Second order two step schemes, based on the numerical flux functions of the first order schemes are also described. The schemes are compared in a numerical experiment, and recommendations on their use are included.

Vanleer, B.↗

Improved Ionospheric Electrodynamic Models and Application to Calculating Joule Heating Rates

Improved techniques have been developed for empirical modeling of the high-latitude electric potentials and magnetic field aligned currents (FAC) as a function of the solar wind parameters. The FAC model is constructed using scalar magnetic Euler potentials, and functions as a twin to the electric potential model. The improved models have more accurate field values as well as more accurate boundary locations. Non-linear saturation effects in the solar wind-magnetosphere coupling are also better reproduced. The models are constructed using a hybrid technique, which has spherical harmonic functions only within a small area at the pole. At lower latitudes the potentials are constructed from multiple Fourier series functions of longitude, at discrete latitudinal steps. It is shown that the two models can be used together in order to calculate the total Poynting flux and Joule heating in the ionosphere. An additional model of the ionospheric conductivity is not required in order to obtain the ionospheric currents and Joule heating, as the conductivity variations as a function of the solar inclination are implicitly contained within the FAC model's data. The models outputs are shown for various input conditions, as well as compared with satellite measurements. The calculations of the total Joule heating are compared with results obtained by the inversion of ground-based magnetometer measurements. Like their predecessors, these empirical models should continue to be a useful research and forecast tools.

Weimer, D. R.↗

The application of the principles of invariance to the radiative transfer equation in plant canopies

Solutions of the radiative transfer equation describing photon interactions with vegetation canopies are important in remote sensing since they provide the canopy reflectance distribution required in the interpretation of satellite acquired information. The general one-dimensional two-angle transport problem for a finite copy of arbitrary leaf angle distribution is considered. Analytical solutions are obtained in terms of generalized Chandrasekhar's X- and Y-functions by invoking the principles of invariance. A critical step in the formulation involves the decomposition of the integral of the scattering phase function into a product of known functions of the incident and scattered photon directions. Several simplified cases previously considered in the literature are derived from the generalized solution. Various symmetries obeyed by the scattering operator and reciprocity relations are formally proved.

Ganapol, B. D.↗

An approach to verification and validation of a reliable multicasting protocol: Extended Abstract

This paper describes the process of implementing a complex communications protocol that provides reliable delivery of data in multicast-capable, packet-switching telecommunication networks. The protocol, called the Reliable Multicasting Protocol (RMP), was developed incrementally using a combination of formal and informal techniques in an attempt to ensure the correctness of its implementation. Our development process involved three concurrent activities: (1) the initial construction and incremental enhancement of a formal state model of the protocol machine; (2) the initial coding and incremental enhancement of the implementation; and (3) model-based testing of iterative implementations of the protocol. These activities were carried out by two separate teams: a design team and a V&V team. The design team built the first version of RMP with limited functionality to handle only nominal requirements of data delivery. This initial version did not handle off-nominal cases such as network partitions or site failures. Meanwhile, the V&V team concurrently developed a formal model of the requirements using a variant of SCR-based state tables. Based on these requirements tables, the V&V team developed test cases to exercise the implementation. In a series of iterative steps, the design team added new functionality to the implementation while the V&V team kept the state model in fidelity with the implementation. This was done by generating test cases based on suspected errant or off-nominal behaviors predicted by the current model. If the execution of a test in the model and implementation agreed, then the test either found a potential problem or verified a required behavior. However, if the execution of a test was different in the model and implementation, then the differences helped identify inconsistencies between the model and implementation. In either case, the dialogue between both teams drove the co-evolution of the model and implementation. We have found that this interactive, iterative approach to development allows software designers to focus on delivery of nominal functionality while the V&V team can focus on analysis of off nominal cases. Testing serves as the vehicle for keeping the model and implementation in fidelity with each other. This paper describes (1) our experiences in developing our process model; and (2) three example problems found during the development of RMP. Although RMP has provided our research effort with a rich set of test cases, it also has practical applications within NASA. For example, RMP is being considered for use in the NASA EOSDIS project due to its significant performance benefits in applications that need to replicate large amounts of data to many network sites.

Callahan, John R.↗

Computations of internal turbulent flow with large separated flow regions

An implicit two-equation turbulence solver in generalized coordinates has been developed and is used in conjunction with the three-dimensional incompressible Navier-Stokes solver, INS3D, to calculate the internal flow in one channel and in an additional channel with a sudden 2:3 expansion. A new and consistent boundary procedure for a low Reynolds number form of the kappa-epsilon turbulence model is chosen to integrate the equations up to the wall. The high Reynolds number form of the equations are integrated using wall functions. The latter approach yields a faster convergence to the steady state solution than the former. For the case of channel flow, both the wall function and wall boundary condition approaches yield results in good agreement with the experimental data. the back-step (sudden expansion) flow is calculated using wall function approach. The predictions are in reasonable agreement with the experimental data.

Kaul, U. K.↗

Enhancing nanoscale charged colloid crystallization near a metastable liquid binodal

Achieving predictive control over crystallization using non-classical nucleation while avoiding kinetic traps would be a step towards designing materials with new functionalities. We address these challenges by inducing the bottom-up assembly of nanocrystals into ordered arrays, or superlattices. Using electrostatics—rather than density—to tune the interactions between particles, we watch self-assembly proceed through a metastable liquid phase. Here, we systematically investigate the phase behaviour as a function of quench conditions in situ and in real time using small-angle X-ray scattering. By fitting to colloid, liquid and superlattice models, we extract the time evolution of each phase and the system phase diagram, which we find to be consistent with short-range attractive interactions. Using the predictive power of the phase diagram, we establish control of the self-assembly rate over three orders of magnitude, and we identify one- and two-step self-assembly regimes, with only the latter implicating the metastable liquid as an intermediate. The presence of the metastable liquid increases the superlattice formation rate relative to the equivalent one-step pathway, and the superlattice order increases with the rate, revealing a generalizable kinetic strategy for promoting and enhancing ordered assembly.

Tanner, Christian P. N. [University of California,↗

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

Upcycling Polyethylene Waste into Hybrid Graphitic Porous Carbon Materials Used in High‐Performance Zinc‐Ion Hybrid Capacitors

Polyethylene (PE) waste is a challenge to upcycle into useful materials because this plastic tends to decompose into volatile compounds when heated at relatively low temperatures. In this work, mixtures of PE wastes into a hybrid graphitic porous carbon (HGPC) by a thermal oxidation pretreatment step, with assistance of an inert solid additive (KCl), to functionalize, crosslink, and stabilize the PE waste followed by carbonization and catalytic graphitization steps with a potassium carbonate catalyst, are upcycled. The PE waste‐derived HGPC (PW‐HGPC) has a hybrid structure composed of graphene‐like carbon nanosheets grown on the surface of carbon particles, high porosity with specific surface area, up to 1,763 m 2 g −1 , and good graphitic degree with average Raman I 2D / I G ratios of 0.53. When used as cathode material for zinc‐ion hybrid capacitors, this PW‐HGPC exhibits an excellent specific capacity, up to 126.7 mAh g −1 , at high mass loading of 10 mg cm −2 . Moreover, PW‐HGPC exhibits remarkable cycling stability with capacity retention of >94% after 10 000 cycles. Additionally, the KCl is recycled and reused over five times. This method provides a new solution for upcycling PE wastes into high value‐added carbon materials, not only for zinc‐ion hybrid capacitors but also for other electrochemical energy storage device applications.

hybrid graphitic porous carbon↗

Constraints on triple oxygen isotope kinetics

Triple oxygen isotope measurements have been used as tracers for classifying meteorites, as indicators of atmospheric chemistry and as a proxy for primary productivity among other uses. These utilities rely on the mass dependency of most chemical reactions: δ17O is proportional to δ18O by roughly a factor of 0.5. However,in detail it is known that mass-dependency factorsvary. Constraints on mass dependency are well established for equilibrium isotope effects, but similar constraints for kinetic isotope effects are lacking.Here, we provide constraints for isotope effects generated by unidirectional chemical reactions (i.e. kinetic isotope effects, KIEs). A simplified Monte Carlo model is used to estimatethe range of kinetic isotope effects. It is found that the range of KIEs exceeds the range of equilibrium isotope effects such that some mass dependent kinetic isotope effects could be mistaken for being mass-independent. For the16O-17O-18O system,Δ′17Ovariations as large 0.8‰ can be generated from a single reaction step. Two case studies from published literature, on the thermal decomposition of calcite (CaCO3) and brucite (Mg(OH)2), are used to test the extreme limits of the modeled results. Isotope fractionationsgenerated fromdensity functional theory models of the surface reaction steps are shownto match the published data. The fractionation models are strictly mass dependent thus this result suggests thermal decomposition of calcite is also mass dependent. Although oxygen is the focus here, our approach applies to any system with more than two isotopes.

J. A. Hayles↗

Dynamic remapping of parallel computations with varying resource demands

A large class of computational problems is characterized by frequent synchronization, and computational requirements which change as a function of time. When such a problem must be solved on a message passing multiprocessor machine, the combination of these characteristics lead to system performance which decreases in time. Performance can be improved with periodic redistribution of computational load; however, redistribution can exact a sometimes large delay cost. We study the issue of deciding when to invoke a global load remapping mechanism. Such a decision policy must effectively weigh the costs of remapping against the performance benefits. We treat this problem by constructing two analytic models which exhibit stochastically decreasing performance. One model is quite tractable; we are able to describe the optimal remapping algorithm, and the optimal decision policy governing when to invoke that algorithm. However, computational complexity prohibits the use of the optimal remapping decision policy. We then study the performance of a general remapping policy on both analytic models. This policy attempts to minimize a statistic W(n) which measures the system degradation (including the cost of remapping) per computation step over a period of n steps. We show that as a function of time, the expected value of W(n) has at most one minimum, and that when this minimum exists it defines the optimal fixed-interval remapping policy. Our decision policy appeals to this result by remapping when it estimates that W(n) is minimized. Our performance data suggests that this policy effectively finds the natural frequency of remapping. We also use the analytic models to express the relationship between performance and remapping cost, number of processors, and the computation's stochastic activity.

Nicol, D. M.↗