Search NASA⌕ Search

SEARCH · Search NASA

Results for “divided differences”

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 37 records · Page 2

An Automated Approach to Very High Order Aeroacoustic Computations in Complex Geometries

Computational aeroacoustics requires efficient, high-resolution simulation tools. And for smooth problems, this is best accomplished with very high order in space and time methods on small stencils. But the complexity of highly accurate numerical methods can inhibit their practical application, especially in irregular geometries. This complexity is reduced by using a special form of Hermite divided-difference spatial interpolation on Cartesian grids, and a Cauchy-Kowalewslci recursion procedure for time advancement. In addition, a stencil constraint tree reduces the complexity of interpolating grid points that are located near wall boundaries. These procedures are used to automatically develop and implement very high order methods (>15) for solving the linearized Euler equations that can achieve less than one grid point per wavelength resolution away from boundaries by including spatial derivatives of the primitive variables at each grid point. The accuracy of stable surface treatments is currently limited to 11th order for grid aligned boundaries and to 2nd order for irregular boundaries.

Dyson, Rodger W.↗

Automated Approach to Very High-Order Aeroacoustic Computations

Computational aeroacoustics requires efficient, high-resolution simulation tools. For smooth problems, this is best accomplished with very high-order in space and time methods on small stencils. However, the complexity of highly accurate numerical methods can inhibit their practical application, especially in irregular geometries. This complexity is reduced by using a special form of Hermite divided-difference spatial interpolation on Cartesian grids, and a Cauchy-Kowalewski recursion procedure for time advancement. In addition, a stencil constraint tree reduces the complexity of interpolating grid points that am located near wall boundaries. These procedures are used to develop automatically and to implement very high-order methods (> 15) for solving the linearized Euler equations that can achieve less than one grid point per wavelength resolution away from boundaries by including spatial derivatives of the primitive variables at each grid point. The accuracy of stable surface treatments is currently limited to 11th order for grid aligned boundaries and to 2nd order for irregular boundaries.

Dyson, Rodger W.↗

Ozone Profiles in the High-latitude Stratosphere and Lower Mesosphere Measured by the Improved Limb Atmospheric Spectrometer (ILAS)-II: Comparison with other Satellite Sensors and Ozonesondes

A solar occultation sensor, the Improved Limb Atmospheric Spectrometer (ILAS)-II, measured 5890 vertical profiles of ozone concentrations in the stratosphere and lower mesosphere and of other species from January to October 2003. The measurement latitude coverage was 54-71degN and 64-88degS, which is similar to the coverage of ILAS (November 1996 to June 1997). One purpose of the ILAS-II measurements was to continue such high-latitude measurements of ozone and its related chemical species in order to help accurately determine their trends. The present paper assesses the quality of ozone data in the version 1.4 retrieval algorithm, through comparisons with results obtained from comprehensive ozonesonde measurements and four satellite-borne solar occultation sensors. In the Northern Hemisphere (NH), the ILAS-II ozone data agree with the other data within +/-10% (in terms of the absolute difference divided by its mean value) at altitudes between 11 and 40 km, with the median coincident ILAS-II profiles being systematically up to 10% higher below 20 km and up to 10% lower between 21 and 40 km after screening possible suspicious retrievals. Above 41 km, the negative bias between the NH ILAS-II ozone data and the other data increases with increasing altitude and reaches 30% at 61-65 km. In the Southern Hemisphere, the ILAS-II ozone data agree with the other data within 10% in the altitude range of 11-60 km, with the median coincident profiles being on average up to 10% higher below 20 km and up to 10% lower above 20 km.

Sugita, T.↗

Comparison of Sigma-Point and Extended Kalman Filters on a Realistic Orbit Determination Scenario

Sigma-point filters have received a lot of attention in recent years as a better alternative to extended Kalman filters for highly nonlinear problems. In this paper, we compare the performance of the additive divided difference sigma-point filter to the extended Kalman filter when applied to orbit determination of a realistic operational scenario based on the Interstellar Boundary Explorer mission. For the scenario studied, both filters provided equivalent results. The performance of each is discussed in detail.

