skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Effective floating volume: a highly parallelizable mesh-free approach for solving transient multiphysics problems in multi-scale geometries with non-linear material properties
A novel numerical method is proposed for the solution of transient multi-physics problems involving heat conduction, electrical current sharing and Joule heating. The innovation consists of a mesh-free Monte Carlo approach that eliminates or drastically reduces the particle scattering requirements typical of conventional Monte-Carlo methods. The proposed algorithm encapsulates a volume around each point that affects the solution at a given point in the domain; the volume includes other points that represent small perturbations along the path of energy transfer. The proposed method is highly parallelizable and amenable for GPU computing, and its computational performance was substantially increased by the elimination of scattered interpolation. The accuracy and simulation time of the proposed method are compared against a finite element solution and also against experimental results from existing literature. The proposed method provides accuracy comparable to that of finite element methods, achieving an order of magnitude reduction in simulation time.  more » « less
Award ID(s):
1827730
PAR ID:
10163441
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Computational mechanics
Volume:
65
ISSN:
1432-0924
Page Range / eLocation ID:
839–852
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Statistical solutions are time-parameterized probability measures on spaces of integrable functions, which have been proposed recently as a framework for global solutions and uncertainty quantification for multi-dimensional hyperbolic system of conservation laws. By combining high-resolution finite volume methods with a Monte Carlo sampling procedure, we present a numerical algorithm to approximate statistical solutions. Under verifiable assumptions on the finite volume method, we prove that the approximations, generated by the proposed algorithm, converge in an appropriate topology to a statistical solution. Numerical experiments illustrating the convergence theory and revealing interesting properties of statistical solutions are also presented. 
    more » « less
  2. The negative impact of measurement time skew on the static state estimation of the power grid has been exacerbated by increasing variation of system operating conditions. To mitigate the time skew problem, this paper proposes a regression model forecasting (RMF) method to forecast the time-skewed measurements, along with a confidence interval estimation (CIE) method to determine the weights associated with the forecasted measurements. The proposed RMF-CIE method is compared against several benchmark methods through Monte-Carlo simulation on the IEEE 16-machine, 68-bus model. It was observed that the proposed RMF-CIE consistently achieved more accurate state estimation on average. In addition, it was found that its estimation accuracy increases with the decrease of the skew time and variation levels. 
    more » « less
  3. This paper describes a novel numerical algorithm for the simulation of the along-wind dynamic response of a prototype of slender towers under turbulent winds, using a Layered Stochastic Approximation Monte Carlo algorithm (LSAMC). The proposed algorithm is applied to derive the statistics of the dynamic response in the presence of uncertainties in the structural properties and in the wind loading. Standard “brute force” Monte-Carlo methods are also used for validating the LSAMC results. The proposed methodology efficiently estimates structural fragility curves under extreme wind loads. The methodology enables a significant speedup in the computing time compared to standard Monte Carlo sampling. Furthermore, it is demonstrated that accuracy in the estimation of structural fragility curves is superior to ordinary reliability methods (e.g. “First-order reliability methods” or FORM). 
    more » « less
  4. Risk assessment of power transmission systems against strong winds requires models that can accurately represent the realistic performance of the physical infrastructure. Capturing material nonlinearity, p-delta effects in towers, buckling of lattice elements, joint slippage, and joint failure requires nonlinear models. For this purpose, this study investigates the reliability of transmission line systems by utilizing a nonlinear model of steel lattice towers, generated in OpenSEES platform. This model is capable of considering various geometric and material nonlinearities mentioned earlier. In order to efficiently estimate the probability of failure of transmission lines, the current study adopts an advanced reliability method through Error rate-based Adaptive Kriging (REAK) proposed by the authors. This method is capable of significantly reducing the number of simulations compared to conventional Monte Carlo methods such that reliability analysis can be done within a reasonable time. Results indicate that REAK efficiently estimates the reliability of transmission lines with a maximum of 150 Finite Element simulations for various wind intensities. 
    more » « less
  5. The diffusive tortuosity factor of a porous media quantifies the material’s resistance to diffusion, an important component of modeling flows in porous structures at the macroscale. Advances in X-ray micro-computed tomography (-CT) imaging provide the geometry of the material at the microscale (microstructure) thus enabling direct numerical simulation (DNS) of transport at the microscale. The data from these DNS are then used to close material’s macroscale transport models, which rely on effective material properties. In this work, we present numerical methods suitable for large scale simulations of diffusive transport through complex microstructures for the full range of Knudsen regimes. These numerical methods include a finite-volume method for continuum conditions, a random walk method for all regimes from continuum to rarefied, and the direct simulation Monte Carlo method. We show that for particle methods, the surface representation significantly affects the accuracy of the simulation for high Knudsen numbers, but not for continuum conditions. We discuss the upscaling of pore-resolved simulations to single species and multi-species volume-averaged models. Finally, diffusive tortuosities of a fibrous material are computed by applying the discussed numerical methods to 3D images of the actual microstructure obtained from X-ray computed micro-tomography. 
    more » « less