Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normalization”

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

Seawater sea-sand engineered/strain-hardening cementitious composites (ECC/SHCC): Assessment and modeling of crack characteristics

Highlights: • A probabilistic model was proposed to model the crack width evolution of SS-ECC at different strain levels. • The Weibull distribution fit the crack width distribution better than the log-normal distribution. • A 5-D representation was proposed to assess the SS-ECC by considering both cracking and mechanical performance. • Larger sea-sand size and lower fiber dosage led to larger crack widths in SS-ECC under the same tensile strain. • 18-mm PE fiber led to larger crack widths in SS-ECC at strain >2%, due to a large fraction of fiber rupture. Seawater sea-sand Engineered Cementitious Composites (SS-ECC) is a new version of ECC for marine constructions facing the scarcity of freshwater and river/manufactured sand. This study aims to assess and model the crack characteristics of SS-ECC, which are critical for its applications with non-corrosive reinforcements. The influence of sea-sand size, fiber length and fiber dosage on the crack characteristics of SS-ECC was explored. A five-dimensional representation was proposed to assess the overall performance of SS-ECC, by comprehensively considering both the crack characteristics (i.e., crack width and its variation) and the mechanical properties (i.e., compressive and tensile properties). A probabilistic model was also proposed to describe the stochastic nature and evolution of crack width, and it can be used to estimate the critical tensile strain on SS-ECC for a given crack-width limit and cumulative probability. The findings and proposed methods can facilitate the design of SS-ECC in marine and coastal structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of a novel starch-based temperature rise inhibitor on cement hydration and microstructure development: The second peak study

Temperature rise inhibitor (TRI) modifies the exothermic process of cement hydration at early age and mitigates the buildup of cumulative heat, which sheds new light on solving the thermal cracking issue. After the main hydration peak, a broader “second peak” appears and brings the cumulative heat back to normal level within 7 days (for additions less than 0.15%). This study shows that the “second peak” can be attributed, almost exclusively, to C{sub 3}S hydration while triggered by C{sub 3}A hydration. In a pure C{sub 3}S system, the heat flow is depressed by TRI during the main hydration peak but never recovers after. In a monophase system composed of C{sub 3}S, C{sub 3}A and gypsum or a real cement paste, the ongoing precipitation of ettringite helps to pump out TRI polymer from the pore solution and eliminate the restriction on C-S-H nucleation. Afterwards, the precipitation of C-S-H accelerates again to form the “second peak”.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strain hardening magnesium-silicate-hydrate composites (SHMSHC) reinforced with short and randomly oriented polyvinyl alcohol microfibers

Magnesium-silicate-hydrate (M-S-H) cement has emerged as an alternative binder with potentially lower energy requirements and emissions. Due to the intrinsically low pH values of its matrix, however, normal steel reinforcement is not appropriate for M-S-H system. To toughen the matrix and to overcome the brittle nature of the material, a new strain hardening magnesium-silicate-hydrate composite (SHMSHC) is developed for the first time by incorporating 2 vol% short and randomly oriented polyvinyl alcohol (PVA) microfibers. The resulting SHMSHC exhibits significant strain hardening with a tensile strain capacity of more than 3%, a compressive strength beyond 50 MPa and a tensile strength of around 3 MPa. Remarkably, saturated multiple cracking with a tight crack width less than 10 μm invisible to the human naked eyes was observed. The newly developed SHMSHC addresses the fundamental challenge of M-S-H binder and is expected to widen possible application areas of M-S-H system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Comparative analysis of temperature dependent properties of commercial nuclear fuel pellet and surrogates undergoing cracking: A review

A variety of normal and accident scenarios can generate thermal stresses enough to cause cracking in light water reactor (LWR) fuel. To better understand this behavior, cracking experiments have been carried out to induce thermal gradients in fuel pellets via induction heating and direct resistance heating. In this study, Ceria (CeO 2 ) and yttria stabilized zirconia (YSZ) pellets were chosen as a surrogate fuel material for uranium dioxide (UO 2 ). Induction heating was done using copper coils and molybdenum susceptors which heated the surrogates to a threshold temperature that is sufficiently high for the fuel material to conduct current. Thereafter, direct resistance heating was achieved by passing current through the specimen using a DC power supply to introduce volumetric heating to replicate reactor operating conditions. Simultaneous real-time dual imaging of the ceria pellet surface has been designed using optical and infra-red camera system to capture images of cracks and full-field temperature gradients on pellet. It was observed that YSZ pellets need a much higher voltage-current for volumetric heating and do not exhibit the cracking pattern as ceria and UO 2 , hence YSZ was ruled out for further consideration. Ceria was found to be an appropriate surrogate for UO 2 . Furthermore, the experimental set up and test conditions demand the need to primarily understand the thermo-physical, mechanical and optical properties of UO 2 and the surrogates. Thus, it is essential to have an in-depth knowledge about the various temperature dependent properties of UO 2 and CeO 2 for establishing a comparative analysis between the surrogate and UO 2 for conducting experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bridging knowledge gaps in liquid- and vapor-fed CO 2 electrolysis through active electrode area

