Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational approach”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Ultrasound Thermal Field Imaging of Opaque Fluids

We have initiated an experimental program to develop an ultrasound system for non-intrusively imaging the thermal field in opaque fluids under an externally imposed temperature gradient. Many industrial processes involve opaque fluids, such as molten metals, semiconductors, and polymers, often in situations in which thermal gradients are important. For example, one may wish to understand semiconductor crystal growth dynamics in a Bridgman apparatus. Destructive testing of the crystal after the process is completed gives only indirect information about the fluid dynamics of the formation process. Knowledge of the coupled thermal and velocity fields during the growth process is then essential. Most techniques for non-intrusive velocity and temperature measurement in fluids are optical in nature, and hence the fluids studied must be transparent. In some cases (for example, LDV (laser Doppler velocimetry) and PIV (particle imaging velocimetry)) the velocities of small neutrally buoyant seed particles suspended in the fluid, are measured. Without particle seeding one can use the variation of the index of refraction of the fluid with temperature to visualize, through interferometric, Schlieren or shadowgraph techniques, the thermal field. The thermal field in turn gives a picture of the pattern existing in the fluid. If the object of study is opaque, non-optical techniques must be used. In this project we focus on the use of ultrasound, which propagates easily through opaque liquids and solids. To date ultrasound measurements have almost exclusively relied on the detection of sound scattered from density discontinuities inside the opaque material of interest. In most cases it has been used to visualize structural properties, but more recently the ultrasound Doppler velocimeter has become available. As in the optical case, it relies on seed particles that scatter Doppler shifted sound back to the detector. Doppler ultrasound techniques are, however, not useful for studying convective fluid flow in crystal growth, because particle seeding is unacceptable and flow velocities are typically too low to be resolved, and may be even lower in microgravity conditions where buoyancy forces are negligible. We will investigate a different use of ultrasound to probe the flows of opaque fluids non-intrusively and without the use of seed particles: our goal is to ultrasonically visualize the thermal field of opaque fluids with relatively high spatial resolution. The proposed technique relies upon the variation of sound speed with temperature of the fluid. A high frequency ultra-sound pulse passing through a fluid-filled chamber will traverse the chamber in a time determined by the relevant chamber dimension and the temperature of the fluid through which the pulse passes. With high time-resolution instrumentation that compares the excitation signal with the received pulse we can detect the influence of the fluid temperature on the pulse travel time. This is effectively an interferometric system, which in its optical form is an extremely sensitive approach to measuring thermal fields in fluids. Moreover, the temperature dependence of sound velocity in liquid metals is comparable to the temperature dependence of the speed of light required for accurate interferometric thermal images in transparent fluids. With an array of transducers scanned electronically a map of the thermal field over the chamber could be produced. An alternative approach would be to use the ultrasound analog of the shadowgraph technique. In the optical version, collimated light passes through the fluid, where it is focused or defocused locally by temperature field induced variations of the index of refraction. The resulting image reveals the thermal field through the spatial pattern of light intensity variations. By analogy, an ultrasound plane wave traversing an opaque fluid sample would be also locally focused or defocused depending on the speed of sound variations, giving rise to spatial variations in sound intensity that will reveal the thermal field pattern. These approaches could be applied to any situation in which temperature differences are expected to occur, and will rapidly provide information about the flow that simply cannot be obtained by any current intrusive or non-intrusive diagnostic technique. As materials processing in microgravity matures it will become increasingly important to have available simple and versatile diagnostic tools, such as we will develop, for studying the flows of opaque fluids under thermal forcing.

Andereck, C. David↗
DOE OSTI2020

Characterizing Peatland Microtopography Using Gradient and Microform-Based Approaches

Abstract Peatlands represent an important component of the global carbon cycle, storing 180–621 Gt of carbon (C). Small-scale spatial variations in elevation, frequently referred to as microtopography, influence ecological processes associated with the peatland C cycle, including Sphagnum photosynthesis and methane flux. Microtopography can be characterized with measures of topographic variability and by using conceptual classes (microforms) linked to function: most commonly hummocks and hollows. However, the criteria used to define these conceptual classes are often poorly described, if at all, and vary between studies. Such inconsistencies compel development of explicit quantitative methods to classify microforms. Furthermore, gradient-based characterizations that describe spatial variability without the use of microforms are lacking in the literature. Therefore, the objectives of this study were to (1) calculate peatland microtopographical elevation gradients and measures of spatial variability, (2) develop three microform classification methods intended for specific purposes, and (3) evaluate and contrast classification methods. Our results suggest that at spatial scales much larger than microforms, elevation distributions are unimodal and are well approximated with parametric probability density functions. Results from classifications were variable between methods and years and exhibited significant differences in mean hollow areal coverages of a raised ombrotrophic bog. Our results suggest that the conceptualization and classification of microforms can significantly influence microtopographic structural metrics. The three explicit methods for microform classification described here may be used and built upon for future applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Atomic Ordering-Induced Ensemble Variation in Alloys Governs Electrocatalyst On/Off States

