skip to main content

Title: An isogeometric collocation method for efficient random field discretization

This paper presents an isogeometric collocation method for a computationally expedient random field discretization by means of the Karhunen‐Loève expansion. The method involves a collocation projection onto a finite‐dimensional subspace of continuous functions over a bounded domain, basis splines (B‐splines) and nonuniform rational B‐splines (NURBS) spanning the subspace, and standard methods of eigensolutions. Similar to the existing Galerkin isogeometric method, the isogeometric collocation method preserves an exact geometrical representation of many commonly used physical or computational domains and exploits the regularity of isogeometric basis functions delivering globally smooth eigensolutions. However, in the collocation method, the construction of the system matrices for ad‐dimensional eigenvalue problem asks for at mostd‐dimensional domain integrations, as compared with 2d‐dimensional integrations required in the Galerkin method. Therefore, the introduction of the collocation method for random field discretization offers a huge computational advantage over the existing Galerkin method. Three numerical examples, including a three‐dimensional random field discretization problem, illustrate the accuracy and convergence properties of the collocation method for obtaining eigensolutions.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical Methods in Engineering
Page Range / eLocation ID:
p. 344-369
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Summary

    This study presents a thermo‐hydro‐mechanical (THM) model of unsaturated soils using isogeometric analysis (IGA). The framework employs Bézier extraction to connect IGA to the conventional finite element analysis (FEA), featuring the current study as one of the first attempts to develop an IGA‐FEA framework for solving THM problems in unsaturated soils. IGA offers higher levels of interelement continuity making it an attractive method for solving highly nonlinear problems. The governing equations of linear momentum, mass, and energy balance are coupled based on the averaging procedure within the hybrid mixture theory. The Drucker‐Prager yield surface is used to limit the modified effective stress where the model follows small strain, quasi‐static loading conditions. Temperature dependency of the surface tension is implemented in the soil‐water retention curve. Nonuniform rational B‐splines (NURBS) basis functions are used in the standard Galerkin method and weak formulations of the balance equations. Displacement, capillary pressure, gas pressure, and temperature are four independent quantities that are approximated by NURBS in spatial discretization. The framework is used to simulate strain localization in an undrained dense sand subjected to plane strain biaxial compression under different temperatures and displacement velocities. Results show that an increase in the displacement rate leads to reduction in the equivalent plastic strain while an increase in the temperature leads to an increase in the equivalent plastic strain. The findings suggest that the proposed IGA‐based framework offers a viable alternative for solving THM problems in unsaturated soils.

    more » « less
  2. Summary

    This paper presents an approach for efficient uncertainty analysis (UA) using an intrusive generalized polynomial chaos (gPC) expansion. The key step of the gPC‐based uncertainty quantification(UQ) is the stochastic Galerkin (SG) projection, which can convert a stochastic model into a set of coupled deterministic models. The SG projection generally yields a high‐dimensional integration problem with respect to the number of random variables used to describe the parametric uncertainties in a model. However, when the number of uncertainties is large and when the governing equation of the system is highly nonlinear, the SG approach‐based gPC can be challenging to derive explicit expressions for the gPC coefficients because of the low convergence in the SG projection. To tackle this challenge, we propose to use a bivariate dimension reduction method (BiDRM) in this work to approximate a high‐dimensional integral in SG projection with a few one‐ and two‐dimensional integrations. The efficiency of the proposed method is demonstrated with three different examples, including chemical reactions and cell signaling. As compared to other UA methods, such as the Monte Carlo simulations and nonintrusive stochastic collocation (SC), the proposed method shows its superior performance in terms of computational efficiency and UA accuracy.

    more » « less
  3. Neurons exhibit remarkably complex geometry in their neurite networks. So far, how materials are transported in the complex geometry for survival and function of neurons remains an unanswered question. Answering this question is fundamental to understanding the physiology and disease of neurons. Here, we have developed an isogeometric analysis (IGA) based platform for material transport simulation in neurite networks. We modeled the transport process by reaction-diffusion-transport equations and represented geometry of the networks using truncated hierarchical tricubic B-splines (THB-spline3D). We solved the Navier-Stokes equations to obtain the velocity field of material transport in the networks. We then solved the transport equations using the streamline upwind/Petrov-Galerkin (SU/PG) method. Using our IGA solver, we simulated material transport in three basic models of the network geometry: a single neurite, a neurite bifurcation, and a neurite tree with three bifurcations. In addition, the robustness of our solver is illustrated by simulating material transport in three representative and complex neurite networks. From the simulation we discovered several spatial patterns of the transport process. Together, our simulation provides key insights into how material transport in neurite networks is mediated by their complex geometry. 
    more » « less
  4. The turbulent channel flow database is produced from a direct numerical simulation (DNS) of wall bounded flow with periodic boundary conditions in the longitudinal and transverse directions, and no-slip conditions at the top and bottom walls. In the simulation, the Navier-Stokes equations are solved using a wall {normal, velocity {vorticity formulation. Solutions to the governing equations are provided using a Fourier-Galerkin pseudo-spectral method for the longitudinal and transverse directions and seventh-order Basis-splines (B-splines) collocation method in the wall normal direction. De-aliasing is performed using the 3/2-rule [3]. Temporal integration is performed using a low-storage, third-order Runge-Kutta method. Initially, the flow is driven using a constant volume flux control (imposing a bulk channel mean velocity of U = 1) until stationary conditions are reached. Then the control is changed to a constant applied mean pressure gradient forcing term equivalent to the shear stress resulting from the prior steps. Additional iterations are then performed to further achieve statistical stationarity before outputting fields. 
    more » « less
  5. This paper presents the enriched Galerkin discretization for modeling fluid flow in fractured porous media using the mixed-dimensional approach. The proposed method has been tested against published benchmarks. Since fracture and porous media discontinuities can significantly influence single- and multi-phase fluid flow, the heterogeneous and anisotropic matrix permeability setting is utilized to assess the enriched Galerkin performance in handling the discontinuity within the matrix domain and between the matrix and fracture domains. Our results illustrate that the enriched Galerkin method has the same advantages as the discontinuous Galerkin method; for example, it conserves local and global fluid mass, captures the pressure discontinuity, and provides the optimal error convergence rate. However, the enriched Galerkin method requires much fewer degrees of freedom than the discontinuous Galerkin method in its classical form. The pressure solutions produced by both methods are similar regardless of the conductive or non-conductive fractures or heterogeneity in matrix permeability. This analysis shows that the enriched Galerkin scheme reduces the computational costs while offering the same accuracy as the discontinuous Galerkin so that it can be applied for large-scale flow problems. Furthermore, the results of a time-dependent problem for a three-dimensional geometry reveal the value of correctly capturing the discontinuities as barriers or highly-conductive fractures. 
    more » « less