Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution regularity”

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 451 records · Page 25

Patterns of brittle deformation under extension on Venus

The development of fractures at regular length scales is a widespread feature of Venusian tectonics. Models of lithospheric deformation under extension based on non-Newtonian viscous flow and brittle-plastic flow develop localized failure at preferred wavelengths that depend on lithospheric thickness and stratification. The characteristic wavelengths seen in rift zones and tessera can therefore provide constraints on crustal and thermal structure. Analytic solutions were obtained for growth rates in infinitesimal perturbations imposed on a one-dimensional, layered rheology. Brittle layers were approximated by perfectly-plastic, uniform strength, overlying ductile layers exhibiting thermally-activated power-law creep. This study investigates the formation of faults under finite amounts of extension, employing a finite-element approach. Our model incorporates non-linear viscous rheology and a Coulomb failure envelope. An initial perturbation in crustal thickness gives rise to necking instabilities. A small amount of velocity weakening serves to localize deformation into planar regions of high strain rate. Such planes are analogous to normal faults seen in terrestrial rift zones. These 'faults' evolve to low angle under finite extension. Fault spacing, orientation and location, and the depth to the brittle-ductile transition, depend in a complex way on lateral variations in crustal thickness. In general, we find that multiple wavelengths of deformation can arise from the interaction of crustal and mantle lithosphere.

Neumann, G. A.↗

ARTEMIS: The First Mission to the Lunar Libration Orbits

The ARTEMIS mission will be the first to navigate to and perform stationkeeping operations around the Earth-Moon L1 and L2 Lagrangian points. The NASA Goddard Space Flight Center (GSFC) has previous mission experience flying in the Sun-Earth L1 (SOHO, ACE, WIND, ISEE-3) and L2 regimes (WMAP) and have maintained these spacecraft in libration point orbits by performing regular orbit stationkeeping maneuvers. The ARTEMIS mission will build on these experiences, but stationkeeping in Earth-Moon libration orbits presents new challenges since the libration point orbit period is on the order of two weeks rather than six months. As a result, stationkeeping maneuvers to maintain the Lissajous orbit will need to be performed frequently, and the orbit determination solutions between maneuvers will need to be quite accurate. The ARTEMIS mission is a collaborative effort between NASA GSFC, the University of California at Berkeley (UCB), and the Jet Propulsion Laboratory (JPL). The ARTEMIS mission is part of the THEMIS extended mission. ARTEMIS comprises two of the five THEMIS spacecraft that will be maneuvered from near-Earth orbits into lunar libration orbits using a sequence of designed orbital maneuvers and Moon & Earth gravity assists. In July 2009, a series of orbit-raising maneuvers began the proper orbit phasing of the two spacecraft for the first lunar flybys. Over subsequent months, additional propulsive maneuvers and gravity assists will be performed to move each spacecraft though the Sun-Earth weak stability regions and eventually into Earth-Moon libration point orbits. We will present the overall orbit designs for the two ARTEMIS spacecraft and provide analysis results of the 3/4-body dynamics, and the sensitivities of the trajectory design to both · maneuver errors and orbit determination errors. We will present results from the. initial orbit-raising maneuvers.

Woodward, Mark↗

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗

Generating local addresses and communication sets for data-parallel programs

Generating local addresses and communication sets is an important issue in distributed-memory implementations of data-parallel languages such as High Performance Fortran. We show that for an array A affinely aligned to a template that is distributed across p processors with a cyclic(k) distribution, and a computation involving the regular section A, the local memory access sequence for any processor is characterized by a finite state machine of at most k states. We present fast algorithms for computing the essential information about these state machines, and extend the framework to handle multidimensional arrays. We also show how to generate communication sets using the state machine approach. Performance results show that this solution requires very little runtime overhead and acceptable preprocessing time.

Chatterjee, Siddhartha↗

Nonlinear Estimation Approach to Real-Time Georegistration from Aerial Images

When taking aerial images, it is important to know locations of specific points of interest in an Earth-centered coordinate system (latitude, longitude, height). The correspondence between a pixel location in the image and its Earth coordinate is known as georegistration. There are two main technical challenges arising in the intended application. The first is that no known features are assumed to be available in any of the images. The second is that the intended applications are real time. Here, images are taken at regular intervals (i.e. once per second), and it is desired to make decisions in real time based on the geolocation of specific objects seen in the images as they arrive. This is in sharp contrast to most current methods for geolocation that operate "after-the-fact" by processing, on the ground, a database of stored images using computationally intensive methods. The solution is a nonlinear estimation algorithm that combines processed realtime camera images with vehicle position and attitude information ob tained from an onboard GPS receiver. This approach provides accurate georegistration estimates (latitude, longitude, height) of arbitrary features and/or points of interest seen in the camera images. This solves the georegistration problem at the modest cost of augmenting the camera information with a GPS receiver carried onboard the vehicle.

