Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element 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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Convergent laser beam shapes: Unveiling the dynamics of Laser-induced elastic waves in composite materials

Overcoming the low signal-to-noise ratio (SNR) in laser ultrasonic testing of composite materials remains a significant challenge. Current efforts focus on enhancing SNR by inserting more energy into the material through temporal and/or spatial modulation of the laser beam. However, potential SNR improvements through wave convergence and wave energy manipulation have been overlooked. This paper addresses this gap by demonstrating the convergence of different wave types to a designated point and by showing the feasibility of directing absorbed laser energy into a specific wave type through spatial modulation of the laser beam. To achieve this, mathematical expressions for the convergent laser beams are derived. Various laser beam profiles are then introduced to the thermoelastic equations and solved using the finite element method. The sample under investigation is a transversely isotropic unidirectional carbon fiber reinforced plastic, characterized by anisotropic thermal expansion coefficients and thermal conductivities. Results reveal pronounced convergence of the intended wave type at the center due to laser beam shaping. This study showcases the ability to direct absorbed laser energy toward a specific wave type through spatial modulation of the laser beam and highlights the role of material anisotropy in energy focusing.

composite materials↗
DOE OSTI2016

A transmissive, spectrum-splitting concentrating photovoltaic module for hybrid photovoltaic-solar thermal energy conversion

Here, a spectrum splitting, transmissive concentrating photovoltaic (tCPV) module is proposed and designed for a hybrid photovoltaic-solar thermal(PV/T) system. This module enables the system to more fully utilize the full spectrum of incoming sunlight. By utilizing III-V triple junction solar cells with bandgaps of 2.1eV/1.7eV/1.4eVin the module, ultraviolet (UV) and visible light will be absorbed and converted to electricity, while infrared (IR) light will pass through and be captured by a solar thermal receiver and stored as heat. The stored heat energy may be dispatched as electricity or process heat, as needed.Numerical analysis based on the Finite Elemental Method (FEM), ray tracing, circuit analysis, device modelling, and more is employed to design this tCPV module. According to these simulations, the tCPV module can perform with overall power conversion efficiency exceeding 43.5% for above bandgap (in-band)light under a standard AM1.5D solar spectrum, under an average concentration ratio of 400 suns.Passive and active cooling methods, keeping cells below 110°C, are also investigated and discussed, indicating that a transparent active cooling design could improve the CPV module efficiency by around 1%(absolute), relative to a passive design,by reducing the maximum cell working temperature by ~16°C. Furthermore, cost analysis shows that installation cost of around 1.9/W–2.2/W could be reached for the tCPV based PV/T system,which shows a competitive economic advantage compared to a more conventional PV with battery system.

14 SOLAR ENERGY↗
DOE OSTI2016

Multilevel sequential Monte Carlo samplers

Here, we study the approximation of expectations w.r.t. probability distributions associated to the solution of partial differential equations (PDEs); this scenario appears routinely in Bayesian inverse problems. In practice, one often has to solve the associated PDE numerically, using, for instance finite element methods and leading to a discretisation bias, with the step-size level h L . In addition, the expectation cannot be computed analytically and one often resorts to Monte Carlo methods. In the context of this problem, it is known that the introduction of the multilevel Monte Carlo (MLMC) method can reduce the amount of computational effort to estimate expectations, for a given level of error. This is achieved via a telescoping identity associated to a Monte Carlo approximation of a sequence of probability distributions with discretisation levels ${\infty}$ >h 0 >h 1 ...>h L . In many practical problems of interest, one cannot achieve an i.i.d. sampling of the associated sequence of probability distributions. A sequential Monte Carlo (SMC) version of the MLMC method is introduced to deal with this problem. In conclusion, it is shown that under appropriate assumptions, the attractive property of a reduction of the amount of computational effort to estimate expectations, for a given level of error, can be maintained within the SMC context.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Bifurcation of equilibrium positions for ellipsoidal particles in inertial shear flows between two walls

