Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization

Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite-difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 12.78 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.

17 WIND ENERGY↗
DOE OSTI2023

Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data

In this work, we present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.

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

Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities

Here, we present a framework for calibration of parameters in elastoplastic constitutive models that is based on the use of automatic differentiation (AD). The model calibration problem is posed as a partial differential equation-constrained optimization problem where a finite element (FE) model of the coupled equilibrium equation and constitutive model evolution equations serves as the constraint. The objective function quantifies the mismatch between the displacement predicted by the FE model and full-field digital image correlation data, and the optimization problem is solved using gradient-based optimization algorithms. Forward and adjoint sensitivities are used to compute the gradient at considerably less cost than its calculation from finite difference approximations. Through the use of AD, we need only to write the constraints in terms of AD objects, where all of the derivatives required for the forward and inverse problems are obtained by appropriately seeding and evaluating these quantities. We present three numerical examples that verify the correctness of the gradient, demonstrate the AD approach's parallel computation capabilities via application to a large-scale FE model, and highlight the formulation's ease of extensibility to other classes of constitutive models.

42 ENGINEERING↗
DOE OSTI2022

An efficient method to propagate model uncertainty when inverting seismic data for time domain seismic moment tensors

SUMMARY We present a computationally efficient method to approximately propagate uncertainty when linearly inverting seismic data for point source, time variable moment tensor components. The method is based on the assumption that the data residual, given by the difference between the observed seismic data and the data predicated by a linear inversion, contains the effects of both data and model uncertainty. Our method uses a distribution of data residuals, added directly to the data, in a pseudo-Monte Carlo scheme. Using the assumption that the data residual is a stochastic process, we use the well-known Karhunen–Loève (KL) theorem to construct a distribution of data residuals, where the required basis functions are constructed using Fourier series. The Fourier series are scaled by a product of a random variable and the real-valued spectral amplitudes of the original data residual’s spectrum. Thus, the Fourier series and spectral amplitudes are eigenfunction-eigenvalue pairs used in the KL-based construction of data residual distribution. Using tests with synthetic data, we show that our method compares closely with a Finite Difference Monte Carlo (FDMC) method that we presented previously. More importantly, the method presented here is computationally several orders of magnitude faster than our previous FDMC method, and requires no a priori assumptions of model and/or data uncertainty.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2025

Solving high-dimensional partial integral differential equations: The finite expression method

Partial integro-differential equations (PIDEs) have broad applications in the sciences, from electro-magnetism to options pricing. Here, in this paper, we introduce a new finite expression method (FEX) to solve PIDEs. This approach builds upon the original FEX and its inherent advantages with new advances: 1) A novel method of parameter grouping is proposed to reduce the number of coefficients in high-dimensional function approximation; 2) A Taylor series approximation method is implemented to significantly improve the computational efficiency and accuracy of the evaluation of the integral terms of PIDEs. The new FEX based method, denoted FEX-PG to indicate the addition of the parameter grouping (PG) step to the algorithm, provides both high accuracy and interpretable numerical solutions, with the outcome being an explicit equation that facilitates intuitive understanding of the underlying solution structures. These features are often absent in traditional methods, such as finite element methods (FEM) and finite difference methods, as well as in deep learning-based approaches. To benchmark our method against recent advances, we apply the new FEX-PG to solve benchmark PIDEs in the literature. In high-dimensional settings, FEX-PG exhibits strong and robust performance, achieving relative errors on the order of single precision machine epsilon, significantly outperforming existing approaches based on neural networks.

Combinatorial optimization↗
DOE OSTI2024

Efficient sensitivity analysis of the thermal profile in powder bed fusion of metals using hypercomplex automatic differentiation finite element method

Rapid cyclic temperature fluctuation occurring in powder bed fusion of metals using a laser beam (PBF-LB/M) influences the formation of flaws in printed parts. Consequently, there is a pressing need to enhance the quality of printed parts by developing innovative methodologies that can predict thermal histories and help uncover the intricate relationships between process parameters and thermal profiles. Sensitivity Analysis (SA) emerges as an essential tool for this, offering the potential for process optimization and enhanced quality control. Nonetheless, conventional SA methodologies often incur in excessive computational costs and potential numerical approximation errors. Here, to address this technical challenge, we present a novel method for SA that integrates the HYPercomplex-based Automatic Differentiation (HYPAD) technique with transient thermal simulations conducted via the finite element method (FEM). Leveraging this methodology, we efficiently and accurately perform SA for PBF-LB/M processes in a post-processing step. Compared to traditional methods like Finite Differences (FD), HYPAD-FEM required 96 % less computational time for obtaining sensitivities for 22 process parameters, under a comparative study conducted within the context of the 2018–02 AM benchmark of the National Institute of Standards and Technology. In summary, HYPAD-FEM offers superior efficiency and accuracy in SA over conventional methods, delivering the best sensitivity of a model without the need for step-size selection and problem or parameter-based implementations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Electrode and Microstructure Dependence of Oxygen Diffusion in Ferroelectric Hafnium Zirconium Oxide Thin Films

