Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Alternating minimization”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Energy efficiency and renewable energy utilization in desalination systems

Desalination technologies and industry have advanced significantly in the past two decades to meet the growing freshwater demands stimulated by the compounding issues of both water quality and quantity in many regions of the world. As desalination processes are energy demanding, there have been many efforts dedicated to improve energy-efficiency of the process units, enhance energy conservation and recovery, and increase renewable energy integration in desalination plants. This review article highlights recent key advances and discusses possible venues for further development in desalination energy portfolio to reduce specific energy consumption and, via integration with solar energy, to minimize the environmental footprint associated with freshwater production. First, an overview of current desalination technologies and their energy requirements are presented followed by a discussion on opportunities for improving energy efficiency and energy recovery in both membrane and thermal desalination technologies. Then, various combinations of renewable energy driven desalination plants are discussed with some recent highlights in solar energy driven membrane, thermal and hybrid desalination processes. Technological readiness levels for novel desalination processes, their perceived impact and expected near-future developments in renewable energy integrated desalination technologies are presented. Lastly, the potential for solar driven desalination as a cost-competitive freshwater supply alternative is discussed

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Distributed Optimization in Distribution Systems: Use Cases, Limitations, and Research Needs

We report electric distribution grid operations typically rely on both centralized optimization and local non-optimal control techniques. As an alternative, distribution system operational practices can consider distributed optimization techniques that leverage communications among various neighboring agents to achieve optimal operation. With the rapidly increasing integration of distributed energy resources (DERs), distributed optimization algorithms are growing in importance due to their potential advantages in scalability, flexibility, privacy, and robustness relative to centralized optimization. Implementation of distributed optimization offers multiple challenges and also opportunities. This paper provides a comprehensive review of the recent advancements in distributed optimization for electric distribution systems and classifications using key attributes. Problem formulations and distributed optimization algorithms are provided for example use cases, including volt/var control, market clearing process, loss minimization, and conservation voltage reduction. Finally, this paper also presents future research needs for the applicability of distributed optimization algorithms in the distribution system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

New mechanism of structuring associated with the quasi-merohedral twinning by an example of Ca{sub 1–x}La{sub x}F{sub 2+x} ordered solid solutions

Merohedry is considered an inseparable property of atomic structures, and uses for the refinement of structural data in a process of correct determination of structure of compounds. Transformation of faulty structures stimulated by decreasing of systemic cumulative energy leads to generation of merohedral twinning type. Ordering is accompanied by origin of antiphase domains. If ordering belongs to the CuAu type, it is accompanied by tetragonal distortions along different (100) directions. If a crystal consists of mosaic of nanodimensional antiphase domains, the conjugation of antiphase domains with different tetragonality leads to monoclinic distortions, at that, conjugated domains are distorted mirrorly. Similar system undergoes further transformation by means of quasi-merohedral twinning. As a result of quasi-merohedry, straight-lines of lattices with different monoclinic distortions are transformed into coherent lattice broken-lines providing minimization of the cumulative energy. Structuring is controlled by regularities of the self-organization. However stochasticity of ordering predetermines the origin areas where few domains with different tetragonality contact which leads to the origin of faulty fields braking regular passage of structuring. Resulting crystal has been found structurally non-uniform, furthermore structural non-uniformity permits identifying elements and stages of a process. However there is no precondition preventing arising the origin of homogenous states. Effect has been revealed in Ca{sub 1–x}La{sub x}F{sub 2+x} solid solution, but it can be expected that distortions of regular alternation of ions similar to antiphase domains can be obtained in non-equilibrium conditions in compounds and similar effect of the quasi-merohedry can falsify results of structural analysis.

36 MATERIALS SCIENCE↗
NASA NTRS1993

An approach to estimate the areal rain-rate distribution from spaceborne radar by the use of multiple thresholds

Estimates of rain rate derived from a spaceborne weather radar will be most reliable over an intermediate range of values. In forming an area average of the rain rate, an alternative to the averaging of the high-resolution estimates, irrespective of their individual accuracies, is a multiple threshold approach. The method is based on the fact that the fractional area above a particular rain-rate threshold R(j) is related to the cumulative distribution of rain rates evaluated at R(j). Varying the threshold over the effective dynamic range of the radar yields the cumulative distribution function over this range. To obtain the distribution at all rain rates, a lognormal or gamma test function is selected such that the mean-square error between the test function and the measured values is minimized. Once the unknown parameters are determined, the first-order statistics of the areawide rain-rate distribution can be found. Tests of the method with data from the SPANDAR radar provide comparisons between it and the single threshold and the direct averaging approaches.