Bayard, David S.↗

Nonperturbative resolution of strong coupling singularities in 4D 𝒩 = 1 heterotic M-theory

We investigate the interior of the moduli space of four-dimensional 𝒩 = 1 theories of gravity arising from compactifications of the 𝐸 8 × 𝐸 8 heterotic string on Calabi-Yau threefolds. By studying the threshold corrections to the coupling of the heterotic gauge groups, we infer the existence of a strong coupling singularity for one of the perturbative heterotic gauge groups, which effectively yields an additional finite distance boundary of the classical scalar field space. In heterotic M-theory, this boundary maps to a domain wall solution for which the gauge coupling and the warp factor on one of the Horava-Witten 9-branes diverge, thus highlighting the gravitational origin of the classical strong coupling singularity. The divergence of the warp factor is, however, regulated once nonperturbative effects are taken into account, as we demonstrate by studying the instanton corrections to the 5D BPS domain wall equations. This regularization implies that the classical strong coupling boundary of the scalar field 4D 𝒩 = 1 heterotic M-theory is resolved, indicating that, at the quantum level, the field space can be extended beyond this classical boundary.

Cvetič, Mirjam↗

Leveraging Multitime Hamilton–Jacobi PDEs for Certain Scientific Machine Learning Problems

Hamilton-Jacobi partial differential equations (HJ PDEs) have deep connections with a wide range of fields, including optimal control, differential games, and imaging sciences. By considering the time variable to be a higher dimensional quantity, HJ PDEs can be extended to the multi-time case. In this paper, we establish a novel theoretical connection between specific optimization problems arising in machine learning and the multi-time Hopf formula, which corresponds to a representation of the solution to certain multi-time HJ PDEs. Through this connection, we increase the interpretability of the training process of certain machine learning applications by showing that when we solve these learning problems, we also solve a multi-time HJ PDE and, by extension, its corresponding optimal control problem. As a first exploration of this connection, we develop the relation between the regularized linear regression problem and the Linear Quadratic Regulator (LQR). We then leverage our theoretical connection to adapt standard LQR solvers (namely, those based on the Riccati ordinary differential equations) to design new training approaches for machine learning. Lastly, we provide some numerical examples that demonstrate the versatility and possible computational advantages of our Riccati-based approach in the context of continual learning, post-training calibration, transfer learning, and sparse dynamics identification.

97 MATHEMATICS AND COMPUTING↗

Clean Grid Vision: A U.S. Perspective - Chapter 2. Distribution Issues and Tools

The continued growth in distributed energy resources is fundamentally changing the way the distribution system is operated and planned. Additionally, an increasing number of grid issues are being regularly experienced by distribution engineers within an increasing number of utilities. Fortunately, these fundamental changes and increasingly regular issues are being addressed via continual improvements in existing distribution system modeling tools—for both tools used by industry and researchers—and via novel modeling methods, which are combining transmission and distribution system modeling to investigate new, often complex, system interactions at the transmission/distribution power system interface. As DERs continue to increase more and more utilities, regulators, and developers are finding usable solutions in the deployment of smart inverters with ever more capable functionalities. This increase in functionality is welcome but comes at the cost of more complex analysis when determining what specific settings should be used for an individual interconnection. With more time, experience, and likely the further development of advanced analysis tools, even higher levels of DER integration will result in a reliable, resilient, low-cost, and safely operating grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Impedance properties of circular microstrip antenna

A moment method solution to the input impedance of a circular microstrip antenna excited by either a microstrip feed or a coaxial probe is presented. Using the exact dyadic Green's function and the Fourier transform the problem is formulated in terms of Richmond's reaction integral equation from which the unknown patch current can be solved for. The patch current is expanded in terms of regular surface patch modes and an attachment mode (for probe excited case) which insures continuity of the current at probe/patch junction, proper polarization and p-dependance of patch current in the vicinity of the probe. The input impedance of a circular microstrip antenna is computed and compared with earlier results. Effect of attachment mode on the input impedance is also discussed.

Deshpande, M. D.↗

