Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiffness”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tuning crosslinking of hybrid preceramic polymers in vat photopolymerization toward controlled ceramic yields

Control of preceramic polymer crosslinking for UV-curable processing is essential for fine 3D printing with high ceramic conversion for sustainable polymer-derived ceramics (PDC) engineering. While various factors influencing ceramic yield have been studied, the systematic exploration of the relationship between crosslinking and ceramic yield, especially when crosslinking increases volatile elements, remains open for further investigation. This study addresses this gap by utilizing vat photopolymerization (VP) additive manufacturing (AM) as a versatile platform for controlling preceramic crosslinking and ceramic yield. By rationally designing and tuning the photochemical crosslinking through digital light processing (DLP), we demonstrate that the ceramic yield can be enhanced from 64% to over 86%, even with added volatile elements. We reveal that the post-pyrolysis ceramic yield can be closely correlated with the pre-pyrolysis crosslinking of the preceramic network represented by its stiffness. This correlation thereby suggests a fast, energy-efficient, non-destructive methodology to predict and improve ceramic yield. Combined with these findings, our one-pot thiol-ene hybridization of polycarbosilane and polycarbosiloxane via DLP offers an exemplary method to generate hybrid preceramic polymers with tailored material properties toward target applications, and potentially even higher ceramic yields. This study thus contributes to achieving better resource- and energy-efficient preceramic polymer and PDC processing routes toward sustainable, advanced organic–inorganic materials manufacturing.

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

Programmable hydrogels by combining persistent and transient dynamic bonds

Out-of-equilibrium chemistry is currently being applied to polymer systems to mimic the autonomous behavior of biological materials. In this study, hydrogels with self-healing properties were developed that combine persistent crosslinks from dynamic metal-ligand coordination with transient crosslinks from dynamic anhydride bonds. Polymers containing terpyridine ligands and carboxylic acid groups were synthesized and crosslinked with divalent metal ions (Fe 2+ , Ni 2+ , Co 2+ , Zn 2+ , Cu 2+ ). The coordination bonds from terpyridine–metal coordination impart persistent stability, while transient anhydrides formed on treatment with 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide hydrochloride (EDC) temporarily increase crosslink density. Distinct behaviors are observed depending on the choice of metal, with Ni 2+ forming robust, stable networks; Zn 2+ creating moderately dynamic gels; and Cu 2+ yielding highly dynamic, soft materials. Treatment with EDC increased storage moduli significantly, with transient effects lasting up to 280 min depending on the metal ion. Self-healing experiments demonstrated independent contributions from metal coordination and transient anhydrides, enabling recovery of stress and strain under varying conditions. Additionally, complex and reversible 2D stiffness patterns were generated by spatially controlled EDC treatment of Zn 2+ and Cu 2+ hydrogel films, demonstrating programmability and reusability.

Rajawasam, Chamoni W. H. [Miami University, Oxford↗
DOE OSTI2022

Ablating Lgr5-expressing prostatic stromal cells activates the ERK-mediated mechanosensory signaling and disrupts prostate tissue homeostasis

The heterogeneity of the prostate stromal cells is widely appreciated but the functional implication remains incompletely understood. Using genetic lineage tracing and light-sheet imaging, we show that some fibroblast cells near the junction of the mouse proximal prostatic ducts and prostatic urethra highly express Lgr5. Genetic ablation of these anatomically restricted stromal cells, but not nonselective ablation of prostatic stromal cells, rapidly induces prostate epithelial turnover and dedifferentiation that are reversed following spontaneous restoration of the Lgr5 + stromal cells. RNA-Seq analysis implies that ablating the Lgr5 + stromal cells activates a mechanosensory response. We show that ablating the Lgr5 + stromal cells impairs the control of prostatic ductal outlet, increases prostate tissue stiffness, and activates the mitogen-activated protein kinase (MAPK). Suppressing MAPK overrides the increased epithelial proliferation. This study demonstrates that the Lgr5 + stromal cells regulate tissue homeostasis in a long-distance manner by maintaining anatomic integrity and implies that the cells near the transitional regions between organs likely control organ homeostasis by sustaining a balanced mechanoforce.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Overall assessment of alkali-silica reaction affected recycled concrete aggregate mixtures derived from construction and demolition waste

A large amount of research has been conducted on the reuse of construction and demolition waste (CDW) in concrete. Yet, CDW may display pre-existing damage which raises concerns on its use. In this work, coarse recycled concrete aggregate (RCA) is reclaimed from distinct members of an alkali-silica reaction (ASR) affected overpass. RCA mixtures incorporating 50 and 100% replacement are then manufactured and stored in conditions enabling further ASR development. Mechanical (Stiffness Damage Test- SDT) and microscopic (Damage Rating Index- DRI) analyses are conducted at a fixed “secondary” induced expansion. Results indicate that the overall damage of ASR-affected RCA mixtures is different from conventional concrete and depends upon the “past” RCA condition. Furthermore, the DRI can capture the “past” and “secondary” expansion while the SDT simply detects the “secondary” distress. Finally, an adapted DRI version is proposed to better evaluate and distinguish “past” and “secondary” damage in ASR-affected recycled concrete.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modelling the development of capillary pressure in freshly 3D-printed concrete elements

3D concrete printing is a promising technology recently developed to automate construction. Since no formwork is used in this technology to support and protect fresh concrete, there are two aspects which considerably accelerate the development of capillary pressure in 3D-printed concrete in comparison to conventionally placed concrete: i) high stiffness of 3D-printed needed to provide sufficient buildability, and ii) very early and fast evaporation of pore water. Accelerated development of capillary pressure may lead to severe plastic shrinkage cracking in 3D-printed elements and, hence, need to be mitigated. This investigation aims at providing a poromechanical model for capillary pressure development in 3D-printed elements. To simulate the development of capillary pressure and plastic shrinkage, environmental factors, material properties, and element geometry need to be considered as a whole. The model inputs – coefficient of permeability, static bulk modulus, air entry pressure and chemical shrinkage rate – were determined experimentally. The model was validated for two fine-grained concretes. Both 3D-printed materials yielded faster capillary pressure increase in comparison to cast concrete, while partial substitution of cement with silica fume further accelerated the capillary pressure development. Furthermore, due to the lower permeability of the mixture containing silica fume, the gradient of capillary pressure between 3D-printed layers increased, as did the gradient of plastic shrinkage.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development and characteristics of ultra high-performance lightweight cementitious composites (UHP-LCCs)