Gaebler, John↗

A Walsh Function Module Users' Manual

The solution of partial differential equations (PDEs) with Walsh functions offers new opportunities to simulate many challenging problems in mathematical physics. The approach was developed to better simulate hypersonic flows with shocks on unstructured grids. It is unique in that integrals and derivatives are computed using simple matrix multiplication of series representations of functions without the need for divided differences. The product of any two Walsh functions is another Walsh function - a feature that radically changes an algorithm for solving PDEs. A FORTRAN module for supporting Walsh function simulations is documented. A FORTRAN code is also documented with options for solving time-dependent problems: an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the usage of the Walsh function module including such features as operator overloading, Fast Walsh Transforms in multi-dimensions, and a Fast Walsh reciprocal.

Gnoffo, Peter A.↗

Conjunction Assessment Late-Notice High-Interest Event Investigation: Space Weather Aspects

Late-notice events usually driven by large changes in primary (protected) object or secondary object state. Main parameter to represent size of state change is component position difference divided by associated standard deviation (epsilon divided by sigma) from covariance. Investigation determined actual frequency of large state changes, in both individual and combined states. Compared them to theoretically expected frequencies. Found that large changes ( (epsilon divided by sigma) is greater than 3) in individual object states occur much more frequently than theory dictates. Effect is less pronounced in radial components and in events with probability of collision (Pc) greater than 1 (sup -5) (1e-5). Found combined state matched much closer to theoretical expectation, especially for radial and cross-track. In-track is expected to be the most vulnerable to modeling errors, so not surprising that non-compliance largest in this component.

Conjunction↗

In situ mechanical characterization of a single cell-cell adhesion interface under large strain

A method of measuring a stress-strain curve in a cell-cell adhesion interface, the method including: providing a structure including a first movable island supported by a first beam, a second movable island supported by a second beam, and a gap therebetween connected by a pair of cells forming a junction, the pair of cells comprising a cell-cell adhesion interface having an initial length defined by a distance between nuclei of the pair of cells; moving the second movable island with a defined displacement; determining a displacement of the first movable island based on moving the second movable island; calculating a difference between the displacement of the first movable island and the defined displacement of the second movable island based on moving the second movable island; determining an applied strain in the cell-cell adhesion interface between the pair of cells based on the difference divided by the initial length of the cell-cell adhesion interface; calculating a force between the cell-cell adhesion interface of the pair of cells based on the displacement of the first movable island; calculating a stress in the cell-cell adhesion interface between the pair of cells based on the force; and determining the stress-strain curve of the cell-cell adhesion interface between the pair of cells by plotting the calculated stress against the applied strain.

Yang, Ruiguo↗

Apparatus and Method to Enable Precision and Fast Laser Frequency Tuning

An apparatus and method is provided to enable precision and fast laser frequency tuning. For instance, a fast tunable slave laser may be dynamically offset-locked to a reference laser line using an optical phase-locked loop. The slave laser is heterodyned against a reference laser line to generate a beatnote that is subsequently frequency divided. The phase difference between the divided beatnote and a reference signal may be detected to generate an error signal proportional to the phase difference. The error signal is converted into appropriate feedback signals to phase lock the divided beatnote to the reference signal. The slave laser frequency target may be rapidly changed based on a combination of a dynamically changing frequency of the reference signal, the frequency dividing factor, and an effective polarity of the error signal. Feed-forward signals may be generated to accelerate the slave laser frequency switching through laser tuning ports.

Chen, Jeffrey R.↗

“It’s Raining Bits”: Patterns in Directional Precipitation Persistence Across the United States

