Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local time stepping”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Seasonal Glacier and Snow Loading in Svalbard Recovered from Geodetic Observations

We processed time series from seven Global Navigation Satellite System (GNSS) stations and one Very Long Baseline Interferometry (VLBI) station in Svalbard. The goal was to capture the seasonal vertical displacements caused by elastic response of variable mass load due to ice and snow accumulation. We found that estimates of the annual signal in different GNSS solutions disagree by more than 3 mm which makes geophysical interpretation of raw GNSS time series problematic. To overcome this problem, we have used an enhanced Common Mode (CM) filtering technique. The time series are differentiated by the time series from remote station BJOS with known mass loading signals removed a priori. Using this technique, we have achieved a substantial reduction of the differences between the GNSS solutions. We have computed mass loading time series from a regional Climatic Mass Balance (CMB) and snow model that provides the amount of water equivalent at a 1 km resolution with a time step of 7 days. We found that the entire vertical loading signal is present in data of two totally independent techniques at a statistically significant level of 95%. This allowed us to conclude that the remaining errors in vertical signal derived from the CMB model are less than 0.2 mm at that significance level. Refining the land water storage loading model with a CMB model resulted in a reduction of the annual amplitude from 2:1 mm to 1:1 mm in the CM filtered time series, while it had only a marginal impact on raw time series. This provides a strong evidence that CM filtering is essential for revealing local periodic signals when a millimetre level of accuracy is required.

H P Kierulf↗
NASA NTRSDate not supplied

Exploring Innovative Methods for Earth Science Communication Through NASA’s Earth Information Center

Environmental challenges stemming from unprecedented changes in our planet have increased the need to communicate our understanding of Earth as a connected system. In June 2023, NASA unveiled the Earth Information Center (EIC), part of NASA’s Earth Action Strategy aimed at amplifying NASA’s Earth science impact by engaging decision-makers beyond the scientific community. Existing in both physical and virtual space, the EIC provides critical information on how the planet is changing and affecting livelihoods across the globe. For over fifty years, NASA has utilized the unique vantage point space offers to monitor Earth in collaboration with Federal partners. Through an expanding fleet of Earth observing satellites and sensors, NASA informs local to global decision making within nine thematic areas: agriculture, biodiversity, disasters, greenhouse gases, health and air quality, sea level rise and coastal impacts, sustainable energy, water resources, and wildfires. Visitors to the inaugural location of the EIC at NASA Headquarters in Washington D.C., can interact with Earth observation data and models through innovative exhibits and storytelling. Designed by an interdisciplinary team, the EIC is designed for multiple modality learning to reach a wide base of end-users. Data driven visualizations, near real-time information, and curated stories highlight the applications of Earth observations to address environmental challenges. Complementary exhibits and immersive storytelling allow visitors to physically step inside environmental data and connect with information on an emotional level. Taken together, the collective EIC provides scientists and non-scientists alike with the ability to see our home planet as NASA sees it, as a system of systems, while delivering easily accessible information to those who need it most.

Earth Information Center↗
DOE OSTI2024

Monitoring water quality in the lower Kansas River using remote sensing

Abstract We demonstrate how to combine remote sensing data from satellite imagery (Sentinel‐2) with in situ water quality gauging (USGS Super Gages and the Gybe hyperspectral radiometer) to create spatially dense maps of water quality parameters (chlorophyll‐a concentration, turbidity, and nitrate plus nitrite concentration) along the lower Kansas River. The water quality maps are created using locally tuned models of the target water quality parameters, and this study describes the steps used to design, calibrate, and validate the empirical correlations. Water quality parameters such as chlorophyll‐a concentration are correlated with well‐studied absorption and scattering features in the visible spectrum (roughly 400–700 nm). Nutrients (such as nitrate plus nitrite concentration) lack strong absorption features in the visible spectrum, and in those cases we describe a novel surrogate data modeling approach that identifies overlapping water parcels between the in situ gauging and the remote sensing imagery. Measurements from the overlapping water parcels yield excellent correlations () for the target water quality parameters for limited windows of time (or limited sections of river reaches). Examples are provided illustrating how the water quality maps can be used to track river inputs from ungauged sources (such as creeks), or reveal the mixing patterns at the confluences.

