Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial optimization”

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 1,099 records · Page 61

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗

A hybrid architecture for volt-var control in active distribution grids

Modern active distribution grids are characterized by the increasing penetration of distributed energy resources (DERs). The proper coordination and scheduling of a large numbers of these small-scale and spatially distributed DERs is necessary, and warrants the use of novel distributed approaches. In this paper, we propose a hybrid volt-var control architecture for the distribution grid, which leverages existing centralized and local approaches to planning, decision making, and control, and augments it with distributed optimization and distributed control for DER management. First, we propose a convex model to describe the power physics of distribution grids of meshed topology and unbalanced structure, based on current injection and McCormick Envelopes. Second, we employ the distributed proximal atomic coordination (PAC) algorithm to coordinate DERs to provide voltage support. We implement volt-var optimization by optimally coordinating DERs including PV smart inverters and demand response. We present results using the IEEE-34 bus network, using real data from a distribution feeder in Hawaii, to model load and PV generation. Different levels of DER penetration and objective functions are simulated. Finally, our results show the need for the coordination of DERs to improve voltage profiles, even in networks with existing voltage control devices. Further, we show the need for flexible reactive power capabilities to achieve desired grid performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Achieving Realistic Runoff in the Western United States with a Land Surface Model Forced by Dynamically Downscaled Meteorology

Abstract In this study, we calibrate a regional climate model’s (RCM) underlying land surface model (LSM). In addition to providing a realistic representation of runoff across the hydroclimatically diverse western United States, this is done to take advantage of the RCM’s ability to physically resolve meteorological forcing data in ungauged regions, and to prepare the calibrated hydrologic model for tight coupling, or the ability to represent land surface–atmosphere interactions, with the RCM. Specifically, we use a 9-km resolution meteorological forcing dataset across the western United States, from the fifth generation ECMWF Reanalysis (ERA5) downscaled by the Weather Research Forecasting (WRF) regional climate model, as an offline forcing for Noah-Multiparameterization (Noah-MP). We detail the steps involved in producing an LSM capable of accurately representing runoff, including physical parameterization selection, parameter calibration, and regionalization to ungauged basins. Based on our model evaluation from 1954 to 2021 for 586 basins with daily natural streamflow, the streamflow bias is reduced from 24.2% to 4.4%, and the median daily Nash–Sutcliffe efficiency (NSE) is improved from 0.12 to 0.36. When validating against basins with monthly natural streamflow data, we obtain a similar reduction in bias and a median monthly NSE improvement from 0.18 to 0.56. In this study, we also discover the optimal setup when using a donor-basin method to regionalize parameters to ungauged basins, which can vary by 0.06 NSE for unique designs of this regionalization method. Significance Statement This study provides useful guidance for improving a land surface model to accurately represent runoff across a spatially extensive and hydroclimatically diverse region (the western United States). The land surface model is updated to represent runoff more accurately at gauged basins, and then additionally updated for basins without observational data using a mathematical approach called the donor-basin method. We make use of a regional climate model’s reanalysis-derived meteorological data and its underlying land surface model to achieve realistic runoff. The calibrated land surface model can thus be tightly coupled in subsequent studies in a manner that should more accurately reflect runoff conditions. Findings from this study will serve as a useful reference for the atmospheric (and hydrologic) modeling communities and their ability to represent large-scale hydrology accurately.

Meteorology & Atmospheric Sciences↗

Resolution improvement of remote sensing data.

Discussion of the theory of a technique for image restoration in the reconstruction of a spatial input scene mapped by a line scanner. Special attention is given to the application of this technique in reducing the data telemetered from the Infrared Scanning Radiometer of the Apollo 17 spacecraft. The technique is based on the theory of splines and uses the results of optimal approximation studies by Colomb (1959), Sard (1967), and Anselone (1968). A recursive algorithm is also developed to improve the reconstruction. The study treats a one-dimensional case, but the treatment can be readily extended to multi-dimensional situations.

Caprihan, A.↗

Spectral methods in time for hyperbolic equations

A pseudospectral numerical scheme for solving linear, periodic, hyperbolic problems is described. It has infinite accuracy both in time and in space. The high accuracy in time is achieved without increasing the computational work and memory space which is needed for a regular, one step explicit scheme. The algorithm is shown to be optimal in the sense that among all the explicit algorithms of a certain class it requires the least amount of work to achieve a certain given resolution. The class of algorithms referred to consists of all explicit schemes which may be represented as a polynomial in the spatial operator.

Tal-Ezer, H.↗

A Comparison of Approximation Modeling Techniques: Polynomial Versus Interpolating Models

