Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sigmoid”

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 19 records

Radial Wall Jet Flow Over Sigmoidal Surfaces

Radial wall jet flows across flat smooth surfaces have been studied for decades. These studies show that the radial velocity of these jets decays inversely with distance from the nozzle with modest contribution from friction (Poreh, et al., 1967; Rajaratnam, 1976). However, the extent to which flat surface results apply to curved surfaces remains unclear. In this article we explore the influence of settled particle bed slope on these velocity profiles. We model the step change in thickness as a sigmoidal curve of variable steepness and use conservation of momentum to evaluate the velocity profile. We show that the velocity profile attenuates because of curvature. Next, we evaluate the influence of surface curvature on radial wall jet flow. Jet flows over particle beds often introduce curvature in the particle bed profile but the influence of the developed curvature on the velocity profile has not been explored. Here we develop a solution for steady fixed beds based on conservation of momentum. We find that surface curvature has a significant influence on the velocity decay coefficients and celerity of radial wall jets, provided there is velocity slip in the vicinity of the particle bed interface, which is strictly true for particle surfaces. Conservation of momentum predicts conditions where the forward momentum of the flow is directed completely upward. The solution identifies two new dimensionless groups that determine whether a curved surface is sufficient to block radial flow and force flow vertically.

Surface curvature, jet flow, multiphase flow↗

Simplifying activations with linear approximations in neural networks

A key step in Neural Networks is activation. Among the different types of activation functions, sigmoid, tanh, and others involve the usage of exponents for calculation. From a hardware perspective, exponential implementation implies the usage of Taylor series or repeated methods involving many addition, multiplication, and division steps, and as a result are power-hungry and consume many clock cycles. We implement a piecewise linear approximation of the sigmoid function as a replacement for standard sigmoid activation libraries. This approach provides a practical alternative by leveraging piecewise segmentation, which simplifies hardware implementation and improves computational efficiency. In this paper, we detail piecewise functions that can be implemented using linear approximations and their implications for overall model accuracy and performance gain. Our results show that for the DenseNet, ResNet, and GoogLeNet architectures, the piecewise linear approximation of the sigmoid function provides faster execution times compared to the standard TensorFlow sigmoid implementation while maintaining comparable accuracy. Specifically, for MNIST with DenseNet, accuracy reaches 99.91% (Piecewise) vs. 99.97% (Base) with up to 1.31x speedup in execution time. For CIFAR-10 with DenseNet, accuracy improves to 98.97% (Piecewise) vs. 99.40% (Base) while achieving 1.24x faster execution. Similarly, for CIFAR-100 with DenseNet, the accuracy is 97.93% (Piecewise) vs. 98.39% (Base), with a 1.18x execution time reduction. These results confirm the proposed method’s capability to efficiently process large-scale datasets and computationally demanding tasks, offering a practical means to accelerate deep learning models, including LSTMs, without compromising accuracy.

Activation function↗

A functional form for liquid scintillator pulse shapes

Digitization of detector signals enables analysis of the original waveform to extract timing, particle identification, and energy deposition information. Here we present the use of analytical functions based on sigmoids to model and fit such pulse shapes from liquid organic scintillators, though the method should also be applicable to other detector systems. Neutron and gamma interactions in NE213 detectors were digitized from the phototube anode and fit using a sigmoid-based function. The acuity of the fit in extracting timing information and performing neutron-gamma pulse-shape discrimination are presented and discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modeling approaches for addressing unrelaxable bound constraints with unconstrained optimization methods

We explore novel approaches for solving nonlinear optimization problems with unrelaxable bound constraints, which must be satisfied before the objective function can be evaluated. Our method reformulates the unrelaxable bound-constrained problem as an unconstrained optimization problem that is amenable to existing unconstrained optimization methods. The reformulation relies on a domain warping to form a merit function; the choice of the warping determines the level of exactness with which the unconstrained problem can be used to find solutions to the bound-constrained problem, as well as key properties of the unconstrained formulation such as smoothness. We develop theory when the domain warping is a multioutput sigmoidal warping, and we explore the practical elements of applying unconstrained optimization methods to the formulation. Here, we develop an algorithm that exploits the structure of the sigmoidal warping to guarantee that unconstrained optimization algorithms applied to the merit function will find a stationary point to the desired tolerance.

97 MATHEMATICS AND COMPUTING↗

Influence of alloy solidification path on melt pool behavior in additive manufacturing