Tufillaro, Nicholas↗
DOE OSTI2024

Generation and Life Cycle of Solar Spicules

The physical mechanism for the creation of solar spicules is proposed with three stages of their life cycle. It is assumed that at stage I the density hump is formed locally in the x-y plane in the lower chromosphere in the presence of temperature gradients of electrons and ions along the z-axis (the vertical direction). In this region, the density structure of quasi-neutral (n i ≃ n e = n) plasma after taking birth is accelerated in the vertical direction owing to the thermal force F th ∝ ∇n(x, y, t) × (∇T e + ∇T i ). The exact time-dependent analytical solution of two-fluid plasma equations is presented assuming that density is maximum at the center of the density structure and decays away from it gradually. The 2D density structure is created as a step function H(t) in time at the bottom of the chromosphere, and consequently, the vertical plasma velocity turns out to be the ramp function of time R(t) = tH(t), whereas the source term S(x, y, t) for the density follows the delta function δ(t) form. The upward acceleration a = a(x, y) $\hat{z}$ produced in this density structure is greater than the downward constant solar acceleration −g ⊙ in the chromosphere. In the transition region, the temperature gradients are steeper; therefore, the upward acceleration increases in magnitude g ⊙ ≪ a and the density hump spends less time there. This is stage II of its life cycle. In stage III, the density structure enters into the corona, where the gradients of temperatures vanish and the structure decelerates to zero velocity under the action of the solar gravitational force.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

A fast, time-accurate unsteady full potential scheme

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Shankar, V.↗
NASA NTRS2010

A Controller-in-the Loop Simulation of Ground-Based Automated Separation Assurance in a NextGen Environment

A controller-in-the-loop simulation was conducted in the Airspace Operations Laboratory (AOL) at the NASA Ames Research Center to investigate the functional allocation aspects associated with ground-based automated separation assurance in a far-term NextGen environment. In this concept, ground-based automation handled the detection and resolution of strategic and tactical conflicts and alerted the controller to deferred situations. The controller was responsible for monitoring the automation and managing situations by exception. This was done in conditions both with and without arrival time constraints across two levels of traffic density. Results showed that although workload increased with an increase in traffic density, it was still manageable in most situations. The number of conflicts increased similarly with a related increase in the issuance of resolution clearances. Although over 99% of conflicts were resolved, operational errors did occur but were tied to local sector complexities. Feedback from the participants revealed that they thought they maintained reasonable situation awareness in this environment, felt that operations were highly acceptable at the lower traffic density level but were less so as it increased, and felt overall that the concept as it was introduced here was a positive step forward to accommodating the more complex environment envisioned as part of NextGen.

Homola, J.↗
DOE OSTI2022

Dynamic energy system modeling using hybrid physics-based and machine learning encoder–decoder models

Three model configurations are presented for multi-step time series predictions of the heat absorbed by the water and steam in a thermal power plant. The models predict over horizons of 2, 4, and 6 steps into the future, where each step is a 5-minute increment. The evaluated models are a pure machine learning model, a novel hybrid machine learning and physics-based model, and the hybrid model with an incomplete dataset. The hybrid model deconstructs the machine learning into individual boiler heat absorption units: economizer, water wall, superheater, and reheater. Each configuration uses a gated recurrent unit (GRU) or a GRU-based encoder–decoder as the deep learning architecture. Mean squared error is used to evaluate the models compared to target values. The encoder–decoder architecture is over 11% more accurate than the GRU only models. The hybrid model with the incomplete dataset highlights the importance of the manipulated variables to the system. The hybrid model, compared to the pure machine learning model, is over 10% more accurate on average over 20 iterations of each model. Automatic differentiation is applied to the hybrid model to perform a local sensitivity analysis to identify the most impactful of the 72 manipulated variables on the heat absorbed in the boiler. The models and sensitivity analyses are used in a discussion about optimizing the thermal power plant.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

One dimensional wormhole corrosion in metals