Characterization of Truncated Tumor-Associated NADH Oxidase (ttNOX)

Bacterial, plant and animal cells possess novel surface proteins that exhibit both NADH oxidation (NOX) or hydroquinone and protein disulfide-thiol interchange. These enzymatic activities alternate to yield oscillating patterns wjth period lengths of approximately 24 minutes. The catalytic period of NOX proteins are temperature compensated and gravity responsive. We report the cloning, expression and characterization of truncated tumor-associated NADH oxidase (ttNOX), in which the membrane spanning region has been deleted. The cDNA (originated from HeLa cells) was cloned into pET-34b and pET-14b (Novagen) vectors for E. coli expression. Optimized expression and purification protocols yielded greater than 300mg per liter of culture with greater than 95% purity. Circular dichroism data was collected from a 2.7mg/ml solution in a 0.1mm cuvette with variable scanning using an Olis RSM CD spectrophotometer. The ellipticity values were scanned from 190 to 260nm. The spectra recorded have characteristics for alpha proteins with band maxima at 216nm and a possible shoulder at 212nm at 12OC and 250 C. Protein crystal screens are in progress and, to date, only small crystals have been observed. The regular periodic oscillatory change in the ttNOX protein is indicative of a possible time-keeping functional role. A single protein possessing alternating catalytic activities, with a potential biological clock function, is unprecedented and structural determination is paramount to understanding this role.

Karr, Laurel J.↗

Evaluating Long-Term Changes in Atmospheric Ozone

Ozonesondes have made inexpensive, accurate measurements of ozone from the ground to 30km for more than 50 years. The data are used extensively for trend analyses and for evaluation of satellite and model data products, and are also part of climatologies that are used as a priori data for satellite retrievals. They are essential as a transfer standard when merging shorter satellite-derived time series, and are the most important source of trend-quality long-term records below about 18 km. The importance of long-term ozonesonde records as a stable reference has led to increased attention to quantifying uncertainties and changes in ozonesonde data. Based on past intercomparison data, ECC sondes show a modest (~1-5%) high bias with respect to UV-absorption measurements (including MOZAIC-IAGOS) in the troposphere, with an uncertainty of 5%, but no evidence of a change with time. Other sonde types show an increase of 5-20% in sensitivity to tropospheric ozone from 1970-1995. Agreement in the stratosphere is much better. In the past 20 years ozonesonde precision has improved by a factor of 2, primarily through the adoption of strict standard operating procedures. In addition, many of the most important ozonesonde records have been re-evaluated and corrected, and detailed estimates made of their biases and uncertainties. Corrections to historical data for known issues may reduce biases but simultaneously introduce additional uncertainties. The uncertainty budget for the ozone partial pressure reading has contributions from stoichiometry, cell background current, pump efficiency and temperature, sensing solution type and volume. The much-discussed “background current” has recently been shown to be primarily related to reaction stoichiometry. Several quality assurance issues remain, but are tractable problems that can be addressed with further research. This will be required if the present goal of better than 5% overall uncertainty throughout the global ozonesonde network is to be achieved. Finally, the importance of regular sonde intercomparisons, employing UV photometers traceable to the modern UV-absorption standard, and of regular comparison of operational station records with multiple satellite sensors, is emphasized.

Ozonesondes↗

QA/QC-ed Groundwater Level Time Series in PLM-1 and PLM-6 Monitoring Wells, East River, Colorado (2016-2022)

