Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse problems”

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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Behavior of an inversion-based precipitation retrieval algorithm with high-resolution AMPR measurements including a low-frequency 10.7-GHz channel

A microwave-based, profile-type precipitation retrieval algorithm has been used to analyze high-resolution passsive microwave measurements over an ocean background, obtained by the Advanced Microwave Precipitation Radiometer (AMPR) flown on a NASA ER-2 aircraft. The analysis is designed to first determine the improvements that can be gained by adding brightness temperature information from the AMPR low-frequency channel (10.7 GHz) to a multispectral retrieval algorithm nominally run with satellite information at 19, 37, and 85 GHz. The impact of spatial resolution degradation of the high-resolution brightness temperature information on the retrieved rain/cloud liquid water contents and ice water contents is then quantified in order to assess the possible biases inherent to satellite-based retrieval. Careful inspection of the high-resolution aircraft dataset reveals five distinctive brightness temperature features associated with cloud structure and scattering effects that are not generally detectable in current passive microwave satellite measurements. Results suggest that the inclusion of 10.7-GHz information overcomes two basic problems associated with three-channel retrieval. Intercomparisons of retrievals carried out at high-resolution and then averaged to a characteristic satellite scale to the corresponding retrievals in which the brightness temperatures are first convolved down to the satellite scale suggest that with the addition of the 10.7-GHz channel, the rain liquid water contents will not be negatively impacted by special resolution degradation. That is not the case with the ice water contents as they appear ti be quite sensitive to the imposed scale, the implication being that as spatial resolution is reduced, ice water contents will become increasingly underestimated.

Smith, E. A.↗
DOE OSTI2024

PHASM: A Toolkit for Creating AI Surrogate Models within Legacy Codebases

PHASM (“Parallel Hardware viA Surrogate Models”) is a software toolkit for creating AI-based surrogate models of scientific code. AI-based surrogate models are widely used for creating fast and inverse simulations. PHASM anticipates an additional future use case: adapting legacy code to modern hardware. While data centers are investing in heterogeneous hardware such as GPUs and FPGAs, many established scientific codebases remain unable to take advantage of the hardware’s higher parallelism without undergoing a costly rewrite. An alternative is to train a AI-based surrogate model to mimic computationally intensive functions in the code, and run the surrogate instead. PHASM formalizes a development lifecycle for such surrogate models, including discovering functions amenable to replacement with a surrogate model, predicting the resulting performance, identifying the function’s space of inputs and outputs, binding the model to the code, and managing model versions. A suite of software tools for facilitating these steps was written and validated against a set of model problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces

Effective numerical treatment of multicomponent viscous flow problems involving the advection of sharp interfaces between materials of differing physical properties requires correction techniques to prevent spurious diffusion and dispersion. We develop a particular algorithm, based on modern shock-capture techniques, employing a two-step nonlinear method. The first step involves the global application of a high-order upwind scheme to a hyperbolic advection equation used to model the distribution of distinct material components in a flow field. The second step is corrective and involves the application of a global filter designed to remove dispersion errors that result from the advection of discontinuities (e.g., material interfaces) by high-order, minimally dissipative schemes. The filter introduces no additional diffusion error. Nonuniform viscosity across a material interface is allowed for by the implementation of a compositionally weighted-inverse interface viscosity scheme. The combined method approaches the optimal accuracy of modern shock-capture techniques with a minimal increase in computational time and memory. A key advantage of this method is its simplicity to incorporate into preexisting codes be they finite difference, element, or volume of two or three dimensions.

Lenardic, A.↗
DOE OSTI2025

CO 2 Adsorption and Hydrogenation on Inverse InO x /Cu(111) Catalysts: Active Role of the Oxide–Metal Interface

The direct conversion of carbon dioxide (CO 2 ) into methanol via hydrogenation is essential for industrial applications. Recent studies on catalysts that contain an inverse oxide/metal configuration have shown very good catalytic performance for the CO 2 hydrogenation to methanol process. Here, in this study, we investigated the behavior of indium oxide-Cu(111) interfaces under pure CO 2 and CO 2 /H 2 mixtures using synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS). Initially, a single layer of copper oxide (Cu x O) was grown on the Cu(111) surface by controlled oxidation. On this surface, indium was deposited at room temperature. Oxygen atoms transferred from Cu x O/Cu(111) to the indium metal upon deposition, forming In-O-Cu bonds and active interfaces. Although Cu(111) is not very active for the binding and activation of CO 2 , the formed InO x -Cu(111) interfaces had no problem adsorbing and dissociating the molecule at room temperature. Reaction of CO 2 with H 2 on InO x -Cu(111) yielded surface-bound H 3 CO, CO 2 δ− , CO 3 , and CH x species that are typical intermediates in the production of methanol and other oxygenates. The InO x -Cu(111) interface underwent dynamic chemical changes under reaction conditions, forming In-Cu alloys at low indium coverages (< 0.05 monolayer), while at higher indium coverages a mixture of an In-Cu and InO x was detected in XPS. These findings indicate that InO x /In-Cu interfaces can play a key role in processes aimed at the trapping and valorization of CO 2 .

