Results for “volume constraint”

DOE OSTI2023

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Analysis of Rig Parameter Data Using Drilling Process Modeling Constraints, Volume 5: Utah FORGE Well 16B(78)-32

Drill rig parameter measurements are routinely used during deep well construction to monitor and guide drilling conditions for improved performance and reduced costs. While insightful into the drilling process, these measurements are of reduced value without a standard to aid in data evaluation and decision making. In the main body of this work (Volume 1), a method is demonstrated whereby rock reduction model constraints are used to interpret drilling response parameters; the method could be applied in real-time to improve decision-making in the field and to further discern technology performance during post-drilling evaluations. Drilling parameters are evaluated using laboratory-validated rock reduction models for predicting the phenomenological response of drag bits (Detournay and Defourny, 1992) in computational algorithms. The method presented has applicability to development of advanced analytics on future geothermal wells using real-time electronic data recording for improved performance and reduced drilling costs. A drilling cost model is also used to show the tradeoff between rate of penetration and bit life and the influence on interval drilling costs. Details of the bit specifications and performance are cataloged in an independent volume, documented under separate cover, for each of the four wells, and include Volume 2: Utah FORGE 16A(78)-32; Volume 3: Utah FORGE 56-32; Volume 4: Utah FORGE 78B-32 and Volume 5: Utah FORGE 16B(78)-32.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Analysis of Rig Parameter Data Using Drilling Process Modeling Constraints, Volume 4: Utah FORGE Well 78B-32

Drill rig parameter measurements are routinely used during deep well construction to monitor and guide drilling conditions for improved performance and reduced costs. While insightful into the drilling process, these measurements are of reduced value without a standard to aid in data evaluation and decision making. In the main body of this work (Volume 1), a method is demonstrated whereby rock reduction model constraints are used to interpret drilling response parameters; the method could be applied in real-time to improve decision-making in the field and to further discern technology performance during post-drilling evaluations. Drilling parameters are evaluated using laboratory-validated rock reduction models for predicting the phenomenological response of drag bits (Detournay and Defourny, 1992) in computational algorithms. The method presented has applicability to development of advanced analytics on future geothermal wells using real-time electronic data recording for improved performance and reduced drilling costs. A drilling cost model is also used to show the tradeoff between rate of penetration and bit life and the influence on interval drilling costs. Details of the bit specifications and performance are cataloged in an independent volume, documented under separate cover, for each of the four wells, and include Volume 2: Utah FORGE 16A(78)-32; Volume 3: Utah FORGE 56-32; Volume 4: Utah FORGE 78B-32 and Volume 5: Utah FORGE 16B(78)-32.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Analysis of Rig Parameter Data Using Drilling Process Modeling Constraints, Volume 3: Utah FORGE Well 56-32

Drill rig parameter measurements are routinely used during deep well construction to monitor and guide drilling conditions for improved performance and reduced costs. While insightful into the drilling process, these measurements are of reduced value without a standard to aid in data evaluation and decision making. In the main body of this work (Volume 1), a method is demonstrated whereby rock reduction model constraints are used to interpret drilling response parameters; the method could be applied in real-time to improve decision-making in the field and to further discern technology performance during post-drilling evaluations. Drilling parameters are evaluated using laboratory-validated rock reduction models for predicting the phenomenological response of drag bits (Detournay and Defourny, 1992) in computational algorithms. The method presented has applicability to development of advanced analytics on future geothermal wells using real-time electronic data recording for improved performance and reduced drilling costs. A drilling cost model is also used to show the tradeoff between rate of penetration and bit life and the influence on interval drilling costs. Details of the bit specifications and performance are cataloged in an independent volume, documented under separate cover, for each of the four wells, and include Volume 2: Utah FORGE 16A(78)-32; Volume 3: Utah FORGE 56-32; Volume 4: Utah FORGE 78B-32 and Volume 5: Utah FORGE 16B(78)-32.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Analysis of Rig Parameter Data Using Drilling Process Modeling Constraints, Volume 1: Summary of Utah FORGE Wells 16A(78)-32, 56-32, 78B-32 and 16B(78)-32

Drill rig parameter measurements are routinely used during deep well construction to monitor and guide drilling conditions for improved performance and reduced costs. While insightful into the drilling process, these measurements are of reduced value without a standard to aid in data evaluation and decision making. In the main body of this work (Volume 1), a method is demonstrated whereby rock reduction model constraints are used to interpret drilling response parameters; the method could be applied in real-time to improve decision-making in the field and to further discern technology performance during post-drilling evaluations. Drilling parameters are evaluated using laboratory-validated rock reduction models for predicting the phenomenological response of drag bits (Detournay and Defourny, 1992) in computational algorithms. The method presented has applicability to development of advanced analytics on future geothermal wells using real-time electronic data recording for improved performance and reduced drilling costs. A drilling cost model is also used to show the tradeoff between rate of penetration and bit life and the influence on interval drilling costs. Details of the bit specifications and performance are cataloged in an independent volume, documented under separate cover, for each of the four wells, and include Volume 2: Utah FORGE 16A(78)-32; Volume 3: Utah FORGE 56-32; Volume 4: Utah FORGE 78B-32 and Volume 5: Utah FORGE 16B(78)-32.

