Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Estimation of Syringyl Units in Wood Lignins by FT-Raman Spectroscopy

Syringyl (S) lignin content and the syringyl-to-guaiacyl (S/G) lignin ratio are important characteristics of wood and lignocellulosic biomass. While numerous methods are available for estimating S lignin units and the S/G ratio, in this work, a new method based on Raman spectroscopy that uses the 370 cm –1 Raman band-area intensity (370-area) was developed. The reliability of the Raman method for determining S content was first tested by the quantitative analysis of three syringyl lignin models by sampling them, separately, in dioxane and in Avicel. Good linear correlations between the 370 cm –1 intensity and model concentrations were obtained. Next, the percent syringyl (%S) lignin units in various woods were measured by correlating the 370 cm –1 Raman intensity data with values of S units in lignin determined by three regularly used methods, namely, thioacidolysis, DFRC, and 2D-HSQC NMR. The former two methods take into account only the monomers cleaved from β–O–4-linked lignin units, whereas the NMR method reports S content on the whole cell wall lignin. When the 370-area intensities and %S values from the regularly used methods were correlated, good linear correlations were obtained (R 2 = 0.767, 0.731, and 0.804, respectively, for the three methods). The correlation with the highest R 2 , i.e., with the 2D NMR approach, is introduced for estimating S units in wood lignins by Raman spectroscopy as, in principle, both represent the whole cell wall lignin and not just the portion of lignin that gets cleaved to release monomers. The Raman analysis method is quick, uses minimal harmful chemicals, is carried out nondestructively, and is insensitive to the wet or dry state of the sample. The only limitations are that the sample of wood contains at least 30% S and not be significantly fluorescent, although the latter can be mitigated in some cases.

09 BIOMASS FUELS↗
DOE OSTI2025

PSCF+: An Extended and Improved Open-Source Software Package for Polymer Self-Consistent Field Calculations

This work introduces PSCF+, a freely available open‑source software tool for modeling how block copolymers self‑assemble into ordered nanoscale structures. The software allows researchers to study a wide range of polymer models and interaction types, making it possible to directly compare theoretical predictions with molecular and simulation‑based studies in a consistent way. PSCF+ incorporates several recently developed numerical methods that significantly reduce computer memory requirements and speed up calculations, enabling efficient use of modern GPU hardware. Here, it also provides automated tools to calculate free‑energy landscapes and phase boundaries, which are essential for understanding and predicting material behavior. Together, these advances make PSCF+ a powerful and flexible platform for studying polymer self‑assembly and support the design of new polymeric materials for scientific and technological applications.

Chemical calculations↗
DOE OSTI2022

Numerical Simulation of Vibrational Sum Frequency Generation Intensity for Non-Centrosymmetric Domains Interspersed in an Amorphous Matrix: A Case Study for Cellulose in Plant Cell Wall

Vibrational sum frequency generation (SFG) spectroscopy can specifically probe molecular species non-centrosymmetrically arranged in a centrosymmetric or isotropic medium. This capability has been extensively utilized to detect and study molecular species present at the two-dimensional (2D) interface at which the centrosymmetry or isotropy of bulk phases is naturally broken. The same principle has been demonstrated to be very effective for the selective detection of non-centrosymmetric crystalline nanodomains interspersed in three-dimensional (3D) amorphous phases. However, the full spectral interpretation of SFG features has been difficult due to the complexity associated with the theoretical calculation of SFG responses of such 3D systems. This paper describes a numerical method to predict the relative SFG intensities of non-centrosymmetric nanodomains in 3D systems as functions of their size and concentration as well as their assembly patterns, i.e., the distributions of tilt, azimuth, and rotation angles with respect to the lab coordinate. We applied the developed method to predict changes in the CH and OH stretch modes characteristic to crystalline cellulose microfibrils distributed with various orders, which are relevant to plant cell wall structures. As a result, the same algorithm can also be applied to any SFG-active nanodomains interspersed in 3D amorphous matrices.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Understanding How Acoustic Vibrations Modulate the Optical Response of Plasmonic Metal Nanoparticles