Highlights: • An ultra high-performance lightweight cementitious composite (UHP-LCC) was designed. • The UHP-LCCs had ultra high strength (>120 MPa) and very low density (3 ). • The UHP-LCCs had high structural efficiency, good functional properties and superior durability. • The strategies and mechanisms on achieving UHP-LCCs with excellent performance were proposed. High strength and light weight are two recent opposite development trends of concrete. This study proposed a design concept of an ultra high-performance lightweight cementitious composite (UHP-LCC), which had a compressive strength of higher than 120 MPa and an air-dried density down to around 1800 kg/m{sup 3}. The UHP-LCCs were innovatively developed by incorporating micro-sized hollow particles with a high strength shell (hollow glass microspheres, HGM) into an ultra-high performance cementitious composite (UHPC). The roles of HGM in the UHP-LCCs were investigated by evaluating the reactivity of the HGM and the mechanisms on achieving the excellent mechanical properties, low density and superior durability were revealed. The Chapelle test results showed that the HGM exhibited some pozzolanic reactivity, which facilitated the reaction between the shell of HGM and the alkali hydration products of the paste matrix. This chemical reaction was conducive to improving the HGM-paste interface and enhancing the mechanical properties. With the use of microspheres with a high stiff shell, the fundamental properties of the UHP-LCCs including thermal insulation, sound absorption, resistance to water ingress and electrical resistivity were improved significantly. The strategies for preparing the UHP-LCCs with high structural efficiency and great performance were proposed. The results of this study provide a new approach for designing and producing a lightweight UHPC, which would be a promising material for long-span structures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Numerical investigation of thermomechanical response of multiscale porous Ultra-High Temperature Ceramics

Recent advances in Ultra-High Temperature Ceramics (UHTC) manufacturing have permitted the development of multiscale porous UHTC microstructures. Within the target application of hypersonic vehicles, dense UHTCs are suitable for thermal protection on leading edges, whereas porous UHTCs may find a role in providing thermally insulated interfaces for temperature-sensitive interior components. Designing vehicles incorporating porous UHTCs requires a characterization of their thermomechanical properties across the full range of expected operating temperatures spanning –20 °C–2500 °C. This research represents a preliminary study in performing this characterization. Several numerical experiments are performed using a coupled thermomechanical implementation of the Material Point Method to determine the temperature dependence of effective material properties both with and without damage. Furthermore, complex time-dependent boundary conditions derived from known hypersonic flight profiles are simulated in order to probe the various couplings between deformation, damage and heat transfer. It is shown that the model reveals the importance of micro-buckling in determining effective material stiffness and thermal conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On differentiable local bounds preserving stabilization for Euler equations

This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.

42 ENGINEERING↗
DOE OSTI2020

Simultaneous material, shape and topology optimization

Using three design fields we develop an optimization environment that can simultaneously optimize material, shape and topology. We use the implicit representation of the boundaries with level-set functions that define the shape and topology. Differentiable R-functions allow us to combine these shapes and topology descriptions with Boolean operations. Additionally, we incorporate design dependent-stiffness materials with another design field. Notably, this framework accommodates design dependent loads, has the ability to introduce holes, and ensures the satisfaction of optimality criteria. It builds upon the fictitious domain, ersatz material, material interpolation and level-set methods. Additionally, it also borrows from parameterized density-based topology optimization methods. Since analytical sensitivities can be computed, we use efficient nonlinear programming algorithms to update the design instead of the Hamilton–Jacobi’s scheme of level-set methods. We illustrate the features of our framework by designing a cantilever beam with octet truss microlattice, a dam with design-dependent loads, and a composite clevis plate.