Increased use of gas diffusion electrodes for CO 2 electroreduction widens the experimental phase space that was previously inaccessible using foil electrodes, raising fundamental questions over the impacts of key variables that translate between liquid- and vapor-fed CO 2 electrolysis systems. This work focuses on studying the interplay of current-potential profiles and electrochemically active surface area (ECSA) by implementing a Cu nanoflower catalyst morphology. The results show decreased overpotentials as much as 460 and 174 mV for foil and gas diffusion electrodes, respectively, while maintaining or improving multi-carbon product current density. Furthermore, these overpotential shifts and product activities normalized by ECSA lead to current-potential relationships akin to those of the Tafel description, which are found through a continuum model to be useful for describing the roughness dependence for both liquid- and vapor-fed systems. This analysis establishes a holistic approach for establishing catalyst design criteria to improve materials development for CO 2 .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Metal oxide-promoted calcium cuprate catalysts for diol oxidative dehydrocyclization to lactones

Here, this work investigates structure-function relationships in electronically tunable, redox-active, basic Cu-Ca mixed metal oxide catalysts for oxidative dehydrocyclization of liquid diols to lactones. Compositional screening identified Ni 2+ and Zn 2+ as effective promoters that increase the surface Cu 2+ population by ∼1.7× and Cu-normalized activity for liquid 1,4-butanediol conversion to γ-butyrolactone by ∼3–4×. In situ Raman spectroscopy, in situ X-ray absorption spectroscopy (XAS), in situ diffuse-reflectance Fourier transform infrared spectroscopy (DRIFTS), ex situ X-ray diffraction (XRD), and H 2 -temperature-programmed reduction (H 2 -TPR) show that Ni 2+ or Zn 2+ incorporation promotes the formation of Ca 0.82 Cu 1.00 O 2 nanoparticles under mild calcination conditions. This cuprate phase features stronger and shorter Cu–O bonds (1.90 Å) than inactive bulk CuO (1.95 Å) and square-planar Cu 2+ O 4 sites with enhanced d z2 electrophilicity, strengthening alkoxy adsorption. Pyridine-DRIFTS confirms the purely basic nature of the catalyst surface, while methanol-DRIFTS indicates Cu 2+ surface enrichment with Ni or Zn promotion, where Cu–O(Ca)–Cu sites can exist as amorphous domains or a truncation layer on crystalline nanoparticles.

09 BIOMASS FUELS↗
DOE OSTI2020

Equilibrium carbon isotope fractionation factors of hydrocarbons: Semi-empirical force-field method

Here, we have calculated the reduced partition function ratios for carbon isotopes (β-factor) of 67 hydrocarbons (alkanes, alkenes, alkynes, cycloalkanes, and aromatics), including their 267 single-substituted isotopomers. The calculations were performed using the harmonic oscillator – rigid rotator model and the Urey / Bigeleisen-Mayer approach. Normal frequencies of molecular vibrations of the isotopologues were calculated from the molecular structures data and valence field force constants stored in the Light-handled Elucidation of Vibrations (LEV) database elaborated by Dr. Gribov and his colleagues in Russia. The LEV database is constructed from experimental spectroscopic and structural data by solving inverse problems of molecular vibrations. The LEV is internally consistent and the most comprehensive database available to date. The β-factor were calculated in the temperature range of 200–800 K with a 10 K step. Our calculations predict that the β-factors increase with an increasing number of C atoms within the same groups of hydrocarbons (e.g., alkanes). Our calculations also show a general descending order of 13 C enrichments among the different groups of hydrocarbons: cycloalkanes, aromatics, alkenes (double bonds) and isoalkanes, alkanes, alkynes (triple bonds). Position-specific, intramolecular isotope effects within hydrocarbons are determined by the β-factors of C in the different functional groups in the order: quaternary (C), tertiary (methine, CH), secondary (methylene, CH 2 ), primary (methyl, CH 3 ) and double bond (C=C), saturated bond (C - C), triple bond (HC≡). Our calculations on bulk and position-specific carbon isotope β-factors of the hydrocarbons, which are generally consistent with very limited ab initio calculations in the literature, are internally consistent and the most comprehensive to date for future applications to position-specific isotope geochemistry of hydrocarbons.