Hafnia-based ferroelectrics hold promise to reduce energy demand for computing by enabling compute-in-memory and as non-volatile memories. The ferroelectric phase in this material system is, in part, stabilized by oxygen vacancies. While oxygen vacancies may be a necessity for phase stability, they limit device endurance through diffusion and accumulation into conducting channels. Herein, it is shown that oxygen diffusion is spatially variable within individual grains of ferroelectric hafnium zirconium oxide (HZO). Using 18 O tracers and finite difference modeling, it is shown that grain boundaries and regions near electrode interfaces allow for relatively rapid oxygen diffusion, with values as much as 10 4 larger than the grain cores. Further, the selection of electrode material affects the diffusion coefficients across all microstructural regions. HZO films in contact with TiN electrodes result in more oxygen-deficient HZO films and higher oxygen diffusion coefficients. Tungsten electrodes result in fewer vacancies and lower diffusion coefficients. Diffusion activation energy differences between the HZO with the two electrodes is reconciled by differing populations of charged and uncharged oxygen vacancies. This insight into the local vacancy populations and diffusion pathways provides a platform for designing hafnia-based films, deposition processes, and integration strategies to reduce vacancy gradients and improve performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structural Color Production in Melanin-Based Disordered Colloidal Nanoparticle Assemblies in Spherical Confinement

Melanin is a ubiquitous natural pigment that exhibits broadband absorption and high refractive index. Despite its widespread use in structural color production, how the absorbing material, melanin, affects the generated color is unknown. Using a combined molecular dynamics and finite-difference time-domain computational approach, here we investigates structural color generation in one-component melanin nanoparticle-based supraparticles (called supraballs) as well as binary mixtures of melanin and silica (nonabsorbing) nanoparticle-based supraballs. Experimentally produced one-component melanin and one-component silica supraballs, with thoroughly characterized primary particle characteristics using neutron scattering, produce reflectance profiles similar to the computational analogs, confirming that the computational approach correctly simulates both absorption and multiple scattering from the self-assembled nanoparticles. These combined approaches demonstrate that melanin's broadband absorption increases the primary reflectance peak wavelength, increases saturation, and decreases lightness factor. Additionally, the dispersity of nanoparticle size more strongly influences the optical properties of supraballs than packing fraction, as evidenced by the production of a larger range of colors when size dispersity is varied versus packing fraction. For binary melanin and silica supraballs, the chemistry-based stratification allows for more diverse color generation and finer saturation tuning than does the degree of mixing/demixing between the two chemistries.

broadband absorption↗
DOE OSTI2025

Imaging Photonic Resonances within an All‐Dielectric Metasurface via Photoelectron Emission Microscopy

Dielectric metasurfaces, through volume‐type photonic resonances, enable precise control of light‐matter interactions for applications including imaging, holography, and sensing. The application space of dielectric metasurfaces has extended from infrared to visible wavelengths by incorporating high refractive index materials, such as titanium dioxide (TiO 2 ). Understanding the fundamental and fabrication limits for these applications requires metrology with nanoscale resolution, sensitivity to electromagnetic fields within the meta‐atom volume, and far‐field excitation. In this work, photoelectron emission microscopy (PEEM) is used to image field distributions of photonic resonances in a TiO 2 metasurface excited with far‐field, visible‐wavelength illumination. The local volumetric field variations within the meta‐atoms are analyzed as a function of illumination angle and polarization by comparing photoelectron images to finite‐difference time‐domain simulations. This study determines the inelastic mean free path of very low‐energy (<1 eV) photoelectrons to be 35 ± 10 nm, which is comparable to the meta‐atom height thereby highlighting PEEM sensitivity to resonances within the volume. Additionally, the simulations reveal high sensitivity of PEEM images to an in‐plane component of the illumination k ‐vector. These results demonstrate that photoelectron imaging with subwavelength resolution offers unique advantages for examining light‐matter interactions in volume‐type (as opposed to surface) photonic modes within dielectric nanophotonic structures.

