Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

srMO-BO-3GP: A sequential regularized multi-objective Bayesian optimization for constrained design applications using an uncertain Pareto classifier

Bayesian optimization (BO) is an efficient and flexible global optimization framework that is applicable to a very wide range of engineering applications. To leverage the capability of the classical BO, many extensions, including multi-objective, multi-fidelity, parallelization, and latent-variable modeling, have been proposed to address the limitations of the classical BO framework. In this work, we propose a novel multi-objective BO formalism, called srMO-BO-3GP, to solve multi-objective optimization problems in a sequential setting. Three different Gaussian processes (GPs) are stacked together, where each of the GPs is assigned with a different task. The first GP is used to approximate a single-objective computed from the multi-objective definition, the second GP is used to learn the unknown constraints, and the third one is used to learn the uncertain Pareto frontier. At each iteration, a multi-objective augmented Tchebycheff function is adopted to convert multi-objective to single-objective, where the regularization with a regularized ridge term is also introduced to smooth the single-objective function. Finally, we couple the third GP along with the classical BO framework to explore the convergence and diversity of the Pareto frontier by the acquisition function for exploitation and exploration. The proposed framework is demonstrated using several numerical benchmark functions, as well as a thermomechanical finite element model for flip-chip package design optimization.

42 ENGINEERING↗
DOE OSTI2022

Continuous conditional generative adversarial networks for data-driven solutions of poroelasticity with heterogeneous material properties

Machine learning-based data-driven modeling can allow computationally efficient time-dependent solutions of PDEs, such as those that describe subsurface multiphysical problems. In this work, our previous approach (Kadeethum et al., 2021d) of conditional generative adversarial networks (cGAN) developed for the solution of steady-state problems involving highly heterogeneous material properties is extended to time-dependent problems by adopting the concept of continuous cGAN (CcGAN). The CcGAN that can condition continuous variables is developed to incorporate the time domain through either element-wise addition or conditional batch normalization. Moreover, this framework can handle training data that contain different timestamps and then predict timestamps that do not exist in the training data. As a numerical example, the transient response of the coupled poroelastic process is studied in two different permeability fields: Zinn & Harvey transformation and a bimodal transformation. The proposed CcGAN uses heterogeneous permeability fields as input parameters while pressure and displacement fields over time are model output. Our results show that the model provides sufficient accuracy with computational speed-up. This robust framework will enable us to perform real-time reservoir management and robust uncertainty quantification in poroelastic problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Tensile behavior of diffusion bonded AA6061 - AA6061 with variation in cooling method

Hot isostatically pressed AA6061 cladding is an important structural component of the high performance, Zr-laminated U-10Mo monolithic fuel system for the application in research and test reactors. In this study, the mechanical behavior of two diffusion bonded aluminum alloy, AA6061, was examined using tensile testing. Solid-to-solid diffusion bonding between two pieces of AA6061 was performed by isothermal annealing at 560 °C for 1.5 h, and diffusion couples were subsequently cooled via three different cooling methods: furnace cooling, air cooling, and water quenching. Dog-bone shaped tensile specimens, with 10 mm in gauge length (with diffusion bonded interface in the middle), and 1.5 × 1.5 mm 2 gauge cross-sections, were fabricated from the diffusion bonded AA6061 by electro-discharge machining. Yield strength (% EL at failure) of furnace cooled, air cooled and water quenched tensile specimens determined was 82 ~ 89 MPa (10 ~ 30%), 112 ~ 116 MPa (10 ~ 14%), and 149 ~ 164 MPa (10 ~ 17%), respectively. This variation in mechanical behavior was examined with cooling-rate dependent, concentrated precipitation of Mg 2 Si at the diffusion bonded interface, with due respect for mechanical properties of the AA6061 alloy that inherently vary as a function of cooling rate from 560 °C. Finite element analysis using ABAQUS was employed to augment experimental findings with the appropriate microstructural constituents and alloy properties. Results suggest that the strength is dominated by matrix/bulk properties of AA6061, while ductility is strongly influenced by the cooling method dependent presence of Mg 2 Si precipitates at the interface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Response of salt intrusion in a tidal estuary to regional climatic forcing

