Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical solution”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2021

Development of data‐driven exponential integrators with application to modeling of delay photocurrents

Abstract Radiation‐induced photocurrent effects represent a threat to microelectronic components operating in space, manmade and terrestrial radiation environments. Analysis of these threats by circuit simulations requires accurate and computationally efficient compact models. Most existing compact models are based on closed form analytic solutions of the governing equations and require empirical assumptions and idealizations that can limit their validity. In this paper we formulate an alternative numerical, data‐driven approach that learns a compact model from data representative of the type of measurements one can obtain in an experimental facility. To develop the model we start from a generic discrete‐time dynamical system and then use physics knowledge to refine its structure. Numerical studies demonstrate the potential of the model and establish some empirical guidelines for its training.

Bochev, Pavel↗
DOE OSTI2022

Feedback between radiation and transport in photovoltaics

By building on the generalized Hovel model, this work develops an analytical solution to unify photon recycling, light trapping, and carrier transport in solar cells. Enhanced transport due to photon recycling diffusion is considered. While photon recycling diffusion is often negligible, it can significantly support transport in cases where charge diffusion is low but radiative efficiency is high. Next, the diffusion equation is solved to determine carrier collection for planar and textured cells. Results elucidate novel features in the external quantum efficiency (EQE) curve that occur for textured cells with high absorption and limited transport. Then, the effect of limited transport on the radiative recombination of solar cells is investigated in detail. While the traditional treatment of detailed balance uses the applied voltage to calculate a cell's radiative recombination, the present approach uses the spatially resolved quasi-Fermi-level splitting (QFLS). Further, the discrepancy between the models is shown to be inconsequential to performance for cells with low radiative recombination (Si) or with high mobilities (lowly doped GaAs). However, for cells with high radiative recombination and low mobility (highly doped GaAs, polycrystalline CdTe, and thin-film perovskites), incorporating variation in QFLS in the determination of radiative recombination can become necessary for accurate photovoltaic modeling.

14 SOLAR ENERGY↗
DOE OSTI2022

Lepto-axiogenesis in minimal $\mathrm{SUSY}$ $\mathrm{KSVZ}$ model

We study the lepto-axiogenesis scenario in the minimal supersymmetric KSVZ axion model. Only one Peccei-Quinn (PQ) field and vector-like fields are introduced besides the MSSM with the type-I see-saw mechanism. The PQ field is stabilized by the radiative correction induced by the Yukawa couplings with the vector-like fields introduced in the KSVZ model. We develop a way to follow the dynamics of the PQ field, in particular we found a semi-analytical solution which describes the rotational motion under the log-arithmic potential with including the thermalization effect via the gluon scattering which preserves the PQ symmetry. Based on the solution, we studied the baryon asymmetry, the effective number of neutrino, and the dark matter density composed of the axion and the neutralino. We found that the baryon asymmetry is successfully explained when the mass of PQ field is $\mathcal{O}$(10 6 GeV) ($\mathcal{O}$(10 5 GeV)) with the power of the PQ breaking term being 10 (8).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

First saturation correction in high energy proton-nucleus collisions. Part II. Single inclusive semi-hard gluon production

Exploiting recently obtained analytic solutions of classical Yang-Mills equations for higher order perturbations in the field of the dilute object (proton), we derive the complete first saturation correction to the single inclusive semi-hard gluon production in high energy proton-nucleus collisions by applying the Lehmann-Symanzik-Zimmermann reduction formula. We thus finalize the program started by Balitsky (see ref. [1]) and independently by Chirilli, Kovchegov and Wertepny (see ref. [2]) albeit using a very different approach to carry out our calculations. We extracted the functional dependence of gluon spectrum on the color charge densities of the colliding objects; thus our results can be used to evaluate complete first saturation correction to the double/multiple inclusive gluon productions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Anticipating gelation and vitrification with medium amplitude parallel superposition (MAPS) rheology and artificial neural networks

