Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

AN ESTIMATE OF SPENT NUCLEAR FUEL MECHANICAL LOADS IN THE GENERAL 30 CM PACKAGE DROP SCENARIO

The US Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is performing research to determine the mechanical loading conditions applied to spent nuclear fuel (SNF) during normal conditions of transport to inform mechanical tests of SNF and close an important knowledge gap related to the practical disposition of SNF in the US. Researchers at Pacific Northwest National Laboratory (PNNL) have completed an extensive finite element study to characterize and estimate the potential mechanical loads on SNF during a hypothetical 30 cm drop of an SNF transportation package. This modeling study is validated with test data collected by the SFWST program during a physical test campaign that included one-third scale package drop tests and full-scale single fuel assembly drop tests. The test campaign was led by Sandia National Laboratories (SNL) and included international collaboration with Equipos Nucleares S.A, S.M.E (ENSA) and Bundesanstalt für Materialforschung und -prüfung (BAM). The key contribution of the modeling study is to go beyond the limitations of the limited number of physical tests to estimate the impact response to variations in impact angle, initial gap conditions, fuel assembly design, burnup and other parameters that affect the mechanical loads. The methodology of this study included a classic parametric study to calculate the impact response of highly detailed fuel assemblies over many combinations of parameters. Models of a 17x17 pressurized water reactor fuel assembly and a generic 10x10 boiling water reactor fuel assembly were both used in this study to cover the major fuel assembly types in the US inventory. Over 2,000 impact responses were calculated. The results of the parametric study were evaluated using traditional methods and basic statistics. The results were also used to construct a damage model using multiple nonlinear regression techniques to predict the mechanical loads over the full range of all input parameters. The damage model was found to work very well for all impact angle cases where the cask came to rest on its side. It was concluded that end drop cases where the cask remained vertical (instead of tipping over onto its side) were not sufficiently characterized by the current set of parametric study cases to include in the damage model, but it was not a priority to fully investigate that range because the highest mechanical loads were observed in the broader range of side impact cases. This modeling work provides sufficient insight into the mechanical loads on SNF during a hypothetical 30 cm package drop that, when considered along with the physical test data collected by the SNL-led team, the SFSWT program can consider the knowledge gap closed.

Klymyshyn, Nicholas A.↗
DOE OSTI2024

Towards exascale for wind energy simulations

We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.

17 WIND ENERGY↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NEAMS-Multiphysics MOOSE End Year Framework Activities FY20

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a general finite element package meant for high performance solution of multiphysics problems in science and engineering. In this report, we present additions and enhancements to MOOSE funded by the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. NEAMS-funded MOOSE framework library improvements include: expansion of support for multi-level multi-application restart; enhancement of coupling between native MOOSE applications and external libraries; addition of a sparse automatic differentiation (AD) container enabling non-local degree of freedom coupling; block-specific quadrature rules; further development of an eigenvalue executioner; faster setup of periodic boundary conditions; and creation of a mesh meta-data system streamlining simulation startup. Besides developing the framework library, NEAMS funds were used to overhaul a swath of important MOOSE infrastructure in order to substantially improve user experience. These critical infrastructure changes included adapting MOOSE to python 3, improving the test harness to catch race conditions, and most importantly transitioning MOOSE to use Conda, a globally known packaging system, that greatly eases software adoption.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

COMPUTER-AIDED LATTICE DESIGN AND ADVANCED MODELING FOR THE DEVELOPMENT OF LIGHTWEIGHT STRUCTURAL MATERIALS