Meneghini, R.↗
NASA NTRS2001

Experiments and Model Development for the Investigation of Sooting and Radiation Effects in Microgravity Droplet Combustion

Today, despite efforts to develop and utilize natural gas and renewable energy sources, nearly 97% of the energy used for transportation is derived from combustion of liquid fuels, principally derived from petroleum. While society continues to rely on liquid petroleum-based fuels as a major energy source in spite of their finite supply, it is of paramount importance to maximize the efficiency and minimize the environmental impact of the devices that burn these fuels. The development of improved energy conversion systems, having higher efficiencies and lower emissions, is central to meeting both local and regional air quality standards. This development requires improvements in computational design tools for applied energy conversion systems, which in turn requires more robust sub-model components for combustion chemistry, transport, energy transport (including radiation), and pollutant emissions (soot formation and burnout). The study of isolated droplet burning as a unidimensional, time dependent model diffusion flame system facilitates extensions of these mechanisms to include fuel molecular sizes and pollutants typical of conventional and alternative liquid fuels used in the transportation sector. Because of the simplified geometry, sub-model components from the most detailed to those reduced to sizes compatible for use in multi-dimensional, time dependent applied models can be developed, compared and validated against experimental diffusion flame processes, and tested against one another. Based on observations in microgravity experiments on droplet combustion, it appears that the formation and lingering presence of soot within the fuel-rich region of isolated droplets can modify the burning rate, flame structure and extinction, soot aerosol properties, and the effective thermophysical properties. These observations led to the belief that perhaps one of the most important outstanding contributions of microgravity droplet combustion is the observation that in the absence of asymmetrical forced and natural convection, a soot shell is formed between the droplet surface and the flame, exerting an influence on the droplet combustion response far greater than previously recognized. The effects of soot on droplet burning parameters, including burning rate, soot shell dynamics, flame structure, and extinction phenomena provide significant testing parameters for studying the structure and coupling of soot models with other sub-model components.

Choi, Mun Young↗
NASA NTRS2003

Software for Training in Pre-College Mathematics

The Intelligent Math Tutor (IMT) is a computer program for training students in pre-college and college-level mathematics courses, including fundamentals, intermediate algebra, college algebra, and trigonometry. The IMT can be executed on a server computer for access by students via the Internet; alternatively, it can be executed on students computers equipped with compact- disk/read-only-memory (CD-ROM) drives. The IMT provides interactive exercises, assessment, tracking, and an on-line graphing calculator with algebraic-manipulation capabilities. The IMT provides an innovative combination of content, delivery mechanism, and artificial intelligence. Careful organization and presentation of the content make it possible to provide intelligent feedback to the student based on performance on exercises and tests. The tracking and feedback mechanisms are implemented within the capabilities of a commercial off-the-shelf development software tool and are written in the Unified Modeling Language to maximize reuse and minimize development cost. The graphical calculator is a standard feature of most college and pre-college algebra and trigonometry courses. Placing this functionality in a Java applet decreases the cost, provides greater capabilities, and provides an opportunity to integrate the calculator with the lessons.

Shelton, Robert O.↗
DOE OSTI2018

Ion movement in thin Nafion films under an applied electric field

In this study, the electromechanical response of Nafion films with and without an ionic liquid (1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide) ([emim]Tf 2 N) additive were characterized under an applied electric field in situ using neutron reflectometry (NR) and voltage modulated atomic force microscopy (VM-AFM). NR showed that pure Nafion films exhibited no response under field strengths of ~0.18 V/μm, while VM-AFM measurements showed a minimal response at higher field strengths (~200 V/μm), which is ascribed to the residual water presence in the films. The addition of ionic liquid resulted in clear electroresponsiveness seen in both NR and VM-AFM. NR results indicated mass migration away from the cathodic interface driven by the [emim] + movement in the direction of the electric field. The lack of ionic liquid accumulation at the electrode interfaces contrasts the bulk electromechanical behavior of similar systems reported in the literature. VM-AFM measurements were able to resolve the relative contributions of the [emim] + cation and Tf 2 N - anion to film deformation by alternating the direction of the applied field and support the [emim] + dominant migration seen in the NR results. Lastly, the findings presented here emphasize the need for the nanoscale analysis of material properties of electroresponsive thin film systems and demonstrate the potential for probing electric field effects using in situ techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2012

