Search

The past few decades have witnessed a growing popularity in Eulerian–Lagrangian solvers due to their significant potential for simulating aerodynamic flows, particularly in cases involving strong body–vortex interactions. In this hybrid approach, the two component solvers are mutually coupled in a two-way fashion. Initially, the Lagrangian solver can supply boundary conditions to the Eulerian solver, while the Eulerian solver functions as a corrector for the Lagrangian solution in regions where the latter cannot achieve high accuracy. To utilize such tools effectively, it is vital for them to be capable of handling dynamic mesh movements. This study builds upon the previous research conducted by our team and extends the capabilities of the hybrid solver to handle dynamic meshes. While OpenFOAM, the Eulerian component of this hybrid code, incorporates built-in dynamic mesh properties, certain modifications are necessary to ensure its compatibility with the Lagrangian solver. More specifically, the evolution algorithm of the pimpleFOAM solver needs to be divided into two discrete steps: first, updating the mesh, and later, evolving the solution. This division enables a proper coupling between pimpleFOAM and the Lagrangian solver as an intermediate step. Therefore, the primary objective of this specific paper is to adapt the OpenFOAM solver to meet the demands of the hybrid solver and subsequently validate that the hybrid solver can effectively address dynamic mesh challenges using this approach. This approach introduces a pioneering method for conducting dynamic mesh simulations within the OpenFOAM framework, showcasing its potential for broader applications. To validate the approach, various test cases involving dynamic mesh movements are employed. Specifically, all these cases employ the Lamb–Oseen diffusing vortex, but each case incorporates different types of mesh movements, including translational, rotational, oscillational, and combinations thereof. The resul, Wind Energy, Aerodynamics

To investigate the effect of force distributions of each turbine component on the power performance of the Darrieus–Savonius combined vertical axis wind turbine (hybrid VAWT), the hybrid VAWT is modeled as idealized turbine under various force distributions. The goal of idealization is to simplify the intricate interactions between the Savonius and Darrieus components. The simulation actuator surfaces with uniform force distributions lead to a cost-effective way to identify the optimal force distribution of each turbine component. The numerical model was validated against momentum theory. The results demonstrated that the numerical and theoretical results yield similar predictions in the low-thrust cases but show differences in the high-thrust cases. The maximum power coefficient (Formula presented.) of an idealized hybrid VAWT with given thrust coefficient (Formula presented.) is lower than that of a single actuator. This is a consequence of the nonoptimal loading on the actuator. The results indicate that an idealized hybrid VAWT does not show a significant power increase compared with an optimal single Darrieus rotor. Therefore, the presence of a Savonius rotor inside a Darrieus rotor leads to a lower power output in any circumstance. The hybrid configuration is primarily advantageous for the start-up performance of the combined rotor, which is not explored in this study., Wind Energy, Aerodynamics

In the field of computational aerodynamics, it is vital to develop tools that can accurately, but also efficiently, simulate the flow around bluff objects and calculate the aerodynamic forces acting on them. When strong body–vortex interactions take place, the simulations become more demanding, since complex phenomena appear. To address this issue, hybrid Eulerian–Lagrangian solvers have been developed and are increasingly used in the field. In this paper, a Vortex Particle Method (VPM) is coupled with the OpenFOAM software. The Eulerian solver (OpenFOAM) resolves the regions close to the solid boundaries, while the vortex particles evolve the wake downstream, significantly reducing artificial diffusion. The coupling strategy and the validation results of a hybrid code based on the domain decomposition technique are presented. This work is the first to couple OpenFOAM with a Lagrangian solver in the framework of a hybrid solver. Our objective is twofold: to verify the capability of OpenFOAM to run with a VPM and to validate the hybrid solver using benchmark cases. We demonstrate the validation of the solver on the Lamb–Oseen vortex case, the dipole case in the unbounded domain, and the flow around a cylinder at Re = 550. Our results show that coupling OpenFOAM with a VPM can be achieved without complications and efficiently reproduces the results of pure Eulerian simulations., Wind Energy, Aerodynamics

In nuclear fuel rod bundles, turbulence-induced pressure fluctuations caused by an axial flow can create small but significant vibrations in the fuel rods, which in turn can cause structural effects such as material fatigue and fretting wear. Fluid-structure interaction simulations can be used to model these vibrations, but for affordable simulations based on the URANS approach, a model for the pressure fluctuations must be utilised. Driven by the goal to improve the current state-of-the-art pressure fluctuation model, AniPFM (Anisotropic Pressure Fluctuation Model) was developed. AniPFM can model velocity fluctuations based on anisotropic Reynolds stress tensors, with temporal correlation through the convection and decorrelation of turbulence. From these velocity fluctuations and the mean flow properties, the pressure fluctuations are calculated. The model was applied to several test cases and shows promising results in terms of reproducing qualitatively similar flow structures, as well as predicting the root-mean-squared pressure fluctuations. While further validation is being performed, the AniPFM has already demonstrated its potential for affordable simulations of turbulence-induced vibrations in industrial nuclear applications., Aerodynamics

This study investigates the implementation of the vortex particle method (VPM) with the goal of efficiently and accurately estimating the power performances and flow characteristics for a Savonius rotor. The accuracy and efficiency of simulation methods are critical for the reliable design of Savonius rotors. Among various approaches, VPM is chosen because it can be flexibly incorporated with self-correction techniques, and the distribution of bound vortex particles can effectively represent complex geometries. In this work, a double-trailing-edge-wake-modeling vortex particle method (DTVPM) is presented to extend the working range of VPM for dealing with large rotating amplitudes and high tip speed ratios (TSRs). DTVPM addresses asymmetrical torque predictions for a Savonius rotor without gap width. However, DTVPM performs poorly at high TSRs due to the absence of viscous effects near the surface. To capture complex wake structures, such as reverse flow structures, the viscous correction for tip vortices is suggested. The current research focuses on the implementation and validation of DTVPM for predicting torque coefficients and wake patterns, as well as comparisons to OpenFOAM results. Two-dimensional and incompressible flow is estimated at (Formula presented.) = 0.2–1.2. For the studied cases, a maximum power coefficient is obtained at (Formula presented.), consistent with published experimental data. In addition, the process of trailing-edge vortices generation and detachment is captured by DTVPM. The comparison results between OpenFOAM and DTVPM show that DTVPM allows to efficiently simulate a Savonius rotor without any empirical parameters. DTVPM will help to improve existing engineering models for wind energy fields., Wind Energy, Aerodynamics