The notion of a “fission battery” conveys a vision focused on realizing very simple “plug-and-play” nuclear systems that can be integrated into a variety of applications requiring affordable, reliable energy in the form of electricity and/or heat and function without operations and maintenance staff. Fission batteries require lightweight structural materials to increase their mobility, and the lightweight materials must demonstrate structural resilience under various conditions. The objective of this work is to develop lattice structured lightweight structural material featuring a good combination of mechanical properties using advanced modeling and simulation together with an advanced additive manufacturing technique such as laser powder bed fusion. The preliminary results show that different lattice structures and types can be successfully meshed using nTopology software, and the lattice structure data can be successfully transformed to Multiphysics Object-Oriented Simulation (MOOSE) Environment input. Finite Element Analysis (FEA) displays that, at macro/engineering scale simulation, the weight saving design has an obvious effect on tensile behavior such as effective elastic modulus and yield stress. The novel approaches of this work are (1) development of lattice structures for improved mechanical properties using advanced simulation and modeling techniques; and (2) model predictions of the mechanical properties (e.g., strength and stress distribution) of macroscopic materials in order to preliminarily select a lattice structure for additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Validation and Parametric Investigations Using a Lumped Thermal Parameter Model of an Internal Permanent Magnet Motor

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multi-physics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a Lumped Parameter Thermal Network (LPTN) model of an Internal Permanent Magnet (IPM) synchronous motor within the commercially available Motor-CAD® modeling environment. The IPM motor’s stator is studied at steady state, and winding losses are generated by a constant DC current. The validation is based on temperature comparison with experimental data and with more detailed Finite Element Analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% — for the highest temperature in the slot windings and the end windings — was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot windings and end windings. As for the LTPN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance. The latter is often ignored in the literature, whereas its impact on temperature rise was found to be more significant than any other contact resistance within the stator.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Validation and Parametric Investigations Using a Lumped Thermal Parameter Model of an Internal Permanent Magnet Motor: Preprint

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multi-physics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a Lumped Parameter Thermal Network (LPTN) model of an Internal Permanent Magnet (IPM) synchronous motor within the commercially available Motor-CAD® modeling environment. The IPM motor’s stator is studied at steady state, and winding losses are generated by a constant DC current. The validation is based on temperature comparison with experimental data and with more detailed Finite Element Analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% - for the highest temperature in the slot windings and the end windings - was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot windings and end windings. As for the LTPN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance. The latter is often ignored in the literature, whereas its impact on temperature rise was found to be more significant than any other contact resistance within the stator.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗
DOE OSTI2024

Development of Lightweight Structural Materials with Improved Properties for Fission Batteries

The notion of a “fission battery” conveys a vision focused on realizing very simple “plug-and-play” nuclear systems that can be integrated into a variety of applications requiring affordable, reliable energy in the form of electricity and/or heat and function without operations and maintenance staff. Fission batteries require lightweight structural materials to increase their mobility, and the lightweight materials must demonstrate structural resilience under various conditions. The objective of this work is to develop lightweight structural material featuring a good combination of mechanical properties using advanced modeling and simulation together with an advanced additive manufacturing technique such as laser powder bed fusion. The preliminary results show that different lattice structures and types can be successfully meshed using nTopology software, and the lattice structure data can be successfully transformed to Multiphysics Object-Oriented Simulation Environment (MOOSE) input. Finite Element Analysis (FEA) displays that, at macro/engineering scale simulation, the weight saving design has an obvious effect on tensile behavior such as effective elastic modulus and yield stress. The novel approaches of this work are (1) development of optimized lattice structures for improved mechanical properties using advanced simulation and modeling techniques; and (2) model predictions of the mechanical properties (e.g., strength and stress distribution) of macroscopic materials in order to preliminarily select a lattice structure for additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The effect of a transducer’s spatial averaging on an elastodynamic guided wave’s wavenumber spectrum

Elastodynamic guided waves propagate in an elastic solid which makes it difficult, if not impractical, to place a receiving transducer in the direct path of the propagating wave (as one would for an acoustic wave in a fluid medium). Thus, to account for this, receiving transducers are often placed on the surface of the solid waveguide such that the transducer surface is parallel to the wave propagation direction. This transducer orientation introduces spatial averaging, which causes the received signal to have an altered signal amplitude and mode bias. We investigate both of these effects and present a simple model from which we derive a scaling-ratio expression that describes the effects of spatial averaging. Furthermore, we then test its performance using finite-element simulations that incorporate “real-world” assumptions (e.g., transient waves, piezoelectric effects, etc.). The results from the simulations demonstrate that the scaling-ratio can characterize the effects resulting from spatial averaging. The scaling-ratio expression will be particularly useful when designing experiments involving high frequency (small wavelength) guided waves. Lastly, the proposed scaling-ratio expression could be applied to other sensing methods, like Laser Doppler Vibrometry (LDV) or piezoelectric waveguides with inter-digitated transducers, because of its generality. However, the authors are careful to note that the scaling-ratio expression is not intended as a replacement of multi-physics analysis or good experimental design, and the effects of spatial averaging should be avoided whenever possible.

2DFFT↗
DOE OSTI2024

Operator product expansion for radial lattice quantization of 3D ϕ 4 theory