We conducted a systematic numerical investigation of spherical, prolate and oblate particles in an inertial shear flow between two parallel walls, using smoothed particle hydrodynamics (SPH). It was previously shown that above a critical Reynolds number, spherical particles experience a supercritical pitchfork bifurcation of the equilibrium position in shear flow between two parallel walls, namely that the central equilibrium position becomes unstable, leading to the emergence of two new off-centre stable positions (Foxet al.,J. Fluid Mech., vol. 915, 2021). This phenomenon was unexpected given the symmetry of the system. In addition to confirming this finding, we found, surprisingly, that ellipsoidal particles can also return to the centre position from the off-centre positions when the particle Reynolds number is further increased, while spherical particles become unstable under this increased Reynolds number. By utilizing both SPH and the finite element method for flow visualization, we explained the underlining mechanism of this reverse of bifurcation by altered streamwise vorticity and symmetry breaking of pressure. Furthermore, we expanded our investigation to include asymmetric particles, a novel aspect that had not been previously modelled, and we observed similar trends in particle dynamics for both symmetric and asymmetric ellipsoidal particles. While further validation through laboratory experiments is necessary, our research paves the road for development of new focusing and separation methods for shaped particles.

Mechanics↗
DOE OSTI2024

The importance of sub-meter-scale snow roughness on conductive heat flux of Arctic sea ice

Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.

Geology↗
DOE OSTI2024

Effect of H + Exchange and Surface Impurities on Bulk and Interfacial Electrochemistry of Garnet Solid Electrolytes

Contact loss and current constriction pose significant challenges at the Li metal interface of solid-state batteries. For garnet-structured Li 7 La 3 Zr 2 O 12 (LLZO), these effects are amplified by Li + /H + exchange and surface contamination reactions, which lead to conductivity losses and poor Li wetting. In this study, we utilize a variety of surface treatment processes across 37 cells to selectively induce proton exchange and contamination reactions in LLZO. The resulting bulk and surface chemistry is systematically characterized and correlated to changes in electrochemical properties. Additionally, we combine impedance analysis and finite element method modeling to deconvolute sources of impedance contributions at the Li metal interface. Specifically, we show that constriction impedance at the Li metal interface arises not solely from voids, but also from ionically-resistive surface contaminants. Further, these findings emphasize the connection between ionic conductivity and constriction, demonstrating that micron-scale ionically-resistive components increase constriction even with identical contact geometries. Finally, we leverage our comprehensive dataset to highlight unstable overpotential growth as a failure mechanism, additionally showing that the phase of a cell’s impedance is a sensitive indicator for the onset of interfacial instability. Overall, this study clarifies the impacts of proton exchange and surface contamination on electrochemical properties at the Li|solid electrolyte interface and elucidates insights that are generalizable to other solid-state battery systems.

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

Tuning Electrogenerated Chemiluminescence Intensity Enhancement Using Hexagonal Lattice Arrays of Gold Nanodisks