Flapping wings display complex flows which can be used to generate large lift forces. Flexibility in wings is widely used by natural flyers to increase the aerodynamic performance. The influence of wing flexibility on the flow can be computed using numerical analysis with Fluid Structure Interaction (FSI). The influence of inertial, elastic and aerodynamic forces is quantified using a 2D wing. A sinusoidal flapping motion is imposed on the leading edge of the vertical wing. The inertial force on the wing dominates for high mass ratios and the wing deflection is rather independent of the flow. For a low mass ratio, the wing deformation scales with the increasing elasticity. The maximum lift and lowest drag were found for the wing with large flexibility and low mass so the passive deformation by aerodynamic forces creates a favourable shape for lift production. Flexible translating and revolving wings at an angle of attack of 45 degrees show that chordwise flexibility decreases both lift and drag, however the lift over drag ratio is increased. The flow around both wings forms a coherent structure with a Root Vortex (RV), Tip Vortex (TV), Leading Edge Vortex (LEV) and Trailing Edge Vortex (TEV). The LEV on the revolving wing is stable for approximately up to half the span because vorticity is transported outward in the vortex core. The flowfield and LEV breakdown are consistent with experimental data of the same wing. The translating wing builds up circulation but the LEV detaches quickly near the centre of the wing. Chordwise bending reduces the angle of attack which decreases the distance to the core of the shed LEVs., Aerodynamics

Flapping wings display complex flows which can be used to generate large lift forces. Flexibility in wings is widely used by natural flyers to increase the aerodynamic performance. The influence of wing flexibility on the flow can be computed using numerical analysis with Fluid Structure Interaction (FSI). The influence of inertial, elastic and aerodynamic forces is quantified using a 2D wing. A sinusoidal flapping motion is imposed on the leading edge of the vertical wing. The inertial force on the wing dominates for high mass ratios and the wing deflection is rather independent of the flow. For a low mass ratio, the wing deformation scales with the increasing elasticity. The maximum lift and lowest drag were found for the wing with large flexibility and low mass so the passive deformation by aerodynamic forces creates a favourable shape for lift production. Flexible translating and revolving wings at an angle of attack of 45 degrees show that chordwise flexibility decreases both lift and drag, however the lift over drag ratio is increased. The flow around both wings forms a coherent structure with a Root Vortex (RV), Tip Vortex (TV), Leading Edge Vortex (LEV) and Trailing Edge Vortex (TEV). The LEV on the revolving wing is stable for approximately up to half the span because vorticity is transported outward in the vortex core. The flowfield and LEV breakdown are consistent with experimental data of the same wing. The translating wing builds up circulation but the LEV detaches quickly near the centre of the wing. Chordwise bending reduces the angle of attack which decreases the distance to the core of the shed LEVs., Aerodynamics

Despite advances in turbulence modelling, the Smagorinsky model remains a popular choice for large eddy simulation (LES) due to its simplicity and ease of use. The dissipation in turbulence energy that the model introduces, is proportional to the Smagorinsky constant, of which many different values have been proposed. These values have been derived for certain simulated test-cases while using a specific set of numerical schemes, to obtain the correct dissipation in energy simply because an incorrect value of the Smagorinsky constant would lead to an incorrect dissipation. However, it is important to bear in mind that numerical codes may suffer from numerical or artificial dissipation, which occurs spuriously through a combination of spatio-temporal and iterative errors. The latter can be controlled through more iterations, the former however, depends on the grid resolution and the time step. Recent research suggests that a complete energy-conserving (EC) spatio-temporal discretisation guarantees zero numerical dissipation for any grid resolution and time step. Therefore, using an EC scheme would ensure that dissipation occurs primarily through the Smagorinsky model (and errors in its implementation) than through the discretisation of the Navier-Stokes (NS) equations. To evaluate the efficacy of these schemes for engineering applications, the article first discusses the use of an EC temporal discretisation as regards to accuracy and computational effort, to ascertain whether EC time advancement is advantageous or not. It was noticed that a simple non-EC explicit method with a smaller time step not only reduces the numerical dissipation to an acceptable level but is computationally cheaper than an implicit-EC scheme for wide range of time steps. Secondly, in terms of spatial discretisation on uniform grids (popular in LES), a simple central-difference scheme is as accurate as an EC spatial discretisation. Finally, following the removal of numerical dissipation with an, Wind Energy, Aerodynamics

Despite advances in turbulence modelling, the Smagorinsky model remains a popular choice for large eddy simulation (LES) due to its simplicity and ease of use. The dissipation in turbulence energy that the model introduces, is proportional to the Smagorinsky constant, of which many different values have been proposed. These values have been derived for certain simulated test-cases while using a specific set of numerical schemes, to obtain the correct dissipation in energy simply because an incorrect value of the Smagorinsky constant would lead to an incorrect dissipation. However, it is important to bear in mind that numerical codes may suffer from numerical or artificial dissipation, which occurs spuriously through a combination of spatio-temporal and iterative errors. The latter can be controlled through more iterations, the former however, depends on the grid resolution and the time step. Recent research suggests that a complete energy-conserving (EC) spatio-temporal discretisation guarantees zero numerical dissipation for any grid resolution and time step. Therefore, using an EC scheme would ensure that dissipation occurs primarily through the Smagorinsky model (and errors in its implementation) than through the discretisation of the Navier-Stokes (NS) equations. To evaluate the efficacy of these schemes for engineering applications, the article first discusses the use of an EC temporal discretisation as regards to accuracy and computational effort, to ascertain whether EC time advancement is advantageous or not. It was noticed that a simple non-EC explicit method with a smaller time step not only reduces the numerical dissipation to an acceptable level but is computationally cheaper than an implicit-EC scheme for wide range of time steps. Secondly, in terms of spatial discretisation on uniform grids (popular in LES), a simple central-difference scheme is as accurate as an EC spatial discretisation. Finally, following the removal of numerical dissipation with an, Wind Energy, Aerodynamics

