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,513 records · Page 84

Planar Investigation of a CobraMRV Reentry Flowfield Using Pulse-Burst, Cross-Correlation DGV

The flowfield around a Co-Optimization Blunt-body Reentry Analysis Mid-lift/drag-ratio Rigid Vehicle (CobraMRV) vehicle is investigated with a combination of high-speed planar laser Mie scattering (PLMS) and pulse-burst cross-correlation Doppler global velocimetry (PB-CC-DGV). Tests were conducted in the NASA Langley 4-foot Supersonic Unitary Plan Wind Tunnel (UPWT) over a range of different tunnel operating conditions and model configurations. Results indicate a complex shock-boundary layer interaction. Scalar information extracted from the PLMS show the evolution of the bow shock structure, while streamwise velocity measurements indicate the spatial evolution of the shock-boundary layer interaction including the growth of the separation shock foot and eventual reacceleration of the flow at farther downstream locations. Assessment of multiple cases show strong Mach and Reynolds numbers driven effects on the character of the shock-boundary layer interactions. Measurement uncertainties ranged from 50 to 150 m/s throughout the region of interest, driven largely by angular uncertainties and instabilities in the laser pointing. The mean accuracy of the freestream measurements was found to be 5.2-percent of tunnel predicted values.

Laser↗

Direct observation of 3D topological spin textures and their interactions using soft x-ray vector ptychography

Magnetic topological defects are energetically stable spin configurations characterized by symmetry breaking. Vortices and skyrmions are two well-known examples of 2D spin textures that have been actively studied for both fundamental interest and practical applications. However, experimental evidence of the 3D spin textures has been largely indirect or qualitative to date, due to the difficulty of quantitively characterizing them within nanoscale volumes. Here, we develop soft x-ray vector ptychography to quantitatively image the 3D magnetization vector field in a frustrated superlattice with 10 nm spatial resolution. By applying homotopy theory to the experimental data, we quantify the topological charge of hedgehogs and anti-hedgehogs as emergent magnetic monopoles and probe their interactions inside the frustrated superlattice. We also directly observe virtual hedgehogs and anti-hedgehogs created by magnetically inert voids. We expect that this new quantitative imaging method will open the door to study 3D topological spin textures in a broad class of magnetic materials. Our work also demonstrates that magnetically frustrated superlattices could be used as a new platform to investigate hedgehog interactions and dynamics and to exploit optimized geometries for information storage and transport applications.

Rana, A↗

Bayesian inference in band excitation scanning probe microscopy for optimal dynamic model selection in imaging

The universal tendency in scanning probe microscopy (SPM) over the last two decades is to transition from simple 2D imaging to complex detection and spectroscopic imaging modes. The emergence of complex SPM engines brings forth the challenge of reliable data interpretation, i.e., conversion from detected signals to descriptors specific to tip–surface interactions and subsequently to material’s properties. In this work, we implemented a Bayesian inference approach for the analysis of the image formation mechanisms in band excitation SPM. Compared to the point estimates in classical functional fit approaches, Bayesian inference allows for the incorporation of extant knowledge of materials and probe behavior in the form of corresponding prior distribution and return the information on the material functionality in the form of readily interpretable posterior distributions. We explore the nonlinear mechanical behaviors spatially in a classical ferroelectric material, PbTiO 3 . We observe the non-trivial evolution of the Duffing stiffness term and the nonlinearity of the sample surface, determine spatial clustering of the nonlinear response, and perform a Landau analysis on predicting the nonlinear coefficient, which indicates that ferroelectric behavior can be a cause of the observed results. These observations suggest that the spectrum of anomalous behaviors at the ferroelectric domain walls may be broader than previously believed and can extend to non-conventional mechanical properties in addition to static and microwave conductance.

36 MATERIALS SCIENCE↗

Photocatalytic Overall Water Splitting at the Integrated Rh–MoRhO x Cluster Heterostructure on InGaN/GaN Nanowires