electron inelastic mean free path↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
DOE OSTI2020

Photonic devices fabricated from (111)‐oriented single crystal diamond

Abstract Diamond is a material of choice in the pursuit of integrated quantum photonic technologies. So far, the majority of photonic devices fabricated from diamond are made from (100)‐oriented crystals. In this work, we demonstrate a methodology for the fabrication of optically active membranes from (111)‐oriented diamond. We use a liftoff technique to generate membranes, followed by chemical vapor deposition of diamond in the presence of silicon to generate homogenous silicon vacancy color centers with emission properties that are superior to those in (100)‐oriented diamond. We further use the diamond membranes to fabricate microring resonators with quality factors exceeding ~ 3000. Supported by finite‐difference time‐domain calculations, we discuss the advantages of (111)‐oriented structures as building blocks for quantum nanophotonic devices. image

Regan, Blake↗
DOE OSTI2023

Determining the twist angle of stacked MoS 2 layers using machine learning-assisted low-frequency interlayer Raman fingerprints

The investigation of twisted stacked few-layer MoS 2 has revealed novel electronic, optical, and vibrational properties over an extended period. For the successful integration of twisted stacked few-layer MoS 2 into a wide range of applications, it is crucial to employ a noninvasive, versatile technique for characterizing the layered architecture of these complex structures. Here, in this work, we introduce a machine learning-assisted low-frequency Raman spectroscopy method to characterize the twist angle of few-layer stacked MoS 2 samples. A feedforward neural network (FNN) is utilized to analyze the low-frequency breathing mode as a function of the twist angle. Moreover, using finite difference method (FDM) and density functional theory (DFT) calculations, we show that the low-frequency Raman spectra of MoS 2 are mainly influenced by the effect of the nearest and second nearest layers. A new improved linear chain model (TA-LCM) with taking the twist angle into the consideration is developed to understand the interlayer breathing modes of stacked few-layer MoS 2 . This approach can be extended to other 2D materials systems and provides an intelligent way to investigate naturally stacked and twisted interlayer interactions.

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

Effect of XC functionals and dispersion corrections on the DFT‐computed structural and vibrational properties of SrCl 2 –NaCl and ZrF 4 –LiF

Density functional theory (DFT) calculations were performed to examine the impact of exchange–correlation (XC) functionals and van der Waals corrections (specifically the D3 method) on the structural and vibrational properties of the SrCl 2 –NaCl and ZrF 4 –LiF salt systems. Multiple XC functionals, including the local density approximation (LDA), the generalized gradient approximation using the Perdew–Burke–Ernzerhof (PBE) model, and its modified form suitable for solids (PBEsol), the dispersion-corrected PBE-D3 and PBEsol-D3, were considered. Of these functionals, LDA was found to exhibit the highest degree of error, while PBEsol and PBE-D3 displayed the least error. Underestimated lattice parameters compared with experimental values were observed to result in higher force constants, leading to an overprediction of vibrational frequencies. Conversely, an overestimation of lattice parameters was associated with lower vibrational frequencies. The methodology presented in this study yielded results that are in good agreement with experiment, irrespective of the method (finite differences vs. density functional perturbation theory) employed for calculating infrared and Raman spectra. It was further demonstrated that for alkali halides with weak Raman scattering, utilizing a supercell constructed from primitive cells better predicts Raman features than does the use of conventional cells.

Raman↗
DOE OSTI2025

Tailored Frequency Conversion in Nonlinear Subwavelength Grating Metaslabs

The use of metaslab waveguides consisting of subwavelength gratings (SWGs) is proposed as a highly flexible approach to obtain phase‐matched frequency conversion in nonlinear integrated devices. Control over the phase matching (PM) condition is achieved through modification of the linear modal propagation constants by the SWG. The method is experimentally demonstrated by fabrication of 0.4 mm‐long GaAs metaslabs designed for second harmonic generation of near‐infrared light. A controlled spectral shift of up to 100 nm in the PM condition, relative to a uniform slab is observed. Excellent agreement is found between the experimental results and finite‐difference eigenmode simulations, which predict a tunable frequency conversion bandwidth of ∼300 nm for the design parameters. The different regimes of PM control supported by the structure are discussed. This work highlights the potential of SWGs as building blocks for tailored frequency conversion in nonlinear integrated devices.

Form birefringence↗