Search NASA⌕ Search

SEARCH · Search NASA

Results for “deformation function”

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,081 records · Page 60

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗

Dislocation drag and its influence on elastic precursor decay

Plastic deformation is mediated by the creation and movement of dislocations, and at high stress the latter is dominated by dislocation drag from phonon wind. Here, by simulating a 1-D shock impact problem we analyze the importance of accurately modeling dislocation drag and dislocation density evolution in the high stress regime. Dislocation drag is modeled according to a first-principles derivation as a function of stress and dislocation character, and its temperature and density dependence are approximated to the extent currently known. Much less is known about dislocation density evolution, leading to far greater uncertainty in these model parameters. In studying anisotropic fcc metals with character dependent dislocations, the present work generalizes similar earlier studies by other authors.

36 MATERIALS SCIENCE↗

Modeling and experimental validation of a W f /W-fabrication by chemical vapor deposition and infiltration

Tungsten (W) has a unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall material in future fusion devices. However, its intrinsic brittleness and its further embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) with extrinsic toughening is being developed. A possible synthesis route is chemical vapor deposition (CVD) using heated W fabrics as substrate. The challenge is that the growing CVD-W can isolate domains from precursor access leading to strength-reducing pores. To deepen the process understanding and to optimize the CVD parameters, models were developed with COMSOL Multiphysics and validated experimentally. W deposition rate equations as function of the temperature and the partial pressures of the precursors H 2 and WF 6 were experimentally validated in previous work. In the present article, these equations are applied to obtain partial pressures within the CVD reactor. The results are taken as input for transient simulations in the microscale, in which W coatings, growing onto multiple adjacent W fibers, were simulated via mesh deformation and remeshing. The surface-to-surface contact of the W coatings and the corresponding potential pore formation were simulated by implementing sophisticated deposition rate stop conditions. Within the measuring uncertainties of ≃ ±1%, the models are validated successfully by experimental comparison regarding the deposition rate, pore structure, and relative densities ranging from 0.6 to 0.9.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms for bacterial gliding motility on soft substrates

The motility mechanism of certain prokaryotes has long been a mystery, since their motion, known as gliding, involves no external appendages. The physical principles behind gliding still remain poorly understood. Using myxobacteria as an example of such organisms, we identify here the physical principles behind gliding motility and develop a theoretical model that predicts a 2-regime behavior of the gliding speed as a function of the substrate stiffness. Our theory describes the elasto-capillary–hydrodynamic interactions between the membrane of the bacteria, the slime it secretes, and the soft substrate underneath. Defining gliding as the horizontal translation under zero net force, we find the 2-regime behavior is due to 2 distinct mechanisms of motility thrust. On mildly soft substrates, the thrust arises from bacterial shape deformations creating a flow of slime that exerts a pressure along the bacterial length. This pressure in conjunction with the bacterial shape provides the necessary thrust for propulsion. On very soft substrates, however, we show that capillary effects must be considered that lead to the formation of a ridge at the slime–substrate–air interface, thereby creating a thrust in the form of a localized pressure gradient at the bacterial leading edge. To test our theory, we perform experiments with isolated cells on agar substrates of varying stiffness and find the measured gliding speeds in good agreement with the predictions from our elasto-capillary–hydrodynamic model. Finally, the mechanisms reported here serve as an important step toward an accurate theory of friction and substrate-mediated interactions between bacteria proliferating in soft media.

59 BASIC BIOLOGICAL SCIENCES↗

Insights into the thermo-hydraulic properties of compacted MX80 bentonite during hydration under elevated temperature