This paper aims to provide an explanation for the overprediction of aerodynamic loads by CFD compared to experiments for the MEXICO wind turbine rotor and improve the CFD prediction by considering laminar-turbulent transition modeling. Large deviations between CFD results and experimental measurements are observed in terms of sectional normal and tangential forces at the blade tip (r/R=0.82 and 0.92) of the MEXICO rotor operating in axial condition at the design tip speed ratio λ=6.7. The first part of this study identifies the effects of ZigZag tape, which is used in the experiment to trigger boundary layer transition, by analyzing the available experimental data of a single, non-rotating MEXICO rotor blade. The analysis indicates that ZigZag tape has a significant impact on sectional aerodynamic tip loads: it alters the boundary layer thickness and additionally reduces the effective airfoil camber besides the expected tripping. These additional effects most likely also occur in the rotating MEXICO experiment, reducing the sectional loads and hence lead to an overprediction by CFD. To eliminate the ZigZag tape interference, experimental data with an untripped blade is preferred to be used as validation case. In the second part of this study, a transitional flow simulation for the MEXICO rotor is performed by using RANS-based transition model k−kL−ω within OpenFOAM-2.1.1. The numerical results are compared against experimental data obtained from the untripped, new MEXICO experiments. The comparison gives that transitional simulation present a very good tip loads prediction for the untripped blade. The measured data also confirms that the ZigZag tape indeed has a significant influence on the blade tip loads in rotating conditions. The transition onset over 3D MEXICO blade is visualized and transition locations are identified. The results shown in the present study can explain the causes of the large differences between CFD and experiment observed in the, Wind Energy, Aerodynamics

An optimization approach is presented that can be used to find the optimal source term distribution in order to represent a high-fidelity vortex-generator (VG) induced flow field on a coarse mesh. The ap- proach employs the continuous adjoint of the problem, from which an exact sensitivity is calculated and used in combination with a trust-region method to find the source term which minimizes the deviation with respect to the reference velocity field. The algorithm is applied to an incompressible flow over a rectangular VG and VG pair on a flat plate and compared to results obtained with the jBAY-model and a body-fitted mesh simulation. The results indicate that a highly accurate flow, yielding only minimal errors with respect to the shape factor, circulation and vortex core, can be obtained on coarse meshes when adding a source term to only a limited number of cells. This approach therefore demonstrates the potential of source-term models to include the effects of VGs in computations of large-scale geometries. It also allows quantification of the achievable accuracy on a particular mesh and the calculation of the source term which is optimal for a specific situation. Furthermore, the optimization approach can be used to diagnose the deficiencies of an existing source-term VG model, in this work the jBAY model., Aerodynamics, Education and Student Affairs

Vortex generators (VGs) are a widely used means of flow control, and predictions of their influence are vital for efficient designs. However, accurate CFD simulations of their effect on the flow field by means of a body fitted mesh are computationally expensive. Therefore the BAY and jBAY models, which represent the effect of VGs on the flow using source terms in the momentum equations, are popular in industry. In this contribution we examine the ability of the BAY and jBAY model to provide accurate flow field results by looking at boundary layer properties close behind VGs. The results are compared with both body fitted mesh and other source term model RANS simulations of 3D incompressible flows, over flat plate and airfoil geometries. We show the influence of mesh resolution and domain of application on the accuracy of the models and investigate the influence of the source term on the generated flow field. Our results demonstrate the grid dependence of the models and indicate the presence of model errors. Furthermore we find that the total applied force has a larger influence on both the intensity and shape of the created vortex than the distribution of the source term over the cells., Aerodynamics, Education and Student Affairs

This paper presents the experimental and numerical study on MEXICO wind turbine blades. Previous work by other researchers shows that large deviations exist in the loads comparison between numerical predictions and experimental data for the rotating MEXICO wind turbine. To reduce complexities and uncertainties, a non-rotating experimental campaign has been carried out on MEXICO blades Delft University of Technology. In this new measurement, quasi-2D aerodynamic characteristics of MEXICO blades on three spanwise sections are measured at different inflow velocities and angles of attack. Additionally, RANS simulations are performed with OpenFOAM-2.1.1 to compare numerical results against measured data. The comparison and analysis of aerodynamic loads on the blade, where three different airfoil families and geometrical transition regions are used, show that for attached flow condition, RANS computation predicts excellent pressure distribution on the NACA airfoil section (r=R D 0.92) and good agreement is observed on the DU (r=R D 0.35) and RISØ (r=R D 0.60) airfoil sections. Unexpected aerodynamic characteristics are observed at the intermediate transition regions connecting the RISØ and DU airfoils, where sudden lift force drop is found at the radial position r=R D 0.55. Through numerical flow visualization, large-scale vortical structures are observed on the suction side of the blade near the mid-span. Moreover, counter-rotating vortices are generated behind the blade at locations where unexpected loads occurs. Consequently, the RISØ airfoil could not give expected 2D aerodynamic characteristics because of upwash/downwash effects induced by these counter-rotating vortices, which make 3D effects play an important role in numerical modeling when calculating the aerodynamic loads for MEXICO rotor., Wind Energy, Aerodynamics

A PhD project is being carried out on the topic of far-wake aerodynamics of Horizontal Axis Wind Turbines (HAWTs) in yawed conditions, which has a large relevance for wind farm design and optimization. Characteristic for a turbine in yaw are the inherent unsteady and non-uniform rotor loading, and the typical wake deflection and strong three-dimensional deformation effects under influence of self-induction (see figure 1). Investigation of HAWTs in yaw is important, as the larce-scale eddies of the turbulent atmosphere dictate that a wind turbine is in practise always operating in unsteady yaw, while the resulting wake effects are already significant for small yaw angles. Despite this relevance, research into the far-wake of yawed wind turbines has been very limited and the symmetry assumptions on which common wake engineering models are based conflict with the physics of the skewed wake of a yawed turbine. Nevertheless, there is an increasing interest into this topic, as it is recognized that the effect of wake deflection can be exploited as a way to optimize the overall wind farm power production through active yaw control. For this purpose, simple two-dimensional models are applied for approximating the wake deflection, but which are unable to capture the typical three-dimensional deformation effects. In summary, there is a large gap of fundamental knowledge on wake physics in yawed conditions, and what the relevance of these phenomena is on the development and issues like the re-energization process of the far-wake. To bridge this gap, the PhD project aims at improving our understanding of the wake physics of HAWTs in yaw and to draft guidelines for reduced-order models that can be applied for wind farm design and optimization. In support of this aim, the objective is to analyze the different physical “modes” that play a role in the yawed wake, through a numerial and experimental investigation of the skewed wakes aft of HAWTs and actuator discs. The results from, Wind Energy, Aerodynamics