Corrosion is a ubiquitous failure mode of materials. Often, the progression of localized corrosion is accompanied by the evolution of porosity in materials previously reported to be either three-dimensional or two-dimensional. However, using new tools and analysis techniques, we have realized that a more localized form of corrosion, which we call 1D wormhole corrosion, has previously been miscategorized in some situations. Using electron tomography, we show multiple examples of this 1D and percolating morphology. To understand the origin of this mechanism in a Ni-Cr alloy corroded by molten salt, we combined energy-filtered four-dimensional scanning transmission electron microscopy and ab initio density functional theory calculations to develop a vacancy mapping method with nanometer-resolution, identifying a remarkably high vacancy concentration in the diffusion-induced grain boundary migration zone, up to 100 times the equilibrium value at the melting point. Deciphering the origins of 1D corrosion is an important step towards designing structural materials with enhanced corrosion resistance.

36 MATERIALS SCIENCE↗
NASA NTRS1991

A 1 MW, 100 kV, less than 100 kg space based dc-dc power converter

A 1 MW dc-dc power converter has been designed which has an input voltage of 5 kV +/-3 percent, an output voltage of 100 kV +/- 0.25 percent, and a run time of 1000 s at full power. The estimated system mass is 83.8 kg, giving a power density of 11.9 kW/kg. The system exceeded the weight goal of 10 kW/kg through the use of innovative components and system concepts. The system volume is approximately 0.1 cu m, and the overall system efficiency is estimated to be 87 percent. Some of the unique system features include a 50-kHz H-bridge inverter using MOS-controlled thyristors as the switching devices, a resonance transformer to step up the voltage, open-cycle cryogenic hydrogen gas cooling, and a nonrigid, inflatable housing which provides on-demand pressurization of the power converter local environment. This system scales very well to higher output powers. The weight of the 10-MW system with the same input and output voltage requirements and overall system configuration is estimated to be 575.3 kg. This gives a power density of 17.4 kW/kg, significantly higher than the 11.9 kW/kg estimated at 1 MW.

Cooper, J. R.↗
NASA NTRS1990

Two- and three-dimensional blade vortex interactions

A three-dimensional time dependent Navier-Stokes analysis was applied to the rotor blade vortex interaction (BVI) problem. The numerical procedure is an iterative implicit procedure using three point central differences to represent spatial derivatives. A series of calculations were made to determine the time steps, pseudo-time steps, iterations, artificial dissipation level, etc. required to maintain a nondissipative vortex. Results show the chosen method to have excellent non-dissipative properties provided the correct parameters are chosen. This study was used to set parameters for both two- and three-dimensional blade vortex interaction studies. The two-dimensional study considered the interaction between a vortex and a NACA0012 airfoil. The results showed the detailed physics during the interaction including the pressure pulse propagating from the blade. The simulated flow physics was qualitatively similar to that experimentally observed. The 2-D BVI phenomena is the result of the buildup and violent collapse of the shock waves and local supersonic pockets on the blade surfaces. The resulting pressure pulse build-up appears to be centered at the blade leading edge. The three-dimensional interaction study considered the case of a vortex at 20 deg incidence to the blade leading edge. Although the qualitative results were similar to that of the two-dimensional interaction, details clearly showed the three-dimensional nature of the interaction process.