Measurements of acoustic vibrations in nanoparticles provide a unique opportunity to study mechanical phenomena at nanometer length scales and picosecond time scales. Vibrations in noble-metal nanoparticles have attracted particular attention, because they couple to plasmon resonances in the nanoparticles, leading to strong modulation of optical absorption and scattering. There are three mechanisms that transduce the mechanical oscillations into changes in the plasmon resonance: (1) changes in polarizability due to changes in the nanoparticle geometry; (2) changes in electron density due to changes in the nanoparticle volume and (3) changes in the interband transition energies due to compression/expansion of the nanoparticle (deformation potential). These mechanisms have been studied in the past to explain the origin of the experimental signals; however, a thorough quantitative connection between the coupling of phonon and plasmon modes and separate contribution of each coupling mechanism has not yet been made. Here, we present a numerical method to quantitatively determine the coupling between vibrational and plasmon modes in noble-metal nanoparticles of arbitrary geometries, and apply it to spheres, shells, rods, and cubes in the context of time resolved measurements. We separately determine the parts of the optical response that are due to shape changes, changes in electron density, and changes in deformation potential (DP). We further show that coupling is in general strongest when the regions of largest electric field (plasmon mode) and largest displacement (phonon mode) overlap. Lastly, these results clarify reported experimental results, and should help guide future experiments and potential applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

The GFDL‐CM4X Climate Model Hierarchy, Part II: Case Studies