A PhD project is being carried out on the topic of far-wake aerodynamics of Horizontal Axis Wind Turbines (HAWTs) in yawed conditions, which has a large relevance for wind farm design and optimization. Characteristic for a turbine in yaw are the inherent unsteady and non-uniform rotor loading, and the typical wake deflection and strong three-dimensional deformation effects under influence of self-induction (see figure 1). Investigation of HAWTs in yaw is important, as the larce-scale eddies of the turbulent atmosphere dictate that a wind turbine is in practise always operating in unsteady yaw, while the resulting wake effects are already significant for small yaw angles. Despite this relevance, research into the far-wake of yawed wind turbines has been very limited and the symmetry assumptions on which common wake engineering models are based conflict with the physics of the skewed wake of a yawed turbine. Nevertheless, there is an increasing interest into this topic, as it is recognized that the effect of wake deflection can be exploited as a way to optimize the overall wind farm power production through active yaw control. For this purpose, simple two-dimensional models are applied for approximating the wake deflection, but which are unable to capture the typical three-dimensional deformation effects. In summary, there is a large gap of fundamental knowledge on wake physics in yawed conditions, and what the relevance of these phenomena is on the development and issues like the re-energization process of the far-wake. To bridge this gap, the PhD project aims at improving our understanding of the wake physics of HAWTs in yaw and to draft guidelines for reduced-order models that can be applied for wind farm design and optimization. In support of this aim, the objective is to analyze the different physical “modes” that play a role in the yawed wake, through a numerial and experimental investigation of the skewed wakes aft of HAWTs and actuator discs. The results from, Wind Energy, Aerodynamics

This paper presents the experimental and numerical study on MEXICO wind turbine blades. Previous work by other researchers shows that large deviations exist in the loads comparison between numerical predictions and experimental data for the rotating MEXICO wind turbine. To reduce complexities and uncertainties, a non-rotating experimental campaign has been carried out on MEXICO blades Delft University of Technology. In this new measurement, quasi-2D aerodynamic characteristics of MEXICO blades on three spanwise sections are measured at different inflow velocities and angles of attack. Additionally, RANS simulations are performed with OpenFOAM-2.1.1 to compare numerical results against measured data. The comparison and analysis of aerodynamic loads on the blade, where three different airfoil families and geometrical transition regions are used, show that for attached flow condition, RANS computation predicts excellent pressure distribution on the NACA airfoil section (r=R D 0.92) and good agreement is observed on the DU (r=R D 0.35) and RISØ (r=R D 0.60) airfoil sections. Unexpected aerodynamic characteristics are observed at the intermediate transition regions connecting the RISØ and DU airfoils, where sudden lift force drop is found at the radial position r=R D 0.55. Through numerical flow visualization, large-scale vortical structures are observed on the suction side of the blade near the mid-span. Moreover, counter-rotating vortices are generated behind the blade at locations where unexpected loads occurs. Consequently, the RISØ airfoil could not give expected 2D aerodynamic characteristics because of upwash/downwash effects induced by these counter-rotating vortices, which make 3D effects play an important role in numerical modeling when calculating the aerodynamic loads for MEXICO rotor., Wind Energy, Aerodynamics

Vortex generators (VGs) are a widely used means of flow control, and predictions of their influence are vital for efficient designs. However, accurate CFD simulations of their effect on the flow field by means of a body fitted mesh are computationally expensive. Therefore the BAY and jBAY models, which represent the effect of VGs on the flow using source terms in the momentum equations, are popular in industry. In this contribution we examine the ability of the BAY and jBAY model to provide accurate flow field results by looking at boundary layer properties close behind VGs. The results are compared with both body fitted mesh and other source term model RANS simulations of 3D incompressible flows, over flat plate and airfoil geometries. We show the influence of mesh resolution and domain of application on the accuracy of the models and investigate the influence of the source term on the generated flow field. Our results demonstrate the grid dependence of the models and indicate the presence of model errors. Furthermore we find that the total applied force has a larger influence on both the intensity and shape of the created vortex than the distribution of the source term over the cells., Aerodynamics, Education and Student Affairs

An optimization approach is presented that can be used to find the optimal source term distribution in order to represent a high-fidelity vortex-generator (VG) induced flow field on a coarse mesh. The ap- proach employs the continuous adjoint of the problem, from which an exact sensitivity is calculated and used in combination with a trust-region method to find the source term which minimizes the deviation with respect to the reference velocity field. The algorithm is applied to an incompressible flow over a rectangular VG and VG pair on a flat plate and compared to results obtained with the jBAY-model and a body-fitted mesh simulation. The results indicate that a highly accurate flow, yielding only minimal errors with respect to the shape factor, circulation and vortex core, can be obtained on coarse meshes when adding a source term to only a limited number of cells. This approach therefore demonstrates the potential of source-term models to include the effects of VGs in computations of large-scale geometries. It also allows quantification of the achievable accuracy on a particular mesh and the calculation of the source term which is optimal for a specific situation. Furthermore, the optimization approach can be used to diagnose the deficiencies of an existing source-term VG model, in this work the jBAY model., Aerodynamics, Education and Student Affairs

Experimental results and complimentary computations for airfoils with vortex generators are compared in this paper, as part of an e_ort within the AVATAR project to develop tools for wind turbine blade control devices. Measurements from two airfoils equipped with passive vortex generators, a 30% thick DU97W300 and an 18% thick NTUA T18 have been used for benchmarking several simulation tools. These tools span low-to-high complexity, ranging from engineering-level integral boundary layer tools to fully-resolved computational uid dynamics codes. Results indicate that with appropriate calibration, engineering-type tools can capture the e_ects of vortex generators and outperform more complex tools. Fully resolved CFD comes at a much higher computational cost and does not necessarily capture the increased lift due to the VGs. However, in lieu of the limited experimental data available for calibration, high _delity tools are still required for assessing the e_ect of vortex generators on airfoil performance., Wind Energy, Aerodynamics, Education and Student Affairs

In this article, we endeavour to find a fast solver for finite volume discretizations for compressible unsteady viscous flows. Thereby, we concentrate on comparing the efficiency of important classes of time integration schemes, namely time adaptive Rosenbrock, singly diagonally implicit (SDIRK) and explicit first stage singly diagonally implicit Runge-Kutta (ESDIRK) methods. To make the comparison fair, efficient equation system solvers need to be chosen and a smart choice of tolerances is needed. This is determined from the tolerance TOL that steers time adaptivity. For implicit Runge-Kutta methods, the solver is given by preconditioned inexact Jacobian-free Newton-Krylov (JFNK) and for Rosenbrock, it is preconditioned Jacobian-free GMRES. To specify the tolerances in there, we suggest a simple strategy of using TOL/100 that is a good compromise between stability and computational effort. Numerical experiments for different test cases show that the fourth order Rosenbrock method RODASP and the fourth order ESDIRK method ESDIRK4 are best for fine tolerances, with RODASP being the most robust scheme., Aerodynamics, Science Centre & Programmering