THE DEPENDENCE OF QUENCHING UPON THE INNER STRUCTURE OF GALAXIES AT 0.5 {<=} z < 0.8 IN THE DEEP2/AEGIS SURVEY

The shutdown of star formation in galaxies is generally termed 'quenching'. Quenching may occur through a variety of processes, e.g., active galactic nucleus (AGN) feedback, stellar feedback, or the shock heating of gas in the dark matter halo. However, which mechanism(s) is, in fact, responsible for quenching is still in question. This paper addresses quenching by searching for traces of possible quenching processes through their effects on galaxy structural parameters such as stellar mass (M {sub *}), M {sub *}/r {sub e}, surface stellar mass density ({approx}M {sub *}/r {sup 2} {sub e}), and Sersic index (n). We analyze the rest-frame U - B color correlations versus these structural parameters using a sample of galaxies in the redshift range 0.5 {<=} z < 0.8 from the DEEP2/AEGIS survey. In addition to global radii, stellar masses, and Sersic parameters, we also use 'bulge' and 'disk' photometric measurements from GIM2D fits to HST/ACS V and I images. We assess the tightness of the color relationships by measuring their 'overlap regions', defined as the area in color-parameter space in which red and blue galaxies overlap; the parameter that minimizes these overlap regions is considered to be the most effective color discriminator. We find that Sersic index (n) has the smallest overlap region among all tested parameters and resembles a step function with a threshold value of n = 2.3. There exists, however, a significant population of outliers with blue colors yet high n values that seem to contradict this behavior; they make up Almost-Equal-To 40% of n > 2.3 galaxies. We hypothesize that their Sersic values may be distorted by bursts of star formation, AGNs, and/or poor fits, leading us to consider central surface stellar mass density, {Sigma}*{sub 1kpc}, as an alternative to Sersic index. Not only does {Sigma}*{sub 1kpc} correct the outliers, but it also forms a tight relationship with color, suggesting that the innermost structure of galaxies is most physically linked with quenching. Furthermore, at z {approx} 0.65, the majority of the blue cloud galaxies cannot simply fade onto the red sequence since their GIM2D bulge masses are only half as large on average as the bulge masses of similar red sequence galaxies, thus demonstrating that stellar mass must absolutely increase at the centers of galaxies as they quench. We discuss a two-stage model for quenching in which galaxy star formation rates are controlled by their dark halos while they are still in the blue cloud and a second quenching process sets in later, associated with the central stellar mass buildup. The mass buildup is naturally explained by any non-axisymmetric features in the potential, such as those induced by mergers and/or disk instabilities. However, the identity of the second quenching agent is still unknown. We have placed our data catalog online.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Inverse problem in the large momentum effective theory framework

One proposal to compute parton distributions from first principles is the large momentum effective theory (LaMET), which requires the Fourier transform of matrix elements computed nonperturbatively. Lattice quantum chromodynamics (QCD) provides calculations of these matrix elements over a finite range of Fourier harmonics that are often noisy or unreliable in the largest computed harmonics. It has been suggested that enforcing an exponential decay of the missing harmonics helps alleviate this issue. Using nonperturbative data, we show that the uncertainty introduced by this inverse problem in a realistic setup remains significant without very restrictive assumptions, and that the importance of the exact asymptotic behavior is minimal for values of 𝑥 where the framework is currently applicable. We show that the crux of the inverse problem lies in harmonics of the order of 𝜆 = 𝑧⁢𝑃 𝑧 ∼ 5–15, where the signal in the lattice data is often barely existent in current studies, and the asymptotic behavior is not firmly established. We stress the need for more sophisticated techniques to account for this inverse problem, whether in the LaMET or related frameworks like the short-distance factorization. We also address a misconception that, with available lattice methods, the LaMET framework allows a “direct” computation of the 𝑥-dependence, whereas the alternative short-distance factorization only gives access to moments or fits of the 𝑥-dependence.