The catalytic behavior of a material is influenced by ensembles—the geometric configuration of atoms. Traditional approaches, mainly utilizing solid-solution alloys in electrocatalysis, have often overlooked the challenges posed by concurrent changes in the electronic structure (i.e. d-band center) when the composition is altered. Here, this study introduces a methodology that distinctly separates the geometric effects (i.e. ensembles) from the electronic structure. We compare the reactivity of compositionally identical, but structurally different Pd 3 Bi ordered intermetallic and solid-solution alloys. Remarkably, we find that Pd 3 Bi intermetallics display nearly no reactivity for the methanol oxidation (MOR), while their solid-solution counterparts have significant reactivity. This highlights a unique case where materials with identical chemical compositions demonstrate drastically different catalytic behavior underscoring the critical importance of ensembles in electrocatalysis. Specifically, Pd 3 Bi intermetallics form smaller ensembles (average coordination number: 4.5 ± 1.6) with almost no measurable MOR activity at room temperature, in contrast to the solid-solution Pd 3 Bi that exhibit larger ensembles (average coordination number: 6.8 ± 0.9) and considerable MOR reactivity (0.5 mA cm −2 Pd ). An ordered Pd 3 Bi alloy, with an intermediate ensemble size (average coordination number: 5.3 ± 1.2), displays moderate MOR activity (0.1 mA cm −2 Pd ), further confirming the direct correlation between ensemble size and catalytic activity. Notably, all Pd 3 Bi alloys maintain similar electronic structures, because the chemical composition of the alloys is fixed, indicating that the differences in reactivity are predominantly from changes to the ensemble size. Our findings offer an approach for precisely controlling catalytic activity through manipulating the geometric configuration of the atoms within an alloy, paving the way for more efficient catalyst design.

alloys↗
NASA NTRS1994

Northern middle-latitude ozone profile features and trends observed by SBUV and Umkehr, 1979-1990

A comparison of Umkehr ozone profile data with the reprocessed solar backscatter ultraviolet (SBUV) ozone data in the northern middle-latitude region, 30 deg to 50 deg N, is reported. Although significant biases exist between the two types of observations, the long-term variations and least squares linear regression trends agree remarkably well over the comparison period of 1979 to 1990. The ozone trend in the upper stratosphere is of the order of -0.9%/yr. Near 25 km, little if any trend appears, but a larger negative trend is seen in the lower stratosphere near 15 km. Comparisons show that the average annual ozone cycles in the profiles also agree well. The upper stratospheric ozone results are consistent with photochemical model predictions of ozone depletion near 40 km that are due to the release of anthropogenically produced chlorofluorocarbons. The lower stratospheric ozone trend results are in reasonable agreement with published ozonesonde data trends. It is shown that the ozone trends in the lower stratospheric layers impact significantly on the total ozone trend of the order of -0.47%/yr. The good agreement now seen between the two types of observations suggest that the combined ground-based and satellite approach could provide a valuable database for long-term monitoring of stratospheric ozone for trends and extraordinary variations.

Deluisi, J. J.↗
NASA NTRS2004

Laser Sounder for Global Measurement of CO2 Concentrations in the Troposphere from Space: Progress