The FP7 AdVanced Aerodynamic Tools for lArge Rotors - Avatar project aims to develop and validate advanced aerodynamic models, to be used in integral design codes for the next generation of large scale wind turbines (10-20MW). One of the approaches towards reaching rotors for 10-20MW size is the application of flow control devices, such as flaps. In Task 3.2: Development of aerodynamic codes for modelling of flow devices on aerofoils and, rotors of the Avatar project, aerodynamic codes are benchmarked and validated against the experimental data of a DU95W180 airfoil in steady and unsteady flow, for different angle of attack and flap settings, including unsteady oscillatory trailing-edge-flap motion, carried out within the framework of WP3: Models for Flow Devices and Flow Control, Task 3.1: CFD and Experimental Database. The aerodynamics codes are: AdaptFoil2D, Foil2W, FLOWer, MaPFlow, OpenFOAM, Q3UIC, ATEFlap. The codes include unsteady Eulerian CFD simulations with grid deformation, panel models and indicial engineering models. The validation cases correspond to 18 steady flow cases, and 42 unsteady flow cases, for varying angle of attack, flap deflection and reduced frequency, with free and forced transition. The validation of the models show varying degrees of agreement, varying between models and flow cases., Wind Energy, Aerodynamics

High computational cost associated with the high fidelity wake models such as RANS or LES serves as a primary bottleneck to perform a direct high fidelity wind farm layout optimization (WFLO) using accurate CFD based wake models. Therefore, a surrogate based multi-fidelity WFLO methodology (SWFLO) is proposed. The surrogate model is built using an SBO method referred as manifold mapping (MM). As a verification, optimization of spacing between two staggered wind turbines was performed using the proposed surrogate based methodology and the performance was compared with that of direct optimization using high fidelity model. Significant reduction in computational cost was achieved using MM: a maximum computational cost reduction of 65%, while arriving at the same optima as that of direct high fidelity optimization. The similarity between the response of models, the number of mapping points and its position, highly influences the computational efficiency of the proposed method. As a proof of concept, realistic WFLO of a small 7-turbine wind farm is performed using the proposed surrogate based methodology. Two variants of Jensen wake model with different decay coefficients were used as the fine and coarse model. The proposed SWFLO method arrived at the same optima as that of the fine model with very less number of fine model simulations., Aerodynamics, Aerospace Engineering

Experimental results and complimentary computations for airfoils with vortex generators are compared in this paper, as part of an e_ort within the AVATAR project to develop tools for wind turbine blade control devices. Measurements from two airfoils equipped with passive vortex generators, a 30% thick DU97W300 and an 18% thick NTUA T18 have been used for benchmarking several simulation tools. These tools span low-to-high complexity, ranging from engineering-level integral boundary layer tools to fully-resolved computational uid dynamics codes. Results indicate that with appropriate calibration, engineering-type tools can capture the e_ects of vortex generators and outperform more complex tools. Fully resolved CFD comes at a much higher computational cost and does not necessarily capture the increased lift due to the VGs. However, in lieu of the limited experimental data available for calibration, high _delity tools are still required for assessing the e_ect of vortex generators on airfoil performance., Wind Energy, Aerodynamics, Education and Student Affairs

In this article, we endeavour to find a fast solver for finite volume discretizations for compressible unsteady viscous flows. Thereby, we concentrate on comparing the efficiency of important classes of time integration schemes, namely time adaptive Rosenbrock, singly diagonally implicit (SDIRK) and explicit first stage singly diagonally implicit Runge-Kutta (ESDIRK) methods. To make the comparison fair, efficient equation system solvers need to be chosen and a smart choice of tolerances is needed. This is determined from the tolerance TOL that steers time adaptivity. For implicit Runge-Kutta methods, the solver is given by preconditioned inexact Jacobian-free Newton-Krylov (JFNK) and for Rosenbrock, it is preconditioned Jacobian-free GMRES. To specify the tolerances in there, we suggest a simple strategy of using TOL/100 that is a good compromise between stability and computational effort. Numerical experiments for different test cases show that the fourth order Rosenbrock method RODASP and the fourth order ESDIRK method ESDIRK4 are best for fine tolerances, with RODASP being the most robust scheme., Aerodynamics, Science Centre & Programmering

The FP7 AdVanced Aerodynamic Tools for lArge Rotors - Avatar project aims to develop and validate advanced aerodynamic models, to be used in integral design codes for the next generation of large scale wind turbines (10-20MW). One of the approaches towards reaching rotors for 10-20MW size is the application of flow control devices, such as flaps. In Task 3.2: Development of aerodynamic codes for modelling of flow devices on aerofoils and, rotors of the Avatar project, aerodynamic codes are benchmarked and validated against the experimental data of a DU95W180 airfoil in steady and unsteady flow, for different angle of attack and flap settings, including unsteady oscillatory trailing-edge-flap motion, carried out within the framework of WP3: Models for Flow Devices and Flow Control, Task 3.1: CFD and Experimental Database. The aerodynamics codes are: AdaptFoil2D, Foil2W, FLOWer, MaPFlow, OpenFOAM, Q3UIC, ATEFlap. The codes include unsteady Eulerian CFD simulations with grid deformation, panel models and indicial engineering models. The validation cases correspond to 18 steady flow cases, and 42 unsteady flow cases, for varying angle of attack, flap deflection and reduced frequency, with free and forced transition. The validation of the models show varying degrees of agreement, varying between models and flow cases., Wind Energy, Aerodynamics

High computational cost associated with the high fidelity wake models such as RANS or LES serves as a primary bottleneck to perform a direct high fidelity wind farm layout optimization (WFLO) using accurate CFD based wake models. Therefore, a surrogate based multi-fidelity WFLO methodology (SWFLO) is proposed. The surrogate model is built using an SBO method referred as manifold mapping (MM). As a verification, optimization of spacing between two staggered wind turbines was performed using the proposed surrogate based methodology and the performance was compared with that of direct optimization using high fidelity model. Significant reduction in computational cost was achieved using MM: a maximum computational cost reduction of 65%, while arriving at the same optima as that of direct high fidelity optimization. The similarity between the response of models, the number of mapping points and its position, highly influences the computational efficiency of the proposed method. As a proof of concept, realistic WFLO of a small 7-turbine wind farm is performed using the proposed surrogate based methodology. Two variants of Jensen wake model with different decay coefficients were used as the fine and coarse model. The proposed SWFLO method arrived at the same optima as that of the fine model with very less number of fine model simulations., Aerodynamics, Aerospace Engineering