58 GEOSCIENCES↗
DOE OSTI2024

Fe Kα XANES, Fe Kβ HERFD XANES and EPMA flank method determinations of the oxidation state of Fe in garnet

The ferric to total iron ratios (Fe 3+ /ΣFe) of garnets can be paired with thermodynamic mineral activity models to quantify the oxygen fugacity of garnet-bearing rocks. However, techniques with a high analytical and spatial resolution are necessary to distinguish differences in garnet Fe 3+ /ΣFe ratios at the percent level and to accurately measure garnets that are zoned or contain inclusions. We acquired conventional Fe Kα and high-resolution energy fluorescence detection (HERFD) Fe Kβ X-ray absorption near edge structure (XANES) spectra and electron microprobe flank method analyses on a suite of 27 peridotitic and eclogitic garnets with Fe 3+ /ΣFe ratios previously determined by Mössbauer spectroscopy to evaluate the precision of each technique. We examined variations in the energy and intensity of three XANES spectral features as a function of Fe 3+ /ΣFe ratios: 1) the intensity ratio of two-post edge features (I-ratio; Fe Kα only); 2) the energy of the Fe edge at 90% normalized intensity (E 0.9 ; Fe Kα only) and 3) the pre-edge centroid energy (Fe Kα and HERFD Fe Kβ). In accordance with previous work, we find the energies of garnet pre-edge centroids are relatively insensitive to Fe 3+ /ΣFe ratios. The I-ratios of peridotitic and eclogitic garnets are offset from each other at low Fe 3+ /ΣFe ratios (≤0.13); I-ratio garnet XANES calibrations are composition-specific. The E 0.9 feature is independent of garnet major element composition in spectra that have been corrected for the effects of self-absorption. We produce two Fe Kα garnet XANES calibrations based on variations in the E0.9 feature; one calibration with all garnet reference materials included (Fe 3+ /ΣFe up to 1.0; “all garnet calibration”) and another calibration specific to garnets with low Fe 3+ /ΣFe ratios (“low ferric calibration”). Fe 3+ /ΣFe ratios calculated from the mean of up to 25 flank method measurements on eight garnet reference materials fall within 4% absolute of a one-to-one correlation with Fe 3+ /ΣFe ratios measured by Mössbauer. The standard error of the mean Fe 3+ /ΣFe ratio calculated from flank method approaches the Mössbauer-determined Fe 3+ /ΣFe ratio within estimated error (3%) after three analyses. Flank method precision is enhanced at higher beam current; however, the precision of the flank method does not approach the precision of XANES under any microprobe analytical condition tested here. Garnet reference materials detailed here are available by request to the Smithsonian Institution.

58 GEOSCIENCES↗
DOE OSTI2022

Weighted nodal domain averages of eigenstates for quantum Monte Carlo and beyond

In this report we study the nodal properties of many-body eigenstates of stationary Schrödinger equation that affect the accuracy of real-space quantum Monte Carlo calculations. In particular, we introduce weighted nodal domain averages that provide a new probe of nodal surfaces beyond the usual expectations. Particular choices for the weight function reveal, for example, that the difference between two arbitrary fermionic eigenvalues is given by the nodal hypersurface integrals normalized by overlaps with the bosonic ground state of the given Hamiltonian. Noninteracting and fully interacting Be atom with corresponding almost exact and approximate wave functions are used to illustrate several aspects of these concepts. Variational formulations that employ different weights are proposed for prospective improvement of nodes in variational and fixed-node diffusion Monte Carlo calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Metallurgical alloy approach to two-dimensional supramolecular materials

Alloys are normally mixtures of metallic elements used to search for superior functional properties compared with their individual components. Although a quintessential approach to better materials in metallurgy, the concept of alloying has been expanded to other inorganic materials and covalent polymers. In this work, we report on alloy formation of small organic molecules structured as amphiphilic chromophores in water. Using X-ray scattering, microscopy, and optical spectroscopy, we unambiguously demonstrate the formation of supramolecular alloys with changes in physical, photophysical, and mechanical properties. Additionally, we identify single-crystalline, two-dimensional supramolecular polymers containing multiple molecular species, analogous to intermetallic compounds. The mechanism of alloy formation suggests the possibility of controlled growth in these supramolecular assemblies to create emergent functional structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A comparative TEA of a two-step process using chemical solvents for producing an ultra-sweet natural gas