We describe progress toward developing a laser-based technique for the remote measurement of the tropospheric CO2 concentrations from orbit. Our goal is to demonstrate a lidar technique and instrument technology that will permit measurements of the CO2 column abundance in the lower troposphere from aircraft at the few ppm level, with a capability of scaling to permit global CO2 measurements from orbit. Accurate measurements of the tropospheric CO2 mixing ratio from space are challenging due to the many potential error sources. These include possible interference from other trace gas species, the effects of temperature, clouds, aerosols & turbulence in the path, changes in surface reflectivity, and variability in dry air density caused by changes in atmospheric pressure, water vapor and topographic height. Some potential instrumental errors include frequency drifts in the transmitter, small transmission and sensitivity drifts in the instrument. High signal-to-noise ratios and measurement stability are needed for mixing ratio estimates at the few ppm level. We have been developing a laser sounder approach as a candidate for a future space mission. It utilizes multiple different laser transmitters to permit simultaneous measurement of CO2 and O2 extinction, and aerosol backscatter in the same measurement path. It directs the narrow co-aligned laser beams from the instrument's fiber lasers toward nadir, and measures the energy of the strong laser echoes reflected from the Earth's land and water surfaces. During the measurement its narrow linewidth lasers are rapidly tuned on- and off- selected CO2 line near 1572 nm and an O2 absorption line near 770 nm. The receiver measures the energies of the laser echoes from the surface and any clouds and aerosols in the path with photon counting detectors. Ratioing the on- to off-line echo pulse energies for each gas permits the column extinction and column densities of CO2 and O2 to be estimated simultaneously via the differential absorption lidar technique. For the on-line wavelengths, the side of the selected absorption lines are used, which due to pressure broadening, weights the measurements to the lower troposphere, where CO2 variations caused by surface sources and sinks are largest. Simultaneous measurements of O2 column abundance are made using an identical approach using an O2 line. The laser backscatter profiles from clouds and aerosols are measured with other lidar channels, which permits identifying measurements influenced by clouds and/or aerosol scattering in the path. For space use, our lidar would continuously measure at nadir in near polar circular orbit. Using dawn and dusk measurements made over the same region will make it possible to sample the diurnal variations in CO2 mixing ratios. A 1-m diameter telescope is used for the receiver for all wavelengths. When averaging over 50 seconds, our calculations show a SNR of approximately 1500 is achievable for each gas at each on- and off-line measurement. Measurements from such a mission can be used to generate monthly global maps of the lower tropospheric CO2 column abundance. Our calculations show global coverage with an accuracy of a few ppm with a spatial resolution of approximately 50,000 sq. km are achievable each month. We have demonstrated some key elements of the laser, detector and receiver approaches in the laboratory and with measurements over a 206 m horizontal path. These include stable measurements of CO2 line shapes in an absorption cell using a fiber laser amplifier seeded by a tunable diode laser, measurement of small amplitude changes at low optical signal levels with the PMT receiver, and comparison of the horizontal path measurements of CO2 against those from an in-situ instrument.

Abshire, J. B.↗
DOE OSTI2025

Genomic Analysis of the Natural Variation of Fatty Acid Composition in Seed Oils of Camelina sativa

Camelina sativa is an oilseed crop that has shown strong promise as a biofuel feedstock. The profile of fatty acids greatly influences the oil quality; however, genetic mechanisms that determine the natural variation of fatty acid composition in camelina are not fully understood. A genome wide association study (GWAS) was performed to uncover genetic loci that may contribute to the contents of major fatty acids such as oleic and linolenic acids in camelina seed. Two approaches were taken to improve the GWAS efficiency. First, growing a diversity panel of 212 accessions in four locations and two nitrogen fertilization conditions revealed great variation in fatty acid contents in seeds. Second, using an improved reference genome, abundant markers, including 203,320 single nucleotide polymorphisms (SNPs) and 99,067 insertions/deletions (indels), were developed, which refined the population structure of the diversity panel. GWAS resulted in 118 genetic markers across 31 trait/treatment conditions. Closely linked markers were determined based on linkage decay and by comparing secondarily associated markers when highly associated ones were removed. Candidate genes were examined by comparing the pangenomes of 12 high-quality reference genomes. This study provides new resources to understand seed lipid metabolism and improve camelina oils through molecular breeding.

Life Sciences & Biomedicine - Other Topics↗
DOE OSTI2024

Weak boundary conditions for Lagrangian shock hydrodynamics: A high-order finite element implementation on curved boundaries

Here, we propose a new Nitsche-type approach for weak enforcement of normal velocity boundary conditions for a Lagrangian discretization of the compressible shock-hydrodynamics equations using high-order finite elements on curved boundaries. Specifically, the variational formulation is appropriately modified to enforce free-slip wall boundary conditions, without perturbing the structure of the function spaces used to represent the solution, with a considerable simplification with respect to traditional approaches. Total energy is conserved and the resulting mass matrices are constant in time. The robustness and accuracy of the proposed method are validated with an extensive set of tests involving nontrivial curved boundaries.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Architecture-Based Unit Testing of the Flight Software Product Line