We present an integrated in-situ visualization approach for partitioned multi-physics simulation of fluid-structure interaction. The simulation itself is treated as a black box and only the information at the fluid-structure interface is considered, and communicated between the fluid and solid solvers with a separate coupling tool. The visualization of the interface data is performed in conjunction with the fluid solver. Furthermore, we present new visualization techniques for the analysis of the interrelation of the two solvers , with emphasis on the involved error due to discretization in space and time and the reconstruction. Our visualization approach also enables the investigation of these errors with respect of their mutual influence on the two simulation codes and their space-time discretization. For efficient interactive visualization, we employ the concept of explorable spatiotemporal images, which also enables finite-time temporal navigation in an in-situ context. We demonstrate our overall approach and its utility by means of a fluid-structure simulation using OpenFOAM that is coupled by the preCICE software layer., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

The flow around a beveled flat plate model with an asymmetric 25 degrees trailing edge with three rounding radii is analyzed using a Navier-Stokes based open source software package OpenFOAM in order to predict the aeroacoustic properties of the models. A Large Eddy Simulation with a dynamic Smagorinsky and implicit model are used as closure model for the flow solver, and are compared regarding their aeroacoustic performance. Velocity coherence and pressure correlation is determined in spanwise direction. The acoustic far field spectrum is obtained by solving Curle’s analogy in frequency domain as a post-processing step., Aerodynamics, Aerospace Engineering

Multi-physics simulations, such as fluid-structure-acoustics interaction (FSA), require a high performance computing environment in order to perform the simulation in a reasonable amount of computation time. Currently used coupling methods use a staggered execution of the fluid and solid solver [6], which leads to inherent load imbalances. In [12] a new coupling scheme based on a quasi-Newton method is proposed for fluidstructure interaction which coupled the fluid and solid solver in parallel. The quasi-Newton method requires approximately the same number of coupling iterations per time step compared to a staggered coupling approach, resulting in a better load balance when running in a parallel environment. This contribution investigates the scalability limit and load-balancing for a strongly coupled fluid-structure interaction problem, and also for a fluid-structure-acoustics interaction problem. The acoustic far field of the fluid-structure-acoustics interaction problem is loosely coupled with the flow field., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

We present an integrated in-situ visualization approach for partitioned multi-physics simulation of fluid-structure interaction. The simulation itself is treated as a black box and only the information at the fluid-structure interface is considered, and communicated between the fluid and solid solvers with a separate coupling tool. The visualization of the interface data is performed in conjunction with the fluid solver. Furthermore, we present new visualization techniques for the analysis of the interrelation of the two solvers , with emphasis on the involved error due to discretization in space and time and the reconstruction. Our visualization approach also enables the investigation of these errors with respect of their mutual influence on the two simulation codes and their space-time discretization. For efficient interactive visualization, we employ the concept of explorable spatiotemporal images, which also enables finite-time temporal navigation in an in-situ context. We demonstrate our overall approach and its utility by means of a fluid-structure simulation using OpenFOAM that is coupled by the preCICE software layer., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

Multi-physics simulations, such as fluid-structure-acoustics interaction (FSA), require a high performance computing environment in order to perform the simulation in a reasonable amount of computation time. Currently used coupling methods use a staggered execution of the fluid and solid solver [6], which leads to inherent load imbalances. In [12] a new coupling scheme based on a quasi-Newton method is proposed for fluidstructure interaction which coupled the fluid and solid solver in parallel. The quasi-Newton method requires approximately the same number of coupling iterations per time step compared to a staggered coupling approach, resulting in a better load balance when running in a parallel environment. This contribution investigates the scalability limit and load-balancing for a strongly coupled fluid-structure interaction problem, and also for a fluid-structure-acoustics interaction problem. The acoustic far field of the fluid-structure-acoustics interaction problem is loosely coupled with the flow field., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

The flow around a beveled flat plate model with an asymmetric 25 degrees trailing edge with three rounding radii is analyzed using a Navier-Stokes based open source software package OpenFOAM in order to predict the aeroacoustic properties of the models. A Large Eddy Simulation with a dynamic Smagorinsky and implicit model are used as closure model for the flow solver, and are compared regarding their aeroacoustic performance. Velocity coherence and pressure correlation is determined in spanwise direction. The acoustic far field spectrum is obtained by solving Curle’s analogy in frequency domain as a post-processing step., Aerodynamics, Aerospace Engineering

This paper is the starting point for the development of a multi-fidelity modeling approach for vortex generators (VG) arrays, where a fully resolved VG model will be coupled with an approximate model in order to improve both accuracy and flexibility without increasing the required computational cost. As a first step thereto, an analysis of the ability of the BAY-model to simulate incompressible flows around a VG on a flat plate has been performed. Results are presented for several turbulence models and using different cell selection approaches, where comparison is made with fully resolved VG results. In addition to a coarse uniform mesh, the BAY-model has been evaluated on a densely gridded mesh in order to distinguish between model and mesh related discrepancies., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

A Large Eddy Simulation (LES) with four different closure models are analyzed in OpenFOAM, an open source Computional Fluid Dynamics (CFD) package and validated for the determination of the streamwise and spanwise coherence length of the pressure field below a turbulent boundary layer at low Reynolds numbers. Matching results are found for outer scaling mean and fluctuating velocity data as well as for the pressure spectrum data. The coherence function shows a similar decay with respect to various literature studies. An exponential fit is applied to determine the coherence length. Agreement within one displacement thickness error in streamwise and spanwise direction is found for the coherence length with semi-empirical data. The spanwise coherence length is considerably smaller than the streamwise coherence length, but indicates a clear peak in the low frequency regime originating from large coherent structures with relatively small amplitudes., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

This paper is the starting point for the development of a multi-fidelity modeling approach for vortex generators (VG) arrays, where a fully resolved VG model will be coupled with an approximate model in order to improve both accuracy and flexibility without increasing the required computational cost. As a first step thereto, an analysis of the ability of the BAY-model to simulate incompressible flows around a VG on a flat plate has been performed. Results are presented for several turbulence models and using different cell selection approaches, where comparison is made with fully resolved VG results. In addition to a coarse uniform mesh, the BAY-model has been evaluated on a densely gridded mesh in order to distinguish between model and mesh related discrepancies., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