The spatial and temporal ordering of precipitation occurrence impacts ecosystems, streamflow, and water availability. For example, both large-scale climate patterns and local landscapes drive weather events, and the typical speeds and directions of these events moving across a basin dictate the timing of flows at its outlet. We address the predictability of precipitation occurrence at a given location, based on the knowledge of past precipitation at surrounding locations. We identify ‘‘dominant directions of precipitation influence’’ across the continental United States based on a gridded daily dataset. Specifically, we apply information theory–based measures that characterize dominant directions and strengths of spatial and temporal precipitation dependencies. On a national average, this dominant direction agrees with the prevalent direction of weather movement from west to east across the country, but regional differences reflect topographic divides, precipitation gradients, and different climatic drivers of precipitation. Trends in these information relationships and their correlations with climate indices over the past 70 years also show seasonal and spatial divides. This study expands upon a framework of information-based predictability to answer questions about spatial connectivity in addition to temporal persistence. The methods presented here are generally useful to understand many aspects of weather and climate variability.

North America↗

Optoelectronic System Measures Distances to Multiple Targets

An optoelectronic metrology apparatus now at the laboratory-prototype stage of development is intended to repeatedly determine distances of as much as several hundred meters, at submillimeter accuracy, to multiple targets in rapid succession. The underlying concept of optoelectronic apparatuses that can measure distances to targets is not new; such apparatuses are commonly used in general surveying and machining. However, until now such apparatuses have been, variously, constrained to (1) a single target or (2) multiple targets with a low update rate and a requirement for some a priori knowledge of target geometry. When fully developed, the present apparatus would enable measurement of distances to more than 50 targets at an update rate greater than 10 Hz, without a requirement for a priori knowledge of target geometry. The apparatus (see figure) includes a laser ranging unit (LRU) that includes an electronic camera (photo receiver), the field of view of which contains all relevant targets. Each target, mounted at a fiducial position on an object of interest, consists of a small lens at the output end of an optical fiber that extends from the object of interest back to the LRU. For each target and its optical fiber, there is a dedicated laser that is used to illuminate the target via the optical fiber. The targets are illuminated, one at a time, with laser light that is modulated at a frequency of 10.01 MHz. The modulated laser light is emitted by the target, from where it returns to the camera (photodetector), where it is detected. Both the outgoing and incoming 10.01-MHz laser signals are mixed with a 10-MHz local-oscillator to obtain beat notes at 10 kHz, and the difference between the phases of the beat notes is measured by a phase meter. This phase difference serves as a measure of the total length of the path traveled by light going out through the optical fiber and returning to the camera (photodetector) through free space. Because the portion of the path length inside the optical fiber is not ordinarily known and can change with temperature, it is also necessary to measure the phase difference associated with this portion and subtract it from the aforementioned overall phase difference to obtain the phase difference proportional to only the free-space path length, which is the distance that one seeks to measure. Therefore, the apparatus includes a photodiode and a circulator that enable measurement of the phase difference associated with propagation from the LRU inside the fiber to the target, reflection from the fiber end, and propagation back inside the fiber to the LRU. Because this phase difference represents twice the optical path length of the fiber, this phase difference is divided in two before subtraction from the aforementioned total-path-length phase difference. Radiation-induced changes in the photodetectors in this apparatus can affect the measurements. To enable calibration for the purpose of compensation for these changes, the apparatus includes an additional target at a known short distance, located inside the camera. If the measured distance to this target changes, then the change is applied to the other targets.

Liebe, Carl Christian↗

Characterization of Transport Errors in Chemical Forecasts from a Global Tropospheric Chemical Transport Model