15 GEOTHERMAL ENERGY↗

DOE OSTI2018

Multi-material continuum topology optimization with arbitrary volume and mass constraints

Not Available

Sanders, Emily D.↗

DOE OSTI2018

PolyMat: an efficient Matlab code for multi-material topology optimization

We introduce a Matlab implementation of topology optimization for compliance minimization on unstructured polygonal finite element meshes that efficiently accommodates many materials and many volume constraints. Leveraging the modular structure of the educational code, PolyTop, we extend it to the multi-material version, PolyMat, with only a few modifications. First, a design variable for each candidate material is defined in each finite element. Next, we couple a Discrete Material Optimization interpolation with the existing penalization and introduce a new parameter such that we can employ continuation and smoothly transition from a convex problem without any penalization to a non-convex problem in which material mixing and intermediate densities are penalized. Mixing that remains due to the density filter operation is eliminated via continuation on the filter radius. To accommodate flexibility in the volume constraint definition, the constraint function is modified to compute multiple volume constraints and the design variable update is modified in accordance with the Zhang-Paulino-Ramos Jr. (ZPR) update scheme, which updates the design variables associated with each constraint independently. The formulation allows for volume constraints controlling any subset of the design variables, i.e., they can be defined globally or locally for any subset of the candidate materials. Borrowing ideas for mesh generation on complex domains from PolyMesher, we conclude which design variables are associated with each local constraint of arbitrary geometry. A number of examples are presented to demonstrate the many material capability, the flexibility of the volume constraint definition, the ease with which we can accommodate passive regions, and how we may use local constraints to break symmetries or achieve graded geometries.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Whole-brain Irradiation Field Design: A Comparison of Parotid Dose

Whole-brain radiation therapy (WBRT) plays an important role in patients with diffusely metastatic intracranial disease. Whether the extent of the radiation field design to C1 or C2 affects parotid dose and risk for developing xerostomia is unknown. The goal of this study is to examine the parotid dose based off of the inferior extent of WBRT field to either C1 or C2. Patients treated with WBRT with either 30 Gy or 37.5 Gy from 2011 to 2014 at a single institution were examined. Parotid dose constraints were compared with Radiation Therapy Oncology Group (RTOG) 0615 nasopharyngeal carcinoma for a 33-fraction treatment: mean <26 Gy, volume constraint at 20 Gy (V20) < 20 cc, and dose at 50% of the parotid volume (D50) < 30 Gy. Biologically effective dose (BED) conversions with an α/β of 3 for normal parotid were performed to compare with 10-fraction and 15-fraction treatments of WBRT. The constraints are as follows: mean < BED 32.83 Gy, V15.76 (for 10-fraction WBRT) or V17.35 (for 15-fraction WBRT) < 20 cc, and D50 < BED 39.09 Gy. Nineteen patients treated to C1 and 26 patients treated to C2 were analyzed. Comparing WBRT to C1 with WBRT to C2, the mean left, right, and both parotids' doses were lower when treated to C1. Converting mean dose to BED{sub 3}, the parotid doses were lower than BED{sub 3} constraint of 32.83 Gy: left (30.12 Gy), right (30.69 Gy), and both parotids (30.32 Gy). V20 to combined parotids was lower in patients treated to C1. When accounting for fractionation of WBRT received, the mean corrected V20 volume was less than 20 cc when treating to C1. D50 for C1 was lower than C2 for the left parotid, right parotid, and both parotids. BED{sub 3} conversion for the mean D50 of the left, right, and both parotids was less than 39.09 Gy. In conclusion, WBRT to C1 limits parotid dose, and parotid dose constraints are achievable compared with inferior border at C2. A possible mean parotid dose constraint with BED{sub 3} should be less than 32.83 Gy.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2018

Comparison of calculated dose distributions reported as dose-to-water and dose-to-medium for intensity-modulated radiotherapy of nasopharyngeal cancer patients