A Large Eddy Simulation (LES) with four different closure models are analyzed in OpenFOAM, an open source Computional Fluid Dynamics (CFD) package and validated for the determination of the streamwise and spanwise coherence length of the pressure field below a turbulent boundary layer at low Reynolds numbers. Matching results are found for outer scaling mean and fluctuating velocity data as well as for the pressure spectrum data. The coherence function shows a similar decay with respect to various literature studies. An exponential fit is applied to determine the coherence length. Agreement within one displacement thickness error in streamwise and spanwise direction is found for the coherence length with semi-empirical data. The spanwise coherence length is considerably smaller than the streamwise coherence length, but indicates a clear peak in the low frequency regime originating from large coherent structures with relatively small amplitudes., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

This contribution compares the efficiency of Rosenbrock time integration schemes with ESDIRK schemes, applicable to unsteady flow and fluid-structure interaction simulations. Compared to non-linear ESDIRK schemes, the linear implicit Rosenbrock-Wanner schemes require subsequent solution of the same linear systems with different right hand sides. By solving the linear systems with the iterative solver GMRES, the preconditioner can be reused for the subsequent stages of the Rosenbrock-Wanner scheme. Unsteady flow simulations show a gain in computational efficiency of approximately factor three to five in comparison with ESDIRK., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

This contribution compares the efficiency of Rosenbrock time integration schemes with ESDIRK schemes, applicable to unsteady flow and fluid-structure interaction simulations. Compared to non-linear ESDIRK schemes, the linear implicit Rosenbrock-Wanner schemes require subsequent solution of the same linear systems with different right hand sides. By solving the linear systems with the iterative solver GMRES, the preconditioner can be reused for the subsequent stages of the Rosenbrock-Wanner scheme. Unsteady flow simulations show a gain in computational efficiency of approximately factor three to five in comparison with ESDIRK., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

Thermal interaction of fluids and solids, or conjugate heat transfer (CHT), is encountered in many engineering applications. Since time-accurate computations of such coupled problems can be computationally expensive, we consider loosely-coupled and strongly- coupled solution algorithms in which higher order multi-stage Runge-Kutta schemes are employed for time integration. The higher order time integration schemes have the potential to improve the computational efficiency at arriving at a certain accuracy relative to the traditionally used 1st and 2nd order implicit schemes. The spatial coupling between the subdomains is realized using Dirichlet-Neumann interface conditions and the coupled domains are solved in a sequential manner at each stage (Block Gauss-Seidel). In this paper, the stability of two partitioned algorithms is analyzed by considering a one dimensional model problem. The model problem consists of two thermally coupled domains where the governing equation within each subdomain is unsteady linear heat conduction. In the loosely-coupled approach, a family of multi-stage IMEX schemes is used for time integration. By observing similarities between the second stage of the IMEX schemes and the ? scheme with ? = 0.5 (Crank-Nicolson), the stability of the partitioned algorithm in which the Crank-Nicolson scheme is used for time integration is first analyzed by applying the stability theory of Godunov-Ryabenkii. The stability of the IMEX schemes is next investigated by numerically solving the model problem and comparing the results to the conclusions of the stability analysis for the Crank-Nicolson scheme. Due to partly explicit nature of the IMEX schemes, the loosely-coupled algorithm becomes unstable for sufficiently large Fourier numbers (similar to the Crank-Nicolson scheme). When the ratio of the thermal effusivities of the coupled domains is much smaller than unity, time step restriction due to stability is sufficiently weak that computations can be perf, Aerdynamics, Wind Energy and Propulsion, Aerospace Engineering

The focus of this paper is on acceleration of strong partitioned coupling algorithms for fluid-structure interaction. Strong partitioned coupling requires the solution of a coupled problem at each time step during the simulation. Hereto, an interface residual is defined such that the kinematic and dynamic interface conditions on the fluid-structure interface are satisfied when it amounts to zero. Subsequently, the coupled problem is formulated as a minimization problem of the interface residual which can efficiently be performed using Newton’s method. However, Newton’s method cannot be applied when the fluid and structure solvers are considered black-boxes since the Jacobian of the interface residual is not available. For this reason, Quasi-Newton methods were developed that approximate either the Jacobian or the inverse Jacobian of the interface residual directly from input/output information. In this contribution we present a new algorithm that uses a technique from multifidelity optimization – called space-mapping – to efficiently perform the minimization of the interface residual. The space-mapping technique exploits a computationally inexpensive lowfidelity model in order to accelerate an expensive high-fidelity model using black-box information only. The space-mapping algorithm is applied to the supersonic panel flutter problem in order to demonstrate its effectiveness. The speedup – defined with respect to a Quasi-Newton algorithm – is found to be 1-1.5 for typical time step sizes. It is expected that higher speedups can be obtained when problems are considered that require strong coupling as the time step decreases, e.g. due to the added mass effect when the structure is in interaction with an incompressible fluid., Aerodynamics, Wind Energy and Propulsion, Aerospace Engineering

Noting that time-accurate computations of the unsteady incompressible Navier-Stokes (INS) equations can be computationally expensive, a family of higher order implicit multi-stage time integration schemes (namely ESDIRK) is used for advancing the solution to the unsteady INS in time. The higher order time integration schemes have the potential to decrease the computational cost of obtaining engineering levels of accuracy relative to the traditionally used 2nd order implicit schemes. The finite volume method is used for spatial discretization, and co-located arrangement of the primitive variables is considered. Furthermore, an iterated PISO algorithm is used to solve the incompressible Navier-Stokes equations. By using a temporally consistent Rhie-Chow interpolation, higher order temporal accuracy in solving the INS equations on co-located grids is achieved. For a two-dimensional lid driven cavity test case, the temporal convergence of the solution is investigated, with the third and fourth order ESDIRK schemes for time integration. The results demonstrate the temporal consistency and temporal order preservation of the algorithm., Aerodynamics, Wind Energy and Propulsion, Aerospace Engineering