Abstract Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Global optimization of harmonic oscillator basis in covariant density functional theory

The present investigation focuses on the improvement of the accuracy of the description of binding energies within moderately sized fermionic basis. Using the solutions corresponding to infinite fermionic basis it was shown that in the case of meson exchange (ME) covariant energy density functionals (CEDFs) the global accuracy of the description of binding energies in the finite $N_F$ = 16 - 20 bases can be drastically (by a factor ranging from ~3 up to ~9 dependent on the functional and $N_F$) improved by a global optimization of oscillator frequency of the basis. This is a consequence of the unique feature of the ME functionals in which with increasing fermionic basis size fermionic and mesonic energies approach the exact (infinite basis) solution from above and below, respectively. As a consequence, an optimal oscillator frequency $\hbar\omega_0$ of the basis can be defined which provides an accurate reproduction of exact total binding energies by the ones calculated in truncated basis. This leads to a very high accuracy of the calculations in moderately sized $N_F=20$ basis when mass dependent oscillator frequency is used: global rms differences $\delta B_{rms}$ between the binding energies calculated in infinite and truncated bases are only 0.025 MeV and 0.031 MeV for the NL5(Z) and DD-MEZ functionals, respectively. Optimized values of the oscillator frequency $\hbar\omega_0$ are provided for three major classes of CEDFs, i.e. for density dependent meson exchange functionals, nonlinear meson exchange ones and point coupling functionals.

Binding energy & masses↗
DOE OSTI2024

Stripe and bubble ratchets on asymmetric substrates

We show that various nonmonotonic ratchet effects can arise when mesophase pattern-forming systems, which exhibit anisotropic crystal, stripe, and bubble regimes, are coupled to one-dimensional asymmetric substrates under ac driving. The patterns emerge in the absence of a substrate when the ratio of attraction to repulsion is varied for particles with competing short-range attraction and long-range repulsion potentials. In the presence of the substrate, we find that the ratchet efficiency varies nonmonotonically with increasing attraction, depending upon how well the mesophase morphology matches the substrate spacing and periodicity. For strong repulsion, there is a weak but finite ratchet effect, while at intermediate attraction, there is a robust ratchet effect in which the system forms stripes aligned with the substrate symmetry direction. For strong attraction, large bubbles appear that have weak or no ratchet effects when the bubble width exceeds the substrate lattice spacing, causing the bubble to be only weakly coupled to the substrate. For very strong attraction, small bubbles form and undergo a strong ratchet effect with an efficiency that oscillates as a function of ac drive amplitude. Although the small bubble regime is a strongly correlated regime in which each bubble contains many particles, the system behaves as if it is in a single-particle regime since the bubble width is much smaller than the substrate lattice spacing. We map out the different rectification phases as a function of the pattern morphology, substrate strength, and ac drive amplitude. The pronounced ratchet effects that we observe in some regimes can be exploited for pattern sorting in hard and soft matter systems. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Analysis of the performance of different culvert boundary conditions in 2D shallow flow models

Culverts allow roads to safely traverse small streams or drainage ditches, and their proper design is critical to ensure a safe and reliable transportation network. A correct modelization of these hydraulic structures becomes crucial in the assessment of flood footprints or discharge peak estimation in a risk evaluation plan. The question of how to include culverts comes up frequently when assembling a hydraulic model that requires the presence of as many singular elements as possible. In this work, three different culvert integrations with the surface domain are studied and compared in the context of a 2D shallow water (SW) model. All of them are based on the Federal Highway Administration (FHWA) formulation for the culvert discharge estimation but differ in complexity and in the interaction with the numerical model for surface flow, some of them as internal boundary conditions. Several steady and unsteady validation test cases are presented and the numerical results are compared with the predictions from HEC-RAS 1D and HY-8 software. Furthermore, the culvert area, shape and their sensitivity to the 2D computational mesh is also analyzed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Response Of Salt Intrusion In A Tidal Estuary To Regional Climatic Forcing (FVCOM Output for the baseline scenario)

Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels. These findings have been summarized into a journal manuscript. The data included in this information release request are FVCOM model output associated with the numerical simulations described in the journal manuscript.

