Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Green's function methods”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A Test Problem for Codes Solving the Discretized Diffusion Equation in Cartesian Geometry Derived Via Discrete Green’s Functions

We obtain here a solution to a zone-centered discretization of the one dimensional time-dependent diffusion equation with arbitrary initial conditions and source, constant absorption and scattering opacity, and constant zone size and time step. The solution is obtained using the discrete Green’s functions of the discretized equation. The solution of the discretized equation is useful in the testing of computer codes, because the code can be expected to agree with the solution to the discrete equation, to within small errors caused by roundoff. This is in contrast to solutions of the differential equation, with which the code results only approximately agree. Finally, the usefulness of the solution for tests of an inertial confinement fusion code is demonstrated.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2007

Simple and Efficient Numerical Evaluation of Near-Hypersingular Integrals

Recently, significant progress has been made in the handling of singular and nearly-singular potential integrals that commonly arise in the Boundary Element Method (BEM). To facilitate object-oriented programming and handling of higher order basis functions, cancellation techniques are favored over techniques involving singularity subtraction. However, gradients of the Newton-type potentials, which produce hypersingular kernels, are also frequently required in BEM formulations. As is the case with the potentials, treatment of the near-hypersingular integrals has proven more challenging than treating the limiting case in which the observation point approaches the surface. Historically, numerical evaluation of these near-hypersingularities has often involved a two-step procedure: a singularity subtraction to reduce the order of the singularity, followed by a boundary contour integral evaluation of the extracted part. Since this evaluation necessarily links basis function, Green s function, and the integration domain (element shape), the approach ill fits object-oriented programming concepts. Thus, there is a need for cancellation-type techniques for efficient numerical evaluation of the gradient of the potential. Progress in the development of efficient cancellation-type procedures for the gradient potentials was recently presented. To the extent possible, a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. However, since the gradient kernel involves singularities of different orders, we also require that the transformation leaves remaining terms that are analytic. The terms "normal" and "tangential" are used herein with reference to the source element. Also, since computational formulations often involve the numerical evaluation of both potentials and their gradients, it is highly desirable that a single integration procedure efficiently handles both.

Fink, Patrick W.↗
DOE OSTI2019

Cluster, classify, regress: A general method for learning discontinuous functions

This paper presents a method for solving the supervised learning problem in which the output is highly nonlinear and discontinuous. Here, it is proposed to solve this problem in three stages: (i) cluster the pairs of input-output data points, resulting in a label for each point; (ii) classify the data, where the corresponding label is the output; and finally (iii) perform one separate regression for each class, where the training data corresponds to the subset of the original input-output pairs which have that label according to the classifier. It has not yet been proposed to combine these 3 fundamental building blocks of machine learning in this simple and powerful fashion. This can be viewed as a form of deep learning, where any of the intermediate layers can itself be deep. The utility and robustness of the methodology is illustrated on some toy problems, including one example problem arising from simulation of plasma fusion in a tokamak.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2018

Characterization of fractional cuts of co-solvent enhanced lignocellulosic fractionation lignin isolated by sequential precipitation

We report that lignin valorization is significantly hindered by the intrinsic heterogeneity of its complex structures and variability of biomass feedstocks. Fractionation of lignin can overcome these challenges by producing functionally distinct lignin cuts that can be further tailored to end products. Herein, lignin was extracted and depolymerized from poplar by the co-solvent enhanced lignocellulosic fractionation method with renewable THF to obtain CELF lignin. Several solvents were screened to separate soluble and insoluble fractions from the parent CELF lignin. The ethanol soluble portion was then fractionated into different molecular weight cuts via sequential precipitation of the lignin by reducing the concentration of THF. The physicochemical structures of different CELF lignin cuts were elucidated by GPC and NMR techniques. In conclusio, these results suggest that CELF lignin cuts with lower molecular weight contain progressively higher phenolic and carboxylic acid OH groups, which can be more suitable as green antioxidants than the parent lignin.

09 BIOMASS FUELS↗
NASA NTRS1993

Double layer propagation in experiments with electron beam injection

Electron beam injection into a plasma is investigated using the analytical inverted Bernstein-Green-Kruskal method. Particle number and momentum conservation laws are applied to evaluate the propagation velocity and potential drop on the leading edge of the beam. Electric potential is supposed to be monotonic, thus the leading front has a double-layer-like structure. For the case of cold particles, analytical expressions for the double layer velocity and potential drop are obtained. It is pointed out that double layer velocity differs from the initial electron speed: even for weak beams a noticeable deceleration takes place. Strong beams are found incapable of penetrating into plasma - their propagation velocity is very small. Ambient electrons undergo a considerable acceleration forming a return current which neutralizes the injector. Possible instability of the distribution functions is discussed.

