We present a fully-coupled, implicit-in-time framework for solving a thermodynamically-consistent Cahn-Hilliard Navier-Stokes system that models two-phase flows. In this work, we extend the block iterative method presented in Khanwale et al. [Simulating two-phase flows with thermodynamically consistent energy stable Cahn-Hilliard Navier-Stokes equations on parallel adaptive octree based meshes, J. Comput. Phys. (2020)], to a fully-coupled, provably second-order accurate scheme in time, while maintaining energy-stability. The new method requires fewer matrix assemblies in each Newton iteration resulting in faster solution time. The method is based on a fully-implicit Crank-Nicolson scheme in time and a pressure stabilization for an equal order Galerkin formulation. That is, we use a conforming continuous Galerkin (cG) finite element method in space equipped with a residual-based variational multiscale (RBVMS) procedure to stabilize the pressure. We deploy this approach on a massively parallel numerical implementation using parallel octree-based adaptive meshes. We present comprehensive numerical experiments showing detailed comparisons with results from the literature for canonical cases, including the single bubble rise, Rayleigh-Taylor instability, and lid-driven cavity flow problems. We analyze in detail the scaling of our numerical implementation.
more »
« less
A dynamic variational multiscale method on unstructured meshes for stationary transport problems
This paper presents a variational multiscale (VMS) based finite element method where the stabilization parameter is computed dynamically. The current dynamic procedure takes in a general structure/form of the stabilization parameter with unknown coefficients and computes them dynamically in a local fashion resulting in a dynamic VMS‐based finite element method. Thus, a static stabilization parameter with pre‐defined coefficients is not needed. A variational Germano identity (VGI) based local procedure suitable for unstructured meshes is developed to perform the dynamic computation in a local fashion. The local VGI based procedure is applied for each interior vertex in the mesh and unknown coefficients are first determined locally at each vertex, and subsequently, for each element a maximum value is taken over the vertices of the element. To make the current procedure practical, a coarser secondary solution is constructed from the primary coarse‐scale solution, which is done locally over a patch of elements around each interior vertex. Further, averaging steps are employed to make the local dynamic procedure robust. Currently, the new dynamic VMS formulation is applied to steady problems governed by the advection‐diffusion and incompressible Navier‐Stokes equations in both 1D and 2D to demonstrate its efficacy and effectiveness.
more » « less- NSF-PAR ID:
- 10418806
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Fluids
- Volume:
- 95
- Issue:
- 7
- ISSN:
- 0271-2091
- Page Range / eLocation ID:
- p. 1117-1147
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.more » « less
-
In this paper, we present a GPU algorithm for finite element hyperelastic simulation. We show that the interior-point method, known to be effective for robust collision resolution, can be coupled with non-Newton procedures and be massively sped up on the GPU. Newton's method has been widely chosen for the interior-point family, which fully solves a linear system at each step. After that, the active set associated with collision/contact constraints is updated. Mimicking this routine using a non-Newton optimization (like gradient descent or ADMM) unfortunately does not deliver expected accelerations. This is because the barrier functions employed in an interior-point method need to be updated at every iteration to strictly confine the search to the feasible region. The associated cost (e.g., per-iteration CCD) quickly overweights the benefit brought by the GPU, and a new parallelism modality is needed. Our algorithm is inspired by the domain decomposition method and designed to move interior-point-related computations to local domains as much as possible. We minimize the size of each domain (i.e., a stencil) by restricting it to a single element, so as to fully exploit the capacity of modern GPUs. The stencil-level results are integrated into a global update using a novel hybrid sweep scheme. Our algorithm is locally second-order offering better convergence. It enables simulation acceleration of up to two orders over its CPU counterpart. We demonstrate the scalability, robustness, efficiency, and quality of our algorithm in a variety of simulation scenarios with complex and detailed collision geometries.more » « less
-
more » « less
Abstract Conformal mapping techniques have been used in many applications in the two-dimensional environments of engineering and physics, especially in the two-dimensional incompressible flow field that was introduced by Prandtl and Tietjens. These methods show reasonable results in the case of comprehensive analysis of the local coefficients of complex airfoils. The mathematical form of conformal mapping always locally preserves angles of the complex functions but it may change the length of the complex model. This research is based on the design of turbine blades as hydrofoils divided into different individual hydrofoils with decreasing thickness from root to tip. The geometric shapes of these hydrofoils come from the original FX77W121 airfoil shape and from interpolating between the FX77W121, FX77W153, and FX77W258 airfoil shapes. The last three digits of this airfoil family approximate the thickness ratio times 1000 (FX77153 => 15.3 % thickness ratio). Of the different airfoil shapes specified for the optimal rotor, there are 23 unique shapes.[15, 16, 17, 18, 19, 20, 21, 22, 24, 25, 28] This study describes the advantage of using at least one complex variable technique of transformation conformal mapping in two dimensions.
Conformal mapping techniques are used to form a database for sectional lift and drag coefficients based on turbine blade design to be used in Blade Element Momentum (BEM) theory to predict the performance of a three bladed single rotor horizontal axis ocean current turbine (1.6-meter diameter) by considering the characteristics of the sea-water. In addition, by considering the fact that in the real ocean, the underwater ocean current turbines encounter different velocities, the maximum brake power will be investigated for different incoming current velocities. The conformal mapping technique is used to calculate the local lift coefficients of different hydrofoils with respect to different angles of attack: −180 ≤ AOA ≤ +180. These results will be compared to those from other methods obtained recently by our research group. This method considers the potential flow analysis module that follows a higher-order panel method based on the geometric properties of each hydrofoil cross section. The velocity and pressure fields are obtained directly by the applications of Bernoulli’s principle, then the lift coefficients are calculated from the results of the integration of the pressure field along the hydrofoil surface for any angle of attack. Ultimately, the results of this research will be used for further investigation of the design and construction of a small-scale experimental ocean current turbine to be tested in the towing tank at the University of New Orleans.
-
more » « less
The solution of compressible flow equations is of interest with many aerospace engineering applications. Past literature has focused primarily on the solution of Computational Fluid Dynamics (CFD) problems with low-order finite element and finite volume methods. High-order methods are more the norm nowadays, in both a finite element and a finite volume setting. In this paper, inviscid compressible flow of an ideal gas is solved with high-order spectral/hp stabilized formulations using uniform high-order spectral element methods. The Euler equations are solved with high-order spectral element methods. Traditional definitions of stabilization parameters used in conjunction with traditional low-order bilinear Lagrange-based polynomials provide diffused results when applied to the high-order context. Thus, a revision of the definitions of the stabilization parameters was needed in a high-order spectral/hp framework. We introduce revised stabilization parameters, τsupg, with low-order finite element solutions. We also reexamine two standard definitions of the shock-capturing parameter, δ: the first is described with entropy variables, and the other is the YZβ parameter. We focus on applications with the above introduced stabilization parameters and analyze an array of problems in the high-speed flow regime. We demonstrate spectral convergence for the Kovasznay flow problem in both L1 and L2 norms. We numerically validate the revised definitions of the stabilization parameter with Sod’s shock and the oblique shock problems and compare the solutions with the exact solutions available in the literature. The high-order formulation is further extended to solve shock reflection and two-dimensional explosion problems. Following, we solve flow past a two-dimensional step at a Mach number of 3.0 and numerically validate the shock standoff distance with results obtained from NASA Overflow 2.2 code. Compressible flow computations with high-order spectral methods are found to perform satisfactorily for this supersonic inflow problem configuration. We extend the formulation to solve the implosion problem. Furthermore, we test the stabilization parameters on a complex flow configuration of AS-202 capsule analyzing the flight envelope. The proposed stabilization parameters have shown robustness, providing excellent results for both simple and complex geometries.