We propose a new methodology to characterize errors in the representation of transport processes in chemical transport models. We constrain the evaluation of a global three-dimensional chemical transport model (GEOS-CHEM) with an extended dataset of carbon monoxide (CO) concentrations obtained during the Transport and Chemical Evolution over the Pacific (TRACE-P) aircraft campaign. The TRACEP mission took place over the western Pacific, a region frequently impacted by continental outflow associated with different synoptic-scale weather systems (such as cold fronts) and deep convection, and thus provides a valuable dataset. for our analysis. Model simulations using both forecast and assimilated meteorology are examined. Background CO concentrations are computed as a function of latitude and altitude and subsequently subtracted from both the observed and the model datasets to focus on the ability of the model to simulate variability on a synoptic scale. Different sampling strategies (i.e., spatial displacement and smoothing) are applied along the flight tracks to search for systematic model biases. Statistical quantities such as correlation coefficient and centered root-mean-square difference are computed between the simulated and the observed fields and are further inter-compared using Taylor diagrams. We find no systematic bias in the model for the TRACE-P region when we consider the entire dataset (i.e., from the surface to 12 km ). This result indicates that the transport error in our model is globally unbiased, which has important implications for using the model to conduct inverse modeling studies. Using the First-Look assimilated meteorology only provides little improvement of the correlation, in comparison with the forecast meteorology. These general statements can be refined when the entire dataset is divided into different vertical domains, i.e., the lower troposphere (less than 2 km), the middle troposphere (2-6 km), and the upper troposphere (greater than 6 km). The best agreement between the observations and the model is found in the lower and middle troposphere. Downward displacements in the lower troposphere provide a better fit with the observed value, which could indicate a problem in the representation of boundary layer height in the model. Significant improvement is also found for downward and southward displacements in the upper troposphere. There are several potential sources of errors in our simulation of the continental outflow in the upper troposphere which could lead to such biases, including the location and/or the strength of deep convective cells as well as that of wildfires in Southeast Asia.

Bey, I.↗

Bringing Planetary Science Mission Outreach to the Deaf and Blind Communities

Introduction: Technology for enhancing outreach, like 3D printing, and science communication products, such as videos and podcasts, can be utilized within the planetary science community, especially for the engagement and excitement of current or upcoming planetary exploration missions. However, these communication products can also be further enhanced for the benefit of the blind and deaf communities. While such products may already be readily available, such projects are not easily accessible to blind and/or deaf certified educators, which often rely on making their own resources or do not have the funds to provide such resources (e.g., cost of 3D printers or cost of braille books). The planetary science community can have better practices to reach these broader audiences. Best practices can include transcripts from podcasts, transcripts in videos, and large-font captions. Images on websites and social media accounts should also include alt-text descriptive captions. 3D printing can also enhance planetary science for the blind community, through tactile posters, maps, and pamphlets. Planetary data can also be augmented by providing different tactile geological maps (e.g., topography or various datasets), and audio-visual videos freely available for educators. Visual Engagement: Visual engagement consists of several avenues to consider, the three main themes includes: 1) swag; 2) videos; 3) interactive exploration. Swag can include the fun visual take-home materials, such as stickers, posters, bookmarks, etc. Videos can include educational-specific videos (available freely via YouTube or by other educational-specific streaming avenues, such as Nebula or Curiosity Stream), provided they have Closed Captioning (CC). The use of QR codes to such videos or websites can also benefit to being added on swag. Interactive exploration can also be sub-divided by different types of engagement. A popular and still fairly new technology for public engagement is the use of virtual reality (VR). While this has been mainly for martian and lunar surface exploration [1], the deaf communities can benefit from VR through a more extensive look at our solar system and beyond (for example, a VR experience of the flight path, or visual map of the heliosphere/dynamics of our Sun). Audio Engagement: Audio tools can also be a useful avenue of communication, especially for the blind communities. Audio archiving can certainly be transcripts from the video engagements, but also the use of podcasts can also be a benefit. Podcasting can take on two forms: 1) interview engagement; and 2) update engagement. For interviews, scientists can communicate with STEM-specific podcast platforms to make other listener-bases aware of what is going on with a specific mission. For update-type communication, missions may opt to have an archived podcast of news, updates, and the teams involved. The most important aspect of such podcasts would be for the need of complimentary transcripts (including descriptive transcripts if sounds are included), and the limited use of jargon. Research Engagement: There have been several examples of involving the blind and low-vision communities in citizen science, such as through the NASA Heliophysics division. Examples include the NASA PUNCH (Polarimeter to Unify the Corona and Heliosphere) mission led by the Southwest Research Institute [2], which include blind and visually-impaired citizens to assist in the Sun’s coronal rhythms.” Another example is the Eclipse Soundscapes: Citizen Science Project (ES:CSP), which documents observations of acoustical changes of nature and ecosystems during solar eclipse events [3]. Inclusivity: A major theme that is necessary for public engagement is inclusivity and the awareness of reaching broader audiences. Outreach to include hearing/seeing impaired communities are still lacking in the sciences. There are several opportunities that the planetary sciences could take. Other projects that have emerged from the space sciences include adding transcripts to visual engagement [4], and the use of 3D printing for the visually-impaired [5]. References: [1] Olgin, J. (2020) 51st LPSC, Abstract 2137. [2] https://scitechdaily.com/outreach-for-nasa-punch-mission-embraces-ancient-and-modern-sun-watching-theme/ [3] https://science.nasa.gov/science-activation-team/eclipse-soundscapes [4] NASA International Observe the Moon Night (Blind and Deaf Accessible), Youtube Video ( https://www.youtube.com/watch?v=neHCfg0S3-Q) [5] Richardson, J., et al. (2018) AGU Fall Meeting, Abstract ED23F-0963.