Thermal interaction of fluids and solids, or conjugate heat transfer (CHT), is encountered in many engineering applications. Since time-accurate computations of such coupled problems can be computationally expensive, we consider loosely-coupled and strongly- coupled solution algorithms in which higher order multi-stage Runge-Kutta schemes are employed for time integration. The higher order time integration schemes have the potential to improve the computational efficiency at arriving at a certain accuracy relative to the traditionally used 1st and 2nd order implicit schemes. The spatial coupling between the subdomains is realized using Dirichlet-Neumann interface conditions and the coupled domains are solved in a sequential manner at each stage (Block Gauss-Seidel). In this paper, the stability of two partitioned algorithms is analyzed by considering a one dimensional model problem. The model problem consists of two thermally coupled domains where the governing equation within each subdomain is unsteady linear heat conduction. In the loosely-coupled approach, a family of multi-stage IMEX schemes is used for time integration. By observing similarities between the second stage of the IMEX schemes and the ? scheme with ? = 0.5 (Crank-Nicolson), the stability of the partitioned algorithm in which the Crank-Nicolson scheme is used for time integration is first analyzed by applying the stability theory of Godunov-Ryabenkii. The stability of the IMEX schemes is next investigated by numerically solving the model problem and comparing the results to the conclusions of the stability analysis for the Crank-Nicolson scheme. Due to partly explicit nature of the IMEX schemes, the loosely-coupled algorithm becomes unstable for sufficiently large Fourier numbers (similar to the Crank-Nicolson scheme). When the ratio of the thermal effusivities of the coupled domains is much smaller than unity, time step restriction due to stability is sufficiently weak that computations can be perf, Aerdynamics, Wind Energy and Propulsion, Aerospace Engineering

The focus of this paper is on acceleration of strong partitioned coupling algorithms for fluid-structure interaction. Strong partitioned coupling requires the solution of a coupled problem at each time step during the simulation. Hereto, an interface residual is defined such that the kinematic and dynamic interface conditions on the fluid-structure interface are satisfied when it amounts to zero. Subsequently, the coupled problem is formulated as a minimization problem of the interface residual which can efficiently be performed using Newton’s method. However, Newton’s method cannot be applied when the fluid and structure solvers are considered black-boxes since the Jacobian of the interface residual is not available. For this reason, Quasi-Newton methods were developed that approximate either the Jacobian or the inverse Jacobian of the interface residual directly from input/output information. In this contribution we present a new algorithm that uses a technique from multifidelity optimization – called space-mapping – to efficiently perform the minimization of the interface residual. The space-mapping technique exploits a computationally inexpensive lowfidelity model in order to accelerate an expensive high-fidelity model using black-box information only. The space-mapping algorithm is applied to the supersonic panel flutter problem in order to demonstrate its effectiveness. The speedup – defined with respect to a Quasi-Newton algorithm – is found to be 1-1.5 for typical time step sizes. It is expected that higher speedups can be obtained when problems are considered that require strong coupling as the time step decreases, e.g. due to the added mass effect when the structure is in interaction with an incompressible fluid., Aerodynamics, Wind Energy and Propulsion, Aerospace Engineering

Noting that time-accurate computations of the unsteady incompressible Navier-Stokes (INS) equations can be computationally expensive, a family of higher order implicit multi-stage time integration schemes (namely ESDIRK) is used for advancing the solution to the unsteady INS in time. The higher order time integration schemes have the potential to decrease the computational cost of obtaining engineering levels of accuracy relative to the traditionally used 2nd order implicit schemes. The finite volume method is used for spatial discretization, and co-located arrangement of the primitive variables is considered. Furthermore, an iterated PISO algorithm is used to solve the incompressible Navier-Stokes equations. By using a temporally consistent Rhie-Chow interpolation, higher order temporal accuracy in solving the INS equations on co-located grids is achieved. For a two-dimensional lid driven cavity test case, the temporal convergence of the solution is investigated, with the third and fourth order ESDIRK schemes for time integration. The results demonstrate the temporal consistency and temporal order preservation of the algorithm., Aerodynamics, Wind Energy and Propulsion, Aerospace Engineering

In this paper, the influence of wing flexibility and the 'clap-and-peel' motion (one versus two wings) on the upward (lifting) force in hover is studied. CFD simulations based on the ALE Navier-Stokes equations together with radial basis function mesh interpolation are used. For the 'clap-and-peel' motion an immersed symmetry plane is introduced into the computational domain. The flexible shapes of the wing are obtained from experiments. From the preliminary results it can be concluded that flexibility enhances the production of the upward force.

High fidelity analysis of fluid-structure interaction systems is often too timeconsuming when a large number of model evaluations are required. The choice for a solution procedure depends often on the efficiency of the method and the possibility of reusing existing field solvers. Aggressive Space-Mapping, a technique originally developed for multi-fidelity optimization, is applied to accelerate the partitioned solution procedure of a high fidelity fluid-structure interaction model. The method supports software modularity. Aggressive Space-Mapping (ASM) is applied to an academic testcase and the results are compared with the corresponding Incremental Quasi-Newton (IQN) method. An efficiency metric is defined to facilitate the comparison. The ASM method is found to be more efficient than the corresponding IQN method for the testcases considered. The efficiency of space-mapping increases with increasing fluid-to-structure mass ratio, indicating that the method is especially useful for strongly coupled problems., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

The application of higher order implicit time integration schemes to conjugate heat transfer problems is analyzed with Dirichlet-Neumann as the decomposition method. In the literature, only up to second order implicit time integration schemes have been reported while there is a potential for gaining computational efficiency using higher orders. For loose coupling of the domains, the IMEX scheme consisting of the ESDIRK scheme for integrating the governing equations within the subdomains and an ERK scheme for explicit integration of the explicit coupling terms is utilized. The IMEX scheme is analyzed for two cases. In one, the material properties of the coupled domains are the same and in the other they are different. While for both cases, the IMEX scheme preserves the design order of the time integration scheme, different stability and accuracy properties are observed for the two. Finally, the computational efficiency of the higher order IMEX schemes relative to the second order scheme is demonstrated using a test case in 2-D involving coupled conduction problem of three domains., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering

Two reduced order models are presented for the simulation of physically strong coupled fluid-structure interaction problem, based on computationally partitioned flow and structure solvers. The reduced order models used are a class of quasi-Newton coupling methods to obtain a stable solution and to reduce the number of subiterations. The second reduced order model is a multi-level acceleration in with coarse grid computations are performed in order to reduce computational costs. Finally an adaptive multi-level strategy is described, that contains an indicator for when to switch from coarse to fine grid level and vice versa., Aerodynamics, Wind Energy & Propulsion, Aerospace Engineering