Davoudzadeh, F.↗
NASA NTRS2001

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. In the absence of any record of the earliest ancestors of contemporary cells, protocells, the most direct way to test our understanding of their characteristics is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs of molecules that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing and storing energy, and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following, questions about these proteins: (1) How do small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides aggregate to form membrane-spannin(y structures (e.g., channels)? (3) By what mechanisms do such aggregates perform their functions? The simulations are performed using the molecular dynamics (MD) method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated directly in real space while long-range forces are evaluated in reciprocal space, usually using a particle-mesh algorithm which is of order O(NlnN). Currently, a time step of 2 femtoseconds is typically used, thereby making studies of problems occurring on multi-nanosecond time scales (10(exp 6) - 10(exp 8) time steps) accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude. Such an extension requires both algorithmic improvements and codes scalable to a large number of parallel processors. Work in this direction is in progress. Two specific series of simulations that demonstrate how peptides self-organize and function in membranes are discussed. In one series of simulations, it was shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and, simultaneously, fold into two different helical structures, which remain in equilibrium. Once in the membrane, the peptides can readily change their orientation, especially in response to local electric fields. This structural and orientational flexibility of peptides with changing conditions may have provided a mechanism of transmitting signals between the environment and the interior of the protocell. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations demonstrate that protons move through the gate by a "shuttle" mechanism, wherein one histidine is protonated on the extracellular side and, subsequently, the proton bound on the opposite side is released.

Pohorille, Andrew↗
NASA NTRS2000

Gravity independence of seed-to-seed cycling in Brassica rapa

Growth of higher plants in the microgravity environment of orbital platforms has been problematic. Plants typically developed more slowly in space and often failed at the reproductive phase. Short-duration experiments on the Space Shuttle showed that early stages in the reproductive process could occur normally in microgravity, so we sought a long-duration opportunity to test gravity's role throughout the complete life cycle. During a 122-d opportunity on the Mir space station, full life cycles were completed in microgravity with Brassica rapa L. in a series of three experiments in the Svet greenhouse. Plant material was preserved in space by chemical fixation, freezing, and drying, and then compared to material preserved in the same way during a high-fidelity ground control. At sampling times 13 d after planting, plants on Mir were the same size and had the same number of flower buds as ground control plants. Following hand-pollination of the flowers by the astronaut, siliques formed. In microgravity, siliques ripened basipetally and contained smaller seeds with less than 20% of the cotyledon cells found in the seeds harvested from the ground control. Cytochemical localization of storage reserves in the mature embryos showed that starch was retained in the spaceflight material, whereas protein and lipid were the primary storage reserves in the ground control seeds. While these successful seed-to-seed cycles show that gravity is not absolutely required for any step in the plant life cycle, seed quality in Brassica is compromised by development in microgravity.

Flight Experiment↗
NASA NTRS1995

An Extended Lagrangian Method

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Liou, Meng-Sing↗
DOE OSTI2024

Control of Quantized Spontaneous Emission from Single GaAs Quantum Dots Embedded in Huygens’ Metasurfaces

Advancements in photonic quantum information systems (QIS) have driven the development of high-brightness, on-demand, and indistinguishable semiconductor epitaxial quantum dots (QDs) as single photon sources. Strain-free, monodisperse, and spatially sparse local-droplet-etched (LDE) QDs have recently been demonstrated as a superior alternative to traditional Stranski–Krastanov QDs. However, integration of LDE QDs into nanophotonic architectures with the ability to scale to many interacting QDs is yet to be demonstrated. Here, we present a potential solution by embedding isolated LDE GaAs QDs within an Al 0.4 Ga 0.6 As Huygens’ metasurface with spectrally overlapping fundamental electric and magnetic dipolar resonances. We demonstrate for the first time a position- and size-independent, 1 order of magnitude increase in the collection efficiency and emission lifetime control for single-photon emission from LDE QDs embedded within the Huygens’ metasurfaces. Our results represent a significant step toward leveraging the advantages of LDE QDs within nanophotonic architectures to meet the scalability demands of photonic QIS.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1998

Parallel Directionally Split Solver Based on Reformulation of Pipelined Thomas Algorithm

In this research an efficient parallel algorithm for 3-D directionally split problems is developed. The proposed algorithm is based on a reformulated version of the pipelined Thomas algorithm that starts the backward step computations immediately after the completion of the forward step computations for the first portion of lines This algorithm has data available for other computational tasks while processors are idle from the Thomas algorithm. The proposed 3-D directionally split solver is based on the static scheduling of processors where local and non-local, data-dependent and data-independent computations are scheduled while processors are idle. A theoretical model of parallelization efficiency is used to define optimal parameters of the algorithm, to show an asymptotic parallelization penalty and to obtain an optimal cover of a global domain with subdomains. It is shown by computational experiments and by the theoretical model that the proposed algorithm reduces the parallelization penalty about two times over the basic algorithm for the range of the number of processors (subdomains) considered and the number of grid nodes per subdomain.

Povitsky, A.↗
NASA NTRS2012

Coronal Magnetic Field Measurement Using CME-Driven Shock Observations

Collisionless shocks form ahead of coronal mass ejections (CMEs) when the CME speed exceeds the Alfven speed of the ambient plasma in the corona and interplanetary medium. The shock stands at a distance from the CME flux rope that depends on the shock Mach number, the geometry of the driver, and the adiabatic index. While the shock ahead of the CME has been observed for a long time in the in situ data, it has been identified recently near the Sun in the coronagraphic and EUV images. Unlike in situ observations, the imaging observations are two dimensional, so one can better discern the CME-shock relationship near the Sun. Gopalswamy and Yashiro demonstrated that the coronal magnetic field can be derived from the shock standoff distance measured in coronagraphic images. The method involves measuring the standoff distance, the radius of curvature of the flux rope, and assuming the value of the adiabatic index and deriving the Alfvenic Mach number. The next step is to derive the Alfvenic Mach number from the measured shock speed and an estimate of the local solar wind speed. The final step involves deriving the magnetic field from the Alfven speed by measuring the local plasma density either from coronagraphic (polarized brightness) images or from the band-splitting of type II radio bursts. In this paper, we derive the combined magnetic field profile from near the Sun to the edge of the LASCO field of view (1.5 to 30 solar radii) and compare it with the current model profiles.

Gopalswarmy, Nat↗
NASA NTRS2012

Methods for Detection and Correction of Sudden Pixel Sensitivity Drops

PDC 8.0 includes implementation of a new algorithm to detect and correct step discontinuities appearing in roughly one of every twenty stellar light curves during a given quarter. An example of such a discontinuity in an actual light curve is shown in fig. 1. The majority of such discontinuities are believed to result from high-energy particles (either cosmic or solar in origin) striking the photometer and causing permanent local changes (typically -0.5% in summed apertures) in quantum efficiency, though a partial exponential recovery is often observed. Since these features, dubbed sudden pixel sensitivity dropouts (SPSDs), are uncorrelated across targets they cannot be properly accounted for by the current detrending algorithm. PDC de-trending is based on the assumption that features in flux time series are due either to intrinsic stellar phenomena or to systematic errors and that systematics will exhibit measurable correlations across targets. SPSD events violate these assumptions and their successful removal not only rectifies the flux values of affected targets, but demonstrably improves the overall performance of PDC de-trending.

Kolodziejczak, Jeffery↗
NASA NTRS2019

Satellite-Based Assessment of Grassland Conversion and Related Fire Disturbance in the Kenai Peninsula, Alaska

Spruce beetle-induced (Dendroctonus rufipennis (Kirby)) mortality on the Kenai Peninsula has been hypothesized by local ecologists to result in the conversion of forest to grassland and subsequent increased fire danger. This hypothesis stands in contrast to empirical studies in the continental US which suggested that beetle mortality has only a negligible effect on fire danger. In response, we conducted a study using Landsat data and modeling techniques to map land cover change in the Kenai Peninsula and to integrate change maps with other geospatial data to predictively map fire danger for the same region. We collected Landsat imagery to map land cover change at roughly five-year intervals following a severe, mid-1990s beetle infestation to the present. Landcover classification was performed at each time step and used to quantify grassland encroachment patterns over time. The maps of land cover change along with digital elevation models (DEMs),temperature, and historical fire data were used to map and assess wildfire danger across the study area. Results indicate the highest wildfire danger tended to occur in herbaceous and black spruce land cover types, suggesting that the relationship between spruce beetle damage and wildfire danger in costal Alaskan forested ecosystems differs from the relationship between the two in the forests of the coterminous United States. These change detection analyses and fire danger predictions provide the Kenai National Wildlife Refuge (KENWR) ecologists and other forest managers a better understanding of the extent and magnitude of grassland conversion and subsequent change in fire danger following the 1990s spruce beetle outbreak.

wildfire↗