Electrogenerated chemiluminescence (ECL) microscopy shows promise as a technique for mapping chemical reactions on single nanoparticles. The technique’s spatial resolution is limited by the quantum yield of the emission and the diffusive nature of the ECL process. To improve signal intensity, ECL dyes have been coupled with plasmonic nanoparticles, which act as nanoantennas. Here, we characterize the optical properties of hexagonal arrays of gold nanodisks and how they impact the enhancement of ECL from the coreaction of tris(2,2'-bipyridyl)dichlororuthenium(II) hexahydrate and tripropylamine. Here, we find that varying the lattice spacing results in a 23-fold enhancement of ECL intensity because of increased dye-array near-field coupling as modeled using finite element method simulations.

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

Spatially resolved probing of electrochemical reactions via energy discovery platforms

The electrochemical reactivity of solid surfaces underpins functionality of a broad spectrum of materials and devices ranging from energy storage and conversion, to sensors and catalytic devices. The surface electrochemistry is, however, a complex process, controlled by the interplay of charge generation, field-controlled and diffusion-controlled transport. Here we explore the fundamental mechanisms of electrochemical reactivity on nanocrystalline ceria, using the synergy of nanofabricated devices and time-resolved Kelvin probe force microscopy (tr-KPFM), an approach we refer to as energy discovery platform. Through tr-KPFM, the surface potential mapping in both the space and time domains and current variation over time are obtained, enabling analysis of local ionic and electronic transport and their dynamic behavior on the 10 ms to 10 s scale. Based on their different responses in the time domain, conduction mechanisms can be separated and identified in a variety of environmental conditions, such as humidity and temperature. The theoretical modeling of ion transport through finite element method allows for creation of a minimal model consistent with observed phenomena, and establishing of the dynamic characteristics of the process, including mobility and diffusivity of charged species. Furthermore, the future potential of the energy discovery platforms is also discussed.

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

High-energy surface and volume plasmons in nanopatterned sub-10 nm aluminum nanostructures

In this paper, we use electron energy-loss spectroscopy to map the complete plasmonic spectrum of aluminum nanodisks with diameters ranging from 3 to 120 nm fabricated by high-resolution electron-beam lithography. Our nanopatterning approach allows us to produce localized surface plasmon resonances across a wide spectral range spanning 2–8 eV. Electromagnetic simulations using the finite element method support the existence of dipolar, quadrupolar, and hexapolar surface plasmon modes as well as centrosymmetric breathing modes depending on the location of the electron-beam excitation. In addition, we have developed an approach using nanolithography that is capable of meV control over the energy and attosecond control over the lifetime of volume plasmons in these nanodisks. The precise measurement of volume plasmon lifetime may also provide an opportunity to probe and control the DC electrical conductivity of highly confined metallic nanostructures. Lastly, we show the strong influence of the nanodisk boundary in determining both the energy and lifetime of surface plasmons and volume plasmons locally across individual aluminum nanodisks, and we have compared these observations to similar effects produced by scaling the nanodisk diameter.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Emergence of Diverse Failure Patterns in Weathering‐Induced Landslides: Insights From Particle Finite Element Simulations

Weathering is a fundamental driver of landslide evolution over geological timescales. Despite its ubiquity and importance, quantifying how weathering drives the progressive destabilization of rock slopes remains challenging. In this work, we develop a unified computational framework based on the particle finite element method to investigate the evolution of weathering‐induced landslides, from long‐term weathering to short‐term slope failure and runout dynamics. The framework integrates key processes, including weathering front propagation, time‐dependent strength degradation, rupture surface development, and post‐failure runout dynamics. Through numerical simulation experiments, we elucidate how interactions among weathering characteristics (type, intensity, and rate law), bedrock strength, fracture distribution, and slope geometry govern the failure modes and kinematics of weathering‐induced landslides. Simulations show that matrix‐dominated weathering leads to shallow translational failures, whereas fracture‐dominated weathering produces deep‐seated rotational and compound landslides. Pre‐existing fractures and slope morphology also strongly influence the movement of destabilized landmasses, affecting the failure pattern (e.g., kinematic mode and rupture surface geometry) and post‐failure behavior (e.g., runout velocity). We further demonstrate that the failure time and volume of weathered slopes are governed by the competition between gravitational driving forces and cohesive resisting forces during progressive destabilization. These findings provide new insights into the fundamental mechanisms that drive the emergence of diverse failure patterns of weathering‐induced landslides with important implications for landslide hazard assessment.

Wang, Liang [Eidgenoessische Technische Hochschule↗
DOE OSTI2021

PRISMS-Fatigue computational framework for fatigue analysis in polycrystalline metals and alloys

Abstract The PRISMS-Fatigue open-source framework for simulation-based analysis of microstructural influences on fatigue resistance for polycrystalline metals and alloys is presented here. The framework uses the crystal plasticity finite element method as its microstructure analysis tool and provides a highly efficient, scalable, flexible, and easy-to-use ICME community platform. The PRISMS-Fatigue framework is linked to different open-source software to instantiate microstructures, compute the material response, and assess fatigue indicator parameters. The performance of PRISMS-Fatigue is benchmarked against a similar framework implemented using ABAQUS. Results indicate that the multilevel parallelism scheme of PRISMS-Fatigue is more efficient and scalable than ABAQUS for large-scale fatigue simulations. The performance and flexibility of this framework is demonstrated with various examples that assess the driving force for fatigue crack formation of microstructures with different crystallographic textures, grain morphologies, and grain numbers, and under different multiaxial strain states, strain magnitudes, and boundary conditions.

Chemistry↗
DOE OSTI2024

Quantitative kinetic rules for plastic strain-induced α - ω phase transformation in Zr under high pressure

Plastic strain-induced phase transformations (PTs) and chemical reactions under high pressure are broadly spread in modern technologies, friction and wear, geophysics, and astrogeology. However, because of very heterogeneous fields of plastic strain $E$ p and stress σ tensors and volume fraction c of phases in a sample compressed in a diamond anvil cell (DAC) and impossibility of measurements of σ and $E$ p , there are no strict kinetic equations for them. Here, we develop a kinetic model, finite element method (FEM) approach, and combined FEM-experimental approaches to determine all fields in strongly plastically predeformed Zr compressed in DAC, and specific kinetic equation for α-ω PT consistent with experimental data for the entire sample. Since all fields in the sample are very heterogeneous, data are obtained for numerous complex 7D paths in the space of 3 components of the plastic strain tensor and 4 components of the stress tensor. Kinetic equation depends on accumulated plastic strain (instead of time) and pressure and is independent of plastic strain and deviatoric stress tensors, i.e., it can be applied for various above processes. Our results initiate kinetic studies of strain-induced PTs and provide efforts toward more comprehensive understanding of material behavior in extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Quantitative correlation between the void morphology of niobium-tin wires and their irreversible critical current degradation upon mechanical loading

Understanding the critical current performance variation of Nb 3 Sn superconducting wires under mechanical loading is a crucial issue for the design of next generation accelerator and fusion magnets. In these applications, the mechanical properties of the conductors may become a limiting factor due to the strong electro-magnetic forces resulting from the combination of large magnets and intense magnetic fields. In particular, the presence of voids in the superconducting filament structure, which are formed during the fabrication and the reaction heat treatment, determines localized stress concentrations and possibly the formation of cracks. In this work, we demonstrate a quantitative correlation between the void morphology and the electro-mechanical limits measured on different Bronze route Nb 3 Sn wires. Hot Isostatic Pressing (HIP) prior to the reaction heat treatment is utilized to partially eliminate the voids. The wires’ void distributions - with and without HIP treatment - are detected and statistically analyzed using high energy X-ray micro tomography. The stress concentration due to the shape and distribution of the voids as well as their impact on the electro-mechanical properties are determined through finite element method modeling. Lastly, the results are quantitatively correlated with the experimentally determined limits of the irreversible critical current degradation upon mechanical loading.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Solving sparse finite element problems on neuromorphic hardware

The finite element method (FEM) is one of the most important and ubiquitous numerical methods for solving partial differential equations (PDEs) on computers for scientific and engineering discovery. Applying the FEM to larger and more detailed scientific models has driven advances in high-performance computing for decades. Here we demonstrate that scalable spiking neuromorphic hardware can directly implement the FEM by constructing a spiking neural network that solves the large, sparse, linear systems of equations at the core of the FEM. We show that for the Poisson equation, a fundamental PDE in science and engineering, our neural circuit achieves meaningful levels of numerical accuracy and close to ideal scaling on modern, inherently parallel and energy-efficient neuromorphic hardware, specifically Intel’s Loihi 2 neuromorphic platform. We illustrate extensions to irregular mesh geometries in both two and three dimensions as well as other PDEs such as linear elasticity. Our spiking neural network is constructed from a recurrent network model of the brain’s motor cortex and, in contrast to black-box deep artificial neural network-based methods for PDEs, directly translates the well-understood and trusted mathematics of the FEM to a natively spiking neuromorphic algorithm.

Applied mathematics↗
DOE OSTI2024

Classifying topology in photonic crystal slabs with radiative environments

Abstract In the recent years, photonic Chern materials have attracted substantial interest as they feature topological edge states that are robust against disorder, promising to realize defect-agnostic integrated photonic crystal slab devices. However, the out-of-plane radiative losses in those photonic Chern slabs has been previously neglected, yielding limited accuracy for predictions of these systems’ topological protection. Here, we develop a general framework for measuring the topological protection in photonic systems, such as in photonic crystal slabs, while accounting for in-plane and out-of-plane radiative losses. Our approach relies on the spectral localizer that combines the position and Hamiltonian matrices of the system to draw a real-picture of the system’s topology. This operator-based approach to topology allows us to use an effective Hamiltonian directly derived from the full-wave Maxwell equations after discretization via finite-elements method (FEM), resulting in the full account of all the system’s physical processes. As the spectral FEM-localizer is constructed solely from FEM discretization of the system’s master equation, the proposed framework is applicable to any physical system and is compatible with commonly used FEM software. Moving forward, we anticipate the generality of the method to aid in the topological classification of a broad range of complex physical systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Topological phase transformations and intrinsic size effects in ferroelectric nanoparticles

Here, composite materials comprised of ferroelectric nanoparticles in a dielectric matrix are being actively investigated for a variety of functional properties attractive for a wide range of novel electronic and energy harvesting devices. However, the dependence of these functionalities on shapes, sizes, orientation and mutual arrangement of ferroelectric particles is currently not fully understood. In this study, we utilize a time-dependent Ginzburg-Landau approach combined with coupled-physics finite-element-method based simulations to elucidate the behavior of polarization in isolated spherical PbTiO 3 or BaTiO 3 nanoparticles embedded in a dielectric medium, including air. The equilibrium polarization topology is strongly affected by particle diameter, as well as the choice of inclusion and matrix materials, with monodomain, vortex-like and multidomain patterns emerging for various combinations of size and materials parameters. This leads to radically different polarization vs electric field responses, resulting in highly tunable size-dependent dielectric properties that should be possible to observe experimentally. Our calculations show that there is a critical particle size below which ferroelectricity vanishes. For the PbTiO 3 particle, this size is 2 and 3.4 nm, respectively, for high- and low-permittivity media. For the BaTiO 3 particle, it is ~3.6 nm regardless of the medium dielectric strength.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Effects of strain-dependent surface stress on the adhesive contact of a rigid sphere to a compliant substrate

Recent experiments have reported that the surface stress of soft elastic solids can increase rapidly with surface strain. For example, when a small hard sphere in adhesive contact with a soft silicone gel is slowly retracted from its rest position, it was found that the retraction force versus displacement relation cannot be explained either by the Johnson–Kendall–Roberts (JKR) theory or a recent indentation theory based on an isotropic surface stress that is independent of surface strain. In this paper, we address this problem using a finite element method to simulate the retraction process. Our numerical model does not have the restrictions of the aforementioned theories; that is, it can handle large nonlinear elastic deformation as well as a surface-strain-dependent surface stress. Our simulation is in good agreement with experimental force versus displacement data with no fitting parameters. Therefore, our results lend further support to the claim that significant strain-dependent surface stresses can occur in simple soft elastic gels. However, significant challenges remain in the reconciliation of theory and experiments, particularly regarding the geometry of the contact and substrate deformation.

36 MATERIALS SCIENCE↗