Numerical models used to study transport phenomena in laser-powder bed fusion processes often rely on assumptions and simplifications to reduce their computational expense. One common simplification is in the description of latent heat evolution during the solid-liquid phase change (i.e., the solidification pathway), justified by the fact that the mushy zone thickness is similar to the numerical grid spacing used for continuum transport models. The lack of resolution of transport phenomena in the mushy zone motivates the use of computationally convenient solidification paths such as linear or sigmoidal relationships over pathways derived from fundamental solidification theory such as equilibrium or Scheil models. In the present work, an uncertainty quantification (UQ) framework is used to analyze the influence of solidification pathway selection on the solidification dynamics and melt pool geometries in laser based additive manufacturing (AM) of IN625. Results show the solidification pathway has a quantifiable influence on the cooling rate at the liquidus isotherm, mushy zone thickness, and solidification time. Due to similarities in the latent heat evolution at the beginning of solidification, the equilibrium and Scheil models predict similar cooling rates near the liquidus isotherm, however the wider freezing range of Scheil leads to a wider mushy zone compared to equilibrium. The non-physical latent heat release profiles of sigmoidal and linear paths lead to significant overpredictions of cooling rates at the liquidus isotherm compared to equilibrium and Scheil. Finally, these results indicate that careful consideration should be given to the choice of solidification pathway to ensure reliable model predictions.

36 MATERIALS SCIENCE↗

Predicting Li-Ion Battery Capacity Fade Using Early-Life Data and a Hybrid Data-Driven Gaussian Process-Bayesian Regression Approach