Bruskin, L. G.↗
NASA NTRS1995

Electromagnetic Scattering From a Rectangular Cavity Recessed in a 3-D Conducting Surface

The problem of electromagnetic (EM) scattering from an aperture backed by a rectangular cavity recessed in a three-dimensional conducting body is analyzed using the coupled field integral equation approach. Using the free-space Green's function, EM fields scattered outside the cavity are determined in terms of (1) an equivalent electric surface current density flowing on the three-dimensional conducting surface of the object including the cavity aperture and (2) an equivalent magnetic surface current density flowing over the aperture only. The EM fields inside the cavity are determined using the waveguide modal expansion functions. Making the total tangential electric and magnetic fields across the aperture continuous and subjecting the total tangential electric field on the outer conducting three-dimensional surface of the object to zero, a set of coupled integral equations is obtained. The equivalent electric and magnetic surface currents are then obtained by solving the coupled integral equation using the Method of Moments. The numerical results on scattering from rectangular cavities embedded in various three-dimensional objects are compared with the results obtained by other numerical techniques.

Deshpande, M. D.↗
NASA NTRS2017

Cyclotron Resonant Scattering Feature Simulations II. Description of the CRSF Simulation Process

Context. Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons o_ quantized plasma electrons in the strong magnetic field (of the order 1012 G) close to the surface of an accreting X-ray pulsar. Due to the complex scattering cross-sections, the line profiles of CRSFs cannot be described by an analytic expression. Numerical methods, such as Monte Carlo (MC) simulations of the scattering processes, are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems.Aims. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time.Methods. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time-consuming simulations possible on convenient time scales. Furthermore, it can generate responses to monoenergetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua.Results. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online (see link in footnote, page 1).

Schwarm, F.-W.↗
DOE OSTI2021

High-Resolution Spatio-Temporal Estimation of Net Ecosystem Exchange in Ice-Wedge Polygon Tundra Using In Situ Sensors and Remote Sensing Data

Land-atmosphere carbon exchange is known to be extremely heterogeneous in arctic ice-wedge polygonal tundra regions. In this study, a Kalman filter-based method was developed to estimate the spatio-temporal dynamics of daytime average net ecosystem exchange (NEEday) at 0.5-m resolution over a 550 m by 700 m study site. We integrated multi-scale, multi-type datasets, including normalized difference vegetation indices (NDVIs) obtained from a novel automated mobile sensor system (or tram system) and a greenness index map obtained from airborne imagery. We took advantage of the significant correlations between NDVI and NEEday identified based on flux chamber measurements. The weighted average of the estimated NEEday within the flux-tower footprint agreed with the flux tower data in term of its seasonal dynamics. We then evaluated the spatial variability of the growing season average NEEday, as a function of polygon geomorphic classes; i.e., the combination of polygon types—which are known to present different degradation stages associated with permafrost thaw—and microtopographic features (i.e., troughs, centers and rims). Our study suggests the importance of considering microtopographic features and their spatial coverage in computing spatially aggregated carbon exchange.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Fabrication and characterization of powder metallurgy tantalum components prepared by high compaction pressure technique

The present study has investigated the consolidation behaviors of tantalum powders during compaction and sintering, and the characteristics of sintered components. For die compaction, the densification behaviors of the powders are simulated by finite element analyses based on the yield function proposed by Shima and Oyane. Accordingly, the green density distribution for coarser particles is predicted to be more uniform because they exhibits higher initial relative tap density owing to lower interparticle friction. It is also found that cold isostatic pressing is capable of producing higher dense compacts compared to the die pressing. However, unlike the compaction behavior, the sintered density of smaller particles is found to be higher than those of coarser ones owing to their higher specific surface area. The maximum sintered density was found to be 0.96 of theoretical density where smaller particles were pressed isostatically at 400 MPa followed by sintering at 2000 °C. Moreover, the effects of processing conditions on grain size and texture were also investigated. The average grain size of the sintered specimen is 30.29 μm and its texture is less than 2 times random intensity. Consequently, it is concluded that the higher pressure compaction technique is beneficial to produce high dense and texture-free tantalum components compared to hot pressing and spark plasma sintering. - Highlights: • Higher Ta density is obtained from higher pressure and sintering temperature. • High compaction method enables P/M Ta to achieve the density of 16.00 g·cm{sup −3}. • A P/M Ta component with fine microstructure and random orientation is developed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Three-dimensional analysis of microbunched electron cooling