Two methods of creating approximation models are compared through the calculation of the modeling accuracy on test problems involving one, five, and ten independent variables. Here, the test problems are representative of the modeling challenges typically encountered in realistic engineering optimization problems. The first approximation model is a quadratic polynomial created using the method of least squares. This type of polynomial model has seen considerable use in recent engineering optimization studies due to its computational simplicity and ease of use. However, quadratic polynomial models may be of limited accuracy when the response data to be modeled have multiple local extrema. The second approximation model employs an interpolation scheme known as kriging developed in the fields of spatial statistics and geostatistics. This class of interpolating model has the flexibility to model response data with multiple local extrema. However, this flexibility is obtained at an increase in computational expense and a decrease in ease of use. The intent of this study is to provide an initial exploration of the accuracy and modeling capabilities of these two approximation methods.

Giunta, Anthony A.↗

The PolCube CubeSat Polarimeter for Earth Science

PolCube is a 12U CubeSat + polarimeter instrument designed by NASA Langley and the Korea Astronomy and Space Science Institute (KASI) for Earth Science. PolCube is based on the PolCam polarimeter onboard the Korean Pathfinder Lunar Observatory (KPLO) that launched in August 2022. The objective of the PolCube instrument is to retrieve detailed fine-mode (pollution and smoke) and coarse-mode (sea-salt and dust) aerosol properties over the ocean for a range of light to heavy aerosol loadings using its polarimetric-imaging capability at multiple angles and wavelengths from 410 − 865 nm. An additional objective is to discriminate aerosols from thin clouds. We quantify the performance of aerosol and ocean remote sensing products from the PolCube polarimeter instrument using the Microphysical Aerosol Properties from Polarimetry (MAPP) remote sensing retrieval algorithm. PolCube’s accurate and high-resolution aerosol-retrieval products will provide unique spatial and temporal coverage of the Earth that can be used synergistically with other instruments, such as the PACE (Plankton, Aerosols, Clouds and Ecosystems) and GEMS (Geostationary Environmental Monitoring Spectrometer) mission to improve air-quality forecasting. We present the PolCube-MAPP retrieval algorithm, which used optimal estimation and artificial intelligence, as well as multiple powerful inherent optical property look-up-tables for the Earth’s aerosol, cloud, and hydrosol particles. We estimate that PolCube can retrieve total aerosol optical depth at 555 nm (AOD555) within ±0.068, fine-mode AOD555 within ±0.078, and fine-mode single-scattering albedo within ±0.036, where all uncertainties are expressed as one standard deviation (1σ).

Snorre Stamnes↗

Near Infrared In-situ Measurements for the Detection of Keyhole Porosity in Autogenous Ti-6AL-4V Welds

Certification of additive manufactured metal parts requires nondestructive evaluation (NDE) to ensure build quality. NDE can be performed during the build process or post build. For large parts with complex geometries, post build NDE can be challenging. In-situ NDE potentially provides a way to perform the inspection layer by layer. This work explores the use of a high speed near infrared (NIR) camera that is focused in-line with a laser to obtain high spatial and temporal resolution thermal imagery of the melt pool and associated cooling areas. The thermal data is obtained during a laser melting process using Ti-6AL-4V plate and of particular interest is the detection of keyhole porosity. Keyhole porosity can result from non-optimal build conditions, such as excessive laser power at a given laser scanning speed, that creates an entrapped bubble. The NIR measured melt pool and cooling areas are processed to detect keyhole porosity. The results are compared to X-ray CT results for validation.

additive manufacturing↗

A multiscale landscape approach for prioritizing river and stream protection and restoration actions

River and stream conservation programs have historically focused on a single spatial scale, for example, a watershed or stream site. Recently, the use of landscape information (e.g., land use and land cover) at multiple spatial scales and over large spatial extents has highlighted the importance of incorporating a landscape perspective into stream protection and restoration activities. Previously, we developed a novel framework that links information about watershed-, catchment-, and reach-scale integrity with stream biological condition using scatterplots and a landscape integrity map. Here we examined an application of this approach for streams in urban and other settings in King County, Washington State, United States, where we related stream macroinvertebrate condition to two indices of landscape integrity, the US Environmental Protection Agency's (USEPA) nationally available Index of Watershed Integrity (IWI) and Index of Catchment Integrity (ICI). We generated a scatterplot of IWI versus ICI for sample sites, where points represented site macroinvertebrate condition from poor to good. The same data were also visualized as a landscape integrity map that displayed catchments of King County according to the level of watershed and catchment integrity (high or low IWI/ICI). Almost three-quarters of poor-condition sites were associated with high-integrity watersheds and catchments (i.e., underperforming sites), which suggested that either one or both national indicators were insufficient for this area, and that sites underperformed because of local-scale factors. In response, we used a catchment-scale indicator related to forest condition (PctForestCat) after examining several GIS-based dispersal indicators from the National Hydrography Dataset and other candidates from the USEPA's StreamCat dataset. We then compared the results of the scatterplots and maps based on the current and original analyses and found that many of the sites previously classified as underperforming now performed as expected, that is, they were poor-condition sites in poor-condition catchments. This analysis demonstrates how results based on a national dataset can be improved by developing an alternative that represents regionally important stressors. The methods used to develop an effective landscape indicator based on StreamCat datasets, and the utility of the multiscale approach, could provide important tools for prioritizing, optimizing, and communicating stream conservation actions.