Liu, Jinliang↗
DOE OSTI2020

Peak Ground Velocity Spatial Variability Revealed by Dense Seismic Array in Southern California

There are an infinite number of locations where seismic waves interact with Earth's surface, yet seismologists can only record a finite set of ground motions. This requires assumptions about how known measurements or predictions can be extrapolated to nearby locations. We test to what extent extrapolation is viable using an array of 1,108 geophones installed 10–20 m apart within a 0.36 km 2 study area for 1 month atop the San Jacinto Fault in Southern California. Measuring ground motions from 38 small local earthquakes, we observe peak ground velocity (PGV) variations up to 22% across the array. We find the basin structure consistently traps energy, amplifying the seismic velocities, and that the damaged region within the fault reduces the PGVs. Comparing surface and borehole measurements (148 m deep), we find PGVs recorded at the surface are 3–10 times greater than the measurements at depth. We attribute these differences to amplification from soft sediments within the basin. We show that by leveraging information from >1,000 measurements per earthquake the ground motion velocities can be tightly constrained. These results suggest the spatial extrapolation of PGVs in regions of known faults and basins, even across a small area, should be done with caution.

58 GEOSCIENCES↗
DOE OSTI2020

Flaw Tolerance of Octet Truss Lattices with Random Flaws

We investigate the behavior of octet truss lattices in which some proportion of the bars are randomly removed from a perfect lattice. Using a lattice specimen of 9 × 9 × 9 octet truss unit cells as the nominal design, we simulate the response of prefect and flawed lattices with a finite element method implementation of rod theory. For each of two simple load cases, we perform a Monte Carlo analysis where we consider many instantiations of flawed lattices where between 1% and 20% of the bars are missing, and compare the axial loads in the remaining bars to those in the perfect lattice. We also compare the overall stiffness of the flawed lattices to the perfect lattice. Our results show progressive deterioration of lattice strength and stiffness with increasing proportion of missing bars, as expected, with the location and magnitude of the highest-loaded (most failure-prone) bars depending more on the load case and bar orientation than on the location and connectivity of missing bars. As such, the results reported herein may not be directly transferable to different lattice geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Biological mechanisms may contribute to soil carbon saturation patterns

Increasing soil organic carbon (SOC) storage is a key strategy to mitigate rising atmospheric CO 2 , yet SOC pools often appear to saturate, or increase at a declining rate, as carbon (C) inputs increase. Soil C saturation is commonly hypothesized to result from the finite amount of reactive mineral surface area available for retaining SOC, and is accordingly represented in SOC models as a physicochemically determined SOC upper limit. However, mineral-associated SOC is largely microbially generated. In this perspective, we present the hypothesis that apparent SOC saturation patterns could emerge as a result of ecological constraints on microbial biomass—for example, via competition or predation—leading to reduced C flow through microbes and a reduced rate of mineral-associated SOC formation as soil C inputs increase. Microbially explicit SOC models offer an opportunity to explore this hypothesis, yet most of these models predict linear microbial biomass increases with C inputs and insensitivity of SOC to input rates. Synthesis of 54 C addition studies revealed constraints on microbial biomass as C inputs increase. Different hypotheses limiting microbial density were embedded in a three-pool SOC model without explicit limits on mineral surface area. As inputs increased, the model demonstrated either no change, linear, or apparently saturating increases in mineral-associated and particulate SOC pools. Taken together, our results suggest that microbial constraints are common and could lead to reduced mineral-associated SOC formation as input rates increase. We conclude that SOC responses to altered C inputs—or any environmental change—are influenced by the ecological factors that limit microbial populations, allowing for a wider range of potential SOC responses to stimuli. Understanding how biotic versus abiotic factors contribute to these patterns will better enable us to predict and manage soil C dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Benders Cut Classification via Support Vector Machines for Solving Two-Stage Stochastic Programs