This paper presents an analysis of the unit testing approach developed and used by the Core Flight Software (CFS) product line team at the NASA GSFC. The goal of the analysis is to understand, review, and reconunend strategies for improving the existing unit testing infrastructure as well as to capture lessons learned and best practices that can be used by other product line teams for their unit testing. The CFS unit testing framework is designed and implemented as a set of variation points, and thus testing support is built into the product line architecture. The analysis found that the CFS unit testing approach has many practical and good solutions that are worth considering when deciding how to design the testing architecture for a product line, which are documented in this paper along with some suggested innprovennents.

Ganesan, Dharmalingam↗
NASA NTRS1975

Time-dependent hydrodynamic phenomena in the envelopes of the Be stars

Temporal and angular variations in the motion and distribution of circumstellar material about Be stars are studied by means of a simple hydrodynamical approach. The flow variables are each written as the sum of a known, steady-state axisymmetric term and a smaller component which contains both temporal and angular dependence. Substitution of these dichotomous expressions for the flow variables into Euler's equation and the continuity equation leads to a set of linearized equations for the terms which depend on time and angle. Solutions can be found to a somewhat approximated version of the linearized equations, and these solutions explain the qualitative features of the observed short-time variations.

Morgan, T. H.↗
NASA NTRS2004

Human-Induced Climate Variations Linked to Urbanization: From Observations to Modeling

The goal of this session is to bring together scientists from interdisciplinary backgrounds to discuss the data, scientific approaches and recent results focusing on the impact of urbanization on the climate. The discussion will highlight current observational and modeling capabilities being employed for investigating the urban environment and its linkage to the change in the Earth's climate system. The goal of the session is to identify our current stand and the future direction on the topic. Urbanization is one of the extreme cases of land use change. Most of population of the world has moved to urban areas. By 1995, more than 70% of population of North America and Europe were living in cities. By 2025, the United Nations estimates that 60% of the worlds population will live in cities. Although currently only 1.2% of the land is urban, better understanding of how the atmosphere-ocean-land-biosphere components interact as a coupled system and the influence of human activities on this system is critical. Our understanding of urbanization effect is incomplete, partly because human activities induce new changes on climate in addition to the original natural variations, and partly because previously few data available for study urban effect globally. Urban construction changes surface roughness, albedo, heat capacity and vegetation coverage. Traffic and industry increase atmospheric aerosol. It is suggested that urbanization may modify rainfall processes through aerosol-cloud interactions or dynamic feedbacks. Because urbanization effect on climate is determined by many factors including land cover, the city's microscale features, population density, and human lifestyle patterns, it is necessary to study urban areas over globe.

Shepherd, J. Marshall↗
DOE OSTI2020

Kensor: Coordinated Intelligence from Co-Located Sensors

Internet of Things (IoT) is becoming more pervasive in many installations, including homes, manufacturing plants, and industrial facilities of all kinds. The data that IoT produces is a reflection of usual behavior such as daily routines and scheduled tasks, but also from unexpected behavior due to unintentional or undesirable abnormalities. Here, we focus on achieving coordinated intelligence about normal and abnormal phenomena from multiple sensors that are geographically colocated in close proximity, monitoring and controlling a set of co-located devices. Given a set of co-located sensors, we seek an intelligent approach that would automatically determine the “normal” patterns of behaviors among the correlated sensors. After normal behavior is extracted, later monitoring should detect any deviant variations over time. An example application is an entry monitoring and alert system for facilities such as nuclear reactors, where badge readers, door locks, lights, weight trackers and other co-located sensors at the entry point are collectively tracked. To address this problem, we identify the possible solution approach that can be used to solve its different variants. The implemented model is developed as a combination of rules and Markov Chain methods.

Kotevska, Olivera↗
DOE OSTI2024

Detecting Anomalies in Time Series Using Kernel Density Approaches