Microbunched electron cooling (MBEC, Ratner, ), a special category of coherent electron cooling, has been under consideration as a candidate method for the strong hadron cooling component of the Electron Ion Collider (EIC) at BNL. Most of the theoretical treatments of this scheme so far have been based on simplified models for the key space charge effect in which electrons and/or hadrons are approximated by charged disks. Although such disk-based models provide a flexible tool that captures a good part of the underlying physics, they lack the fidelity and robustness of a rigorous analysis that explicitly treats the hadron and cooler electron beams as collections of point charges. In this paper, we present such a rigorous, three-dimensional (3D) model, both from a theory-based and a simulation perspective. Specifically, we calculate the generalized wakefield (or Green’s function) of the cooling system—a crucial quantity that determines the overall performance—paying special attention to the dependence of the wake on the hadron transverse positions. Additionally, our treatment takes into account previously neglected 3D effects that can cause damping of the plasma oscillations in the amplification section of the cooler, such as focusing and the angular spread of the electron beam. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Electrode-molecule energy level offsets in a gold-benzene diamine-gold single molecule tunnel junction

One means for describing electron transport across single molecule tunnel junctions (MTJs) is to use density functional theory (DFT) in conjunction with a nonequilibrium Green’s function formalism. This description relies on interpreting solutions to the Kohn–Sham (KS) equations used to solve the DFT problem as quasiparticle (QP) states. Many practical DFT implementations suffer from electron self-interaction errors and an inability to treat charge image potentials for molecules near metal surfaces. For MTJs, the overall effect of these errors is typically manifested as an overestimation of electronic currents. Correcting KS energies for self-interaction and image potential errors results in MTJ current–voltage characteristics in close agreement with measured currents. An alternative transport approach foregoes a QP picture and solves for a many-electron wavefunction on the MTJ subject to open system boundary conditions. It is demonstrated that this many-electron method provides similar results to the corrected QP picture for electronic current. The analysis of these two distinct approaches is related through corrections to a junction’s electronic structure beyond the KS energies for the case of a benzene diamine molecule bonded between two gold electrodes.

Szepieniec, Mark S.↗
DOE OSTI2025

Archaeal protein containing domain of unknown function 2193 undergoes oligomeric reconfiguration upon iron–sulfur cluster binding

Methanogenic archaea are particularly rich in iron–sulfur proteins, yet their roles remain largely enigmatic. Here, we characterized aMethanococcus voltae(Mvo) protein from the domain of unknown function (DUF) 2193 family, a group of proteins present primarily in archaea and characterized by a conserved cysteine‐rich C‐terminal motif.MvoDUF2193 was heterologously expressed and characterized by a range of spectroscopic and analytical methods. The results demonstrate thatMvoDUF2193 binds a single [4Fe–4S] cluster per subunit and that cluster occupancy regulates the transition from an apo tetramer to a [4Fe–4S] monomeric form. We hypothesize thatMvoDUF2193 serves a regulatory role in the cell, mediated by [Fe–S] cluster binding and changes in oligomeric state.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Millimeter-Scale Exfoliation of hBN with Tunable Flake Thickness for Scalable Encapsulation

As a two-dimensional (2D) dielectric material, hexagonal boron nitride (hBN) is in high demand for applications in photonics, nonlinear optics, and nanoelectronics. Unfortunately, the high-throughput preparation of macroscopic-scale, high-quality hBN flakes with controlled thickness is an ongoing challenge, limiting device fabrication and technological integration. In this report we present a metal thin-film exfoliation method to prepare hBN flakes with millimeter-scale dimension, near-unity yields, and tunable flake thickness distribution from 1-7 layers, a substantial improvement over scotch tape exfoliation. The single crystallinity and high quality of the exfoliated hBN are demonstrated with optical microscopy, atomic force microscopy, Raman spectroscopy, and second harmonic generation. We further explore a possible mechanism for the effectiveness and selectivity based on thin-film residual stress measurements, density functional theory calculations, and transmission electron microscopy imaging of the deposited metal films. We find that the magnitude of the residual tensile stress induced by thin film deposition plays a key role in determining exfoliated flake thickness in a manner which closely resembles 3D semiconductor spalling. Lastly, we demonstrate that our exfoliated, large-area hBN flakes can be readily incorporated as encapsulating layers for other 2D monolayers. Altogether, this method brings us one step closer to the high throughput, mass production of hBN-based 2D photonic, optoelectronic, and quantum devices.

2D materials↗
DOE OSTI2026

Engineering Escherichia coli for Urease-Driven Synthesis of Metal Oxide Nanomaterials

The development of functional nanomaterials with controlled morphologies is essential for advancements in medicine, electronics and computing, energy, catalysis, and environmental applications. However, conventional synthesis methods often demand high energy input and pose significant environmental challenges. Urease-based biomineralization presents an efficient, eco-friendly alternative for nanomaterial production under mild conditions. In this study, we engineered Escherichia coli ( E. coli ) to express a urease gene cluster from Sporosarcina pasteurii using CRAGE-Duet technology. The engineered strain successfully synthesized calcium carbonate and calcium phosphate crystals. Expanding the approach, we synthesized metal oxide nanoparticles, including hematite (Fe 2 O 3 ), and nanocrystalline anatase titanium dioxide (TiO 2 ). These nanomaterials were characterized by electron microscopy, demonstrating the potential of E. coli as a sustainable and versatile platform for green nanomaterial synthesis.

