Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compression techniques”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measurement of needle and armature dynamics in a gasoline direct injector by high-speed neutron imaging

In modern spark ignition engines, precise delivery of fuel with gasoline direct injection has become increasingly important in the effort to meet ever stricter efficiency and emissions regulations. Use of multiple small close-coupled injections has become more common in attempt to precisely control fuel distribution in the cylinder, but these strategies are hindered by nonlinear injection effects due to operation in the ballistic region of the solenoid-operated valve and due to armature bounce at the end of injection. Understanding the internal dynamics of the injector is crucial to minimizing and controlling non-linearity and shot-to-shot variation, but the experimental techniques available to date are capable only of tracking the position of either the top of the needle (via laser sensors or by monitoring current and voltage in the solenoid coil) or the bottom of the needle (via transparent nozzles or high-speed x-ray imaging). Furthermore, a complete picture of the axial and radial motion of valve needle has until now remained elusive. In this work, we present high-speed ensemble neutron transmission imaging of an entire 8-hole gasoline direct injector operating at 200 bar, allowing for both visualization and quantification of the actuation dynamics including lift and wobble of the valve ball, oscillation and bending of the valve needle, lift, rocking, and bounce of the armature, compression of the springs, and radial swelling of the solenoid during energization. Because neutrons offer high penetration through the metal injector while also being sensitive to the 1 H nuclei in fuel molecules, it is also possible to simultaneously see the fluid dynamics of the injection process, including filling of the sac volume, emanation of the spray through the nozzle holes and into the downstream gas, and the formation and evolution of fuel films on the tip of the injector and the walls of the spray container.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Thermodynamic properties and enhancement of diamagnetism in nitrogen doped lutetium hydride synthesized at high pressure

Nitrogen doped lutetium hydride has drawn global attention in the pursuit of room-temperature superconductivity near ambient pressure and temperature. However, variable synthesis techniques and uncertainty surrounding nitrogen concentration have contributed to extensive debate within the scientific community about this material and its properties. We used a solid-state approach to synthesize nitrogen doped lutetium hydride at high pressure and temperature (HPT) and analyzed the residual starting materials to determine its nitrogen content. High temperature oxide melt solution calorimetry determined the formation enthalpy of LuH 1.96 N 0.02 (LHN) from LuH 2 and LuN to be −28.4 ± 11.4 kJ/mol. Magnetic measurements indicated diamagnetism which increased with nitrogen content. Ambient pressure conductivity measurements observed metallic behavior from 5 to 350 K, and the constant and parabolic magnetoresistance changed with increasing temperature. High pressure conductivity measurements revealed that LHN does not exhibit superconductivity up to 26.6 GPa. We compressed LHN in a diamond anvil cell to 13.7 GPa and measured the Raman signal at each step, with no evidence of any phase transition. Despite the absence of superconductivity, a color change from blue to purple to red was observed with increasing pressure. Thus, our findings confirm the thermodynamic stability of LHN, do not support superconductivity, and provide insights into the origins of its diamagnetism.

Science & Technology - Other Topics↗
DOE OSTI2026

Dark Energy Survey Year 3 results: Simulation-based 𝑤CDM inference from weak lensing and galaxy clustering maps with deep learning: Analysis design

Data-driven approaches using deep learning are emerging as powerful techniques to extract non-Gaussian information from cosmological large-scale structure. Here, this work presents the first simulation-based inference (SBI) pipeline that combines weak lensing and galaxy clustering maps in a realistic Dark Energy Survey Year 3 (DES Y3) configuration and serves as preparation for a forthcoming analysis of the survey data. We develop a scalable forward model based on the CosmoGridV1 suite of N-body simulations to generate over one million self-consistent mock realizations of DES Y3 at the map level. Leveraging this large dataset, we train deep graph convolutional neural networks on the full survey footprint in spherical geometry to learn low-dimensional features that approximately maximize mutual information with target parameters. These learned compressions enable neural density estimation of the implicit likelihood via normalizing flows in a ten-dimensional parameter space spanning cosmological 𝑤CDM, intrinsic alignment, and linear galaxy bias parameters, while marginalizing over baryonic, photometric redshift, and shear bias nuisances. To ensure robustness, we extensively validate our inference pipeline using synthetic observations derived from both systematic contaminations in our forward model and independent Buzzard galaxy catalogs. Our forecasts yield significant improvements in cosmological parameter constraints, achieving 2−3× higher figures of merit in the 𝛺 𝑚 − 𝑆 8 plane relative to our implementation of baseline two-point statistics and effectively breaking parameter degeneracies through probe combination. These results demonstrate the potential of SBI analyses powered by deep learning for upcoming Stage-IV wide-field imaging surveys.