C J Ahrens↗

Carbon Transport and Storage Planning and Viability Support Tools

The EDX disCO2ver Carbon Transport and Storage Planning and Viability Support Tools are made up of the Carbon Storage Planning Inquiry Tool (CS PlanIT, Justman et al. 2024) and the Carbon Storage Technical Viability Approach Support Tool (CS TVA). Together, these tools support data access to support understanding data availability to support planning efforts for carbon transport and storage. The Carbon Storage Planning Inquiry Tool (CS PlanIT) is an online web mapping application designed to help users explore, query, and evaluate multiple data layers to support and accelerate carbon storage resource and feasibility assessments and planning efforts. CS PlanIT currently contains a range of datasets associated with geologic, technical, and infrastructure factors. The data sets can be filtered geographically for an area of interest to update statistics and charts within the dashboard. The dashboard is divided into different sections called widgets, relating to different steps in the carbon storage planning process. The resources in this submission include a link to PlanIT, as well as a data catalog and link to user documentation. The original citation for the CS PlanIT tool, which has now been integrated into the toolset here, was: - Devin Justman, Scott Pantaleone, Maneesh Sharma, Lucy Romeo, Paige Morkner, CS PlanIT (Carbon Storage Planning Inquiry Tool) , 6/28/2024, https://edx.netl.doe.gov/dataset/cs-planit-carbon-storage-planning-inquiry-tool, DOI: 10.18141/2377953 The Carbon Storage Technical Viability Approach Support (CS TVA) Tool displays spatial data availability for the many components of Geologic Carbon Storage (GCS) technical viability assessment (Creason al 2025). Identifying sites suitable for GCS requires evaluating the intersection of myriad factors, including reservoir conditions, subsurface and surface hazards, infrastructure requirements, and energy community metrics. The technical viability of a site can only be confirmed for instances where all these factors have data available, and where those data support viability. Additional Resources related to the Technical Viability Assessment Tool: - Julia Mulhern, Casey White, Araceli Lara, Neyda Cordero Rodriguez, Zachary Jackson, Jacob Shay, Gabriel Creason, MacKenzie Mark-Moser, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Database, 3/26/2025, https://edx.netl.doe.gov/dataset/edx4ccs-carbon-storage-technical-viability-approach-database , DOI:10.18141/1984655 - Julia Mulhern, MacKenzie Mark-Moser, Gabriel Creason, Casey White, Araceli Lara, Neyda Cordero Rodriguez, Zach Jackson, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Matrix, 3/27/2025, https://edx.netl.doe.gov/dataset/carbon-storage-technical-viability-approach-cs-tva-matrix , DOI: 10.18141/2539979 - Gabriel Creason, Zach Jackson, Neyda Cordero Rodriguez, Julia Mulhern, Casey White, Araceli Lara, MacKenzie Mark-Moser, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Data Availability Results Database, 3/27/2025, https://edx.netl.doe.gov/dataset/carbon-storage-technical-viability-approach-cs-tva-data-availability-results-database, DOI:10.18141/2538557

Carbon storage↗

Precision voltage division without precision parts

Active voltage divider applies different logic signals to field-effect transistor switches to generate various voltage levels at output. Circuit is compact and reduces errors from drifting component values.

Strange, M. G.↗

Automatic differentiation of advanced CFD codes for multidisciplinary design