54 ENVIRONMENTAL SCIENCES↗

Modeling Electrolytic Conversion of Metabolic CO2 and Optimizing a Macrofluidic Electrochemical Reactor for Advanced Closed Loop Life Support Systems

The International Space Station (ISS) is currently equipped with a complex, heavy, and power consuming system that recovers approximately 50% of O2 from metabolic CO2. Future long duration missions will require a sustainable and highly efficient system capable of yielding a minimum of 75% O2 recovery. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to significantly increase current O2 recovery efficiency and reduce complexity of the system. This paper presents a comprehensive multi-physic 3D model developed at MSFC on CO2 conversion to O2 and C2H4 at standard conditions via MFECR. The 3D spatial domain of the model is a replica of the actual MFECR’s 3D drawing generated for the MFECR fabrication and operated to recover O2 from CO2 yielding C2H4 as byproduct. Electrochemical (EC) physics that includes EC multicomponent reaction mechanisms, mass transport, and current density distributions is coupled in the model with all the other physics phenomena involved in the process, such as free and porous fluid flow, multicomponent mass transfer, heat transfer, and DC electrical current generation along with Joule heating effect. The authors plan to use experimental results to validate this comprehensive and rigorous model and build a reliable simulator that will not only assist the authors on the MFECR design but also optimize its operation.

Dominguez, Jesus A.↗

Spatial and temporal traction response in human airway smooth muscle cells

Tractions that cells exert on their substrates are essential in cell spreading, migration, and contraction. These tractions can be determined by plating the cells on a flexible gel and measuring the deformation of the gel by using fluorescent beads embedded just below the surface of the gel. In this article we describe the image correlation method (ICM) optimized for determining the displacement field of the gel under a contracting cell. For the calculation of the traction field from the displacement field we use the recently developed method of Fourier transform traction cytometry (FTTC). The ICM and FTTC methods are applied to human airway smooth muscle cells during stimulation with the contractile agonist histamine or the relaxing agonist isoproterenol. The overall intensity of the cell contraction (the median traction magnitude, the energy transferred from the cell to the gel, and the net contractile moment) increased after activation with histamine, and decreased after treatment with isoproterenol. Cells exhibited regional differences in the time course of traction during the treatment. Both temporal evolution and magnitude of traction increase induced by histamine varied markedly among different cell protrusions, whereas the nuclear region showed the smallest response. These results suggest that intracellular mediators of cell adhesion and contraction respond to contractile stimuli with different rates and intensities in different regions of the cell.

NASA Discipline Cell Biology↗

An explicit, energy-conserving particle-in-cell scheme

We present an explicit temporal discretization of particle-in-cell schemes for the non-relativistic Vlasov equation that results in exact energy conservation when combined with an appropriate spatial discretization. The scheme is inspired by a simple, second-order explicit scheme that conserves energy exactly in the Eulerian context. We show that direct translation to particle-in-cell does not result in strict conservation, but derive a simple correction based on an analytically solvable optimization problem that recovers conservation. While this optimization problem is not guaranteed to have a real solution for every particle, we provide a correction that makes imaginary values extremely rare and still admits $\mathcal{O}$(10 –12 ) fractional errors in energy for practical simulation parameters. We present the scheme in both electrostatic – where we use the Ampère formulation – and electromagnetic contexts. With an electromagnetic field solve, the field update is most naturally linearly implicit, but the more computationally intensive particle update remains fully explicit. Here, we also show how the scheme can be extended to use the fully explicit leapfrog and pseudospectral analytic time-domain (PSATD) field solvers. The scheme is tested on standard kinetic plasma problems, confirming its conservation properties.

Energy conservation↗

Pilot microbial electrolysis cell closes the hydrogen loop for hydrothermal wet waste conversion to jet fuel