In high-level radioactive waste geological repositories, compacted bentonite undergoes coupled thermo-hydraulic processes due to heat released from a central waste canister and groundwater imbibition from the surrounding host rock. An understanding of these processes is essential for long-term simulations of radionuclide migration and canister corrosion, which requires an understanding of temperature effects on the coupled thermo-hydraulic properties governing these processes. In this study, a tank-scale radial infiltration test was used to investigate water imbibition processes in compacted bentonite under a central heater temperature of 200 °C that simulates high thermal gradients in a repository. Interpretation of this test focuses on evaluation of the liquid water wetting front during hydration and the interpretation of the temperature-dependent transient soil water retention curve (SWRC), thermal conductivity function (TCF), and hydraulic conductivity function (HCF). The SWRC during imbibition follows a temperature-dependent wetting path. Temperature effects on the hydraulic conductivity of bentonite in saturated conditions had the greatest effect on the shape of the HCF, with minimal temperature effects at higher suctions. The transient thermal conductivity data matched well with a new TCF linked with the SWRC shape. Although the bentonite layer was restrained, local deformations during hydration may have affected the shapes of the TCF and HCF.

Engineering↗

Sierra/SolidMechanics 4.56 User's Guide

Siera/SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra/SM.

97 MATHEMATICS AND COMPUTING↗

SSME main combustion chamber life prediction

Typically, low cycle fatigue life is a function of the cyclic strain range, the material properties, and the operating temperature. The reusable life is normally defined by the number of strain cycles that can be accrued before severe material degradation occurs. Reusable life is normally signified by the initiation or propagation of surface cracks. Hot-fire testing of channel wall combustors has shown significant mid-channel wall thinning or deformation during accrued cyclic testing. This phenomenon is termed cyclic-creep and appears to be significantly accelerated at elevated surface temperatures. This failure mode was analytically modelled. The cyclic life of the baseline SSME-MCC based on measured calorimeter heat transfer data, and the life sensitivity of local hot spots caused by injector effects were determined. Four life enhanced designs were assessed.

Cook, R. T.↗

AC Electric Field Activated Shape Memory Polymer Composite

Shape memory materials have drawn interest for applications like intelligent medical devices, deployable space structures and morphing structures. Compared to other shape memory materials like shape memory alloys (SMAs) or shape memory ceramics (SMCs), shape memory polymers (SMPs) have high elastic deformation that is amenable to tailored of mechanical properties, have lower density, and are easily processed. However, SMPs have low recovery stress and long response times. A new shape memory thermosetting polymer nanocomposite (LaRC-SMPC) was synthesized with conductive fillers to enhance its thermo-mechanical characteristics. A new composition of shape memory thermosetting polymer nanocomposite (LaRC-SMPC) was synthesized with conductive functionalized graphene sheets (FGS) to enhance its thermo-mechanical characteristics. The elastic modulus of LaRC-SMPC is approximately 2.7 GPa at room temperature and 4.3 MPa above its glass transition temperature. Conductive FGSs-doped LaRC-SMPC exhibited higher conductivity compared to pristine LaRC SMP. Applying an electric field at between 0.1 Hz and 1 kHz induced faster heating to activate the LaRC-SMPC s shape memory effect relative to applying DC electric field or AC electric field at frequencies exceeding1 kHz.

Kang, Jin Ho↗

Jig Twist Optimization of Mach 0.8 Transonic Truss-Braced Wing Aircraft

This paper presents a jig twist optimization study of Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft using FUN3D. The solver has been validated against wind tunnel test data. An aeroelastic model deformation correction method is developed for the wind tunnel model to improve the validation. The flight condition corresponding to Mach 0.8 and design lift coefficient 0.695 is selected for the TTBW aircraft jig twist optimization to reduce the drag coefficient. The static aeroelastic simulation module in FUN3D is used for the aeroelastic simulation. Three different orders of the Chebyshev polynomials are used as shape function to represent the variations of the jig twist distribution along the wing span. The preliminary optimization results show that the Mach 0.8 TTBW aircraft with the optimized jig twist achieves a drag reduction of about 5.7 counts at the design condition when 5 𝑡ℎ order Chebyshev polynomials are used as shape function.

TTBW↗

Aeroelastic Modeling of Elastically Shaped Aircraft Concept via Wing Shaping Control for Drag Reduction

