Search NASA⌕ Search

Engineering topics

Chen, Jiun-Shyan

Publications and source records attributed to Chen, Jiun-Shyan.

13 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗
DOE OSTI2024

A Comprehensive Review of Latent Space Dynamics Identification Algorithms for Intrusive and Non-Intrusive Reduced-Order-Modeling

Numerical solvers of partial differential equations (PDEs) have been widely employed for simulating physical systems. However, the computational cost remains a major bottleneck in various scientific and engineering applications, which has motivated the development of reduced-order models (ROMs). Recently, machine-learning-based ROMs have gained significant popularity and are promising for addressing some limitations of traditional ROM methods, especially for advection dominated systems. In this chapter, we focus on a particular framework known as Latent Space Dynamics Identification (LaSDI), which transforms the high-fidelity data, governed by a PDE, to simpler and low-dimensional latent-space data, governed by ordinary differential equations (ODEs). These ODEs can be learned and subsequently interpolated to make ROM predictions. Each building block of LaSDI can be easily modulated depending on the application, which makes the LaSDI framework highly flexible. In particular, we present strategies to enforce the laws of thermodynamics into LaSDI models (tLaSDI), enhance robustness in the presence of noise through the weak form (WLaSDI), select high-fidelity training data efficiently through active learning (gLaSDI, GPLaSDI), and quantify the ROM prediction uncertainty through Gaussian processes (GPLaSDI). We demonstrate the performance of different LaSDI approaches on Burgers equation, a non-linear heat conduction problem, and a plasma physics problem, showing that LaSDI algorithms can achieve relative errors of less than a few percent and up to thousands of times speed-ups.

Computational Engineering, Finance, and Science (c↗
DOE OSTI2023

gLaSDI: Parametric physics-informed greedy latent space dynamics identification

A parametric adaptive physics-informed greedy Latent Space Dynamics Identification (gLaSDI) method is proposed for accurate, efficient, and robust data-driven reduced-order modeling of high-dimensional nonlinear dynamical systems. In the proposed gLaSDI framework, an autoencoder discovers intrinsic nonlinear latent representations of high-dimensional data, while dynamics identification (DI) models capture local latent-space dynamics. Here, an interactive training algorithm is adopted for the autoencoder and local DI models, which enables identification of simple latent-space dynamics and enhances accuracy and efficiency of data-driven reduced-order modeling. To maximize and accelerate the exploration of the parameter space for the optimal model performance, an adaptive greedy sampling algorithm integrated with a physics-informed residual-based error indicator and random-subset evaluation is introduced to search for the optimal training samples on the fly. Further, to exploit local latent-space dynamics captured by the local DI models for an improved modeling accuracy with a minimum number of local DI models in the parameter space, a -nearest neighbor convex interpolation scheme is employed. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including Burgers equations, nonlinear heat conduction, and radial advection. The proposed adaptive greedy sampling outperforms the conventional predefined uniform sampling in terms of accuracy. Compared with the high-fidelity models, gLaSDI achieves 17 to 2,658× speed-up with 1 to 5% relative errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Initial Simulations of Empty Room Collapse and Reconsolidation at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

A MUSCL-SCNI approach for meshfree modeling of shock waves in fluids

Here, a stable and nodally integrated meshfree formulation for modeling shock waves in fluids is developed. The reproducing kernel approximation is employed to discretize the conservation equations for compressible flow, and a flux vector splitting approach is applied to allow proper numerical treatments for the advection and pressure parts, respectively, based on the characteristics of each flux term. To capture the essential shock physics in fluids, including the Rankine–Hugoniot jump conditions and the entropy condition, local Riemann enrichment is introduced under the stabilized conforming nodal integration (SCNI) framework. Meanwhile, numerical instabilities associated with the advection flux are eliminated by adopting a modified upwind scheme. To further enhance accuracy, a MUSCL-type method is introduced in conjunction with an oscillation limiter to avoid Gibbs phenomenon and ensure monotonic piecewise linear reconstruction in the smooth region. The present meshfree formulation is free from tunable artificial parameters and is capable of capturing shock and rarefaction waves without over/undershoots. Finally, several numerical examples are analyzed to demonstrate the effectiveness of the proposed MUSCL-SCNI approach in meshfree modeling of complex shock phenomena, including shock diffraction, shock–vortex interaction, and high energy explosion processes.

42 ENGINEERING↗
DOE OSTI2019

Treatment of near-incompressibility in meshfree and immersed-particle methods

We propose here new projection methods for treating near-incompressibility in small and large deformation elasticity and plasticity within the framework of particle and meshfree methods. Using the $\bar{B}$ and $\bar{F}$ techniques as our point of departure, we develop projection methods for the conforming reproducing kernel method and the immersed-particle or material point-like methods. The methods are based on the projection of the dilatational part of the appropriate measure of deformation onto lower-dimensional approximation spaces, according to the traditional $\bar{B}$ and $\bar{F}$ approaches, but tailored to meshfree and particle methods. The presented numerical examples exhibit reduced stress oscillations and are free of volumetric locking and hourglassing phenomena.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Conforming window functions for meshfree methods

Window functions provide a base for the construction of approximation functions in many meshfree methods. They control the smoothness and extent of the approximation functions and are commonly defined using Euclidean distances which helps eliminate the need for a meshed discretization, simplifying model development for some classes of problems. However, for problems with complicated geometries such as nonconvex or multi-body domains, poor solution accuracy and convergence can occur unless the extents of the window functions, and thus approximation functions, are carefully controlled, often a time consuming or intractable task. Here in this paper, we present a method to provide more control in window function design, allowing efficient and systematic handling of complex geometries. “Conforming” window functions are constructed using Bernstein–Bézier splines defined on local triangulations with constraints imposed to control smoothness. Graph distances are used in conjunction with Euclidean metrics to provide adequate information for shaping the window functions. The conforming window functions are demonstrated using the Reproducing Kernel Particle Method showing improved accuracy and convergence rates for problems with challenging geometries. Conforming window functions are also demonstrated as a means to simplify the imposition of essential boundary conditions.

, rapid design-to-analysis↗
NASA NTRS2012

Determining Wheel-Soil Interaction Loads Using a Meshfree Finite Element Approach Assisting Future Missions with Rover Wheel Design

A wheel experiencing sinkage and slippage events poses a high risk to rover missions as evidenced by recent mobility challenges on the Mars Exploration Rover (MER) project. Because several factors contribute to wheel sinkage and slippage conditions such as soil composition, large deformation soil behavior, wheel geometry, nonlinear contact forces, terrain irregularity, etc., there are significant benefits to modeling these events to a sufficient degree of complexity. For the purposes of modeling wheel sinkage and slippage at an engineering scale, meshfree finite element approaches enable simulations that capture sufficient detail of wheel-soil interaction while remaining computationally feasible. This study demonstrates some of the large deformation modeling capability of meshfree methods and the realistic solutions obtained by accounting for the soil material properties. A benchmark wheel-soil interaction problem is developed and analyzed using a specific class of meshfree methods called Reproducing Kernel Particle Method (RKPM). The benchmark problem is also analyzed using a commercially available finite element approach with Lagrangian meshing for comparison. RKPM results are comparable to classical pressure-sinkage terramechanics relationships proposed by Bekker-Wong. Pending experimental calibration by future work, the meshfree modeling technique will be a viable simulation tool for trade studies assisting rover wheel design.

Contreras, Michael T.↗