This paper introduces a novel anomaly detection approach tailored for time series data with exclusive reliance on normal events during training. Our key innovation lies in the application of kernel-density estimation (KDE) to scrutinize reconstruction errors, providing an empirically derived probability distribution for normal events post-reconstruction. This non-parametric density estimation technique offers a nuanced understanding of anomaly detection, differentiating it from prevalent threshold-based mechanisms in existing methodologies. In post-training, events are encoded, decoded, and evaluated against the estimated density, providing a comprehensive notion of normality. In addition, we propose a data augmentation strategy involving variational autoencoder-generated events and a smoothing step for enhanced model robustness. The significance of our autoencoder-based approach is evident in its capacity to learn normal representation without prior anomaly knowledge. Through the KDE step on reconstruction errors, our method addresses the versatility of anomalies, departing from assumptions tied to larger reconstruction errors for anomalous events. Our proposed likelihood measure then distinguishes normal from anomalous events, providing a concise yet comprehensive anomaly detection solution. The extensive experimental results support the feasibility of our proposed method, yielding significantly improved classification performance by nearly 10% on the UCR benchmark data.

Frehner, Robin↗
DOE OSTI2020

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors

Abstract Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum , ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Reviews and syntheses: Opportunities for robust use of peak intensities from high-resolution mass spectrometry in organic matter studies

Abstract. Earth's biogeochemical cycles are intimately tied to the biotic and abiotic processing of organic matter (OM). Spatial and temporal variations in OM chemistry are often studied using direct infusion, high-resolution Fourier transform mass spectrometry (FTMS). An increasingly common approach is to use ecological metrics (e.g., within-sample diversity) to summarize high-dimensional FTMS data, notably Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS). However, problems can arise when FTMS peak-intensity data are used in a way that is analogous to abundances in ecological analyses (e.g., species abundance distributions). Using peak-intensity data in this way requires the assumption that intensities act as direct proxies for concentrations. Here, we show that comparisons of the same peak across samples (within-peak) may carry information regarding variations in relative concentration, but comparing different peaks (between-peak) within or between samples does not. We further developed a simulation model to study the quantitative implications of using peak intensities to compute ecological metrics (e.g., intensity-weighted mean properties and diversity) that rely on information about both within-peak and between-peak shifts in relative abundance. We found that, despite analytical limitations in linking concentration to intensity, ecological metrics often perform well in terms of providing robust qualitative inferences and sometimes quantitatively accurate estimates of diversity and mean molecular characteristics. We conclude with recommendations for the robust use of peak intensities for natural organic matter studies. A primary recommendation is the use and extension of the simulation model to provide objective guidance on the degree to which conceptual and quantitative inferences can be made for a given analysis of a given dataset. Broad use of this approach can help ensure rigorous scientific outcomes from the use of FTMS peak intensities in environmental applications.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2007

Laser Sounder for Global Measurement of CO2 Concentrations in the Troposphere from Space

Measurements of tropospheric CO2 abundance with global-coverage, a few hundred km spatial and monthly temporal resolution are needed to quantify processes that regulate CO2 storage by the land and oceans. The Orbiting Carbon Observatory (OCO) is the first space mission focused on atmospheric CO2 for measuring total column CO, and O2 by detecting the spectral absorption in reflected sunlight. The OCO mission is an essential step, and will yield important new information about atmospheric CO2 distributions. However there are unavoidable limitations imposed by its measurement approach. These include best accuracy only during daytime at moderate to high sun angles, interference by cloud and aerosol scattering, and limited signal from CO2 variability in the lower tropospheric CO2 column. We have been developing a new laser-based technique for the remote measurement of the tropospheric CO2 concentrations from orbit. Our initial goal is to demonstrate a lidar technique and instrument technology that will permit measurements of the CO2 column abundance in the lower troposphere from aircraft. Our final goal is to develop a space instrument and mission approach for active measurements of the CO2 mixing ratio at the 1-2 ppmv level. Our technique is much less sensitive to cloud and atmospheric scattering conditions and would allow continuous measurements of CO2 mixing ratio in the lower troposphere from orbit over land and ocean surfaces during day and night. Our approach is to use the 1570nm CO2 band and a 3-channel laser absorption spectrometer (i.e. lidar used an altimeter mode), which continuously measures at nadir from a near polar circular orbit. The approach directs the narrow co-aligned laser beams from the instrument's lasers toward nadir, and measures the energy of the laser echoes reflected from land and water surfaces. It uses several tunable fiber laser transmitters which allowing measurement of the extinction from a single selected CO2 absorption line in the 1570 nm band. This band is free from interference from other gases and has temperature insensitive absorption lines. During the measurement the lasers are tuned on- and off- a selected CO2 line near 1572 nm and a selected O2 line near 768 nm in the Oxygen A band at kHz rates. The lasers use tunable diode seed lasers followed by fiber amplifiers, and have spectral widths much narrower than the gas absorption lines. The receiver uses a 1-m diameter telescope and photon counting detectors and measures the background light and energies of the laser echoes from the surface. The extinction and column densities for the CO2 and O2 gases are estimated from the ratio of the on and offline surface echo via the differential optical absorption technique. Our technique rapidly alternates between several on-line wavelengths set to the sides of the selected gas absorption lines. It exploits the atmospheric pressure broadening of the lines to weight the measurement sensitivity to the atmospheric column below 5 km. This maximizes sensitivity to CO2 in the boundary layer, where variations caused by surface sources and sinks are largest. Simultaneous measurements of O2 column will use an identical approach with an O2 line. Thee laser frequencies are tunable and have narrow (MHz) line widths. In combination with sensitive photon counting detectors these enables much higher spectral resolution and precision than is possible with passive spectrometer. 1aser backscatter profiles are also measured, which permits identifying measurements made to cloud tops and through aerosol layers. The measurement approach using lasers in common-nadir-zenith path allows retrieving CO2 column mixing ratios in the lower troposphere irrespective of sun angle. Pulsed laser signals, time gated receiver and a narrow receiver field-of-view are used to isolate the surface laser echo signals and to exclude photons scattered from clouds and aerosols. Nonetheless, the optical absorption change due to a change of a few ppO2 is small, <1 % which makes achieving the needed measurement sensitivities and stabilities quite challenging. Measurement SNRs and stabilities of >600:1 are needed to estimate CO2 mixing ratio at the 1-2 ppm level. We have calculated characteristics of the technique and have demonstrated aspects of the laser, detector and receiver approaches in th e laboratory We have also measured O2 in an absorption cell, and made C02 measurements over a 400 m long (one way) horizontal path using a sensor breadboard. We will describe these and more details of our approach in the paper.