Lightweight aircraft design has received a considerable attention in recent years as a means for improving cruise efficiency. Reducing aircraft weight results in lower lift requirements which directly translate into lower drag, hence reduced engine thrust requirements during cruise. The use of lightweight materials such as advanced composite materials has been adopted by airframe manufacturers in current and future aircraft. Modern lightweight materials can provide less structural rigidity while maintaining load-carrying capacity. As structural flexibility increases, aeroelastic interactions with aerodynamic forces and moments become an increasingly important consideration in aircraft design and aerodynamic performance. Furthermore, aeroelastic interactions with flight dynamics can result in issues with vehicle stability and control. Abstract This paper describes a recent aeroelastic modeling effort for an elastically shaped aircraft concept (ESAC). The aircraft model is based on the rigid-body generic transport model (GTM) originally developed at NASA Langley Research Center. The ESAC distinguishes itself from the GTM in that it is equipped with highly flexible wing structures as a weight reduction design feature. More significantly, the wings are outfitted with a novel control effector concept called variable camber continuous trailing edge (VCCTE) flap system for active control of wing aeroelastic deflections to optimize the local angle of attack of wing sections for improved aerodynamic efficiency through cruise drag reduction and lift enhancement during take-off and landing. The VCCTE flap is a multi-functional and aerodynamically efficient device capable of achieving high lift-to-drag ratios. The flap system is comprised of three chordwise segments that form the variable camber feature of the flap and multiple spanwise segments that form a piecewise continuous trailing edge. By configuring the flap camber and trailing edge shape, drag reduction could be achieved. Moreover, some parts of the flap system can be made to have a high frequency response for roll control, gust load alleviation, and aeroservoelastic (ASE) modal suppression control. Abstract The aeroelastic model of the ESAC is based on one-dimensional structural dynamic theory that captures the aeroelastic deformation of a wing structure in a combined motion that involves flapwise bending, chordwise bending, and torsion. The model includes the effect of aircraft propulsion due to wing flexibility which causes the propulsive forces and moments to couple with the wing elastic motion. Engine mass is also accounted in the model. A fuel management model is developed to describe the wing mass change due to fuel usage in the main tank and wing tanks during cruise. Abstract The model computes both static and dynamic responses of the wing structures. The static aeroelastic deflections are used to estimate the effect of wing flexibility on induced drag and the potential drag reduction by the VCCTE flap system. A flutter analysis is conducted to estimate the flutter speed boundary. Gust load alleviation via adaptive control has been recently investigated to address flexibility of aircraft structures. A multi-objective flight control approach is presented for drag reduction control. The approach is based on an optimal control framework using a multi-objective cost function. Future studies will demonstrate the potential benefits of the approach.

Aeroelasticity↗

Elucidating Lithium Transport Mechanisms in Disordered LiF from Machine-Learning Molecular Dynamics Simulations

Lithium fluoride (LiF) is a ubiquitous component of solid- and cathode–electrolyte interphases, yet its functional role remains unclear under the structural and chemical heterogeneity typical of cycling batteries. Here, we systematically quantify how structural disorder, off-stoichiometry, and strain govern Li-ion transport in LiF. Using a machine-learning potential to enable extensive molecular-dynamics sampling, we compare crystalline and amorphous LiF, Li 0.95 F, and LiF 0.95 , and evaluate the impact of small homogeneous deformations. Defect-free crystalline LiF is effectively ion-blocking at 300−500 K, whereas amorphization generates free-volume–assisted percolation pathways that facilitates Li-ion diffusion. At elevated temperatures, thermodynamically driven crystallization disrupts these pathways, leading to non-Arrhenius behavior. In crystalline phases, Li deficiency activates vacancy-mediated diffusion, while in amorphous LiF, transport is governed primarily by network connectivity. Strain is found to have only a marginal effect on Li mobility in both crystalline and amorphous structures.

Batteries↗

Pinning dislocations in colloidal crystals with active particles that seek stacking faults