Here, a comprehensive Techno-Economic Analysis (TEA) was performed to evaluate the economic feasibility of a novel two-step process (TSP) developed in Aspen Plus V12.1 to desulfurize and decarbonize a raw natural gas containing (2 mol% H 2 S and 5 mol% CO 2 ) into an ultra-sweet natural gas containing (1.72 ppmv H 2 S and 4.19 ppmv CO 2 ). The raw natural gas flow rate used in the TSP was 117.74 kg/s at 60 °C and 50 bar. The TSP combines an H 2 S desulfurization step using potassium carbonate (K 2 CO 3 ) and a CO 2 capture step using 3 different chemical solvents, monoethanolamine (MEA), sodium glycinate (SGS), and potassium glycinate (PGS). Both steps employ fixed-bed absorbers packed with Mellapak 250Y structured packing. The hydraulics and mass transfer characteristics for the TSP were calculated, indicating normal operation with higher gas-side (k G ) than liquid-side (k L ) mass transfer coefficients. The TEA of TSP indicated that PGS had the most promising economic feasibility among the 3 solvents as it exhibited the lowest Levelized Cost of CO 2 capture (LCOC) of $\$$47.54/ton.CO 2 at a Capital Expenditure (CAPEX) of $\$$24.98 million, and an Operating Expenditure (OPEX) of $\$$12.20 million/year. Also, the TSP could produce one MMSCF of ultra-sweet natural gas at a total cost of $\$$339.55.

03 NATURAL GAS↗
DOE OSTI2022

Comparing and contrasting In-Vial and full-scale systems for sparging volatile analytes

Here, in-vial sparging was demonstrated as an effective, practical alternative to a full-scale sparging system for supporting the analysis of volatile constituents. Using elemental mercury (Hg 0 ) and toluene as representative purgeable analytes, the mass removal for various sparge configurations was measured and a reduced order model was developed and validated. In the primary experiments, Hg 0 in the sparge gas was trapped on activated carbon or gold, thermally desorbed, and quantified using atomic absorption or atomic fluorescence spectroscopy. Toluene experiments using the same in-vial sparge apparatus and sparge parameters were performed to demonstrate the applicability of the reduced order model to a broad range of compounds. Toluene removal was tracked by measuring the remaining toluene in sparged aliquots using Ultraviolet-visible (UV–Vis) spectroscopy. For the sparging, flow rates varied from 25 to 75 mL/min for periods from 0 to 30 min. Sparge performance, mass removal as a function of time, and sparge gas volume were measured for both in-vial and full-scale systems. A model based on dimensionless Henry's Law coefficient, normalized sparge gas volume, and fractional extent of equilibrium matched the experimental data for both compounds and provides a practical tool for future applications. For the conditions tested in this study, the calibrated model indicated that the sparge gas in the in-vial system reached approximately 33% of its equilibrium value before exiting the water surface, while a full-scale system reached approximately 100%. The tests validated the quality, reproducibility, and predictability of sparging performance for both full scale and in-vial sparge systems. Related factors such as waste generation, worker risk, and labor were also assessed. Full scale sparge systems provide the advantage of lower detection levels due to larger sample volume, while the in-vial sparge systems provide advantages for most other factors; including automatability, reducing secondary wastes, lessening the need to clean and check the sparge apparatus, and lowering labor and costs. The data and associated reduced order model support continued development and deployment of in-vial sparge platforms as a practical option for analysis of purgeable analytes such as volatile organic compounds and volatile metals/organometallics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Adsorption of terbium (III) on DGA and LN resins: Thermodynamics, isotherms, and kinetics