Abshire, James B.↗
NASA NTRS2005

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.↗
DOE OSTI2019

Toward the evaluation of intersystem crossing rates with variational relativistic methods

The change in electronic state from one spin multiplicity to another, known as intersystem crossing, occurs in molecules via the relativistic phenomenon of spin-orbit coupling. Current means of estimating intersystem crossing rates rely on the perturbative evaluation of spin-orbit coupling effects. This perturbative approach, valid in lighter atoms where spin-orbit coupling is weaker, is expected to break down for heavier elements where relativistic effects become dominant. Methods which incorporate spin-orbit effects variationally, such as the exact-two-component (X2C) method, will be necessary to treat this strong-coupling regime. In this paper, we present a novel procedure which produces a diabatic basis of spin-pure electronic states coupled by spin-orbit terms, generated from fully variational relativistic calculations. This method is implemented within X2C using time-dependent density-functional theory and is compared to results from a perturbative relativistic study in the weak spin-orbit coupling regime. Additional calculations on a more strongly spin-orbit-coupled [UO 2 Cl 4 ] 2- complex further illustrate the strengths of this method. This procedure will be valuable in the estimation of intersystem crossing rates within strongly spin-coupled species.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2004

Mars Express: The Exploration of Phobos

The origin, structure and composition of Phobos are essentially still unresolved issue, with major outcomes for understanding the origin and evolution of the Solar system, and of the Mars system. Several missions have been designed with this unique body the only target. Is Phobos a captured small body (and in this case, can we identify the degree of differentiation such a small object has undergone ?), or is there any co-generic process with Mars involved ? Are there still volatile species trapped, and if so where, and in what form ? What are the reciprocal contributions of Mars and Phobos material accreted in the other body? The number of key clues for deciphering the early solar system processes, the Mars-Phobos binary evolution, and the role Phobos could play in the future of Mars exploration (including human expeditions), are numerous, and exciting for a large community. Part of the answers are in the high resolution coverage of this object (optical and IR, for composition variation), and in its full spectral analysis, from UV to thermal IR : Mars Express is the unique mission to provide this opportunity, and orbit 756, August 22, a totally unique opportunity to approach Phobos so closely. OMEGA should provide all across Stickney, the deepest accessible material, and over more than 20 km along, a spatial resolution less than 200 m, and assess the composition (and possibly their variations) wrt silicates (mafic materials), oxides, and eventually hydrated minerals and organics. By no means we should miss this opportunity: no one would understand the rationale not to turn the instruments ON while so close, and pointing them towards the Phobos surface.

Phobos↗