36 MATERIALS SCIENCE↗
NASA NTRS1993

Research and development of high thermal stability fuels

Increases in aircraft performance are leading to increases in the thermal stress on the primary aircraft coolant--the fuel. Fuel thermal stability limitations may offset future aircraft performance gains. The Air Force's Wright Laboratory is sponsoring several research programs to address this problem. The development of an additive package for JP-8 to improve its thermal stability is the primary focus of this paper. This program involves extensive testing of fuels and additives in a variety of test devices, culminating in tests in a fuel system simulator and engine tests. These tests involve Air Force personnel, on-site contractors (University of Dayton Research Institute, Systems Research Laboratories), Pratt and Whitney Aircraft Co., additive manufacturers, and Sandia National Laboratory. The test devices include several flowing and static tests, where the behavior of a fuel is investigated in a wide variety of environments. The study of several baseline fuels in these devices has led to some new insights into the mechanisms of fuel thermal (in)stability. It is becoming clear that a fuel's tendency to oxidize (to form peroxides, for example) is often inversely proportional to its tendency to form insoluble deposits.

T Edwards↗
NASA NTRS1993

Cerebral blood velocity and other cardiovascular responses to 2 days of head-down tilt

Spaceflight induces a cephalad redistribution of fluid volume and blood flow within the human body, and space motion sickness, which is a problem during the first few days of space flight, could be related to these changes in fluid status and in blood flow of the cerebrum and vestibular system. To evaluate possible changes in cerebral blood flow during simulated weightlessness, we measured blood velocity in the middle cerebral artery (MCA) along with retinal vascular diameters, intraocular pressure, impedance cardiography, and sphygmomanometry on nine men (26.2 +/- 6.6 yr) morning and evening for 2 days during continuous 10 deg head-down tilt (HDT). When subjects went from seated to head-down bed rest, their heart rate and retinal diameters decreased, and intraocular pressures increased. After 48 h of HDT, blood flow velocity in the MCA was decreased and thoracic impedance was increased, indicating less fluid in the thorax. Percent changes in blood flow velocities in the MCA after 48 h of HDT were inversely correlated with percent changes in retinal vascular diameters. Blood flow velocities in the MCA were inversely correlated (intersubject) with arterial pressures and retinal vascular diameters. Heart rate, stroke volume, cardiac output, systolic arterial pressure, and at times pulse pressure and blood flow velocities in the MCA were greater in the evening. Total peripheral resistance was higher in the morning. Although cerebral blood velocity is reduced after subjects are head down for 2 days, the inverse relationship with retinal vessel diameters, which have control analogous to that of cerebral vessels, indicates cerebral blood flow is not reduced.

Frey, Mary A. B.↗
DOE OSTI2021

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

Mapping From an Instrumented Glove to a Robot Hand

An algorithm has been developed to solve the problem of mapping from (1) a glove instrumented with joint-angle sensors to (2) an anthropomorphic robot hand. Such a mapping is needed to generate control signals to make the robot hand mimic the configuration of the hand of a human attempting to control the robot. The mapping problem is complicated by uncertainties in sensor locations caused by variations in sizes and shapes of hands and variations in the fit of the glove. The present mapping algorithm is robust in the face of these uncertainties, largely because it includes a calibration sub-algorithm that inherently adapts the mapping to the specific hand and glove, without need for measuring the hand and without regard for goodness of fit. The algorithm utilizes a forward-kinematics model of the glove derived from documentation provided by the manufacturer of the glove. In this case, forward-kinematics model signifies a mathematical model of the glove fingertip positions as functions of the sensor readings. More specifically, given the sensor readings, the forward-kinematics model calculates the glove fingertip positions in a Cartesian reference frame nominally attached to the palm. The algorithm also utilizes an inverse-kinematics model of the robot hand. In this case, inverse-kinematics model signifies a mathematical model of the robot finger-joint angles as functions of the robot fingertip positions. Again, more specifically, the inverse-kinematics model calculates the finger-joint commands needed to place the fingertips at specified positions in a Cartesian reference frame that is attached to the palm of the robot hand and that nominally corresponds to the Cartesian reference frame attached to the palm of the glove. Initially, because of the aforementioned uncertainties, the glove fingertip positions calculated by the forwardkinematics model in the glove Cartesian reference frame cannot be expected to match the robot fingertip positions in the robot-hand Cartesian reference frame. A calibration must be performed to make the glove and robot-hand fingertip positions correspond more precisely. The calibration procedure involves a few simple hand poses designed to provide well-defined fingertip positions. One of the poses is a fist. In each of the other poses, a finger touches the thumb. The calibration subalgorithm uses the sensor readings from these poses to modify the kinematical models to make the two sets of fingertip positions agree more closely.