Advanced dose calculation algorithms for radiation therapy treatment planning can report external beam photon dose 2-sided, in terms of dose-to-medium (D{sub m}) and dose-to-water (D{sub w}). The purpose of our study was to determinate the effect of D{sub w} and D{sub m} reporting modes built in Elekta Monaco treatment planning system on intensity-modulated radiotherapy dose distributions for patients with nasopharyngeal cancer. For 13 patients involved in this retrospective study, 2 plans were created: 1 using D{sub w} and another according to D{sub m} reporting mode. Treatment plans were normalized such that 100% planning target volume should be covered by 95% of prescribed dose. Dose-volume constraints were assigned according to international standards. The comparison between dose distributions was performed evaluating quantities important for respective volumes of interest. For target volumes, heterogeneity index and conformity index methodology were used along with the maximum dose concept. Also, for the comparisons over particular organ at risk, maximum dose or mean dose as well as dose-volume concepts were used. For all target volumes and majority of organs at risk, the differences between 2 reporting modes are statistically insignificant, but this is not the case for bony structured organs at risks: mandible and cochlea. It was observed that D{sub w} is higher than D{sub m} with mean difference of 9.91% (p = 0.000009) of the mandible volume covered with 70 Gy. The same trend was observed for left and right cochlea with difference in mean dose of 8.74% (p = 0.037) and 6.87% (p = 0.029), respectively. The comparative analysis of dosimetric parameters in this study shows that the selection of reporting modes in Monaco treatment planning system can produce dose differences up to 15% in high-density volumes such as mandible and cochlea, which might have clinical consequences.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2018

Achieving stress-constrained topological design via length scale control

A new suite of computational procedures for stress-constrained continuum topology optimization is presented. In contrast to common approaches for imposing stress constraints, herein it is proposed to limit the maximum stress by controlling the length scale of the optimized design. Several procedures are formulated based on the treatment of the filter radius as a design variable. This enables to automatically manipulate the minimum length scale such that stresses are constrained to the allowable value, while the optimization is driven to minimizing compliance under a volume constraint – without any direct constraints on stresses. Numerical experiments are presented that incorporate the following : 1) Global control over the filter radius that leads to a uniform minimum length scale throughout the design; 2) Spatial variation of the filter radius that leads to local manipulation of the minimum length according to stress concentrations; and 3) Combinations of the two above. The optimized designs provide high-quality trade-offs between compliance, stress and volume. From a computational perspective, the proposed procedures are efficient and simple to implement: essentially, stress-constrained topology optimization is posed as a minimum compliance problem with additional treatment of the length scale.

42 ENGINEERING↗

DOE OSTI2024

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗

DOE OSTI2025

Response tailoring of elasto-plastic trusses

Abstract In this work we tailor the response of trusses loaded beyond their yield limit. The truss structures are modeled using finite strain theory and rate-independent elasto-plasticity. We design trusses with a tailored mechanical response that is between “elastic” and “elastic-ideal-plastic” subject to the volume constraint. The design updates are generated by the gradient-based Method of Moving Asymptotes (MMA) solver and the sensitivities of the response functions are computed using a path-dependent adjoint sensitivity analysis. The computations are performed in Matlab.

42 ENGINEERING↗

DOE OSTI2018

Comparison of box-wing and conventional aircraft mission performance using multidisciplinary analysis and optimization

Box-wing aircraft designs have the potential to achieve significant reductions in fuel consumption. Closed non-planar wing designs have been shown to reduce induced drag and the statically indeterminate wing structure can lead to reduced wing weight. In addition, the streamwise separation of the two main wings can provide the moments necessary for static stability and control, eliminating the weight and aerodynamic drag of a horizontal tail. Proper assessment of the disciplinary interactions in box-wing designs is essential to determine any realistic performance benefits arising from the use of such a configuration. This study analyzes both box-wing and conventional aircraft designed for representative regional-jet missions. A preliminary parametric investigation shows a lift-to-drag ratio advantage for box-wing designs, while a more detailed multidisciplinary study indicates that the requirement to carry the mission fuel in the wings leads to an increase of between 5% and 1% in total fuel burn compared to conventional designs. Furthermore, the multidisciplinary study identified operating conditions where the box-wing can have superior performance to conventional aircraft despite the fuel volume constraint.

42 ENGINEERING↗

DOE OSTI2023

An extended three-field principle to scale-bridge the granular micromechanics of polymer-bonded particulate materials

Observed phenomena of damage, viscoelasticity, and viscoplasticity in polymer-bonded particulate materials (PBPM), like polymer-bonded explosive (PBX), are accommodated in a scale-bridging constitutive model based on granular micromechanics. An extended three-field variational principle with min - max optimization problem is introduced to homogenize the bonded aggregates’ grain-scale contact mechanics. Similarly to prior applications of a three-field principle, volume constraints are invoked to embed microscale-to-macroscale correspondences for kinematic and static volumetric measures, like mean strain and its energy-conjugate stress invariant. Unlike prior applications, our three-field principle realizes a strain-driven model for effective constitutive behavior due to bonded intergranular contact. In conclusion, the resultant micro-sphere integration-based homogenization is demonstrated in comparison to existing publicly-available measurements to capture the extreme tension/compression constitutive asymmetry exhibited by PBPMs (specifically PBX-9501) during material softening.