The global shift toward net-zero emissions necessitates resource recovery from wet waste. In this study, we demonstrate the first feasibility of combining pilot-scale microbial electrolytic cells (MECs) with hydrothermal liquefaction (HTL) for simultaneous post-hydrothermal liquefaction wastewater (PHW) treatment and efficient hydrogen (H₂) production to meet biocrude upgrading requirements. Long-term single reactor operation revealed that fixed anode potential enabled rapid startup, and low catholyte pH and high salinity were effective in suppression of cathodic methanogenesis and acetogenesis – resulting in high current density of 16.6 A m –2 and 9.3 A m –2 when feeding synthetic wastewater and PHW respectively. Additionally, the anode biofilm exhibited spatial variations in response to local environmental conditions. In conclusion, onsite parallel or serial operations of multiple MECs showed good performance using actual PHW with a record-high H 2 production rate of 0.5 L L R day –1 for MEC over 10 liters scale, and the optimal chemical oxygen demand (COD)-to-H 2 yield reached 0.127 kg-H 2 per kg-COD, supporting a self-sufficient, closed-loop upgrade to jet fuel.

08 HYDROGEN↗

Inverse-Designed Photonics for Semiconductor Foundries

Silicon photonics is becoming a leading technology in photonics, displacing traditional fiber optic transceivers and enabling new applications. Further improving the density and performance of silicon photonics, however, has been challenging due to the large size and limited performance of traditional semianalytically designed components. Automated optimization of photonic devices using inverse design is a promising path forward but has, until now, faced difficulties in producing designs that can be fabricated reliably at scale. In this paper we experimentally demonstrate four inverse-designed devices made successfully in a commercial silicon photonics foundry: a spatial mode multiplexer, wavelength demultiplexer, 50–50 directional coupler, and 3-way power splitter. These devices are efficient, robust to fabrication variability, and compact, with footprints only a few micrometers across. They pave the way forward for the widespread practical use of inverse design.

36 MATERIALS SCIENCE↗

Visualizing heterogeneous dipole fields by terahertz light coupling in individual nano-junctions

The challenge underlying superconducting quantum computing is to remove materials bottleneck for highly coherent quantum devices. The nonuniformity and complex structural components in the underlying quantum circuits often lead to local electric field concentration, charge scattering, dissipation and ultimately decoherence. Here we visualize interface dipole heterogeneous distribution of individual Al/AlO$_{x}$/Al junctions employed in transmon qubits by broadband terahertz scanning near-field microscopy that enables the non-destructive and contactless identification of defective boundaries in nano-junctions at an extremely precise nanoscale level. Our THz nano-imaging tool reveals an asymmetry across the junction in electromagnetic wave-junction coupling response that manifests as hot (high intensity) vs cold (low intensity) spots in the spatial electrical field structures and correlates with defected boundaries from the multi-angle deposition processes in Josephson junction fabrication inside qubit devices. The demonstrated local electromagnetic scattering method offers high sensitivity, allowing for reliable device defect detection in the pursuit of improved quantum circuit fabrication for ultimately optimizing coherence times.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The preferred-mode coherent structure in the near field of an axisymmetric jet with and without excitation

Visualization and phase-locked hot-wire measurements have been used to investigate the preferred-mode coherent structure in the near field of an axisymmetric jet for different Reynolds numbers and for the two limiting asymptotic states of the exit boundary layer: laminar and fully developed turbulent. Distributions of phase-averaged longitudinal and lateral velocities, coherent vorticity, coherent and incoherent turbulence Reynolds stresses, and coherent strain rate and production have been obtained over the spatial extent of the coherent structure for a Strouhal number of 0.3. It is found that these properties are independent of the initial conditions, but weakly dependent on the Reynolds number up to about 100,000. The naturally occurring coherent structure in the unexcited jet, educed via an optimized conditional sampling scheme, shows no clear dependence on the jet Reynolds number up to 800,000.

Hussain, A. K. M. F.↗

Polarimetric Multispectral Imaging Technology

The Jet Propulsion Laboratory is developing a remote sensing technology on which a new generation of compact, lightweight, high-resolution, low-power, reliable, versatile, programmable scientific polarimetric multispectral imaging instruments can be built to meet the challenge of future planetary exploration missions. The instrument is based on the fast programmable acousto-optic tunable filter (AOTF) of tellurium dioxide (TeO2) that operates in the wavelength range of 0.4-5 microns. Basically, the AOTF multispectral imaging instrument measures incoming light intensity as a function of spatial coordinates, wavelength, and polarization. Its operation can be in either sequential, random access, or multiwavelength mode as required. This provides observation flexibility, allowing real-time alternation among desired observations, collecting needed data only, minimizing data transmission, and permitting implementation of new experiments. These will result in optimization of the mission performance with minimal resources. Recently we completed a polarimetric multispectral imaging prototype instrument and performed outdoor field experiments for evaluating application potentials of the technology. We also investigated potential improvements on AOTF performance to strengthen technology readiness for applications. This paper will give a status report on the technology and a prospect toward future planetary exploration.

Cheng, L.-J.↗