Abstract Anticipating qualitative changes in the rheological response of complex fluids (e.g., a gelation or vitrification transition) is an important capability for processing operations that utilize such materials in real-world environments. One class of complex fluids that exhibits distinct rheological states are soft glassy materials such as colloidal gels and clay dispersions, which can be well characterized by the soft glassy rheology (SGR) model. We first solve the model equations for the time-dependent, weakly nonlinear response of the SGR model. With this analytical solution, we show that the weak nonlinearities measured via medium amplitude parallel superposition (MAPS) rheology can be used to anticipate the rheological aging transitions in the linear response of soft glassy materials. This is a rheological version of a technique called structural health monitoring used widely in civil and aerospace engineering. We design and train artificial neural networks (ANNs) that are capable of quickly inferring the parameters of the SGR model from the results of sequential MAPS experiments. The combination of these data-rich experiments and machine learning tools to provide a surrogate for computationally expensive viscoelastic constitutive equations allows for rapid experimental characterization of the rheological state of soft glassy materials. We apply this technique to an aging dispersion of Laponite ® clay particles approaching the gel point and demonstrate that a trained ANN can provide real-time detection of transitions in the nonlinear response well in advance of incipient changes in the linear viscoelastic response of the system.

Lennon, Kyle R. (ORCID:0000000212515461)↗
DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗
DOE OSTI2023

A helical actuator driven by biased SMA: design, model, and experiment

A helical actuator driven by biased shape memory alloy (SMA) patterns embedded into a soft composite ribbon base is presented in this work. Instead of common U-shape SMA wires, a single SMA wire is woven into planar patterns, which enable helical deformation of the composite ribbon from an initially flat geometry. An analytical static model is established for accurate and rapid prediction of the helical reconfiguration arising from the shape memory effect of woven SMA patterns, followed by validation of the static model using the finite element method (FEM). The finite element results are compared with the analytical solutions given by this static model, which show a high agreement. Parametric study of the influences of eight independent design variables on the dependent helical parameters, such as combined curvatures, pitches, and helical angles, is completed. It is found that the helically deformed geometry is mainly dominated by diameters, biased positions, inclined angles, and numbers of skewed segments of the SMA wire. Fabrication and in-situ experimental test of a prototype of such helical actuators qualitatively demonstrate its dramatic three-dimensional (3D) spiral reconfiguration from a two-dimensional (2D) flat ribbon. In conclusion, such SMA patterns will allow more diverse designs of soft actuators for a wider range of robotic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Package of Momentum and Heat Transfer Coefficients for the Stable Surface Layer Extended by New Coefficients over Sea Ice

The bulk transfer coefficients of momentum, heat, and humidity belong to the main ingredients of numerical weather prediction and climate models. They are needed for the calculation of turbulent fluxes in the surface layer and often rely on the Monin–Obukhov similarity theory requiring universal stability functions. The problem of a derivation of transfer coefficients based on different stability functions has been considered by many researchers over the years but it remains to this day. In this work, dedicated to the memory of S.S. Zilitinkevich, we also address this task, and obtain transfer coefficients from three pairs of theoretically derived stability functions suggested by Zilitinkevich and co-authors for stable conditions. Additionally, we construct non-iterative parametrizations of these transfer coefficients based on earlier work. Results are compared with state-of-the-art coefficients for land, ocean, and sea ice. The combined parametrizations form a package in a universal framework relying on a semi-analytical solution of the Monin-Obukhov similarity theory equations. A comparison with data of the Surface Heat Budget of the Arctic Ocean campaign (SHEBA) over sea ice reveals large differences between the coefficients for land conditions and the measurements over sea ice. However, two schemes of Zilitinkevich and co-authors show, after slight modification, good agreement with SHEBA although they had not been especially developed for sea ice. One pair of the modified transfer coefficients is superior and is compatible to earlier SHEBA-based parametrizations. Finally, an algorithm for practical use of all transfer coefficients in climate models is given.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

TOUGH3-FLAC3D: a modeling approach for parallel computing of fluid flow and geomechanics

The recent development of the TOUGH3 code allows for a faster and more reliable fluid flow simulator. At the same time, new versions of FLAC3D are released periodically, allowing for new features and faster execution. In this paper, we present the first implementation of the coupling between TOUGH3 and FLAC3Dv6/7, maintaining parallel computing capabilities for the coupled fluid flow and geomechanical codes. We compare the newly developed version with analytical solutions and with the previous approach, and provide some performance analysis on different meshes and varying the number of running processors. Finally, we present two case studies related to fault reactivation during CO 2 sequestration and nuclear waste disposal. The use of parallel computing allows for meshes with a larger number of elements, and hence more detailed understanding of thermo-hydro-mechanical processes occurring at depth.