Goza, Michael↗
DOE OSTI2022

Quasielastic least-squares reverse time migration of PS reflections

Conventional analysis of amplitude variation with offset for elastic PS reflections is based on analytical reflection coefficients in a layered medium, and wave-equation-based PS migration is mainly used to produce a structural image. To overcome this problem, we have developed a least-squares reverse time migration (LSRTM) method for elastic PS reflections based on a quasielastic wave equation. Here, the quasielastic wave equation can accurately model PS reflections with elastic amplitudes under the first-order Born approximation. Our LSRTM method inverts for perturbations of the S-wave velocity and density by minimizing the L 2 norm of the difference between recorded and predicted PS reflections modeled using a quasielastic wave equation. We refer to our new method as quasielastic LSRTM of PS reflections. Numerical tests on synthetic and field data indicate that our method can properly handle the amplitudes of elastic PS reflections and provides an accurate estimate of the perturbations of S-wave velocity and density. Extending the method to the 3D case is not straightforward and might require incorporating certain PS data processing techniques into the inversion itself.

58 GEOSCIENCES↗
NASA NTRS1998

Lunar Basins: New Evidence from Gravity for Impact-Formed Mascons

The prominent gravity highs (mascons) associated with uncompensated mass anomalies in lunar mare basins are a dramatic expression of the present-day rigidity of the lunar lithosphere. First discovered in Lunar Orbiter tracking data, these about 350-mGal gravity highs have been redetermined from the analysis of Clementine and historical tracking. These highs coincide with topographic lows, indicating nonisostatic support. One of the rediscoveries of this analysis is the encirclement of the highs by substantial negative anomalies over topographic highs. Recent gravity fields are providing the increased resolution necessary to determine the causes of this unique mascon signature. The compensation of the basin anomalies remains controversial. The mascon highs have long been interpreted as the result of mare loading, subsequent to the decay of residual stresses resulting from the impact. Substantially more mare fill is required to produce mascon highs than has been inferred on geological grounds, and the amount of near-surface mass deficit required to produce a gravity most exceeds bounds inferred from terrestrial examples. This problem is most acute for the youngest basin, Orientale. Recent gravity fields from Lunar Prospector have suggested mascon highs associated with nonmare basins such as Mendel-Rydberg, or minimally filled basins like Humboldtianum, further calling this explanation into question. We suggest that the mascon gravity signal is produced by a combination of crustal thickness changes, manifested by central mantle uplift, outward displacement of crust, and downward flexure of the lithosphere under mare loading. The mantle uplift is superisostatic, maintained by residual stresses resulting from the process of impact cratering and modification. In particular, the process of crater collapse and mantle rebound terminates abruptly, leaving the mantle plug in a non-equilibrium state, surrounded by a ring of thickened crust. Viscous relaxation over geological timescales has erased some but not all of the signature of the impact process. Mantle uplift inferred from gravity modeling is inversely correlated with age. While the oldest basins such as South Pole Aitken are mainly compensated isotatically, the younger basins appear to have been in a state of superisostatic loading prior to mare emplacement. If this is true, this places an important constraint on the impact process at basin scales. The idea that rebound of the transient crater via acoustic fluidization may freeze substantial stresses imcomplete to this day, may be tested by examining the gravity signatures of major basins on terrestrial bodies. The moon provides the clearest resolved examples to date, but uncertainty in gravity knowledge remains problematic.

Neumann, Gregory A.↗
DOE OSTI2022

Neumann Series in MGS-GMRES and Inner-Outer Iterations: Preprint