This data set contains QA/QC-ed (Quality Assurance and Quality Control) water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location at the East River Watershed, Colorado, USA. These wells are used to monitor subsurface water and carbon inventories and fluxes, and to determine the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry (see related references) can be used to estimate rates of solute export from the hillslope to the floodplain and river. QA/QC analysis of measured groundwater levels in monitoring wells PLM-1 and PLM-6 included identification and flagging of duplicated values of timestamps, gap filling of missing timestamps and water levels, removal of abnormal/bad and outliers of measured water levels. The QA/QC analysis also tested the application of different QA/QC methods and the development of regular (5-minute, 1-hour, and 1-day) time series datasets, which can serve as a benchmark for testing other QA/QC techniques, and will be applicable for ecohydrological modeling. The package includes a Readme file, one R code file used to perform QA/QC, a series of 8 data csv files (six QA/QC-ed regular time series datasets of varying intervals (5-min, 1-hr, 1-day) and two files with QA/QC flagging of original data), and three files for the reporting format adoption of this dataset (InstallationMethods, file level metadata (flmd), and data dictionary (dd) files).QA/QC-ed data herein were derived from the original/raw data publication available at Williams et al., 2020 (DOI: 10.15485/1818367). For more information about running R code file (10.15485_1866836_QAQC_PLM1_PLM6.R) to reproduce QA/QC output files, see README (QAQC_PLM_readme.docx). This dataset replaces the previously published raw data time series, and is the final groundwater data product for the PLM wells in the East River. Complete metadata information on the PLM1 and PLM6 wells are available in a related dataset on ESS-DIVE: Varadharajan C, et al (2022). https://doi.org/10.15485/1660962. These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales. 2022/09/09 Update: Converted data files using ESS-DIVE’s Hydrological Monitoring Reporting Format. With the adoption of this reporting format, the addition of three new files (v1_20220909_flmd.csv, V1_20220909_dd.csv, and InstallationMethods.csv) were added. The file-level metadata file (v1_20220909_flmd.csv) contains information specific to the files contained within the dataset. The data dictionary file (v1_20220909_dd.csv) contains definitions of column headers and other terms across the dataset. The installation methods file (InstallationMethods.csv) contains a description of methods associated with installation and deployment at PLM1 and PLM6 wells. Additionally, eight data files were re-formatted to follow the reporting format guidance (er_plm1_waterlevel_2016-2020.csv, er_plm1_waterlevel_1-hour_2016-2020.csv, er_plm1_waterlevel_daily_2016-2020.csv, QA_PLM1_Flagging.csv, er_plm6_waterlevel_2016-2020.csv, er_plm6_waterlevel_1-hour_2016-2020.csv, er_plm6_waterlevel_daily_2016-2020.csv, QA_PLM6_Flagging.csv). The major changes to the data files include the addition of header_rows above the data containing metadata about the particular well, units, and sensor description. 2023/01/18 Update: Dataset updated to include additional QA/QC-ed water level data up until 2022-10-12 for ER-PLM1 and 2022-10-13 for ER-PLM6. Reporting format specific files (v2_20230118_flmd.csv, v2_20230118_dd.csv, v2_20230118_InstallationMethods.csv) were updated to reflect the additional data. R code file (QAQC_PLM1_PLM6.R) was added to replace the previously uploaded HTML files to enable execution of the associated code. R code file (QAQC_PLM1_PLM6.R) and ReadMe file (QAQC_PLM_readme.docx) were revised to clarify where original data was retrieved from and to remove local file paths.

54 ENVIRONMENTAL SCIENCES↗

Generating local addresses and communication sets for data-parallel programs

Generating local addresses and communication sets is an important issue in distributed-memory implementations of data-parallel languages such as High Performance FORTRAN. We show that, for an array A affinely aligned to a template that is distributed across p processors with a cyclic(k) distribution and a computation involving the regular section A(l:h:s), the local memory access sequence for any processor is characterized by a finite state machine of at most k states. We present fast algorithms for computing the essential information about these state machines, and extend the framework to handle multidimensional arrays. We also show how to generate communication sets using the state machine approach. Performance results show that this solution requires very little run-time overhead and acceptable preprocessing time.

Chatterjee, Siddhartha↗

Limited-angle computed tomography with deep image and physics priors

Computed tomography is a well-established x-ray imaging technique to reconstruct the three-dimensional structure of objects. It has been used extensively in a variety of fields, from diagnostic imaging to materials and biological sciences. One major challenge in some applications, such as in electron or x-ray tomography systems, is that the projections cannot be gathered over all the angles due to the sample holder setup or shape of the sample. This results in an ill-posed problem called the limited angle reconstruction problem. Typical image reconstruction in this setup leads to distortion and artifacts, thereby hindering a quantitative evaluation of the results. To address this challenge, we use a generative model to effectively constrain the solution of a physics-based approach. Our approach is self-training that can iteratively learn the nonlinear mapping from partial projections to the scanned object. Because our approach combines the data likelihood and image prior terms into a single deep network, it is computationally tractable and improves performance through an end-to-end training. We also complement our approach with total-variation regularization to handle high-frequency noise in reconstructions and implement a solver based on alternating direction method of multipliers. We present numerical results for various degrees of missing angle range and noise levels, which demonstrate the effectiveness of the proposed approach.

47 OTHER INSTRUMENTATION↗

Deep Space navigation for the BioSentinel spacecraft science orbit