58 GEOSCIENCES↗
DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗
DOE OSTI2020

An Experimental Method to Measure Gaseous Diffusivity in Tight and Partially Saturated Porous Media via Continuously Monitored Mass Spectrometry

Detection of radioxenon and radioargon produced by underground nuclear explosions is one of the primary methods by which the Comprehensive Nuclear-Test–Ban Treaty (CTBT) monitors for nuclear activities. However, transport of these noble gases to the surface via barometric pumping is a complex process relying on advective and diffusive processes in a fractured porous medium to bring detectable levels to the surface. To better understand this process, experimental measurements of noble gas and chemical surrogate diffusivity in relevant lithologies are necessary. However, measurement of noble gas diffusivity in tight or partially saturated porous media is challenging due to the transparent nature of noble gases, the lengthy diffusion times, and difficulty maintaining consistent water saturation. Here, the quasi-steady-state Ney–Armistead method is modified to accommodate continuous gas sampling via effusive flow to a mass spectrometer. An analytical solution accounting for the cumulative sampling losses and induced advective flow is then derived. Experimental results appear in good agreement with the proposed theory, suggesting the presence of retained groundwater reduces the effective diffusivity of the gas tracers by 10–1000 times. Finally, by using a mass spectrometer, the method described herein is applicable to a broad range of gas species and porous media.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Direct Fabrication of Vertically Stacked Double Barrier Tunnel Junctions Based on Graphene and h-BN

Direct manufacturing of two-dimensional material-based double barrier (DB) tunnel junctions, based on a lithography-free approach was developed. Graphene/h-BN/Graphene/h-BN/Graphene devices were deposited on Si/SiO2 substrates by employing a plasma enhanced chemical vapor deposition technique in a sequential manner. DB tunneling junctions with varying barrier widths (by varying the thickness of the second graphene layer) were studied. Samples were characterized using Raman, Atomic Force Microscopy and X-ray photoemission spectroscopy. The I–V characteristics of tunneling current showed resonant tunneling behavior at room temperature with a negative differential conductance. The behavior could be explained with quantum mechanical double barrier tunneling model in which analytic solutions to Schrödinger’s equation were obtained in each region of the system. Resonances in transmission probability coefficient for varying barrier widths were evaluated and compared with the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Lax-Oleinik-Type Formulas and Efficient Algorithms for Certain High-Dimensional Optimal Control Problems

Two of the main challenges in optimal control are solving problems with state-dependent running costs and developing efficient numerical solvers that are computationally tractable in high dimension. In this paper, we provide analytical solutions to certain optimal control problems whose running cost depends on the state variable and with constraints on the control. We also provide Lax-Oleinik-type representation formulas for the corresponding Hamilton-Jacobi partial differential equations with state-dependent Hamiltonians. Additionally, we present an efficient, grid-free numerical solver based on our representation formulas, which is shown to scale linearly with the state dimension, and thus, to overcome the curse of dimensionality. Using existing optimization methods and the min-plus technique, we extend our numerical solvers to address more general classes of convex and nonconvex initial costs. We demonstrate the capabilities of our numerical solvers using implementations on a central processing unit (CPU) and a field-programmable gate array (FPGA). In several cases, our FPGA implementation obtains over a 10 times speedup compared to the CPU, which demonstrates the promising performance boosts FPGAs can achieve. Furthermore, our numerical results show that our solvers have the potential to serve as a building block for solving broader classes of high-dimensional optimal control problems in real-time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Modeling the influence of the solid electrolyte interphase on the sand’s time and dendrite formation on lithium metal electrodes

Lithium metal is a sought after battery material for its high energy density due to the low electrochemical potential and density. However, lithium metal is also highly reactive, which results in a strong propensity for dendrite formation. The Sand’s time has previously been used to predict the time of dendrite initiation on metals that do not form a solid-electrolyte interphase (SEI), but it has been shown that the Sand’s time is not accurate for lithium electrodes when using transport parameters associated with the electrolyte. Thus, we built a numerical model to simulate lithium ion transport through a growing SEI to predict the Sand’s time. The numerical model is shown to be more accurate than previous analytical solutions, especially for low current densities. We then analyze the sensitivity of the Sand’s time to different SEI properties and the chemical potential gradients present in the SEI, driving lithium transport. The results showed that high lithium concentration has a greater impact at high current density, while fast diffusivity is more important at low current density. Lastly, we modeled the influence of surface roughness on the plating evolution and chemical potential gradients when an SEI is present in comparison to the electrolyte. As a result, we demonstrate that the SEI plays a critical role in lithium electrode stability, and that improved characterization techniques are needed to better understand transport through the SEI and increase lithium metal utilization in energy storage devices.