42 ENGINEERING↗
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 OSTI2021

Galerkin formulations of isogeometric shell analysis: Alleviating locking with Greville quadratures and higher-order elements

Here we propose new quadrature schemes that asymptotically require only four in-plane points for Reissner–Mindlin shell elements and nine in-plane points for Kirchhoff–Love shell elements in B-spline and NURBS-based isogeometric shell analysis, independent of the polynomial degree p of the elements. The quadrature points are Greville abscissae associated with pth-order B-spline basis functions whose continuities depend on the specific Galerkin formulations, and the quadrature weights are calculated by solving a linear moment fitting problem in each parametric direction. The proposed shell element formulations are shown through numerical studies to be rank sufficient and to be free of spurious modes. The studies reveal comparable accuracy, in terms of both displacement and stress, compared with fully integrated spline-based shell elements, while at the same time reducing storage and computational cost associated with forming element stiffness and mass matrices and force vectors. The high accuracy with low computational cost makes the proposed quadratures along with higher-order spline bases, in particular polynomial orders, p = 5 and 6, good choices for alleviating membrane and shear locking in shells.

42 ENGINEERING↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Isogeometric Bézier dual mortaring: The Kirchhoff–Love shell problem

Here in this paper we develop an isogeometric Bézier dual mortar method for coupling multi-patch Kirchhoff–Love shell structures. The proposed approach weakly enforces the continuity of the solution at patch interfaces through a dual mortar method and can be applied to both conforming and non-conforming discretizations. As the employed dual basis functions have local supports and satisfy the biorthogonality property, the resulting stiffness matrix is sparse. In addition, the coupling accuracy is optimal because the dual basis possesses the polynomial reproduction property. We also formulate the continuity constraints through the Rodrigues’ rotation operator which gives a unified framework for coupling patches that are intersected with G 1 continuity as well as patches that meet at a kink. Several linear and nonlinear examples demonstrated the performance and robustness of the proposed coupling techniques.

42 ENGINEERING↗
DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗
DOE OSTI2022

Stress-constrained topology optimization of lattice-like structures using component-wise reduced order models

We report lattice-like structures can provide a combination of high stiffness with light weight that is useful in many applications, but a resolved finite element mesh of such structures results in a computationally expensive discretization. This computational expense may be particularly burdensome in many-query applications, such as optimization. We develop a stress-constrained topology optimization method for lattice-like structures that uses component-wise reduced order models as a cheap surrogate, providing accurate computation of stress fields while greatly reducing run time relative to a full order model. We demonstrate the ability of our method to produce large reductions in mass while respecting a constraint on the maximum stress in a pair of test problems. The ROM methodology provides a speedup of about 150x in forward solves compared to full order static condensation and provides a relative error of less than 5% in the relaxed stress.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
DOE OSTI2020

Slow manifold reduction for plasma science

The classical Chapman–Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma physics audience in mind. In particular we illustrate (a) how the slow manifold concept may be used to understand breakdown of a reduced model over sufficiently-long time intervals, and (b) how a discrete-time analogue of slow manifold theory provides a useful framework for developing implicit integrators for temporally-stiff plasma models. For readers with more advanced mathematical training we also use slow manifold reduction to explain the phenomenon of inheritance of Hamiltonian structure in dissipation-free reduced plasma models. Various facets of the theory are illustrated in the context of the Abraham–Lorentz model of a single charged particle experiencing its own radiation drag. As a culminating example we derive the slow manifold underlying kinetic quasineutral plasma dynamics up to first-order in perturbation theory. This first-order result incorporates several physical effects associated with small deviations from exact charge neutrality that lead to slow drift away from predictions based on the leading-order approximation n e = Z i n i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Strategies for design and synthesis of porous liquids toward carbon capture and separation

Porous liquids (PLs) represent a promising category of sorbents in carbon capture and separation capable of integrating the advantages of flowing liquid and porous solid systems. Well-defined pores were engineered into liquid sorbents via liquifying molecules with stiff interior voids, dissolving rigid porous hosts in flowing liquids, or dispersing porous frameworks in high steric hindrance solvents, producing type I, II, or III PLs, respectively. Unique features of PLs have triggered broad interest in exploring their applications in carbon capture and separation, in which diverse design strategies, synthesis approaches, and enhanced performance have been reported. Here, in this minireview, recent progress in the design, synthesis, and structural engineering of PLs and efforts towards the optimization of their carbon capture and separation behavior will be summarized, including the comparison between PLs with varied types. Porosity engineering into liquid sorbents provides opportunities to resolve challenging issues in conventional sorption and separation systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