BioSentinel is an astrobiology small spacecraft mission. The payload consists of two parts, the first has optical and microfluidics sensors, and the second is a Linear Energy Transfer spectrometer that has the objective to measure deep space radiation from events such as coronal mass ejections. The goal of the mission is to observe potential DNA damage due to the radiation in heliocentric space on the living organism Saccharomyces cerevisiae, which is a budding yeast. Two types of this living organism are included in the payload. The first is a natural type that is more radiation tolerant, while the second is a mutant strain that has a deficiency in a gene that allows DNA repair once damage occurs. The impact caused by the radiation on the DNA is compared to an identical sample aboard the International Space Station, as well as another identical sample at a laboratory on the ground. The BioSentinel mission consists of a 6U CubeSat currently ,as of January 2024, active in heliocentric orbit. The spacecraft was launched aboard the first SLS flight as part of the Artemis-I campaign in November 2022. After successful deployment from the launch vehicle, it performed a lunar flyby with an altitude of 406 km. The delta-V imparted by the flyby provided the necessary energy to achieve a heliocentric orbit, in an Earth-trailing pattern. The navigation analysis consisted of a Kalman-filter that utilized data from the Deep Space Network and the ESA Estrack network. All those antennas were needed since the Artemis-1 campaign included the deployment of several other cubesats, therefore the scheduling process required more antenna assets than usual due to simultaneous demands from various missions. The processed tracking data was later also refined with a smoother in order to obtain a more accurate solution. The type of tracking data included TCP, Sequential Range, Doppler and Range formats. The solar radiation pressure coefficient, as well as the delta-V from the deployment and the flyby were modeled to obtain suitable solutions that could decrease the position and velocity uncertainties at several steps along the mission concept of operations. The final product each time resulted in updated ephemeris files that were used by the mission and the antenna networks as the mission progressed. Once in the final science orbit, the utilized antennas are only from the DSN network and the data format is bounded to just TCP. Regular orbit determination is performed, every two weeks. The spacecraft is in a nominal well-known orbit, performing regular operations. This paper includes an analysis of the final science orbit, the techniques and procedures utilized to perform orbit determination and a description of the overall navigation campaign produced during the mission and, more specifically, during the final science operations in Deep Space.

BioSentinel↗

Human-system Interaction Issues and Proposed Solutions To Promote Successful Maturation of the UTM System

Over five years, NASA, together with partnering organizations, has been developing and successfully demonstrating the maturing capabilities of the Unmanned Aircraft Systems (UAS) Traffic Management (UTM) system and its ability to support communication and coordination among small UAS operations through a series of flight tests. During these flight tests, human-system interaction (HSI) elements were also explored in order to identify the barriers to implementation as human operators transitionally fulfill roles that will be ultimately tasked to future automation. Throughout the tests, similar issues were regularly documented and are expected to persist if not formally addressed by consistent procedures or regulation. Documented here, along with suggested mitigations, are the most frequently noted HSI items, which include operator training, data standardization, and information quantity.

UTM↗

Human-system interaction issues and proposed solutions to promote successful maturation of the UTM system

Over five years, NASA, together with partnering organizations, has been developing and successfully demonstrating the maturing capabilities of the Unmanned Aircraft Systems (UAS) Traffic Management (UTM) system and its ability to support communication and coordination among small UAS operations through a series of flight tests. During these flight tests, human-system interaction (HSI) elements were also explored in order to identify the barriers to implementation as human operators transitionally fulfill roles that will be ultimately tasked to future automation. Throughout the tests, similar issues were regularly documented and are expected to persist if not formally addressed by consistent procedures or regulation. Documented here, along with suggested mitigations, are the most frequently noted HSI items, which include operator training, data standardization, and information quantity.

UTM↗

Regularizing the linearly extrapolated BDF2 scheme for incompressible flows with time relaxation

This paper presents a highly-efficient finite element scheme for the time relaxation model (TRM). The efficiency is achieved through the second-order BDF2 time-stepping scheme with linear extrapolation (BDF2LE). The accuracy of the scheme is also greatly enhanced through the use of the divergence-free Scott-Vogeulis finite elements, and van Cittert approximate deconvolution. A complete finite element analysis is provided, which includes rigorous proofs for the stability, well-possessedness, and convergence of both velocity and pressure solutions. Furthermore, we also demonstrate that the inclusion of the linear time relaxation term preserves the long-time stability of the unregularized BDF2LE scheme. Finally, numerical experiments are presented that demonstrate the added stability and accuracy that time relaxation can provide.

97 MATHEMATICS AND COMPUTING↗