42 ENGINEERING↗

DOE OSTI2017

A hybrid, coupled approach for modeling charged fluids from the nano to the mesoscale

Here, we develop and demonstrate a new, hybrid simulation approach for charged fluids, which combines the accuracy of the nonlocal, classical density functional theory (cDFT) with the efficiency of the Poisson–Nernst–Planck (PNP) equations. The approach is motivated by the fact that the more accurate description of the physics in the cDFT model is required only near the charged surfaces, while away from these regions the PNP equations provide an acceptable representation of the ionic system. We formulate the hybrid approach in two stages. The first stage defines a coupled hybrid model in which the PNP and cDFT equations act independently on two overlapping domains, subject to suitable interface coupling conditions. At the second stage we apply the principles of the alternating Schwarz method to the hybrid model by using the interface conditions to define the appropriate boundary conditions and volume constraints exchanged between the PNP and the cDFT subdomains. Numerical examples with two representative examples of ionic systems demonstrate the numerical properties of the method and its potential to reduce the computational cost of a full cDFT calculation, while retaining the accuracy of the latter near the charged surfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Joining technique for in-oven/autoclave molds manufactured by large scale polymer additive manufacturing

Large scale additive manufacturing (AM) of fiber-reinforced polymer composites has gained traction in the mold and die community over the last decade, largely due to reduced material cost and lead time. However, limitations in mold size remain present due to printer volume constraints and costs associated with large format machines. Here, this paper describes a method for producing additively manufactured polymer molds in segments with an O-ring joint design, which can be machined and subsequently assembled for autoclave applications. The developed joining mechanism allows for maintaining the vacuum integrity during the cure cycle while accounting for the anisotropic nature of the additively manufactured polymer molds. Over a range of tested temperatures, 20–200 °C, the multi-part mold joining design was able to achieve full vacuum of ∼73.7 cmHg (29 inchHg) and maintain a vacuum integrity with vacuum losses less than 1.5 cmHg/5min (0.6 inchHg).

Joining↗

DOE OSTI2020

Rate Enhancement of Acid-Catalyzed Alcohol Dehydration by Supramolecular Organic Capsules

The presence of a (sub)nanometric environment has been demonstrated to enhance the rate of homogeneous acid-catalyzed reactions. Conventionally, such catalysts are constructed by directly embedding a catalytically active acid site into the structure of a porous material, linking the acid properties with the properties of the embedding solid. Here, a holoenzyme mimicking approach was used, in which a self-assembled hexameric resorcinarene capsule (1.4 nm 3 cavity) acts like an apoenzyme that transforms into a catalytically active site by hosting hydrated HCl as a cofactor. This capsule binds the hydronium ion together with cyclohexanol (CyOH) and catalyzes the dehydration to cyclohexene with enhanced rates that are 2 orders of magnitude higher than those in the unconstrained environment. Kinetic analysis shows that cyclohexanol dehydration proceeds via an E2 mechanism on hydrated HCl in the unconstrained environment, while the constrained environment of the active capsule induces the stabilization of an ionic intermediate, allowing an E1 mechanism to dominate. By comparison with the reaction in zeolite environments, we show that the direct chemical environment of the capsule or a zeolite wall exerts only a minor direct influence. The volume constraint increases the reaction rate by inducing a larger activation entropy, indicating the stabilization of a later transition state in the E1 mechanism, as well as a higher reaction space of the carbenium ion constituting the transition state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Densified vertically lamellar electrode architectures for compact energy storage

As one of the most compact electrochemical energy storage systems, lithium-ion batteries (LIBs) are playing an indispensable role in the process of vehicle electrification to accelerate the shift to sustainable mobility. Making battery electrodes thicker is a promising strategy for improving the energy density of LIBs which is essential for applications with weight or volume constraints, such as electric-powered transportation; however, their power densities are often significantly restricted due to elongated and tortuous charge traveling distances. Here, we propose an effective methodology that couples bidirectional freeze-casting and compression-induced densification to create densified vertically lamellar electrode architectures for compact energy storage. The vertically lamellar architectures not only overcome the critical thickness limit for conventional electrodes but also facilitate and redistribute the lithium-ion flux enabling both high rate capability and stable cyclability. Furthermore, this proposed methodology is universal as demonstrated in various electrochemical active material systems. This study offers a facile approach that realizes simultaneous high energy and high power in high-loading battery electrodes and provides useful rationales in designing electrode architectures for scalable energy storage systems.

42 ENGINEERING↗