At its critical point, the three-dimensional lattice Ising model is described by a conformal field theory (CFT), the 3D Ising CFT. Instead of carrying out simulations on Euclidean lattices, we use the quantum finite elements method to implement radially quantized critical ϕ 4 theory on simplicial lattices approaching R × S 2 . Computing the four-point function of identical scalars, we demonstrate the power of radial quantization by the accurate determination of the scaling dimensions Δ ε and Δ T as well as ratios of the operator product expansion coefficients f σ σ ε and f σ σ T of the first spin-0 and spin-2 primary operators ε and T of the 3D Ising CFT. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Computationally Efficient Control Co-Design Optimization Framework with Mixed-Fidelity Fluid and Structure Analysis: OpenTurbineCoDe (OTCD) - ARPA-E ATLANTIS Phase I Project (Final Technical Report)

The goal of the Department of Energy (DOE) Advanced Research Projects Agency-Energy (ARPA-E) Aerodynamic Turbines Lighter and Afloat with Nautical Technologies and Integrated Servo-control (ATLANTIS) Program is to develop new technologies for floating offshore wind turbines, or wind farms, using the discipline of control co-design (CCD.) In this context, our goal is to develop a computationally efficient optimization framework for design of floating offshore wind turbines. Our specific aim is to utilize high-fidelity structural, aerodynamic, aero-structural tools, and to derive control-oriented reduced- or low-order models directly from the high-fidelity tools. We proposed a mixed-fidelity modeling approach which means that we are also using low- and mid-fidelity tools when necessary. This research is conducted by a multidisciplinary team consisting of Rutgers University, University of Michigan, Brigham Young University, and the National Renewable Energy Laboratory (NREL). The computational framework, called OpenTurbineCoDe, is designed to integrate, where possible, traditional structural, aerodynamic, aeroelastic models (e.g., OpenFAST) and advanced control algorithms with higher fidelity simulation tools including Reynolds-averaged Navier–Stokes (RANS) solvers, and three-dimensional structural finite element solvers. All the high-fidelity tools used in this research provide numerically exact gradients to facilitate both efficient optimization and local linearization for control implementation.

17 WIND ENERGY↗
DOE OSTI2020

Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture

Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.

Mechanics↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Effective field theory for radiative corrections to charged-current processes. II. Axial-vector coupling

Here we discuss the hadronic structure-dependent radiative corrections to the axial-vector coupling that controls single-nucleon weak charged-current processes—commonly denoted by 𝑔 𝐴 . We match the Standard Model at the GeV scale onto chiral perturbation theory at next-to-leading order in the one-nucleon sector, in the presence of electromagnetic and weak interactions. As a result, we provide a representation for the corrections to 𝑔 𝐴 in terms of infrared finite convolutions of simple kernels with the single-nucleon matrix elements of time-ordered products of two and three quark bilinears (vector, axial-vector, and pseudoscalar). We discuss strategies to determine the required nonperturbative input from data, lattice-QCD (+QED), and possibly hadronic models. This work paves the way for a precise comparison of the values of the ratio 𝑔 𝐴 ⁡/𝑔 𝑉 extracted from experiment and from lattice QCD, which constrain physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Predicting Mechanical Properties from Microstructure Images in Fiber-Reinforced Polymers Using Convolutional Neural Networks

Evaluating the mechanical response of fiber-reinforced composites can be extremely time-consuming and expensive. Machine learning (ML) techniques offer a means for faster predictions via models trained on existing input–output pairs and have exhibited success in composite research. This paper explores a fully convolutional neural network modified from StressNet, which was originally used for linear elastic materials, and extended here for a non-linear finite element (FE) simulation to predict the stress field in 2D slices of segmented tomography images of a fiber-reinforced polymer specimen. The network was trained and evaluated on data generated from the FE simulations of the exact microstructure. The testing results show that the trained network accurately captures the characteristics of the stress distribution, especially on fibers, solely from the segmented microstructure images. The trained model can make predictions within seconds in a single forward pass on an ordinary laptop, given the input microstructure, compared to 92.5 h to run the full FE simulation on a high-performance computing cluster. These results show promise in using ML techniques to conduct fast structural analysis for fiber-reinforced composites and suggest a corollary that the trained model can be used to identify the location of potential damage sites in fiber-reinforced polymers.

Sun, Yixuan (ORCID:0000000311093380)↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