Dutrieux, Hervé [Aix-Marseille Université, Marseil↗
DOE OSTI2023

Technical Report on Waveform Fit Metrics for Global Models

The new WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes an increased emphasis on the development of Earth models and methods to predict entire seismic and acoustic waveforms more accurately. In general, this increased emphasis is predicated on the need to better characterize seismic events and provide improved model-based discrimination between event types including earthquakes and explosions. More specifically, while current moment tensor inversion methods tend to work well for larger events (M>~4) using tuned 1-D Earth models, the development of state-of-the-art 3-D models and methods is required for the prediction of shorter period waves over large areas for discrimination of smaller events. There is no standard metric for model-based waveform prediction accuracy used in the waveform modeling/inversion community. However, there are several popular waveform misfit definitions; and minimizing the corresponding objective functions is the goal of waveform inversion. Some example misfit definitions employed for adjoint waveform tomography include measures of simple travel time differences (e.g. Tape et al., 2010), cross-correlation travel time differences (e.g. Luo and Schuster, 1991), multi-taper frequency dependent methods (e.g. Lei et al., 2020), time-frequency phase misfit functions (e.g. Fichtner 2010; Rodgers et al., 2022), normalized cross-correlation methods (e.g. Tao et al., 2018), and others. In some cases, these misfit definitions also involve complicated weighting schemes and summations over multiple frequency bands making it difficult to duplicate the misfit measurement with alternative models and datasets. Although each of the misfit definitions mentioned above are useful for developing waveform models, the actual misfit values are not usually meaningful outside of a given project, model, and/or dataset. Therefore, it is difficult to understand and communicate model performance for predicting waveforms and comparing to other models and/or new model iterations with a different dataset. Therefore, there is a need for a generalized method for evaluating overall model performance that is independent from the specific misfit chosen to develop the waveform models that is also intuitive and meaningful. In this report, we describe a new metric we refer to as ‘Percent of Correlated Signal’. The following sections describe and demonstrate the metric with a case study event and a more rigorous test using a random selection of globally distributed events. While the focus here is on global tomography models, the metric is meant to applicable to regional ‘wiggle-for-wiggle’ waveform models/studies as well.

58 GEOSCIENCES↗
NASA NTRS2011

Saturn Neutron Exosphere as Source for Inner and Innermost Radiation Belts

Energetic proton and electron measurements by the ongoing Cassini orbiter mission are expanding our knowledge of the highest energy components of the Saturn magnetosphere in the inner radiation belt region after the initial discoveries of these belts by the Pioneer 11 and Voyager 2 missions. Saturn has a neutron exosphere that extends throughout the magnetosphere from the cosmic ray albedo neutron source at the planetary main rings and atmosphere. The neutrons emitted from these sources at energies respectively above 4 and 8 eV escape the Saturn system, while those at lower energies are gravitationally bound. The neutrons undergo beta decay in average times of about 1000 seconds to provide distributed sources of protons and electrons throughout Saturn's magnetosphere with highest injection rates close to the Saturn and ring sources. The competing radiation belt source for energetic electrons is rapid inward diffusion and acceleration of electrons from the middle magnetosphere and beyond. Minimal losses during diffusive transport across the moon orbits, e.g. of Mimas and Enceladus, and local time asymmetries in electron intensity, suggest that drift resonance effects preferentially boost the diffusion rates of electrons from both sources. Energy dependences of longitudinal gradient-curvature drift speeds relative to the icy moons are likely responsible for hemispheric differences (e.g., Mimas, Tethys) in composition and thermal properties as at least partly produced by radiolytic processes. A continuing mystery is the similar radial profiles of lower energy (<10 MeV) protons in the inner belt region. Either the source of these lower energy protons is also neutron decay, but perhaps alternatively from atmospheric albedo, or else all protons from diverse distributed sources are similarly affected by losses at the moon' orbits, e.g. because the proton diffusion rates are extremely low. Enceladus cryovolcanism, and radiolytic processing elsewhere on the icy moon and ring surfaces, are additional sources of protons via ionization and charge exchange from breakup of water molecules. But one must then account somehow for local acceleration to the observed keV-MeV energies, since moon sweeping and E-ring absorption would remove protons diffusing inward from the middle magnetosphere. Although the main rings block further inward diffusion from the inner radiation belts, the exospheric neutron-decay source, combined with much slower diffusion of protons relative to electrons, may produce an innermost radiation belt in the gap between the upper atmosphere and the D-ring. This innermost belt will first be explored in-situ during the final proximal orbits of the Cassini mission.

Cooper, John↗
NASA NTRS2019

Silver Foam as Long-Term Passive Biocide for Potable Water Systems

A spacecraft water disinfection system, suitable for extended length space exploration missions, should prevent or control the growth of microbes, prevent or limit biofilm formation, and prevent microbiologically-influenced corrosion. In addition, the system should have minimal maintenance requirements, the effluent should be chemically compatible with all materials in contact with the water, be safe for human consumption, and suitable to be shared across international spacecraft platforms and mission architectures. Silver ions are a proven broad spectrum biocide. Silver is also the potable water biocide of choice for future exploration missions. Currently, the proposed method for implementing silver biocide in spacecraft systems relies on silver electrode technology to produce a controlled amount of silver ions. Unfortunately, electrolytic-based silver dosing presents multiple inherent challenges that affect performance and increase maintenance requirements over time. To decrease the risk of failure, an alternative silver biocide delivery method is needed. Control-release technology is an attractive option for developing a passive high-reliability silver dosing device. The concept of a nanoparticle/polyurethane (PU) composite foam for the controlled release of silver was prompted by the controlled release technology developed by NASA for the delivery of corrosion inhibitors and indicators. This paper presents the technical background and results from the synthesis and properties testing of the silver nanoparticles (AgNPs)/PU composite foam that is being developed for use in spacecraft potable water systems.

Irwin, Tesia D.↗
DOE OSTI2022

Opportunities for Hybrid Wind and Solar PV Plants in India

The share of variable renewable energy (VRE) on India's grid has surpassed 100 GW, and the government has ambitious plans reach 450 GW by 2030. One strategy to increase wind and solar PV deployment is through the co-location of wind and solar PV plants to form a single hybrid power plant. Hybrid plants have the potential to reduce transmission infrastructure costs and variability in the output power profile compared to a standalone plant with a single technology, and this resource analysis aims to take a first step towards quantifying the potential savings from hybridizing wind and solar PV plants in India and the size of this opportunity. We utilize a brute-force optimization to minimize the levelized cost of energy (LCOE) for standalone wind, standalone solar PV, and hybrid wind/solar PV plants across all of India. By comparing these LCOEs, we determine that locations where hybrid plants exhibit potential cost savings and grid benefits exhibit both; a high interconnection cost and; a wind capacity factor between roughly 34% and 38%. However, because our work does not capture the value of the electricity generated by looking at energy prices, nor does it quantify the potential of hybrids to provide other value streams such as firm capacity and reserves. Further, because the work does not compare solar PV and wind hybrids to alternative generation technologies or storage systems, it cannot be considered a holistic cost-benefit analysis.

14 SOLAR ENERGY↗
DOE OSTI2016

A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems

In this paper a new, simple and universal formulation of the HLLEM Riemann solver (RS) is proposed that works for general conservative and non-conservative systems of hyperbolic equations. For non-conservative PDE, a path-conservative formulation of the HLLEM RS is presented for the first time in this paper. The HLLEM Riemann solver is built on top of a novel and very robust path-conservative HLL method. It thus naturally inherits the positivity properties and the entropy enforcement of the underlying HLL scheme. However, with just the slight additional cost of evaluating eigenvectors and eigenvalues of intermediate characteristic fields, we can represent linearly degenerate intermediate waves with a minimum of smearing. For conservative systems, our paper provides the easiest and most seamless path for taking a pre-existing HLL RS and quickly and effortlessly converting it to a RS that provides improved results, comparable with those of an HLLC, HLLD, Osher or Roe-type RS. This is done with minimal additional computational complexity, making our variant of the HLLEM RS also a very fast RS that can accurately represent linearly degenerate discontinuities. Our present HLLEM RS also transparently extends these advantages to non-conservative systems. For shallow water-type systems, the resulting method is proven to be well-balanced. Several test problems are presented for shallow water-type equations and two-phase flow models, as well as for gas dynamics with real equation of state, magnetohydrodynamics (MHD & RMHD), and nonlinear elasticity. Since our new formulation accommodates multiple intermediate waves and has a broader applicability than the original HLLEM method, it could alternatively be called the HLLI Riemann solver, where the “I” stands for the intermediate characteristic fields that can be accounted for. -- Highlights: •New simple and general path-conservative formulation of the HLLEM Riemann solver. •Application to general conservative and non-conservative hyperbolic systems. •Inclusion of sub-structure and resolution of intermediate characteristic fields. •Well-balanced for single- and two-layer shallow water equations and multi-phase flows. •Euler equations with real equation of state, MHD equations, nonlinear elasticity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2013

Mission Data System Java Edition Version 7

The Mission Data System framework defines closed-loop control system abstractions from State Analysis including interfaces for state variables, goals, estimators, and controllers that can be adapted to implement a goal-oriented control system. The framework further provides an execution environment that includes a goal scheduler, execution engine, and fault monitor that support the expression of goal network activity plans. Using these frameworks, adapters can build a goal-oriented control system where activity coordination is verified before execution begins (plan time), and continually during execution. Plan failures including violations of safety constraints expressed in the plan can be handled through automatic re-planning. This version optimizes a number of key interfaces and features to minimize dependencies, performance overhead, and improve reliability. Fault diagnosis and real-time projection capabilities are incorporated. This version enhances earlier versions primarily through optimizations and quality improvements that raise the technology readiness level. Goals explicitly constrain system states over explicit time intervals to eliminate ambiguity about intent, as compared to command-oriented control that only implies persistent intent until another command is sent. A goal network scheduling and verification process ensures that all goals in the plan are achievable before starting execution. Goal failures at runtime can be detected (including predicted failures) and handled by adapted response logic. Responses can include plan repairs (try an alternate tactic to achieve the same goal), goal shedding, ignoring the fault, cancelling the plan, or safing the system.

Reinholtz, William K.↗
DOE OSTI2026

Constrained variational optimization of counting-time allocation in sequential scattering measurements: Application to Bonse–Hart USANS

Sequential scattering measurements are often performed under a fixed experimental-time budget, even though the expected count rate varies strongly across the measured coordinate. When the dwell time at each measurement position can be controlled independently, this variation creates a general resource-allocation problem: how should the available time be distributed to minimize the uncertainty of the reconstructed profile? We formulate this problem as a constrained variational optimization for measurements governed by Poisson counting statistics. When each measurement is treated independently, minimizing the averaged squared relative uncertainty yields an inverse-square-root intensity allocation. The formulation is then generalized to include correlations between neighboring measurements and an instrumental resolution operator, leading to an allocation criterion that equalizes the marginal reduction in posterior uncertainty per unit measurement time. Bonse–Hart ultra-small-angle neutron scattering (USANS), in which reciprocal space is sampled sequentially through analyzer-angle stepping, provides an experimentally grounded application. Computational benchmarking shows that the optimized allocation outperforms uniform-time and constant-relative-error strategies, while application to an experimentally measured graphite USANS profile from the Spallation Neutron Source, using Poisson resampling under alternative schedules, demonstrates how counting time should be redistributed toward weak-intensity regions under an identical total duration. The resulting framework applies to sequential scattering and related scanning measurements whenever local dwell times are adjustable and directly determine the measurement uncertainties, and when the relevant correlation and instrumental-response models are available.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗
DOE OSTI2020

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗
DOE OSTI2011

Generalized variational principle for excited states using nodes of trial functions

The familiar variational principle provides an upper bound to the ground-state energy of a given Hamiltonian. This allows one to optimize a trial wave function by minimizing the expectation value of the energy. This approach is also trivially generalized to excited states, so that given a trial wave function of a certain symmetry, one can compute an upper bound to the lowest-energy level of that symmetry. In order to generalize further and build an upper bound of an arbitrary excited state of the desired symmetry, a linear combination of basis functions is generally used to generate an orthogonal set of trial functions, all bounding their respective states. However, sometimes a compact wave-function form is sought, and a basis-set expansion is not desirable or possible. Here we present an alternative generalization of the variational principle to excited states that does not require explicit orthogonalization to lower-energy states. It is valid for one-dimensional systems and, with additional information, to at least some n-dimensional systems. This generalized variational principle exploits information about the nodal structure of the trial wave function, giving an upper bound to the exact energy without the need to build a linear combination of basis functions. To illustrate the theorem we apply it to a nontrivial example: the 1s2s {sup 1} S excited state of the helium atom.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