The quest for efficient solar-driven water splitting, a promising avenue for clean fuel production, faces challenges due to limited solar energy conversion efficiency. Traditional approaches study the overall water splitting as two spatially separate half reactions on two unrelated sites, hindering full utilization of photogenerated charge and water molecules. To overcome these limitations, an integrated cluster heterostructure catalyst on InGaN/GaN semiconductor nanowires is proposed for the effective utilization of photogenerated charge carriers and water molecules on the same redox localization. By establishing the fast charge extraction kinetics based on InGaN/GaN nanowires, the integration of Rh and MoRhOx clusters on the nanowire surface enables simultaneous and fast hydrogen/oxygen evolution reactions at the cluster heterostructure. Furthermore, the integrated strategy can enhance the charge redistribution across the heterostructure between the two clusters, further optimizing adsorption of reaction intermediates on each cluster for boosted photocatalytic water splitting activity. Consequently, the integrated heterostructure triggers a 40-fold increased hydrogen production efficiency in an artificial leaf system. This study provides valuable insights for the rational design of advanced heterostructured photocatalysts for water splitting and beyond.

GaN nanowire↗

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗

Characterizing Stress Orientations in Southern Kansas

We repot that Induced seismicity predominantly occurs along faults that are optimally oriented to the local principal compressive stress direction, and the characterization of these stress orientations is an important component of understanding seismic hazards. The seismicity rate in southern Kansas rapidly increased in 2013 primarily due to the disposal of large volumes of wastewater into the Arbuckle Group. Previously, local stress orientations in this area were poorly constrained, which limited our understanding of the complex faulting and diverse earthquake mechanisms in this region. We use shear-wave splitting and focal mechanism inversion techniques to create multiple, independent estimates of maximum horizontal stress directions (S Hmax ) and their spatial variation across the study area. We then create an integrated model of stress orientations for southern Kansas and northern Oklahoma using our local results in conjunction with previous, regional stress orientation estimates. We find that S Hmax in both southern Kansas and central Oklahoma exhibits an east-northeast (~N78°E) orientation, and these regions bound a northeast (~N59°E) rotation within a ~20km area in northern Oklahoma near the Nemaha ridge.

58 GEOSCIENCES↗

Spatial filters for shape control

Recently there has emerged a new class of sensors, called spatial filters, for structures which respond over a significant gauge length. Examples include piezoelectric laminate PVDF film, modal domain optical fiber sensors, and holographic sensors. These sensors have a unique capability in that they can be fabricated to locally alter their sensitivity to the measurand. In this paper we discuss how these sensors can be used for the implementation of control algorithms for the suppression of acoustic radiation from flexible structures. Based on this relationship between the total power radiated to the far field to the modal velocities of the structure, we show how the sensor placement to optimize the control algorithm to suppress the radiated power.

Lindner, Douglas K.↗

Romie: A Domain-Independent Tool for Computer-Aided Robust Operations Management

Romie is a decision support tool based on AI's latest advances in the domain of robust scheduling. Unlike all its predecessors, the tool allows to (i) visually model the operational problem and context entirely (ii) optimize to find near-optimal schedules while taking uncertainty into account and deals with (iii) a combination of various {key performance indicators (KPIs). It comes with a web user interface. Part or all of the modelled activities may be associated to random variables describing their stochastic durations, in order to produce schedules that are robust w.r.t. temporal uncertainty. Hence, depending on the pursued KPIs, the schedules maximize a combination of the following terms: the probability of satisfying the problem constraints, the expected return/efficiency, the expected outcome quality, and even the operators' wellness by minimizing its expected extra-hours. Initially developed for spatial exploration and demonstration in the context of Mars analog missions, this versatile tool is here applied to operations management in both biotechnology manufacturing and robots parametrization in a cave exploration context.

Chien, Steve A.↗

Orbit Optimization for the Geospace Electrodynamics Connections (GEC) Mission

The Geospace Electrodynamics Connections (GEC) mission plan is to launch multiple spacecraft to perform in-situ atmospheric science in the lower ionosphere. There is limited experience in this low altitude region with the Atmospheric Explorer-C (AE-C) being the last spacecraft to explore this region in 1973. AE-C flew an eccentric orbit using maneuvers to lower its perigee to near 130 km at various times during its mission. GEC will advance the science performed by AE-C by performing multiple low-perigee, atmospheric dipping campaigns for extended durations. AE-C kept its perigee near 130 km for only a total of roughly 1 day. Furthermore, GEC plans to carry a more diverse suite of instruments and will be able to capture different temporal and spatial phenomena through the use of multiple spacecraft flying in a string of pearls formation. The mission analysis for GEC has been broken into two parts: the analysis of the parking orbit with the dipping campaigns and the examination of the multi-satellite dynamics of the GEC constellation. The analysis described in this paper examines the capability to meet the requirements necessary to support the 10 dipping campaigns using a single spacecraft as a representative of all three in the constellation. Further analysis is being performed to analyze the multi-satellite nature of the GEC mission.