Chemistry↗
DOE OSTI2020

Slip band characteristics in the presence of grain boundaries in nickel-based superalloy

Shear strain profiles along slip bands in a modified Rolls-Royce nickel superalloy (RR1000) were analyzed for a tensile sample deformed by 2%. The strain increased with distance away from a grain boundary (GB), with maximum shear strain towards the center of the grain, indicating that dislocation nucleation generally occurred in the grain interior. The strain gradients in the neighborhood of the GBs were quantified and generally correlated with rotation about the active slip system line direction. This leads to an ability to determine the active slip system in these regions. The dislocation spacing and pileup stresses were inferred. The dislocation spacing closely follows an Eshelby analytical solution for a single ended pileup of dislocations under an applied stress. The distribution of pileup stress values for GBs of a given misorientation angle follows a log-normal distribution, with no correlation between the pileup stress and the GB misorientation angle. Furthermore, there is no observed correlation between various transmissivity factors and slip band pileup stress. Hence it appears that the obstacle strength of any of the observed GBs is adequate to facilitate the dislocation pileups present in the slip bands. However, slip band transmission does correlate with transmissivity factors, with the current study focusing on the Luster and Morris m’-factor. Observation of strain profiles of transmitted bands indicate dislocation nucleation locations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

GW-PINN: A deep learning algorithm for solving groundwater flow equations

Machine learning methods provide new perspective for more convenient and efficient prediction of groundwater flow. In this study, a deep learning method “GW-PINN” without labeled data for solving groundwater flow equations with wells was proposed. GW-PINN takes the physics inform neural network (PINN) as the backbone and uses either the hard or soft constraint in the loss function for training. A locally refined sampling strategy (LRS) is adopted to generate the consistent spatial sampling points for problems with strong hydraulic head change, and then combined with an appropriate temporal sampling scheme to obtain the final spatial-temporal sampling points. A snowball-style two-stage training strategy by dividing the temporal domain into two subdomains is designed to decrease the sampling points. Five cases were designed to test the training performance of GW-PINN under different sampling strategies and two constraints. The predicted results of GW-PINN were compared with MODFLOW and the analytical solution. The results demonstrate that GW-PINN possesses strong ability in capturing the hydraulic head change for both confined and un-confined aquifers. The hard constraint owns more robust learning ability than the soft constraint. The LRS strategy can generate more accurate results with much fewer sampling points than traditional sampling strategies, and the snowball-style two-stage training strategy is significantly efficient for problems with the drastic change of hydraulic head. Additionally, the application of GW-PINN as a surrogate model for parameterized groundwater flow equations is illustrated. This study provides an option tool for efficient groundwater flow simulation, especially for those with local refinements are needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A compact mechanical power take-off for wave energy converters: Design, analysis, and test verification

Power Take-off (PTO) is one of the most critical components in the ocean wave conversion. This paper introduces an innovative compact mechanical PTO using a ball screw mechanism and a Mechanical Motion Rectifier (MMR) to improve the energy conversion performance. The proposed MMR PTO converts bidirectional oscillatory motion of the ocean waves into unidirectional rotation. The detailed design and working principle of the MMR PTO are explicitly explained. Dynamic analysis is presented for the three most important factors that influence the performance, including the equivalent mass, the equivalent damping, and the excitation frequency. The analytical relationship between these factors and the disengaged freewheeling motion caused by the MMR mechanism is derived, and the analytical solution is verified through dry lab testing. Experimental results in bench tests show that compared with the traditional linear PTO which use the ball screw to drive the generator directly, the MMR PTO is more efficient because of the unique freewheeling motion, achieving a maximum 81.2% energy transfer efficiency. After integrating the prototype PTO with a single-body point absorber type of wave energy converter, testing in a water tank yields total PTO efficiency up to 62.4%, showing a very promising potential for real applications.

16 TIDAL AND WAVE POWER↗