Automated multidisciplinary design of aircraft and other flight vehicles requires the optimization of complex performance objectives with respect to a number of design parameters and constraints. The effect of these independent design variables on the system performance criteria can be quantified in terms of sensitivity derivatives which must be calculated and propagated by the individual discipline simulation codes. Typical advanced CFD analysis codes do not provide such derivatives as part of a flow solution; these derivatives are very expensive to obtain by divided (finite) differences from perturbed solutions. It is shown that sensitivity derivatives can be obtained accurately and efficiently using the ADIFOR source translator for automatic differentiation. In particular, it is demonstrated that the 3-D, thin-layer Navier-Stokes, multigrid flow solver called TLNS3D is amenable to automatic differentiation in the forward mode even with its implicit iterative solution algorithm and complex turbulence modeling. It is significant that by using computational differentiation, consistent discrete nongeometric sensitivity derivatives have been obtained from an aerodynamic 3-D CFD code in a relatively short time, e.g., O(man-week) not O(man-year).

Bischof, C.↗

Tracking capability for entry, descent and landing and its support to NASA Mars Exploration Rovers

This paper described a new capability recently deployed in the NASA Deep Space Network to support the entry, descent and landing of the two Mars Exploration Rovers - Spirit and Opportunity. This final segment of a 7-month journey was critical step in enabling successful mission operation. The EDL Data Analysis (EDA) equipment processed the open loop recording samples from the existing radio science receivers, extracted the MFSK-modulated tones, and provided information on the occurrence of associated events happened onboard spacecraft. The EDA software supported real-time as well as post pass processing. The design flexibility allowed the track to be divided into different segments, each could be processed with a configuration optimized for its signal conditions and uncertainties. The use of dynamic lock state concept enabled signal detection be done with minimum required computing power - less was spent when signal was steadily detected and more when signal fluctuated. With the server-client design, the EDA supported connections to different remote operation centers for monitor purpose. The paper also presented the field results from both Spirit and Opportunity landings. The detection were much better than anticipated. All critical event tones were seen in real time. Surprising detection of LCP signal was found with Spirit landing, resulted in a reconfiguration for Opportunity landing to maximize signal detection.

tracking↗

Black Hole Spectral States and Physical Connections

The dramatic changes seen in the X-ray spectral and timing properties of accreting black hole candidates (BHCs) provide important clues about the accretion and jet formation processes that occur in these systems. Dividing the different source behaviors into spectral states provides a framework for studying BHCs. To date, there have been three main classification schemes with Luminosity-based, Component-based, or Transition-based criteria. The canonical, Luminosity-based criteria and physical models that are based on this concept do not provide clear explanations for several phenomena, including hysteresis of spectral states and the presence of jets. I discuss the re-definitions of states, focusing on an application of the Component-based states to more than 400 RXTE observations of the recurrent BHC 4U 1630^17. We compare the X-ray properties for the recent 2002-2004 outburst to those of an earlier (1998) outburst, during which radio jets were observed. The results suggest a connection between hysteresis of states and major jet ejections, and it is possible that both of these are related to the evolution of the inner radius of the optically thick accretion disk.

Tomsick, John A.↗

Black Hole Spectral States and Physical Connections

The dramatic changes seen in the X-ray spectral and timing properties of accreting black hole candidates (BHCs) provide important clues about the accretion and jet formation processes that occur in these systems. Dividing the different source behaviors into spectral states provides a framework for studying BHCs. To date, there have been three main classification schemes with Luminosity-based, Component-based, or Transition-based criteria. The canonical, Luminosity-based criteria and physical models that are based on this concept do not provide clear explanations for several phenomena, including hysteresis of spectral states and the presence of jets. I discuss the re-definitions of states, focusing on an application of the Component-based states to more than 400 RXTE observations of the recurrent BHC 4U 1630-47. We compare the X-ray properties for the recent 2002-2004 outburst to those of an earlier (1998) outburst, during which radio jets were observed. The results suggest a connection between hysteresis of states and major jet ejections, and it is possible that both of these are related to the evolution of the inner radius of the optically thick accretion disk.

Tomsick, John A.↗