In this work, we consider Benders decomposition for solving two-stage stochastic programs with complete recourse based on finite samples of the uncertain parameters. We define the Benders cuts binding at the final optimal solution or the ones significantly improving bounds over iterations as valuable cuts. We propose a learning-enhanced Benders decomposition (LearnBD) algorithm, which adds a cut classification step in each iteration to selectively generate cuts that are more likely to be valuable cuts. The LearnBD algorithm includes two phases: (i) sampling cuts and collecting information from training problems and (ii) solving testing problems with a support vector machine (SVM) cut classifier. We run the LearnBD algorithm on instances of capacitated facility location and multicommodity network design under uncertain demand. Our results show that SVM cut classifier works effectively for identifying valuable cuts, and the LearnBD algorithm reduces the total solving time of all instances for different problems with various sizes and complexities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

[Presentation Slides] Code-to-Code Benchmark Study for Thermal Stress Modeling and Preliminary Analysis of the High-temperature Single Heat-Pipe Experiment

In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios. This presentation slides discuss finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL), using two commercial software packages, Abaqus and Star-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain preliminary insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the inner surface of the heat pipe hole surrounded by electric heaters where the largest temperature gradients appear. The temperature fields have good agreements between Abaqus and Star-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes. It is found that the local maximum thermal stresses may reach close to the ultimate tensile strength and yield strength of structural material depending on the heater power. Ultimately, the coupled thermal-structural analysis will help guide the current experimental plan and ensure the facility safety for future experimental study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗
DOE OSTI2022

Symmetry Breaking in an Extended-O(2) Model

Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we consider an extended-O(2) model that differs from the ordinary O(2) model by an explicit symmetry breaking term. Its coupling allows to smoothly interpolate between the O(2) model (zero coupling) and a $q$-state clock model (infinite coupling). In the latter case, a $q$-state clock model can also be defined for non-integer values of $q$. Thus, such a limit can also be considered as an analytic continuation of an ordinary $q$-state clock model to non-integer $q$. The phase diagram of the extended-O(2) model in the infinite coupling limit was established in our previous work, where it was shown that for non-integer $q$, there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we investigate the model at finite values of the coupling using Monte Carlo and tensor methods. The results may be relevant for configurable Rydberg-atom arrays.

Hostetler, Leon↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Differential analysis of incompressibility in neutron-rich nuclei

Both the incompressibility K A of a finite nucleus of mass A and that (K ∞ ) of infinite nuclear matter are fundamentally important for many critical issues in nuclear physics and astrophysics. While some consensus has been reached about K ∞ , accurate theoretical predictions and experimental extractions of K τ characterizing the isospin dependence of K A have been very difficult. We propose a differential approach to extract K τ and K ∞ independently from the K A data of any two nuclei in a given isotope chain. Applying this method to the K A data from isoscalar giant monopole resonances (ISGMR) in even-even Pb, Sn, Cd, and Ca isotopes taken by Garg et al. at the Research Center for Nuclear Physics (RCNP), Osaka University, Japan, we find that the 106 Cd– 116 Cd and 112 Sn– 124 Sn pairs having the largest differences in isospin asymmetries in their respective isotope chains measured so far provide consistently the most accurate up-to-date K τ value of K τ = –616 ± 59 MeV and K τ =–623 ± 86 MeV, respectively, largely independent of the remaining uncertainties of the surface and Coulomb terms in expanding K A , while the K ∞ values extracted from different isotopes chains are all well within the current uncertainty range of the community consensus for K ∞ . Moreover, the size and origin of the “soft Sn puzzle” is studied with respect to the “stiff Pb phenomenon.” Furthermore, it is found that the latter is favored due to a much larger (by ≈ 380 MeV) K τ for Pb isotopes than for Sn isotopes, while K ∞ from analyzing the K A data of Sn isotopes is only about 5 MeV less than that from analyzing the Pb data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