Mesarch, Michael A.↗

Recursive inverse kinematics for robot arms via Kalman filtering and Bryson-Frazier smoothing

This paper applies linear filtering and smoothing theory to solve recursively the inverse kinematics problem for serial multilink manipulators. This problem is to find a set of joint angles that achieve a prescribed tip position and/or orientation. A widely applicable numerical search solution is presented. The approach finds the minimum of a generalized distance between the desired and the actual manipulator tip position and/or orientation. Both a first-order steepest-descent gradient search and a second-order Newton-Raphson search are developed. The optimal relaxation factor required for the steepest descent method is computed recursively using an outward/inward procedure similar to those used typically for recursive inverse dynamics calculations. The second-order search requires evaluation of a gradient and an approximate Hessian. A Gauss-Markov approach is used to approximate the Hessian matrix in terms of products of first-order derivatives. This matrix is inverted recursively using a two-stage process of inward Kalman filtering followed by outward smoothing. This two-stage process is analogous to that recently developed by the author to solve by means of spatial filtering and smoothing the forward dynamics problem for serial manipulators.

Rodriguez, G.↗

Leonid Storm Flux Analysis From One Leonid MAC Video AL50R

A detailed meteor flux analysis is presented of a seventeen-minute portion of one videotape, collected on November 18, 1999, during the Leonid Multi-instrument Aircraft Campaign. The data was recorded around the peak of the Leonid meteor storm using an intensified CCD camera pointed towards the low southern horizon. Positions of meteors on the sky were measured. These measured meteor distributions were compared to a Monte Carlo simulation, which is a new approach to parameter estimation for mass ratio and flux. Comparison of simulated flux versus observed flux levels, seen between 1:50:00 and 2:06:41 UT, indicate a magnitude population index of r = 1.8 +/- 0.1 and mass ratio of s = 1.64 +/- 0.06. The average spatial density of the material contributing to the Leonid storm peak is measured at 0.82 +/- 0.19 particles per square kilometer per hour for particles of at least absolute visual magnitude +6.5. Clustering analysis of the arrival times of Leonids impacting the earth's atmosphere over the total observing interval shows no enhancement or clumping down to time scales of the video frame rate. This indicates a uniformly random temporal distribution of particles in the stream encountered during the 1999 epoch. Based on the observed distribution of meteors on the sky and the model distribution, recommendations am made for the optimal pointing directions for video camera meteor counts during future ground and airborne missions.

Gural, Peter S.↗

The Aquarius ocean salinity mission high stability L-band radiometer

The NASA Earth Science System Pathfinder (ESSP) mission Aquarius, will measure global ocean surface salinity with ~120 km spatial resolution every 7-days with an average monthly salinity accuracy of 0.2 psu (parts per thousand) [1]. This requires an L-band low-noise radiometer with the long-term calibration stability of ≤0.15 K over 7 days. The instrument utilizes a push-broom configuration which makes it impractical to use a traditional warm load and cold plate in front of the feedhorns. Therefore, to achieve the necessary performance Aquarius utilizes a Dicke radiometer with noise injection to perform a warm – hot calibration. The radiometer sequence between antenna, Dicke load, and noise diode has been optimized to maximize antenna observations and therefore minimize NEDT. This is possible due the ability to thermally control the radiometer electronics and front-end components to 0.1 °Crms over 7 days.

Freedman, Adam↗

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗

On the theoretical limitations in estimating thickness of a plate-like structure from a full-field single-tone response Lamb wave measurement

A persistent question in wavenumber analysis in the estimation of thickness from a steady-state wave field is identifying the theoretical sensitivity of the system. Here, a well-known trade-off between spatial frequency/wavenumber resolution and thickness resolution exists. The current work presents a calculation of the Cramer-Rao Lower bound (CRLB), specifically as applied to thickness estimates, for a 2-dimensional multi-mode waveform in Gaussian noise. Cases of near-field and far-field excitation are considered, and transducer position with respect to the scan area is also varied. Additionally, we consider the CRLB in a plate with multiple sources, simulated as sources placed on the boundary of the plate. We conclude by presenting the CRLB values in terms of frequency for various thicknesses, and by presenting optimal excitation frequencies for a nominal thickness, based on the CRLB.