Accurately predicting Li-ion battery capacity trajectories using early-life data can dramatically improve battery-life understandings and be used to rapidly evaluate design/cost/performance trade-offs when developing new battery materials. Accurate early-life predictions enable researchers to quickly iterate over cell designs and material precursor properties without consistently cycling cells to failure. To this end, we present a toolbox that uses a combined Gaussian Process and Bayesian regression approach that capitalizes on signals other than just capacity (e.g., dQ/dV, voltage drops) to rapidly predict capacity-fade trajectories. The prediction tool uses Bayesian regression to fit functional forms, e.g., power law, sigmoids, etc., to predict capacity-fade dynamics. By fitting functional forms, the capacity fade can be interrogated at any point in the future, allowing for early cell-failure prediction. Additionally, Bayesian regression allows for accurate uncertainty estimates that account for cell-to-cell variability (aleatoric uncertainty) and the lack of observation data (epistemic uncertainty). By only using early cycle data to predict the capacity fade trajectory, uncertainty bounds at end-of-life can be extremely large. The large uncertainty bounds are further exacerbated because there is no systematic way to define the prior distribution of the functional forms' parameters. We improve our the predicted trajectory confidence interval of our predicted trajectory using two methods. First, we shows that a small amount of held-out cycling data is sufficientuse some train cells, that have been cycled to failure to derive information regarding the appropriate prior distributions for the functional forms' parameters of the functional form, effectively leading to data-driven priors.. We propose constructing the data-driven priors by first running a Bayesian regression starting with uninformed priors to generate intermediate cell-specific posterior parameter distributions. These posterior distributions are combined using a Ggaussian mixture model for each parameter to create the data-driven priors. These mixture models serve as the data-driven prior distributions for the parameters for. Second, we derive multiple features, e.g., C_dchg 0.5 DoD 0.5, log (|mean(dQ/dV_(w_3-w_0 ) (V)|), etc., from the train cellsheld-out cycling data, identify which the features are that best predicting capacity at early/mid-life cycles, and then create Ggaussian process regression models that are used for predicting capacity at early/mid-life cycles for the test cells (see blue dots with error bars in Fig 1b). Finally, these predicted data-points are used in addition to the actual early cycle data capacity fade to construct the Bayesian regression trajectory for the test cell s. Notably. We note that these two methods are complementary and can be combined with each other. We evaluate the performance of our proposed method on an testing open-source dataset from Iowa State University and Iowa Lakes Community College (ISU-ILCC). This dataset comprises of 251 nickel-manganese-cobalt/graphite Lithium-ion cells that are cycled under 63 different conditions. We compute the mean average percentage error (MAPE) and negative log predictive density (NLPD) to quantify the efficacy of our method. Our initial findings suggest that, when only few observations are available, for test cells, when using only Bayesian regression with uninformed priors, a power law functional provides the most accurate predictions. with very few data points. However, asHowever, a the number of data points increases, a twin sigmoidal function becomes more accurate as the number of observations further increases. We also find that using as little as 10% of the data set towards generating data-driven priors can lead to significant improvement in prediction accuracy when using early cycle data. Lastly, we found that augmenting early-cycle data with Gaussian process-predicted capacity data for Bayesian regression greatly improves the prediction accuracy. We will present a comprehensive comparison of our methods to other methods available in the literature and apply this method to additional battery datasets.

42 ENGINEERING↗

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion

Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.

97 - MATHEMATICS AND COMPUTING↗

Automated Image Segmentation and Processing Pipeline Applied to X–Ray Computed Tomography Studies of Pitting Corrosion in Aluminum Wires

Understanding pitting corrosion is critical, yet its kinetics and morphology remain challenging to study from X-ray computed tomography (XCT) due to manual segmentation barriers. To address this, an automated pipeline leveraging deep learning for efficient large-scale XCT analysis is developed, revealing new corrosion insights. The pipeline enables pit segmentation, 3D reconstruction, statistical characterization, and a topological transformation for visualization. Here, the pipeline is applied to 87 648 XCT images capturing commercial purity aluminum (1100 Al) wire exposed to sodium chloride (NaCl) salt particles over a period of 122 h. The pipeline achieves complete feature extraction and statistical quantification across the entire XCT dataset, leveraging distributed computing environment for high efficiency. Global growth kinetics such as high-level stepwise sigmoidal volume loss patterns and granular individual pit developments are both captured for 36 detected pits. By combining automation, computer vision, and extensive XCT datasets, this research accelerates precise corrosion assessment to enable materials science discoveries at scale.

36 MATERIALS SCIENCE↗

Evaluation of the rheological and electrical percolation of high‐density polyethylene/carbon black composites using mathematical models

Abstract In this work, conductive polymer composites (CPCs) of bio‐based polyethylene (BioPe) containing different concentrations of carbon black (CB) were developed. By using oscillatory rheology analysis, a Newtonian plateau was observed in BioPe, and all BioPe/CB composites had a behavior of a pseudo‐solid and that composites with volume fractions ranging from 0.24 to 0.56 presented higher viscosity, storage, and loss modulus. This suggests the formation of a percolated network and by using the power‐law models, it was observed that the electrical percolation threshold was higher than the rheological percolation threshold. The electrical conductivity was measured using the four‐point probe method and a sigmoid model was used to predict the CPCs' electrical conductivity percolation threshold. The results indicated that the four‐point probe method presented satisfactory results according to the calculated standard deviations and voltage–current characteristics for each round of measurements considering the same ranging as used in rheology analysis. The analytical model used showed a coefficient of determination ( R 2 ) higher than 95%, allowing the prediction of the electrical conductivity of the CPC and the percolation threshold as a function of the volumetric fraction of the CB.

da SiIva, Moacy P.↗

Towards a heat- and mass-balanced kinetic model of TATB decomposition

We report TATB (1,3,5-triamino-2,4,6-trinitrobenzene) was thermally degraded by two small-scale analytical methods – simultaneous differential scanning calorimetry and thermogravimetric analysis (SDT) and a hot-stage microscope with Fourier Transform Infrared (FTIR) analysis capabilities. SDT used ramped heating, isothermal soaking, and thermal pretreatment at various conditions. The heat flow and mass loss were monitored during various treatment conditions to derive chemical decomposition kinetics and Arrhenius parameters. FTIR experiments used isothermal heating, and changes were monitored spectroscopically. Solid samples generated at specific conditions were collected from both methods and were analyzed by DMSO extraction followed by chemical speciation by optical and mass spectrometric methods. Characterization provided the following reaction insights: TATB decreases in a sigmoidal pattern in isothermally heated samples. Other soluble products gradually increase in concentration and then abruptly decline in concentration during the second exotherm, such as diamino-dinitro-benzofurazan and amino-nitro-benzodifurazan. FTIR showed gradual changes in the amino and nitro functionality, shifting positions and decreasing intensity for the first 40 min. Then the solid gradually appeared more like an amorphous C with N incorporated, similar to previous studies on thermally degraded TATB-type materials. Extracted residues (DMSO-soluble components removed) examined by FTIR showed an abrupt change in chemical composition between 40- and 45-min isothermal treatment, indicating early forming solids are different than later forming residues. A reliable mass- and energy-balanced global reaction network must include at least two autocatalytic reactions, either in parallel or series, and at least one must have an explicit initiation reaction having a low activation energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A simple introduction to the SiMPL method for density-based topology optimization

We introduce a novel method for solving density-based topology optimization problems: Sigmoidal Mirror descent with a Projected Latent variable (SiMPL). The SiMPL method (pronounced as “the simple method”) optimizes a design using only first-order derivative information of the objective function. The bound constraints on the density field are enforced with the help of the (negative) Fermi–Dirac entropy, which is also used to define a non-symmetric distance function called a Bregman divergence on the set of admissible designs. This Bregman divergence leads to a simple update rule that is further simplified with the help of a so-called latent variable. Because the SiMPL method involves discretizing the latent variable, it produces a sequence of pointwise-feasible iterates, even when high-order finite elements are used in the discretization. Numerical experiments demonstrate that the method outperforms other popular first-order optimization algorithms. In conclusion, to outline the general applicability of the technique, we include examples with (self-load) compliance minimization and compliant mechanism optimization problems.

Calculus of Variations and Optimization↗

Attention-based 3D – convolutional neural network model for mechanical property predictions using visible light images in metal additive manufacturing

Additive manufacturing (AM), while commonly used for rapid prototyping and creating components with complex geometries, has not been widely adopted for critical applications across the aerospace, automotive, defense, energy, and medical industries. This is, in part, due to the challenges of controlling flaws and uncertainty in the mechanical behavior of additively manufactured components. In recent years, there has been an increase in research aimed at predicting the final mechanical properties of additively manufactured components during the printing process. To address these issues, a 3D-CNN model was trained using low-cost in situ visible-light camera data, anomaly classifications, and the chosen process parameters to predict the ultimate tensile strength (UTS), yield strength (YS), total elongation (TE), and uniform elongation (UE). The 3D-CNN layers of the model employed attention mechanisms to prioritize features in the data, thereby improving prediction accuracy. Furthermore, the effect of each process parameter and anomaly class is investigated using attention-based dynamic sigmoid weighted gates to interpret the influence each class has on the final prediction. Different combinations of the in situ data were fed into the 3D-CNN, with varying amounts of image layers, to determine the ideal combination for predicting mechanical properties in situ. Here, the 3D-CNN model achieved mean absolute percentage errors (MAPE) below 5% for both UTS and YS while using only a single camera input and under half of the available image layers.

36 MATERIALS SCIENCE↗

A hierarchical stochastic modeling approach for representing point bar geometries and petrophysical property variations

The flow of fluids in point bars is affected by the existence of heterogeneities like shale drapes that are found on the surfaces of inclined heterolithic stratifications. In fact, these shale drapes can act as fluid flow baffles; therefore, developing a framework for modeling point bars and their associated heterogeneities is vital. In this study, a stochastic process-based modeling approach is presented for capturing the main point bar heterogeneities: accretion surfaces (i.e., the aerial heterogeneity) and inclined heterolithic stratifications (i.e., the vertical heterogeneity). The former was modeled using a sine-generation function and the latter, with a sigmoidal function after which they were combined into a 3D point bar model. To ensure proper modeling of petrophysical properties, we developed a more representative gridding scheme which generates curvilinear grids representative of the point bar geometry. This grid was then transformed into a rectilinear grid to allow for geostatistical simulation after which all petrophysical properties were mapped back into the original curvilinear grid. An essential element of this modeling approach is the stochastic representation of shale drapes at the interface between successive accretion surfaces. The workflow was tested using a real field dataset for the Cranfield field, Mississippi. The constructed model was then subjected to a flow simulation study mimicking a CO 2 storage scenario. Various sensitivities were simulated to evaluate the effect of heterogeneities on CO 2 flow within the point-bar. Results demonstrate the importance of representing point-bar related heterogeneity and the spatial distribution of shale drapes on CO 2 plume migration and storage.

58 GEOSCIENCES↗

Comparison of three state-of-the-art crystal plasticity based deformation twinning models for magnesium alloys

In magnesium alloys, deformation twinning and its interactions with dislocation slip are responsible for a sigmoidal shape stress–strain behavior and an asymmetrical tension–compression yield strength in magnesium alloys. The sensitivity of twinning to the underlying microstructure renders the crystal plasticity method the most commonly adopted modeling approach for magnesium-twinning. This paper compares three state-of-the-art crystal plasticity-based twinning models from the literature, namely the elastic-viscoplastic self-consistent twinning-detwinning (EVPSC-TDT) model, crystal plasticity finite element model based on enhanced predominate twin reorientation approach (CPFE-ePTR), and the crystal plasticity finite element model based on “discrete twinning” approach (CPFE-DT). A polycrystalline microstructure is simulated with all three methods to compare the resulting stress–strain curves and lattice strains to those from the experimentally measured data. All three methods showed the capability of capturing the experimental results with different levels of accuracy. Additionally, the EVPSC-TDT method avoids solving the finite element matrices and showed the highest computational efficiency. The CPFE-ePTR model shows a higher accuracy in capturing the lattice strain. The CPFE-DT relies on high-resolution finite element mesh and is much slower than the other two methods, but it captured the local deformation concentration and stress reversal phenomena near the twin band, which was not possible with the other two methods. Based on the comparison, guidance for the selection of the appropriate model based on the specific modeling target is provided in this paper.

36 MATERIALS SCIENCE↗

Synchrotron X-ray-induced Synthesis of Copper Hydroxide Nitrate Nanoplates on Cu Thin Films in an Ambient Atmosphere

Synchrotron X-rays are widely used for material characterizations. However, they can also ionize atoms and molecules to damage and manipulate probed materials. We report here an X-ray-induced growth of copper hydroxide nitrate, Cu 2 (OH) 3 NO 3 , on copper thin films in the ambient atmosphere without solvents and thermal treatment. In situ synchrotron X-ray diffraction measurements showed that the time-dependent growth process of the Cu 2 (OH) 3 NO 3 is accompanied by the consumption of Cu metal and can be described by a sigmoidal model. The growth rate was reduced after the initial fast growth period. Scanning electron microscopy (SEM) images show that the isolated islands of Cu 2 (OH) 3 NO 3 nanoplates formed in the beginning, which grew together with new nanoplates formed under continued X-ray irradiation. The result demonstrated that high-flux synchrotron X-rays may provide an unconventional approach to synthesizing and manipulating materials, which will inspire future investigation both experimentally and theoretically.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comprehensive Insights into Nucleation, Autocatalytic Growth, and Stripping Efficiency for Lithium Plating in Full Cells

Synchrotron high-energy X-ray diffraction is used to enable nondestructive detection and quantification of heterogeneous lithium plating in working batteries. In this study a LiNi 0.5 Mn 0.3 Co 0.2 O 2 /graphite pouch cell was operated under 6C fast-charge rate for greater than 1200 cycles. The magnitude and spatial distribution of lithium plating, lithium stripping, and the effect of metallic lithium deposition on lithium intercalation into graphite were quantified. Fully intercalated graphite (LiC 6 ) was detected after discharge with a lateral distribution closely correlated with lithium plating, which can be used as a higher-sensitivity indicator for lithium plating. Over an extended cycle life, the overall metallic lithium concentration followed a sigmoidal curve indicating two-stage continuous nucleation and autocatalytic growth. The lithium stripping efficiency underwent an exponential decay as a function of cycle life as the buildup of metallic lithium hindered the efficient dissolution back into the electrolyte. The findings provide direct insights into the characteristics of lithium plating and stripping under realistic fast-charge conditions.

25 ENERGY STORAGE↗

Decadal-scale recovery of carbon stocks after wildfires throughout the boreal forests

Boreal forests store 30% of the world's terrestrial carbon (C). Consequently, climate change mediated alterations in the boreal forest fire regime can have a significant impact on the global C budget. Here, we synthesize the effects of forest fires on the stocks and recovery rates of C in boreal forests from 368 plots from 16 long-term (= 100 years) fire chronosequences spread throughout the boreal zone. The mean annual temperature varied from -9.5°C to +5.2°C and annual precipitation from 295 mm to 824 mm in the study sites. Total C stocks decreased on average by 60% (range from <10% to >80%) because of forest fires, primarily due to large declines in the living trees and O-layer C stocks. C emissions of fires ranged from 0.1 to 11.4 kg m-2, and in most cases, the emissions were < 4 kg m-2. Total C stocks increased with time since fire following a sigmoidal shape Gompertz function, and it took 125 years until an average asymptote of 8.1 kg C m-2 was reached. Total C stocks accumulated at a rate of 2–60 g m-2 yr-1 during the first 100 years. We found that potential evapotranspiration (PET) was significantly correlated with C stocks and their post-fire recovery, likely because it integrates temperature, radiation and the length of the growing season. If the fire return interval shortens to under 100 years in the future, this suggests that many boreal forests will be prevented from reaching their full C storage potential. However, our results suggest that climate warming-induced increases in PET may also speed up the post-fire recovery of C stocks.

carbon accumulation, carbon pools, chronosequence,↗