Thomsen, A. [Zurich, ETH] (ORCID:0000000203099021)↗
DOE OSTI2020

A transported Livengood-Wu integral model for knock prediction in CFD simulation

This work describes the development of a transported Livengood-Wu (L-W) integral model for computational fluid dynamics (CFD) simulation to predict auto-ignition and engine knock tendency. The currently employed L-W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition combustion for three different fuels, showing good accuracy in prediction of auto-ignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L-W integral model is coupled with G-equation model to indicate end-gas auto-ignition and knock tendency in CFD simulations of a direct injection spark ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely Co-Optima Alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multi-variable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency indicated by the fuel energy contained in the auto-ignited region of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

Yue, Zongyu↗
DOE OSTI2025

Compressing Hamiltonians with ab initio downfolding for simulating strongly-correlated materials on quantum computers

The accurate first-principles description of strongly correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving compressed many-body Hamiltonians that maintain the essential physics of strongly correlated materials. The solution of these material-specific models is still exponentially difficult to generate on classical computers, but quantum algorithms allow for a significant speed-up in obtaining the ground states of these compressed Hamiltonians. Here, we demonstrate that using quantum algorithms to obtain the properties of downfolded Hamiltonians can indeed yield high-fidelity solutions. By combining ab initio downfolding and variational quantum eigensolvers, we correctly predict the antiferromagnetic state of one-dimensional cuprate Ca 2 Cu O 3 , the excitonic ground state of monolayer W Te 2 , and the charge-ordered state of correlated metal Sr VO 3 . Numerical simulations using a classical tensor network implementation of variational quantum eigensolvers allow us to simulate large models with up to 54 qubits and encompassing up to four bands in the correlated subspace, which is indicative of the complexity that our framework can address. Through these methods we demonstrate the potential of classical preoptimization and downfolding techniques for enabling efficient materials simulation using quantum algorithms.

Alvertis, Antonios M. [NASA, Ames; LBNL, Berkeley]↗
DOE OSTI2023

Annual report for DOE VTO

Carbon fiber (CF)/polymer composites are a transformative class of high-performance, lightweight material, where high aspect-ratio CFs reinforce a polymer matrix and exceed the strength of steel alloys at a fraction of the density. Despite the advantages of such a class of material, the broader implementation of CF composites in a range of automotive, aerospace, and energy applications is hindered by limitations of current manufacturing methods. These current techniques (e.g., hand lay-up, wet filament winding) are costly and impose severe limitations on fiber placement, orientation, and angle, and thus a composite’s ultimate properties. Today’s CF composites are expensive to manufacture, limited in form factor, and utilize costly and sub-optimal continuous filament CF. Advanced additive manufacturing (AM) processes, combined with computational design optimization and new approaches to resin development, offer alternative design and manufacturing paradigms that have the realistic potential to lift these constraints. Such integrated AM approaches could thus help to realize the full potential of CF composite materials. One relevant application of CF composite materials where manufacturing constraints limit the cost-benefit ratio is in the manufacture of high-performance composite pressure vessels for onboard compressed natural gas (CNG) storage. Current CNG storage vessels (Types 3–5) are made from load-bearing filament-wound carbon-fiber composite and are ~3.5 times as expensive as an all-metallic Type-1 vessel. This cost is invariably tied to the complex and labor-intensive nature of conventional filament winding processes and the large volumes of expensive high tensile-strength CF tow feedstock required in manufacture. Our proposed approach to CNG storage vessel manufacture is based on a combination of AM technologies for CF composite printing and design optimization tools that were pioneered at Lawrence Livermore National Laboratory (LLNL), with advances in resin/composite formulation enabled by chemical and nano-material modification. Through the successful development of this technology, LLNL seeks to demonstrate the capability for advanced CNG storage vessel manufacture at reduced cost with no reduction in performance versus the most advanced, extant Type-5 designs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Carbide Formation in Refractory Mo 15 Nb 20 Re 15 Ta 30 W 20 Alloy under a Combined High-Pressure and High-Temperature Condition

In this work, the formation of carbide with the concertation of carbon at 0.1 at.% in refractory high-entropy alloy (RHEA) Mo 15 Nb 20 Re 15 Ta 30 W 20 was studied under both ambient and high-pressure high-temperature conditions. The x-ray diffraction of dilute carbon (C)-doped RHEA under ambient pressure showed that the phases and lattice constant of RHEA were not influenced by the addition of 0.1 at.% C. In contrast, C-doped RHEA showed unexpected phase formation and transformation under combined high-pressure and high-temperature conditions by resistively employing the heated diamond anvil cell (DAC) technique. The new FCC_L1 2 phase appeared at 6 GPa and 809 °C and preserved the ambient temperature and pressure. High-pressure and high-temperature promoted the formation of carbides Ta 3 C and Nb 3 C, which are stable and may further improve the mechanical performance of the dilute C-doped alloy Mo 15 Nb 20 Re 15 Ta 30 W 20 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

