Sensitivity-function calculation in linear systems using time-shared analog integration.
Sensitivity functions relating parameter variations effect to system variables calculated for linear systems using time shared analog integration
Kahne, S. J.↗
SEARCH · Search NASA
Results for “time integration”
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.
Spectral Element Method for the Simulation of Unsteady Compressible Flows
This work uses a discontinuous-Galerkin spectral-element method (DGSEM) to solve the compressible Navier-Stokes equations [1{3]. The inviscid ux is computed using the approximate Riemann solver of Roe [4]. The viscous fluxes are computed using the second form of Bassi and Rebay (BR2) [5] in a manner consistent with the spectral-element approximation. The method of lines with the classical 4th-order explicit Runge-Kutta scheme is used for time integration. Results for polynomial orders up to p = 15 (16th order) are presented. The code is parallelized using the Message Passing Interface (MPI). The computations presented in this work are performed using the Sandy Bridge nodes of the NASA Pleiades supercomputer at NASA Ames Research Center. Each Sandy Bridge node consists of 2 eight-core Intel Xeon E5-2670 processors with a clock speed of 2.6Ghz and 2GB per core memory. On a Sandy Bridge node the Tau Benchmark [6] runs in a time of 7.6s.
Diosady, Laslo Tibor↗
Deep Spectroscopy with DESI for Photometric Redshift Training and Calibration
Deep spectroscopic samples can improve photometric redshift (photo-z) estimates and reduce uncertainties on redshift distributions. Such improvements can increase the cosmological constraining power of large imaging-based experiments such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) and mitigate what may be a limiting systematic effect. We present results from the “DESI-Deep pilot” program, which was designed to assess the capability of the Dark Energy Spectroscopic Instrument (DESI) on the 4m Mayall telescope to measure redshifts of galaxies as faint as expected lensing samples for early LSST data (m i ≤ 24.5). We find that DESI is remarkably efficient at this task, with redshift success rates comparable to the results of observations from 10 m class telescopes with only ∼2 × longer integration time (rather than ∼8 × longer as would be expected from aperture-area scaling), while simultaneously achieving ∼30 times larger multiplexing. We also find that the signal-to-noise ratio of the spectra scales as expected for background-limited observations even for the longest exposure times (∼7 hr) and faintest targets in the program. These results demonstrate that DESI could provide the definitive redshift sample for the early years of LSST with a modest investment of observing time. Based upon the results of this program, we provide updated predictions for the time required to collect benchmark samples for photo-z training and calibration using a variety of spectroscopic facilities. Finally, we describe a potential “DESI-Deep” survey designed to train and calibrate photo-z’s for imaging experiments, and provide forecasts of its impact on cosmological inference.
Dey, Biprateep [University of Toronto; University ↗
Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems
The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.
Advection-dominated dynamical systems↗
Indium Single-Ion Frequency Standard
A single laser-cooled indium ion is a promising candidate for an ultimate resolution optical time or frequency standard. It can be shown that single ions from group IIIA of the periodic table (indium, thallium, etc.) can have extremely small systematic errors. In addition to being free from Doppler, transit-time and collisional shifts, these ions are also quite insensitive to perturbations from ambient magnetic and electric fields (mainly due to the use of a J=0-0 transition for spectroscopy). Of all group IIIA ions, indium seems to be the most practical, since it is heavy enough to have a tolerable intercombination cooling transition rate and (unlike thallium) has transitions which are easily accessible with frequency multiplied continuous-wave lasers. A single indium ion standard has a potential inaccuracy of one part in 10(exp 18) for integration times of 10(exp 6) seconds. We have made substantial progress during the grant period in constructing a frequency standard based upon a single indium ion. At the beginning of the grant period, single indium ions were being successfully trapped, but the lasers and optical systems were inadequate to achieve the desired goal. We have considerably improved the stability of the dye laser used to cool the ions and locked it to a molecular resonance line, making it possible to observe stable cooling-line fluorescence from a single indium ion for reasonable periods of time, as required by the demands of precision spectroscopy. We have substantially improved the single-ion fluorescence signal with significant benefits for the detection efficiency of forbidden transitions using the 'shelving' technique. Finally, we have constructed a compact, efficient UV 'clock' laser and observed 'clock' transitions in single indium ions using this laser system. We will elaborate on these accomplishments.
Nagourney, Warren↗
Signal-to-noise ratio for the wide field-planetary camera of the Space Telescope
Signal-to-noise ratios for the Wide Field Camera and Planetary Camera of the Space Telescope were calculated as a function of integration time. Models of the optical systems and CCD detector arrays were used with a 27th visual magnitude point source and a 25th visual magnitude per arc-sq. second extended source. A 23rd visual magnitude per arc-sq. second background was assumed. The models predicted signal-to-noise ratios of 10 within 4 hours for the point source centered on a signal pixel. Signal-to-noise ratios approaching 10 are estimated for approximately 0.25 x 0.25 arc-second areas within the extended source after 10 hours integration.
Zissa, D. E.↗
Efficient lattice QCD computation of radiative-leptonic-decay form factors at multiple positive and negative photon virtualities
In previous work [D. Giusti, Methods for high-precision determinations of radiative-leptonic decay form factors using lattice QCD, Phys. Rev. D 107, 074507 (2023)], we showed that form factors for radiative leptonic decays of pseudoscalar mesons can be determined efficiently and with high precision from lattice QCD using the “three-dimensional (3D) method,” in which three-point functions are computed for all values of the current insertion time and the time integral is performed at the data-analysis stage. Here, we demonstrate another benefit of the 3D method: the form factors can be extracted for any number of nonzero photon virtualities from the same three-point functions at no extra cost. We present results for the $D_s → ℓνγ*$ vector form factor as a function of photon energy and photon virtuality, for both positive and negative virtuality, for a single ensemble with 340 MeV pion mass and 0.11 fm lattice spacing. In our analysis, we separately consider the two different time orderings and the different quark flavors in the electromagnetic current. We discuss in detail the behavior of the unwanted exponentials contributing to the three-point functions, as well as the choice of fit models and fit ranges used to remove them for various values of the virtuality. While positive photon virtuality is relevant for decays to multiple charged leptons, negative photon virtuality suppresses soft contributions and is of interest in QCD-factorization studies of the form factors.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Stormtime transport of ring current and radiation belt ions
This is an investigation of stormtime particle transport that leads to formation of the ring current. Our method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants mu) in response to model substorm-associated impulses in the convection electric field. We compare our simulation results qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For mu approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of our model storm, and we find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For mu approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For mu approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in this case) 12 representative ions equally spaced in drift time around the steady-state drift shell of interest (L equivalent to 3), we have estimated (from both our forward and our time-reversed simulations) the time-integrated radial-diffusion coefficients D(sup sim)(sub LL) for particles having selected values of mu approximately greater than 15 MeV/G. The results agree surprisingly well with the predictions (D(sup ql)(sub LL)) of quasilinear radial diffusion theory, despite the rather brief duration (approximately 3 hrs) of our model storm and despite the extreme variability (with frequency) of the spectral-density function that characterizes the applied electric field during our model storm. As expected, the values of D(sup sim)(sub LL) deduced (respectively) from our forward and time-reversed simulations agree even better with each other and with D(sup sim)(sub LL) when the impulse amplitudes which characterize the individual substorms of our model storm are systematically reduced.
Chen, Margaret W.↗
Integrity of the Plasma Magnetic Nozzle
This report examines the physics governing certain aspects of plasma propellant flow through a magnetic nozzle, specifically the integrity of the interface between the plasma and the nozzle s magnetic field. The injection of 100s of eV plasma into a magnetic flux nozzle that converts thermal energy into directed thrust is fundamental to enabling 10 000s of seconds specific impulse and 10s of kW/kg specific power piloted interplanetary propulsion. An expression for the initial thickness of the interface is derived and found to be approx.10(exp -2) m. An algorithm is reviewed and applied to compare classical resistivity to gradient-driven microturbulent (anomalous) resistivity, in terms of the spatial rate and time integral of resistive interface broadening, which can then be related to the geometry of the nozzle. An algorithm characterizing plasma temperature, density, and velocity dependencies is derived and found to be comparable to classical resistivity at local plasma temperatures of approx. 200 eV. Macroscopic flute-mode instabilities in regions of "adverse magnetic curvature" are discussed; a growth rate formula is derived and found to be one to two e-foldings of the most unstable Rayleigh-Taylor (RT) mode. After establishing the necessity of incorporating the Hall effect into Ohm s law (allowing full Hall current to flow and concomitant plasma rotation), a critical nozzle length expression is derived in which the interface thickness is limited to about 1 ion gyroradius.
Gerwin, Richard A.↗
Stability of Supersonic Boundary Layers Over Blunt Wedges
Receptivity and stability of supersonic boundary layers over blunt flat plates and wedges are numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 10(exp 6)/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. Computations are performed for a flat plate with leading edge thicknesses of 0.0001, 0.001, 0.005 and 0.01 inches that give Reynolds numbers based on the leading edge thickness ranging from 1000 to 10000. Calculations are also performed for a wedge of 10 degrees half angle with different leading edge radii 0.001 and 0.01 inches. The linear stability results showed that the bluntness has a strong stabilizing effect on the stability of two-dimensional boundary layers. The transition Reynolds number for a flat plate with a leading edge thickness of 0.01 inches is about 3.5 times larger than it is for the Blasius boundary layer. It was also revealed that boundary layers on blunt wedges are far more stable than on blunt flat plates.
Balakumar, Ponnampalam↗
Modular approach for satellite communication ground terminals
The trend in satellite communications is toward completely digital, time division multiple access (TDMA) systems with uplink and downlink data rates dictated by the type of service offered. Trunking terminals will operate in the 550 MBPS (megabit per second) region uplink and downlink, whereas customer premise service (CPS) terminals will operate in the 25 to 10 MBPS region uplink and in the 200 MBPS region downlink. Additional criteria for the ground terminals will be to maintain clock sychronization with the system and burst time integrity to within a matter of nanoseconds, to process required order-fire information, to provide adaptive data scrambing, and to compensate for variations in the user input output data rates, and for changes in range in the satellite communications links resulting from satellite perturbations in orbit. To achieve the required adaptability of a ground terminal to the above mentioned variables, programmable building blocks can be developed that will meet all of these requirements. To maintain system synchronization, i.e., all bursted data arriving at the satellite within assigned TDMA windows, ground terminal transmit data rates and burst timing must be maintained within tight tolerances. With a programmable synchronizer as the heart of the terminal timing generation, variable data rates and burst timing tolerances are achievable. In essence, the unit inputs microprocessor generated timing words and outputs discrete timing pulses.
Gould, G. R.↗
Sensitivity of the terrestrial planet finder
A key long-term goal of NASA's Origins program is the detection and characterization of habitable planets orbiting stars within the solar neighborhood. A cold, space-borne interferometer operating in the mid-infrared with a approx. 75 m baseline can null the light of a parent star and detect the million-times fainter radiation from an Earth-like planet located in the "habitable zone" around stars as far as 15 pc away. Such an interferometer, designated the Terrestrial Planet Finder (TPF) by NASA, could even detect atmospheric signatures of species such as CO2, O3, and H2O indicative of either the possibility or presence of primitive life. This talk highlights some of the sensitivity issues affecting the detectability of terrestrial planets. Sensitivity calculations show that a system consisting of 2 m apertures operating at 5 AU or 4 m apertures operating at 1 AU can detect terrestrial planets in reasonable integration times for levels of exo-zodiacal emission up to 10 times that seen in our solar system (hereafter denoted as 10xSS). Additionally, simulations show that confusion noise from structures in the exo-zodiacal cloud should not impede planet detection until the exo-zodiacal emission reaches the 10xSS level.
Beichman, Charles↗
Observability of the neutrino flux from the inner region of the galactic disk
The observability of galactic neutrinos in a detector of 10 to the 10th tons of water with an observing time of a few years is explored. Although the atmospheric flux exceeds the galactic flux considerably at energies greater than or equal to 1 TeV, the latter may still provide a marginally observable signal owing to its directionality. Galactic muon neutrinos with energy greater than approximately 1 TeV will produce a signal approximately two standard deviations above the atmospheric background over a four year period. If electron neutrinos can also be studied with DUMAND, then galactic electron neutrinos above 1 TeV would give an approximately four to five standard deviations signal above the electron neutrino background over a four-year integration time.
Silberberg, R.↗
RHESSI Investigations of the Neupert Effect in Solar Flares
In many solar flares the time integral of the hard X-ray emission closely matches the temporal variation of the soft X-ray emission. This is known as the Neupert Effect and is believed to result when the same accelerated electrons that produce the bremsstrahlung hard X-rays also heat the plasma that produces the thermal bremsstrahlung soft X-rays. The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) is uniquely capable of investigating this effect since is covers the soft (down to 3 keV) and hard X-ray energy ranges with keV energy resolution, arcsecond-class angular resolution, and sub-second time resolution. Several GOES M-class flares have already been detected that show this general effect, but the detailed time histories suggest that there must be other heating mechanisms besides the thermalization of accelerated electrons that are heating the soft X-ray emitting plasma, even during the impulsive phase. Spatially-resolved hard and soft X-ray observations of these flares are also being investigated to determine what fraction of the heating is produced by the accelerated electrons.
Dennis, Brian R.↗
The Search for Habitable Worlds. 1. The Viability of a Starshade Mission
As part of NASA's mission to explore habitable planets orbiting nearby stars, this article explores the detection and characterization capabilities of a 4 m space telescope plus 50 m starshade located at the Earth-Sun L2 point, known as the New Worlds Observer (NWO). Our calculations include the true spectral types and distribution of stars on the sky, an iterative target selection protocol designed to maximize efficiency based on prior detections, and realistic mission constraints. We conduct simulated observing runs for a wide range in exozodiacal background levels (epsilon = 1-100 times the local zodi brightness) and overall prevalence of Earth-like terrestrial planets (eta(sub solar halo))0.1-1). We find that even without any return visits, the NWO baseline architecture (IWA = 65 mas, limiting FPB = 4 x 10(exp -11) can achieve a 95% probability of detecting and spectrally characterizing at least one habitable Earth-like planet and an expectation value of approximately 3 planets found, within the mission lifetime and delta V budgets, even in the worst-case scenario (eta(sub solar halo) = 0.1 and = epsilon = 100 zodis for every target). This achievement requires about 1 yr of integration time spread over the 5 yr mission, leaving the remainder of the telescope time for UV-NIR general astrophysics. Cost and technical feasibility considerations point to a "sweet spot" in starshade design near a 50 m starshade effective diameter. with 12 or 16 petals, at a distance of 70,000-100,000 km from the telescope.
Turnbull, Margaret C.↗
A fast implicit upwind solution algorithm for three-dimensional unstructured dynamic meshes
A fast implicit upwind algorithm for the solution of the time-dependent Euler equations is presented for aerodynamic analysis involving unstructured dynamic meshes. The spatial discretization of the scheme is based on the upwind approach of Roe, referred to as flux-difference splitting (FDS). The FDS approach is naturally dissipative and captures shock waves and contact discontinuities sharply. The temporal discretization of the scheme involves an implicit time-integration using a two-sweep Gauss-Seidel relaxation procedure. The procedure is computationally efficient for either steady or unsteady flow problems. A detailed description is given of the implicit upwind solution algorithm along with results which assess the capability. The results are presented for the NACA 0012 airfoil and for the Boeing 747 aircraft. The 747 geometry includes the fuselage, wing, horizontal and vertical tails, under-wing pylons, and flow-through engine nacelles. Euler solutions for the 747 aircraft on an unstructured tetrahedral mesh containing approximately 100,000 cells were obtained to engineering accuracy in less than one hour CPU time on a Cray-2 computer.
Batina, John T.↗
A fast implicit upwind solution algorithm for three-dimensional unstructured dynamic meshes
A fast implicit upwind algorithm for the solution of the time-dependent Euler equations is presented for aerodynamic analysis involving unstructured dynamic meshes. The spatial discretization of the scheme is based on the upwind approach of Roe, referred to as flux-difference splitting (FDS). The FDS approach is naturally dissipative and captures shock waves and contact discontinuities sharply. The temporal discretization of the scheme involves an implicit time-integration using a two-sweep Gauss-Seidel relaxation procedure. The procedure is computationally efficient for either steady or unsteady flow problems. A detailed description is given of the implicit upwind solution algorithm along with results which assess the capability. The results are presented for the NACA 0012 airfoil and for the Boeing 747 aircraft. The 747 geometry includes the fuselage, wing, horizontal and vertical tails, under-wing pylons, and flow-through engine nacelles. Euler solutions for the 747 aircraft on an unstructured tetrahedral mesh containing approximately 100,000 cells were obtained to engineering accuracy in less than one hour CPU time on a Cray-2 computer.
Batina, John T.↗
Observability of the neutrino flux from the inner region of the galactic disk
The observability of galactic neutrinos in a detector of 10 billion tons of water with an observing time of a few years is explored. Although the atmospheric flux exceeds the galactic flux considerably at energies greater than or equal to 1 TeV, the latter may still provide a marginally observable signal owing to its directionality. Galactic muon neutrinos with energy greater than or equal to 1 TeV will produce a signal approximately 2 sigma above the atmospheric background over a four year period. If electron neutrinos can also be studied with the deep underwater muon and neutrino detector, then galactic electron neutrinos above 1 TeV would give an approximate 4 to 5 sigma signal above the electron neutrino background over a four year integration time.
Silberberg, R.↗