A low-synchronization MGS-GMRES Krylov solver employing a truncated Neumann series for the inverse compact WY MGS correction matrix T is presented. A corollary to the backward stability result of Paige et al. [1] establishes that T = I - Lk is sufficient for convergence of GMRES when kLkp F = O("p)_p F (B), where the strictly lower triangular matrix L is defined by the inner products of Krylov vectors V T 1:k-2 vk-1. The preconditioner is the classical Ruge-Stuben AMG algorithm with compatible relaxation and inner-outer Gauss-Seidel smoother. This smoother may also be expressed as a truncated Neumann series. Drop tolerances are applied to the lower triangular matrices arising in the smoother in order to reduce the number of non-zeros and accelerate the time to solution. The number of small matrix elements are found to increase from fine to coarse levels and thus the effciency gains are greater for large problems with many levels in the V -cycle. The solver is applied to the pressure continuity equation for the incompressible Navier-Stokes equations. Unlike the inner-outer iteration, the solver convergence rate with the standard Gauss-Seidel smoother deteriorates with dropping. The solver compute time is reduced by up to 50% without a change in the convergence rate.

Gauss-Seidel smoother↗
DOE OSTI2020

Singlet fission

Singlet fission (SF) is a process whereby a photo-excited singlet exciton spontaneously splits into lower energy triplet excitons. This process has been known to occur in select organic materials for over 50 years but in the past was viewed as a curiosity and not as a problem worthy of focused attention. Over the past decade, this viewpoint has changed as the promising use of SF to improve photovoltaic efficiency has spurred vast activity, which has revealed the richness of the phenomenon. Currently, our understanding of the role played by various bright and dark states and their couplings to vibrations and to each other and by chemical composition and packing, among other factors, has enabled a large increase in the number of viable SF systems as well as a far deeper understanding of the process itself. These advances have allowed for the demonstration of proof-of-concept increases in the solar cell efficiency via the use of SF materials in simple tandem cells and as sensitizers in conjunction with conventional photovoltaics. Indeed, concepts surrounding the field have now invaded other fields, with a prominent example being the use of photon upconversion associated with triplet–triplet fusion (the inverse process to SF) in organic chemistry and photocatalysis. The 23 research articles in this special topic showcase recent cutting-edge advances in the field and highlight the diverse new directions being explored in the field of SF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Reducing measurement costs by recycling the Hessian in adaptive variational quantum algorithms

Abstract Adaptive protocols enable the construction of more efficient state preparation circuits in variational quantum algorithms (VQAs) by utilizing data obtained from the quantum processor during the execution of the algorithm. This idea originated with Adaptive Derivative-Assembled Problem-Tailored variational quantum eigensolver (ADAPT-VQE), an algorithm that iteratively grows the state preparation circuit operator by operator, with each new operator accompanied by a new variational parameter, and where all parameters acquired thus far are optimized in each iteration. In ADAPT-VQE and other adaptive VQAs that followed it, it has been shown that initializing parameters to their optimal values from the previous iteration speeds up convergence and avoids shallow local traps in the parameter landscape. However, no other data from the optimization performed at one iteration is carried over to the next. In this work, we propose an improved quasi-Newton optimization protocol specifically tailored to adaptive VQAs. The distinctive feature in our proposal is that approximate second derivatives of the cost function are recycled across iterations in addition to optimal parameter values. We implement a quasi-Newton optimizer where an approximation to the inverse Hessian matrix is continuously built and grown across the iterations of an adaptive VQA. The resulting algorithm has the flavor of a continuous optimization where the dimension of the search space is augmented when the gradient norm falls below a given threshold. We show that this inter-optimization exchange of second-order information leads the approximate Hessian in the state of the optimizer to be consistently closer to the exact Hessian. As a result, our method achieves a superlinear convergence rate even in situations where the typical implementation of a quasi-Newton optimizer converges only linearly. Our protocol decreases the measurement costs in implementing adaptive VQAs on quantum hardware as well as the runtime of their classical simulation.

Ramôa, Mafalda (ORCID:0000000302187801)↗
NASA NTRS2021

Modeling a Coronal Mass Ejection from an Extended Filament Channel. I. Eruption and Early Evolution

We present observations and modeling of the magnetic field configuration, morphology, and dynamics of a large-scale, high-latitude filament eruption observed by the Solar Dynamics Observatory. We analyze the 2015 July 9–10 filament eruption and the evolution of the resulting coronal mass ejection (CME) through the solar corona. The slow streamer-blowout CME leaves behind an elongated post-eruption arcade above the extended polarity inversion line that is only poorly visible in extreme ultraviolet (EUV) disk observations and does not resemble a typical bright flare-loop system. Magnetohydrodynamic (MHD) simulation results from our data-inspired modeling of this eruption compare favorably with the EUV and white-light coronagraph observations. We estimate the reconnection flux from the simulation’s flare-arcade growth and examine the magnetic-field orientation and evolution of the erupting prominence, highlighting the transition from an erupting sheared-arcade filament channel into a streamer-blowout flux-rope CME. Our results represent the first numerical modeling of a global-scale filament eruption where multiple ambiguous and complex observational signatures in EUV and white light can be fully understood and explained with the MHD simulation. In this context, our findings also suggest that the so-called “stealth CME” classification, as a driver of unexpected or “problem” geomagnetic storms, belongs more to a continuum of observable/non-observable signatures than to separate or distinct eruption processes.

Benjamin J Lynch↗
DOE OSTI2026

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator↗
NASA NTRS1991

A semi-automated process for the production of custom-made shoes

A more efficient, cost-effective and timely way of designing and manufacturing custom footware is needed. A potential solution to this problem lies in the use of computer-aided design and manufacturing (CAD/CAM) techniques in the production of custom shoes. A prototype computer-based system was developed, and the system is primarily a software entity which directs and controls a 3-D scanner, a lathe or milling machine, and a pattern-cutting machine to produce the shoe last and the components to be assembled into a shoe. The steps in this process are: (1) scan the surface of the foot to obtain a 3-D image; (2) thin the foot surface data and create a tiled wire model of the foot; (3) interactively modify the wire model of the foot to produce a model of the shoe last; (4) machine the last; (5) scan the surface of the last and verify that it correctly represents the last model; (6) design cutting patterns for shoe uppers; (7) cut uppers; (8) machine an inverse mold for the shoe innersole/sole combination; (9) mold the innersole/sole; and (10) assemble the shoe. For all its capabilities, this system still requires the direction and assistance of skilled operators, and shoemakers to assemble the shoes. Currently, the system is running on a SUN3/260 workstation with TAAC application accelerator. The software elements of the system are written in either Fortran or C and run under a UNIX operator system.

Farmer, Franklin H.↗
NASA NTRSDate not supplied

Multiobjective Multidisciplinary Optimization of Low-Boom Supersonic Transports Using Multifidelity Models

A multidisciplinary optimization (MDO) method has been developed to design a computational fluid dynamics (CFD) based low-boom configuration that can be obtained from a Pareto solution of a low-fidelity multiobjective MDO problem with mission constraints. This paper refines the developed MDO method using multifidelity models for CFD-based multiobjective MDO. The refined MDO method can generate a low-boom configuration that satisfies the mission requirements, has the lowest takeoff gross weight and the longest range for the low-boom mission, trims the low-boom cruise flight with fuel redistributions, and has a reversed equivalent area distribution closely matching a low-boom target with ground noise level below 70 PLdB. The validity of the refined MDO method is demonstrated by a design study of a low-boom supersonic transport that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.7 and range of 3500 nm, and cruises overwater at Mach 1.8 with range of 3882 nm. Moreover, the refined MDO method eliminates the difference between the assumed cruise weight for CFD-based low-boom inverse design optimization and the estimated cruise weight of the optimal inverse design solution with respect to the mission requirements.

multidisciplinary optimization↗
DOE OSTI2022

SREMI: Super-resolution electromagnetic imaging with single-channel ground-penetrating radar

High-resolution near-surface imaging has important applications in civil engineering, infrastructure inspection, military threat detection, geological characterization, and lunar and planetary exploration. Zero-offset, singlechannel ground penetrating radar (GPR) imaging is an established technique for near-surface target imaging and sensing but often suffers from low spatial resolution and imaging artifacts, especially of deep structures. In response, we formulate the GPR imaging as a dual-sparsity optimization problem, and develop a super-resolution electromagnetic imaging method based on a fast iterative shrinkage-thresholding algorithm. We develop our GPR imaging method in the framework of electromagnetic exploding-reflectors simulation theory, therefore the imaging method is computationally efficient. In this work, we demonstrate through synthetic and field data examples that our method can produce sharper, more reliable images with fewer artifacts compared with single-pass reverse-time migration GPR method, thus leading to improved near-surface interpretation and object identification.

58 GEOSCIENCES↗