In situ x-ray diffraction study of dynamically compressed α -cristobalite using a dynamic diamond anvil cell

Here we present results of the dynamic compression of α -cristobalite up to a pressure of 106 GPa with the use of the dynamic diamond anvil cell. X-ray diffraction images were recorded at different ramp compression and decompression rates to investigate in situ the high-pressure phase transitions of α -cristobalite. Our results suggest that the pressure onset of the phase transformation of α -cristobalite to cristobalite II, cristobalite X-I, and ultimately to seifertite ( α – PbO 2 type SiO 2 ) is dependent on the applied compression rates and stress conditions of the experiment. Increasing compression rates in general shift the studied phase transitions to higher pressures. Furthermore, our results indicate for single crystals under hydrostatic conditions a suppression of a phase transition from cristobalite X-I to seifertite at pressures of up to 82 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The dependence of aluminum lattice orientation on semiconductor lattice parameter in planar InAs/Al hybrid heterostructures

We present a detailed study of 7-nm thick Al layers epitaxially grown on In x Al 1-x As or Al 1-x Ga x Sb y As 1-y (0 0 1) semiconductors using high-resolution transmission electron microscopy (HRTEM) and synchrotron X-ray diffraction (XRD). These techniques provide high spatial and high reciprocal space resolution information about Al lattice orientations and interfacial structure, highlighting the existence of a strong correlation between the surface lattice parameter of the semiconductor and the normal orientation of the aluminum film. The normal orientation of Al layers is found to be [1 1 1] when the in-plane semiconductor surface lattice parameter, a S , is smaller than 5.98 Å, and is [1 1 0] when aS is larger than this value. In the transition between these lattice parameters where one orientation is preferred, both orientations may be observed. Concomitant with the change in normal orientation is an inversion in the strain state of the aluminum film. When the normal orientation is [1 1 1], the Al film is compressively strained; while for [1 1 0] normal orientations, the Al in-plane strain is tensile. An energy balance model, accounting for surface, interfacial, side-face, and strain energy, and the registry of lattice planes across the semiconductor/aluminum interface provide explanations for the observed normal orientation and strain inversion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Intern Poster

Digital Image Correlation (DIC) is an optical technique that combines image registration and tracking methods for accurate two-dimensional and three-dimensional changes in images. DIC software can be used to track the contour, deformation, and strain of a sample. In the Advanced Test Reactor (ATR) at INL (Idaho National Laboratory) there exists a small working window of samples that can become irradiated. Hundreds of graphite disks called piggybacks have undergone this irradiation as part of the Advanced Reactor Technologies (ART) program. After irradiation, it is desirable to understand the change in tensile strength as a function of dose. Due to the limited space in the ATR, typical dog bone tensile tests reduce the number of graphite samples from hundreds to tens. However, there does exist an ASTM standard, D8289, which uses disc compression of graphite to estimate the tensile strength of the specimen with the Brazilian Disk test fixture. While only used as an estimate, which is typically off by a third, it is the purpose of this study to identify how to amend D8289 to remove the word "estimate" with the help of DIC.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Physics-informed graph neural networks for predicting cetane number with systematic data quality analysis

Designing alternative fuels for advanced compression ignition engines necessitates a predictive model for cetane number (CN). In this study, the physics-informed graph neural networks are introduced for a reliable CN prediction by considering molecular features pertinent to the physical properties of molecules that affect CN. The reliability of measured data is another key factor to consider for improving the predictive model. Various experimental instruments for measuring CN exist, including standard and non-standard methods. In this regard, a systematic data quality analysis was carried out for the total 630 CNs collected from literature and new measurements in this study using Advanced Fuel Ignition Delay Analyzer (AFIDA). The results from this data curation process were reflected in the model by imposing lower sample weights on the data coming from less reliable measurement techniques. This approach effectively maximized the prediction accuracy while incorporating data from all available sources. Using the sample weights decreased the mean absolute error (MAE) up to 0.8 CN units. The accuracy was also improved by introducing the CN-related physical properties (the number of hydrogen bond donors and acceptors); the test set MAE is 5.74 and 7.01 for the model with and without such properties, respectively. Investigating molecular structural effects on CN was also carried out to gain chemical insights into factors used to design new fuel candidates. The dimensionality reduction analysis of feature vectors showed a clear clustering in terms of functional groups and CN and the structural effect derived from the model was consistent with the physicochemical insights. Finally, this physics-informed model and data curation would be helpful for accurate CN prediction and inform rational fuel design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Measurement of mix at the fuel–ablator interface in indirectly driven capsule implosions on the National Ignition Facility