This paper is Part II of a two‐part paper that documents the Climate Model version 4X (CM4X) hierarchy of coupled climate models developed at the Geophysical Fluid Dynamics Laboratory. Part I of this paper is presented in Griffies et al. (2025a, https://doi.org/10.1029/2024MS004861 ). Here we present a suite of case studies that examine ocean and sea ice features that are targeted for further research, which include sea level, eastern boundary upwelling, Arctic and Southern Ocean sea ice, Southern Ocean circulation, and North Atlantic circulation. The case studies are based on experiments that follow the protocol of version 6 from the Coupled Model Intercomparison Project. The analysis reveals a systematic improvement in the simulation fidelity of CM4X relative to its CM4.0 predecessor, as well as an improvement when refining the ocean/sea ice horizontal grid spacing from the 0.25° of CM4X‐p25 to the 0.125° of CM4X‐p125. Even so, there remain many outstanding biases, thus pointing to the need for further grid refinements, enhancements to numerical methods, and/or advances in parameterizations, each of which target long‐standing model biases and limitations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Vibrational hierarchy leads to dual-phonon transport in low thermal conductivity crystals

Many low-thermal-conductivity ( κ L ) crystals show intriguing temperature ( T ) dependence of κ L : κ L ∝ T -1 (crystal-like) at intermediate temperatures whereas weak T -dependence (glass-like) at high temperatures. It has been in debate whether thermal transport can still be described by phonons at the Ioffe-Regel limit. In this work, we propose that most phonons are still well defined for thermal transport, whereas they carry heat via dual channels: normal phonons described by the Boltzmann transport equation theory, and diffuson-like phonons described by the diffusion theory. Three physics-based criteria are incorporated into first-principles calculations to judge mode-by-mode between the two phonon channels. Case studies on La 2 Zr 2 O 7 and Tl 3 VSe 4 show that normal phonons dominate low temperatures while diffuson-like phonons dominate high temperatures. Our present dual-phonon theory enlightens the physics of hierarchical phonon transport as approaching the Ioffe-Regel limit and provides a numerical method that should be practically applicable to many materials with vibrational hierarchy.

42 ENGINEERING↗
DOE OSTI2025

One-shot learning for solution operators of partial differential equations

Learning and solving governing equations of a physical system, represented by partial differential equations (PDEs), from data is a central challenge in many areas of science and engineering. Traditional numerical methods can be computationally expensive for complex systems and require complete governing equations. Existing data-driven machine learning methods require large datasets to learn a surrogate solution operator, which could be impractical. Here, we propose a solution operator learning method that requires only one PDE solution, i.e., one-shot learning, along with suitable initial and boundary conditions. Leveraging the locality of derivatives, we define a local solution operator in small local domains, train it using a neural network, and use it to predict solutions of new input functions via mesh-based fixed-point iteration or meshfree neural-network based approaches. We test our method on various PDEs, complex geometries, and a practical spatial infection spread application, demonstrating its effectiveness and generalization capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

A numerical analysis of a magnetocaloric refrigerator with a 16-layer regenerator

A numerical analysis was conducted to study a room temperature magnetocaloric refrigerator with a 16-layer parallel plates active magnetic regenerator (AMR). Sixteen layers of LaFeMnSiH having different Curie temperatures were employed as magnetocaloric material (MCM) in the regenerator. Measured properties data was used. A transient one dimensional (1D) model was employed, in which a unique numerical method was developed to significantly accelerate the simulation speed of the multi-layer AMR system. As a result, the computation speed of a multi-layer AMR case was very close to the single-layer configuration. The performance of the 16-layer AMR system in different frequencies and utilizations has been investigated using this model. To optimize the layer length distribution of the 16-layer MCMs in the regenerator, a set of 137 simulations with different MCM distributions based on the Design of Experiments (DoE) method was conducted and the results were analyzed. The results show that the 16-layer AMR system can operate up to 84% of Carnot cycle COP at a temperature span of 41 K, which cannot be obtained using an AMR with fewer layers. Here, the DoE results indicate that for a 16-layer AMR system, the uniform distribution is very close to the optimized design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Micro-continuum approach for mineral precipitation

Abstract Rates and extents of mineral precipitation in porous media are difficult to predict, in part because laboratory experiments are problematic. It is similarly challenging to implement numerical methods that model this process due to the need to dynamically evolve the interface of solid material. We developed a multiphase solver that implements a micro-continuum simulation approach based on the Darcy–Brinkman–Stokes equation to study mineral precipitation. We used the volume-of-fluid technique in sharp interface implementation to capture the propagation of the solid mineral surface. Additionally, we utilize an adaptive mesh refinement method to improve the resolution of near interface simulation domain dynamically. The developed solver was validated against both analytical solution and Arbitrary Lagrangian–Eulerian approach to ensure its accuracy on simulating the propagation of the solid interface. The precipitation of barite (BaSO 4 ) was chosen as a model system to test the solver using variety of simulation parameters: different geometrical constraints, flow conditions, reaction rate and ion diffusion. The growth of a single barite crystal was simulated to demonstrate the solver’s capability to capture the crystal face specific directional growth.

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

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Stripe correlations in the two-dimensional Hubbard-Holstein model

Several state-of-the-art numerical methods have observed static or fluctuating spin and charge stripes in doped two-dimensional Hubbard models, suggesting that these orders play a significant role in shaping the cuprate phase diagram. Many experiments, however, also indicate that the cuprates have strong electron-phonon (e-ph) coupling, and it is unclear how this interaction influences stripe correlations. We study static and fluctuating stripe orders in the doped single-band Hubbard-Holstein model using zero temperature variational Monte Carlo and finite temperature determinant quantum Monte Carlo. We find that the lattice couples more strongly with the charge component of the stripes, leading to an enhancement or suppression of stripe correlations, depending on model parameters like the next-nearest-neighbor hopping t' or phonon energy Ω. Our results help elucidate how the e-ph interaction can tip the delicate balance between stripe and superconducting correlations in the Hubbard-Holstein model with implications for our understanding of the high-T c cuprates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Magnetic anisotropy and GGG substrate stray field in YIG films down to millikelvin temperatures

Abstract Quantum magnonics investigates the quantum-mechanical properties of magnons, such as quantum coherence or entanglement for solid-state quantum information technologies at the nanoscale. The most promising material for quantum magnonics is the ferrimagnetic yttrium iron garnet (YIG), which hosts magnons with the longest lifetimes. YIG films of the highest quality are grown on a paramagnetic gadolinium gallium garnet (GGG) substrate. The literature has reported that ferromagnetic resonance (FMR) frequencies of YIG/GGG decrease at temperatures below 50 K despite the increase in YIG magnetization. We investigated a 97 nm-thick YIG film grown on 500 μm-thick GGG substrate through a series of experiments conducted at temperatures as low as 30 mK, and using both analytical and numerical methods. Our findings suggest that the primary factor contributing to the FMR frequency shift is the stray magnetic field created by the partially magnetized GGG substrate. This stray field is antiparallel to the applied external field and is highly inhomogeneous, reaching up to 40 mT in the center of the sample. At temperatures below 500 mK, the GGG field exhibits a saturation that cannot be described by the standard Brillouin function for a paramagnet. Including the calculated GGG field in the analysis of the FMR frequency versus temperature dependence allowed the determination of the cubic and uniaxial anisotropies. We find that the total crystallographic anisotropy increases more than three times with the decrease in temperature down to 2 K. Our findings enable accurate predictions of the YIG/GGG magnetic systems behavior at low and ultralow millikelvin temperatures, crucial for developing quantum magnonic devices.

Serha, Rostyslav O.↗
DOE OSTI2015

Conductivity-limiting bipolar thermal conductivity in semiconductors

Intriguing experimental results raised the question about the fundamental mechanisms governing the electron-hole coupling induced bipolar thermal conduction in semiconductors. Our combined theoretical analysis and experimental measurements show that in semiconductors bipolar thermal transport is in general a “conductivity-limiting” phenomenon and it is thus controlled by the carrier mobility ratio and by the minority carrier partial electrical conductivity for the intrinsic and extrinsic cases, respectively. Our numerical method quantifies the role of electronic band structure and carrier scattering mechanisms. We have successfully demonstrated bipolar thermal conductivity reduction in doped semiconductors via electronic band structure modulation and/or preferential minority carrier scatterings. We expect this study to be beneficial to the current interests in optimizing thermoelectric properties of narrow gap semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Extraction of pore-morphology and capillary pressure curves of porous media from synchrotron-based tomography data

The elevated level of atmospheric carbon dioxide (CO 2 ) has caused serious concern of the progression of global warming. Geological sequestration is considered as one of the most promising techniques for mitigating the damaging effect of global climate change. Investigations over wide range of length-scales are important for systematic evaluation of the underground formations from prospective CO 2 reservoir. Understanding the relationship between the micro morphology and the observed macro phenomena is even more crucial. Here we show Synchrotron based X-ray micro tomographic study of the morphological buildup of Sandstones. We present a numerical method to extract the pore sizes distribution of the porous structure directly, without approximation or complex calculation. We have also demonstrated its capability in predicting the capillary pressure curve in a mercury intrusion porosimetry (MIP) measurement. The method presented in this work can be directly applied to the morphological studies of heterogeneous systems in various research fields, ranging from Carbon Capture and Storage, and Enhanced Oil Recovery to environmental remediation in the vadose zone.

2-phase flow↗
DOE OSTI2013

Activity-based protein profiling of secreted cellulolytic enzyme activity dynamics in Trichoderma reesei QM6a, NG14, and RUT-C30

Development of alternative, non-petroleum based sources of bioenergy that can be applied in the short-term find great promise in the use of highly abundant and renewable lignocellulosic plant biomass.1 This material obtained from different feedstocks, such as forest litter or agricultural residues, can yield liquid fuels and other chemical products through biorefinery processes.2 Biofuels are obtained from lignocellulosic materials by chemical pretreatment of the biomass, followed by enzymatic decomposition of cellulosic and hemicellulosic compounds into soluble sugars that are converted to desired chemical products via microbial metabolism and fermentation.3, 4 To release soluble sugars from polymeric cellulose multiple enzymes are required, including endoglucanase, exoglucanase, and β-glucosidase.5, 6 However, the enzymatic hydrolysis of cellulose into soluble sugars remains a significant limiting factor to the efficient and economically viable utilization of lignocellulosic biomass for transport fuels.7, 8 The primary industrial source of cellulose and hemicellulases is the mesophilic soft-rot fungus Trichoderma reesei,9 having widespread applications in food, feed, textile, pulp, and paper industries.10 The genome encodes 200 glycoside hydrolases, including 10 cellulolytic and 16 hemicellulolytic enzymes.11 The hypercellulolytic catabolite derepressed strain RUT-C30 was obtained through a three-step UV and chemical mutagenesis of the original T. reesei strain QM6a,12, 13 in which strains M7 and NG14 were intermediate, having higher cellulolytic activity than the parent strain but less activity and higher catabolite repression than RUT-C30.14 Numerous methods have been employed to optimize the secreted enzyme cocktail of T. reesei including cultivation conditions, operational parameters, and mutagenesis.3 However, creating an optimal and economical enzyme mixture for production-scale biofuels synthesis may take thousands of experiments to identify.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2018

Weak anti-localization of two-dimensional holes in germanium beyond the diffusive regime

Gate-controllable spin-orbit coupling is often one requisite for spintronic devices. For practical spin field-effect transistors, another essential requirement is ballistic spin transport, where the spin precession length is shorter than the mean free path such that the gate-controlled spin precession is not randomized by disorder. In this letter, we report the observation of a gate-induced crossover from weak localization to weak anti-localization in the magneto-resistance of a high-mobility twodimensional hole gas in a strained germanium quantum well. From the magneto-resistance, we extract the phase-coherence time, spin-orbit precession time, spin-orbit energy splitting, and cubic Rashba coefficient over a wide density range. The mobility and the mean free path increase with increasing hole density, while the spin precession length decreases due to increasingly stronger spin-orbit coupling. As the density becomes larger than , 6 x 10 11 cm -2 , the spin precession length becomes shorter than the mean free path, and the system enters the ballistic spin transport regime. We also report here the numerical methods and code developed for calculating the magnetoresistance in the ballistic regime, where the commonly used HLN and ILP models for analyzing weak localization and anti-localization are not valid. These results pave the way toward siliconcompatible spintronic devices.

42 ENGINEERING↗
DOE OSTI2016

Planck 2015 results: XII. Full focal plane simulations

In this paper, we present the 8th full focal plane simulation set (FFP8), deployed in support of the Planck 2015 results. FFP8 consists of 10 fiducial mission realizations reduced to 18 144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising 10 4 mission realizations reduced to about 10 6 maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects, and the remaining subdominant effects will be included in future updates. Finally, generated at a cost of some 25 million CPU-hours spread across multiple high-performance-computing (HPC) platforms, FFP8 is used to validate and verify analysis algorithms and their implementations, and to remove biases from and quantify uncertainties in the results of analyses of the real data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Unexpected topology of the temperature fluctuations in the cosmic microwave background

We study the topology generated by the temperature fluctuations of the cosmic microwave background (CMB) radiation, as quantified by the number of components and holes, formally given by the Betti numbers, in the growing excursion sets. We compare CMB maps observed by the Planck satellite with a thousand simulated maps generated according to the $Λ$CDM paradigm with Gaussian distributed fluctuations. The comparison is multi-scale, being performed on a sequence of degraded maps with mean pixel separation ranging from 0.05 to 7.33°. The survey of the CMB over S 2 is incomplete due to obfuscation effects by bright point sources and other extended foreground objects like our own galaxy. To deal with such situations, where analysis in the presence of “masks” is of importance, we introduce the concept of relative homology. The parametric χ 2 -test shows differences between observations and simulations, yielding p -values at percent to less than permil levels roughly between 2 and 7°, with the difference in the number of components and holes peaking at more than 3 σ sporadically at these scales. The highest observed deviation between the observations and simulations for b 0 and b 1 is approximately between 3 σ and 4 σ at scales of 3–7°. There are reports of mildly unusual behaviour of the Euler characteristic at 3.66° in the literature, computed from independent measurements of the CMB temperature fluctuations by Planck ’s predecessor, the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. The mildly anomalous behaviour of the Euler characteristic is phenomenologically related to the strongly anomalous behaviour of components and holes, or the zeroth and first Betti numbers, respectively. Further, since these topological descriptors show consistent anomalous behaviour over independent measurements of Planck and WMAP, instrumental and systematic errors may be an unlikely source. These are also the scales at which the observed maps exhibit low variance compared to the simulations, and approximately the range of scales at which the power spectrum exhibits a dip with respect to the theoretical model. Non-parametric tests show even stronger differences at almost all scales. Crucially, Gaussian simulations based on power-spectrum matching the characteristics of the observed dipped power spectrum are not able to resolve the anomaly. Understanding the origin of the anomalies in the CMB, whether cosmological in nature or arising due to late-time effects, is an extremely challenging task. Regardless, beyond the trivial possibility that this may still be a manifestation of an extreme Gaussian case, these observations, along with the super-horizon scales involved, may motivate the study of primordial non-Gaussianity. Alternative scenarios worth exploring may be models with non-trivial topology, including topological defect models.

79 ASTRONOMY AND ASTROPHYSICS↗