47 OTHER INSTRUMENTATION↗

Toward the fabrication of a 5- μm -resolution Wolter microscope for the National Ignition Facility (invited)

Advancements in computer-controlled polishing, metrology, and replication have led to an x-ray mirror fabrication process that is capable of producing high-resolution Wolter microscopes. Here, we present the fabrication and test of a nickel–cobalt replicated full-shell x-ray mirror that was electroformed from a finely figured and polished mandrel. This mandrel was designed for an 8-m source-to-detector-distance microscope, with 10× magnification, and was optimized to reduce shell distortions that occur within 20 mm of the shell ends. This, in combination with an improved replication tooling design and refined bath parameters informed by a detailed COMSOL Multiphysics® model, has led to reductions in replication errors in the mirrors. Mandrel surface fabrication was improved by implementing a computer-controlled polishing process that corrected the low-frequency mandrel figure error and achieved <2.0 nm RMS convergence error. X-ray tests performed on a pair of mirror shells replicated from the mandrel have demonstrated <10 μm full-width at half-maximum (FWHM) spatial resolution. Here, we discuss the development process, highlight results from metrology and x-ray testing, and define a path for achieving a program goal of 5 μm FWHM resolution.

47 OTHER INSTRUMENTATION↗

Measurement of Deformation Heterogeneity During Shale Swelling Using Digital Image Correlation

Abstract Rock-fluid interactions in shale formations are one of the main sources of wellbore instability issues and inadequate stimulation performance. For better planning of fracturing operations and optimizing production from these reservoirs, it is important to understand the mechanisms behind these interactions. These issues are especially prevalent in clay-rich shales, such as Tuscaloosa Marine Shale, which is the subject of this study. Conventional techniques to quantify these shale-fluid interactions comprise of measuring swelling in powdered rock grains or measurement of deformation in the whole pieces of the core using a linear variable differential transformer and strain gages. However, the contribution from individual laminae to overall deformation cannot be evaluated using these methods. In this study, we developed an experimental setup to evaluate the spatial deformation in shale during interaction with water using digital image correlation (DIC). Deformation of two shale samples, with 34 wt% to 51 wt% clay content, was studied. White paint was used to generate a random speckle pattern on the specimen and then immersed in deionized water. The deformation process was captured using a digital camera and images were analyzed using DIC to quantify the deformation. The implementation of the DIC technique enables the visualization and quantification of spatial deformation in the specimen during interacting with water. The results show the localization of large strains in select laminations. The results provide a better understanding of shale deformation when interacting with water in comparison to traditional measurements that can provide only an average strain value.

Energy & Fuels↗

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗

An a priori evaluation of a principal component and artificial neural network based combustion model in diesel engine conditions

A principal component analysis (PCA) and artificial neural network (ANN) based chemistry tabulation approach is presented. ANNs are used to map the thermochemical state onto a low-dimensional manifold consisting of five control variables that have been identified using PCA. Three canonical configurations are considered to train the PCA-ANN model: a series of homogeneous reactors, a nonpremixed flamelet, and a two-dimensional lifted flame. The performance of the model in predicting the thermochemical manifold of a spatially-developing turbulent jet flame in diesel engine thermochemical conditions is a priori evaluated using direct numerical simulation (DNS) data. The PCA-ANN approach is compared with a conventional tabulation approach (tabulation using ad hoc defined control variables and linear interpolation). The PCA-ANN model provides higher accuracy and requires several orders of magnitude less memory. Here, these observations indicate that the PCA-ANN model is superior for chemistry tabulation, especially for modelling complex chemistries that present multiple combustion modes as observed in diesel combustion. The performance of the PCA-ANN model is then compared to the optimal estimator, i.e. the conditional mean from the DNS. The results indicate that the PCA-ANN model gives high prediction accuracy, comparable to the optimal estimator, especially for major species and the thermophysical properties. Higher errors are observed for the minor species and reaction rate predictions when compared to the optimal estimator. It is shown that the prediction of minor species and reaction rates can be improved by using training data that exhibits a variation of parameters as observed in the turbulent flame. The output of the ANN is analysed to assess mass conservation. It is observed that the ANN incurs a mean absolute error of 0.05% in mass conservation. Furthermore, it is demonstrated that this error can be reduced by modifying the cost function of the ANN to penalise for deviation from mass conservation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