The interface between the capsule ablator and fuel ice layer is susceptible to hydrodynamic instabilities. The subsequent mixing of hot ablator material into the ice reduces fuel compression at stagnation and is a candidate for reduced capsule performance. The ability to diagnose ice–ablator mix is critical to understanding and improving stability at this interface. Combining the crystal backlighter imager with the single line of sight camera on the National Ignition Facility (NIF) allows direct measurement of ice–ablator mix by providing multiple quasi-monochromatic radiographs of layered capsule implosions per experiment with high spatial (∼12 μm) and temporal (∼35 ps) resolution. The narrow bandwidth of this diagnostic platform allows radiography of the inner edge of the capsule limb close to stagnation without capsule self-emission contaminating the data and removes opacity uncertainties typically associated with the spectral content of the radiograph. Analysis of radiographic data via a parameterized forward-fitting Abel inversion technique provides measurements of the distribution of mix mass inwards from the ice–ablator interface. The sensitivity of this mix measurement technique was demonstrated by applying it to layered experiments in which the stability of the ice–ablator interface was expected to vary significantly. Additional experiments suggest that high-density carbon capsules that employ a buried-layer dopant profile suffer from mixing at the innermost doped–undoped interface. Data from these experiments suggest that opacity models used in hydrodynamic simulations of NIF experiments can potentially over-predict the opacity of doped capsules. LLNL-JRNL-850535-DRAFT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fundamental Investigations into the Properties and Performance of Advanced Materials

Intelligent engineering designs require an accurate understanding of material behavior, since any uncertainties or gaps in knowledge must be counterbalanced with heightened factors of safety, leading to overdesign. Therefore, building better structures and pushing the performance of new components requires an improved understanding of the thermomechanical response of advanced materials under service conditions. This dissertation provides fundamental investigations of several advanced materials: thermoset polymers, a common matrix material for fiber-reinforced composites and nanocomposites; aluminum alloy 7075-T6 (AA7075-T6), a high-performance aerospace material; and ceramic matrix composites (CMCs), an advanced composite for extreme-temperature applications. To understand matrix interactions with various interfaces and nanoinclusions at their fundamental scale, the properties of thermoset polymers are studied at the atomistic scale. An improved proximity-based molecular dynamics (MD) technique for modeling the crosslinking of thermoset polymers is carefully established, enabling realistic curing simulations through its ability to dynamically and probabilistically perform complex topology transformations. The proximity-based MD curing methodology is then used to explore damage initiation and the local anisotropic evolution of mechanical properties in thermoset polymers under uniaxial tension with an emphasis on changes in stiffness through a series of tensile loading, unloading, and reloading experiments. Aluminum alloys in aerospace applications often require a fatigue life of over 109 cycles, which is well over ii the number of cycles that can be practically tested using conventional fatigue testing equipment. In order to study these high-life regimes, a detailed ultrasonic cycle fatigue study is presented for AA7075-T6 under fully reversed tension-compression loading. The geometric sensitivity, frequency effects, size effects, surface roughness effects, and the corresponding failure mechanisms for ultrasonic fatigue across different fatigue regimes are investigated. Finally, because CMCs are utilized in extreme environments, oxidation plays an important role in their degradation. A multiphysics modeling methodology is thus developed to address the complex coupling between oxidation, mechanical stress, and oxygen diffusion in heterogeneous carbon fiber-reinforced CMC microstructures.

Schichtel, Jacob J.↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Considering computational speed vs. accuracy: Choosing appropriate mesoscale RVE boundary conditions