There is growing interest in functional, adaptive devices built from colloidal subunits of micron size or smaller. A colloidal material with dynamic mechanical properties could facilitate such microrobotic machines. Here we study via computer simulation how active interstitial particles in small quantities can be used to modify the bulk mechanical properties of a colloidal crystal. Passive interstitial particles are known to pin dislocations in metals, thereby increasing resistance to plastic deformation. We extend this tactic by employing anisotropic active interstitials that travel super-diffusively and bind strongly to stacking faults associated with partial dislocations. We find that: (1) interstitials that are effective at reducing plasticity compromise between strong binding to stacking faults and high mobility in the crystal bulk. (2) Reorientation of active interstitials in the crystal depends upon rotational transitions between high-symmetry crystal directions. (3) The addition of certain active interstitial shapes at concentrations as low as 60 per million host particles (0.006%) can create a shear threshold for dislocation migration. This work demonstrates how active materials in a dense matrix can locally sense their environment and lead to bulk property changes.

36 MATERIALS SCIENCE↗

Analysis of Size Effects during Martensitic Transitions in Epitaxial Films and Microparticles of the Ni–Mn–Sn Alloy

Abstract—: The available literature data on the effect of the thickness of the epitaxial films and size of the crystallites in the microparticles of the Ni–Mn–Sn alloy powder on the parameters of the martensitic transitions in this alloy have been analyzed in the framework of the theory of defuse martensitic transitions (DMT), which is based on thermodynamic and kinetic ratios. The purpose of the analysis is to establish an explicit (functional) dependence of these parameters on the film thickness D and the size of nanocrystallites d in the film and microparticles. The results of the R-diffraction analysis show that internal elastic microdeformations and stresses arise in films and microparticles of the alloy powder due to the coherent coupling of the epitaxial film with a solid substrate or as a result of severe plastic deformation of the alloy during its grinding in a ball mill. The analysis shows that the local microstresses significantly affect the type of the dependence of the temperature interval (diffusity) ΔT of the transition on the size factors D or d. In the absence of microstresses, these dependences have the form ΔT ∼ 1/D{sup 2} or ΔT ∼ 1/d{sup 2}. In the presence of microstresses, the temperature transition interval depends on the size factors as ΔT ∼ 1/D or ΔT ∼ 1/d.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sierra/SolidMechanics 4.58 User's Guide

Sierra / SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, an d a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra / SM.

36 MATERIALS SCIENCE↗

Equilibrium shape and location of a liquid drop acoustically positioned in a resonant rectangular chamber

The effect of a standing wave field in a rectangular chamber on the shape and location of an acoustically positioned drop or bubble is calculated. The sample deformation and equilibrium position are obtained from an analysis of the spherical harmonic projections of the total surface stress tensor. The method of calculation relies on the assumed condition that the sample is only slightly distorted from a spherical form. The equilibrium location of a levitated drop is combined with a formula introduced by Hasegawa (1979) to calcualte the ka dependence of the radiation force function. The present theory is valid for large as well as small ka values. Calculations in the small ka limit agree with previous theories and experimental results. Examples are presented for nonplane-wave modes as well as plane-wave rectangular modes.

Jackson, H. W.↗

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗

The Exact Solution for Linear Thermoelastic Axisymmetric Deformations of Generally Laminated Circular Cylindrical Shells

A detailed exact solution is presented for laminated-composite circular cylinders with general wall construction and that undergo axisymmetric deformations. The overall solution is formulated in a general, systematic way and is based on the solution of a single fourth-order, nonhomogeneous ordinary differential equation with constant coefficients in which the radial displacement is the dependent variable. Moreover, the effects of general anisotropy are included and positive-definiteness of the strain energy is used to define uniquely the form of the basis functions spanning the solution space of the ordinary differential equation. Loading conditions are considered that include axisymmetric edge loads, surface tractions, and temperature fields. Likewise, all possible axisymmetric boundary conditions are considered. Results are presented for five examples that demonstrate a wide range of behavior for specially orthotropic and fully anisotropic cylinders.

Nemeth, Michael P.↗

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