Two commercially available extraction chromatography (EXC) resins containing N,N,N’,N’-tetra-n-octyldiglycolamide (DGA Resin, Normal, 50 – 100 μm) and Bis(2-ethylhexyl) phosphate (LN Resin, 100 – 150 μm) were used as adsorbents to study fundamental adsorption properties such as thermodynamic values, equilibrium isotherms, and kinetic uptake models for terbium(III) adsorption. Weight distribution ratios (D w ) for terbium on DGA and LN resins were measured using a [ 160 Tb]Tb 3+ radiometric tracer in nitric acid as a function of acidity, temperature, initial analyte concentration, and equilibrium time. The D w values showed increasing binding affinity for DGA resin at high nitric acid concentrations and decreasing binding affinity for LN resins. Thermodynamic studies for DGA and LN resins revealed that the Gibbs free energy (ΔG) increased consistently with temperature. To model equilibrium data, increasingly higher parameter equilibrium isotherm models (Henry (1) < Langmuir, Freundlich (2) < Redlich-Peterson (3) < Fritz-Schluender (4)) were compared on their root mean squared errors (RMSE) and adjusted determination coefficients to determine the most applicable model. In all cases, the empirical four-parameter Fritz-Schluender isotherm demonstrated a superior fit. Similar comparisons for reaction-based kinetic models (Pseudo-first-order < Pseudo-second-order < Pseudo-n-order) revealed that the higher-order PNO model yielded a superior fit of kinetic data for both resins. Furthermore, in some cases, adsorption isotherms and kinetic models could also be modeled by a lower-order model with minimal change in error parameters. Weber-Morris plots revealed that two linear sections are observed for each resin, where the first linear segment is attributed to fast (film diffusion) adsorption of terbium, followed by slower intraparticle diffusion of terbium through the pores as the rate-limiting step. Based on the Weber-Morris plot, both film and intraparticle diffusion are involved in controlling the kinetic rate of adsorption for DGA and LN resins.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Development of a SnO 2 -based 44 Ti/ 44 Sc generator for medical applications

Towards application of 44 Sc for diagnostic nuclear medicine, a 44 Ti/ 44 Sc generator based on an inorganic resin has been evaluated. Unlike other radionuclide generators used for medical applications, the long-term retention of the parent 44 Ti is vital due to its long half life. In this work, tin dioxide (SnO 2 ), a robust inorganic-based resin, has been synthesized and used as the stationary phase for a 44 Ti/ 44 Sc generator. The sorption behavior of 44 Ti/ 44 Sc was tested on SnO 2 with varying acids, concentrations, and times. Preliminary batch study results showed >88 % 44 Ti retention to the resin at lower acid concentrations (0.05 M HNO 3 and 0.05 M HCl). A pilot generator was evaluated for a year, demonstrating 85.3 ± 2.8 % 44 Sc elution yields and 0.71 ± 0.14 % 44 Ti breakthrough in 5 M HNO 3 . Based on capacity studies, a 7.4 MBq (200 µCi) upscaled generator system was constructed for further evaluation of the SnO 2 resin stability and the efficacy of the eluted 44 Sc for radiolabeling. 44 Sc could be regularly eluted from this generator in 5 M HNO 3 with an overall average radiochemical yield 84.7 ± 9.5 %. Post-elution processing of the 44 Sc with DGA-normal resin removed all 44 Ti present and allowed for high 44 Sc-DOTA labeling yields of 94.2 ± 0.5 %. Overall, SnO 2 has been shown to be a viable material for a 44 Ti/ 44 Sc generator.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗
DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗
DOE OSTI2022

Imposition of local boundary conditions in peridynamics without a fictitious layer and unphysical stress concentrations

Here, this study introduces a general approach for the imposition of local boundary conditions in non- ordinary state-based peridynamics (NOSB PD) to eliminate the displacement kinks near the boundary without a fictitious layer under quasi-static loading conditions. It identifies the underlying reason for the unphysical displacement kinks. Under an imposed linear displacement field, the NOSB PD equilibrium equation is not satisfied near the boundary due to the unsymmetric horizon of material points. However, the equilibrium equation derived by using the PD differential operator is satisfied at such material points. Therefore, the material body is divided into three regions to satisfy the equilibrium equations and to impose displacement and tractions boundary conditions. This approach does not deviate from the original NOSB PD; however, it provides a simple solution to eliminate the displacement kink near the boundary, which leads to unphysical stress concentrations. Its efficacy is demonstrated by considering elastic rectangular and square plates subjected to various types of boundary conditions leading to homogeneous as well as nonhomogeneous deformations. The creep response of a rectangular plate further proves the robustness of the present approach. Also, a quasi-static crack propagation from a pre-existing crack in a square plate under mode-I, mode-II and mixed-mode loading conditions demonstrates its capability for failure prediction based on the critical stretch criteria. Finally, its applicability for 3D analysis is demonstrated by considering a rectangular prism under applied stretch and normal stress.

42 ENGINEERING↗
DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