Modeling a material’s microstructure using continuum theories allows for inspection of the relationship between coarse scale and fine scale behaviors. Computational limits generally require selection of a sub-volume from a bulk sample in order to directly model the microstructure. Boundary conditions are applied to the sub-volume to mimic the excluded bulk material. Appropriate selection of boundary conditions helps effectively determine the appropriate spatial scale required of the sub-volume. Applicable boundary conditions include direct displacement, periodic, and uniform traction. While direct displacement and periodic boundary conditions are commonly used, uniform traction boundary conditions have seen limited use due to rigid body stability issues in simulations of compression or shear deformation. A new application of uniform traction boundary conditions was developed through linear constraint equations, similar to approaches employed by direct displacement and periodic boundary conditions, to quench rigid body motions with minimal interference of the relative deformation of the model. These boundary conditions were tested by compressing several synthetically generated periodic microstructures using the finite element method. Evaluating the effective stiffness along the compression axis, the direct displacement boundary condition produced the stiffest response, whereas the uniform traction boundary condition produced the most compliant. Periodic boundary conditions produced the same response for all volumes analyzed and both the direct displacement and uniform traction boundary conditions trended toward the periodic response as the domain volume increased. Computational performance was also evaluated for each boundary condition using implicit and explicit solvers. Direct displacement boundary conditions presented the lowest computational cost of all of the boundary conditions followed by periodic then uniform traction. The computational expense of periodic and uniform traction boundary conditions limited the viable spatial scale and mesh resolutions able to be simulated. Selection of appropriate boundary conditions for specific uses need to be a balance between allowable computational expense and accuracy of the method. Techniques for evaluating which boundary conditions to use are discussed.

42 ENGINEERING↗
DOE OSTI2022

Study of Stark broadening of krypton helium-$\beta$ lines and estimation of electron density and temperature in NIF compressed capsules

Here the National Ignition Facility (NIF) diagnostic instrument manipulator (DIM) - based high resolution (dHIRES) x-ray spectrometer was used to measure the time evolution of the electron density (n e ) and temperature (T e ) in the hot spot of four NIF compressed capsules with 25 ps time resolution during the 'stagnation' phase. The electron density was inferred by comparing the measured Stark broadening of the krypton (Kr) Heβ spectral complex with theoretical calculations that include ion dynamic effects, and the electron temperature was inferred by comparing the measured ratio of the intensity of a dielectronically excited Li-like Kr line to the intensity of the Kr Heβ resonance line with calculations using the spectroscopic collisional radiative atomic model (SCRAM) and CRETIN collisional-radiative models. The inferred, time averaged n e values mainly agree with n e values from neutron diagnostics within uncertainties, but the neutron time-of-flight values of T ion are consistently higher than dHIRES T e values by 200–700 eV. The dHIRES measurements and measurement techniques, method of uncertainty analysis, and discussion of comparisons with measurements from neutron diagnostics are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The Adjoint Petrov–Galerkin method for non-linear model reduction

Here, we formulate a new projection-based reduced-order modeling technique for non-linear dynamical systems. The proposed technique, which we refer to as the Adjoint Petrov–Galerkin (APG) method, is derived by decomposing the generalized coordinates of a dynamical system into a resolved coarse-scale set and an unresolved fine-scale set. A Markovian finite memory assumption within the Mori–Zwanzig formalism is then used to develop a reduced-order representation of the coarse scales. This procedure leads to a closed reduced-order model that displays commonalities with the adjoint stabilization method used in finite elements. The formulation is shown to be equivalent to a Petrov–Galerkin method with a non-linear, time-varying test basis, thus sharing some similarities with the Least-Squares Petrov–Galerkin method. Theoretical analysis examining a priori error bounds and computational cost is presented. Numerical experiments on the compressible Navier–Stokes equations demonstrate that the proposed method can lead to improvements in numerical accuracy, robustness, and computational efficiency over the Galerkin method on problems of practical interest. Improvements in numerical accuracy and computational efficiency over the Least-Squares Petrov–Galerkin method are observed in most cases.

42 ENGINEERING↗
DOE OSTI2025

Hybrid data-driven cement-stabilized soil design: An integration of machine learning, multi-objective optimization, and life cycle assessment

Soil stabilization is crucial in geotechnical engineering, yet conventional methods are often time-consuming, resource-intensive, and environmentally unsustainable. Despite growing interest in Machine Learning (ML) and optimization tools for mix design, few studies integrate these methods with decision-making techniques and environmental assessment to support practical implementation. This study proposes a hybrid data-driven framework for predicting strength, optimizing mix compositions, and evaluating environmental impacts via life cycle assessment of cement-stabilized soft soils. Six ML models were evaluated, and the top-performing eXtreme Gradient Boosting (XGB) model was further improved using the Grey Wolf Optimizer (GWO). The optimized XGB-GWO model, integrated with a polynomial cost function, served as the objective function in a multi-objective optimization problem solved via the Non-Dominated Sorting Genetic Algorithm II (NSGA-II), with final mix selection guided by the entropy-weighted TOPSIS method. Validation through a case study produced mix designs offering superior strength-cost trade-offs, with the optimal mix achieving 2243.2 kPa unconfined compressive strength and a 16.07 % reduction in carbon emissions compared to the highest-cost design. In conclusion, this study offers a sustainable, scalable approach to soil stabilization and supports informed decision-making in construction.

Life cycle assessment↗