bacteria↗
DOE OSTI2023

Defects go green: using defects in nanomaterials for renewable energy and environmental sustainability

Induction of point defects in nanomaterials can bestow upon them entirely new physics or augment their pre-existing physical properties, thereby expanding their potential use in green energy technology. Predicting structure-property relationships for defects a priori is challenging, and developing methods for precise control of defect type, density, or structural distribution during synthesis is an even more formidable task. Hence, tuning the defect structure to tailor nanomaterials for enhanced device performance remains an underutilized tool in materials design. We review here the state of nanomaterial design through the lens of computational prediction of defect properties for green energy technology, and synthesis methods to control defect formation for optimal performance. We illustrate the efficacy of defect-focused approaches for refining nanomaterial physics by describing several specific applications where these techniques hold potential. Most notably, we focus on quantum dots for reabsorption-free solar windows and net-zero emission buildings, oxide cathodes for high energy density lithium-ion batteries and electric vehicles, and transition metal dichalcogenides for electrocatalytic green hydrogen production and carbon-free fuels.

14 SOLAR ENERGY↗
DOE OSTI2021

Chemometric optimisation of pressurised liquid extraction for the determination of alliin and S‐allyl‐cysteine in giant garlic ( Allium ampeloprasum L.) by liquid chromatography tandem mass spectrometry

Abstract Introduction Giant garlic is a functional food that contains different kinds of bioactive molecules with beneficial effects on chronic noncommunicable diseases like diabetes and cardiovascular conditions. Considering biosynthesis pathways, abundance, and biological activity, alliin and S‐allyl‐cysteine were used as chemical markers of organosulphur compounds present in giant garlic. Objective To establish a chemometric optimisation of pressurised liquid extraction for the determination of alliin and S‐allyl‐cysteine in giant garlic by liquid chromatography tandem mass spectrometry (LC–MS/MS). Methodology Samples were blanched (ca. 90°C for 10 min) to inactivate alliinase and γ‐glutamyl transpeptidase enzymes and then freeze‐dried. Chemometric optimisation was performed via response surface methodology based on central composite design (CCD). Organosulphur compound yields were determined applying a validated LC–MS/MS method in multiple reaction monitoring (MRM) mode using the following transitions: for alliin m/z 178 → 74 and for S‐allyl‐cysteine m/z 162 → 41. Results According to CCD results, under constant conditions of pressure (1500 psi) and time (20 min), the optimal conditions for pressurised liquid extraction of alliin and S‐allyl‐cysteine were 70.75 and 68.97% v/v of ethanol in extraction solvent and 76.45 and 98.88°C as extraction temperature, respectively. Multiple response optimisation for the simultaneous extraction of both organosulphur compounds was established via desirability function. Under these conditions, 2.70 ± 0.27 mg g −1 dry weight (DW) of alliin and 2.79 ± 0.22 mg g −1 DW of S‐allyl‐cysteine were extracted. Conclusions These results clearly demonstrated that pressurised liquid extraction is an efficient green technique to extract bioactive organosulphur compounds from giant garlic. Extraction yields were significantly ( p < 0.05) higher than those obtained with conventional ultra‐turrax extraction.

Peterssen‐Fonseca, Darlene↗
DOE OSTI2017

Design of Reduction Process of SnO 2 by CH 4 for Efficient Sn Recovery

We design a novel method for the CH 4 reduction of SnO 2 for the efficient recovery of Sn from SnO 2 through a study combining theory and experiment. The atomic-level process of CH 4 -SnO 2 interaction and temperature-dependent reduction behavior of SnO 2 were studied with a combination of a multi-scale computational method of thermodynamic simulations and density functional theory (DFT) calculations. We found that CH 4 was a highly efficient and a versatile reducing agent, as the total reducing power of CH 4 originates from the carbon and hydrogen of CH 4 , which sequentially reduce SnO 2 . Moreover, as a result of the CH 4 reduction of SnO 2 , a mixture of CO and H 2 was produced as a gas-phase product (syngas). The relative molar ratio of the produced gas-phase product was controllable by the reduction temperature and the amount of supplied CH 4 . The laboratory-scale experimental study confirmed that CH 4 actively reduces SnO 2 , producing 99.34% high-purity Sn and H 2 and CO. Our results present a novel method for an efficient, green, and economical recycling strategy for Sn with economic value added that is held by the co-produced clean energy source (syngas).

Science & Technology - Other Topics↗
DOE OSTI2024

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