A simple non-linear global-local finite element methodology for composite structures

A simple non-linear global-local finite element methodology is presented. A global coarse model, using 2-D shell elements, is solved non-linearly and the displacements and rotations around a region of interest are applied, as displacement boundary conditions, to a refined local 3-D model using Kirchhoff plate assumptions. The global elements' shape functions are used to interpolate between nodes. The local model is then solved non-linearly with an incremental scheme independent of that used for the global model.

Numerical Study of Film Cooling Scheme on a Blunt-Nosed Body in Hypersonic Flow

In hypersonic flight, the prediction of aerodynamic heating and the construction of a proper thermal protection system (TPS) are significantly important. In this study, the method of a film cooling technique, which is already the state of the art in cooling of gas turbine engines, is proposed for a fully reusable and active TPS. Effectiveness of the film cooling scheme to reduce convective heating rates for a blunt-nosed spacecraft flying at Mach number 6.56 and 40 deg angle of attack is investigated numerically. The inflow boundary conditions used the standard values at an altitude of 30 km. The computational domain consists of infinite rows of film cooling holes on the bottom of a blunt-nosed slab. Laminar and several turbulent calculations have been performed and compared. The influence of blowing ratios on the film cooling effectiveness is investigated. The results exhibit that the film cooling technique could be an effective method for an active cooling of blunt-nosed bodies in hypersonic flows.

Substrate Clamping Effects on Irreversible Domain Wall Dynamics in Lead Zirconate Titanate Thin Films

The role of long-range strain interactions on domain wall dynamics is explored through macroscopic and local measurements of nonlinear behavior in mechanically clamped and released polycrystalline lead zirconate-titanate (PZT) films. Released films show a dramatic change in the global dielectric nonlinearity and its frequency dependence as a function of mechanical clamping. Furthermore, we observe a transition from strong clustering of the nonlinear response for the clamped case to almost uniform nonlinearity for the released film. This behavior is ascribed to increased mobility of domain walls. These results suggest the dominant role of collective strain interactions mediated by the local and global mechanical boundary conditions on the domain wall dynamics. The work presented in this Letter demonstrates that measurements on clamped films may considerably underestimate the piezoelectric coefficients and coupling constants of released structures used in microelectromechanical systems, energy harvesting systems, and microrobots.

TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE NONRESONANT, COSMIC-RAY-DRIVEN INSTABILITY IN SUPERNOVA REMNANT SHOCKS

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the "knee" at 10(15.5) eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized nonpropagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times, we observe saturation of the instability in the simulation, but the mechanism responsible is an artifact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

Development and Validation of a Compact Thermal Model for an Aircraft Compartment

The development of accurate structural/thermal numerical models of complex systems, such as aircraft fuselage barrels, is often limited and determined by the smallest scales that need to be modelled. The development of reduced order models of the smallest scales and consequently their integration with higher level models can be a way to minimise the bottle neck present, while still having efficient, robust and accurate numerical models. In this paper a methodology on how to develop compact thermal fluid models (CTFMs) for compartments where mixed convection regimes are present is demonstrated. Detailed numerical simulations (CFD) have been developed for an aircraft crown compartment and validated against experimental data obtained from a 1:1 scale compartment rig. The crown compartment is defined as the confined area between the upper fuselage and the passenger cabin in a single aisle commercial aircraft. CFD results were utilised to extract average quantities (temperature and heat fluxes) and characteristic parameters (heat transfer coefficients) to generate CTFMs. The CTFMs have then been compared with the results obtained from the detailed models showing average errors for temperature predictions lower than 5%. This error can be deemed acceptable when compared to the nominal experimental error associated with the thermocouple measurements.

The CTFMs methodology developed allows to generate accurate reduced order models where accuracy is restricted to the region of Boundary Conditions applied. This limitation arises from the sensitivity of the internal flow structures to the applied boundary condition set. CTFMs thus generated can be then integrated in complex numerical modelling of whole fuselage sections.

Further steps in the development of an exhaustive methodology would be the implementation of a logic ruled based approach to extract directly from the CFD simulations numbers and positions of the nodes for the CTFM.

Formal calibration methodology for CFD models of naturally ventilated indoor environments

Well planned natural ventilation strategies and systems in the built environments may provide healthy and comfortable indoor conditions, while contributing to a significant reduction in the energy consumed by buildings. Computational Fluid Dynamics (CFD) is particularly suited for modelling indoor conditions in naturally ventilated spaces, which are difficult to predict using other types of building simulation tools. Hence, accurate and reliable CFD models of naturally ventilated indoor spaces are necessary to support the effective design and operation of indoor environments in buildings. This paper presents a formal calibration methodology for the development of CFD models of naturally ventilated indoor environments. The methodology explains how to qualitatively and quantitatively verify and validate CFD models, including parametric analysis utilising the response surface technique to support a robust calibration process. The proposed methodology is demonstrated on a naturally ventilated study zone in the library building at the National University of Ireland in Galway. The calibration process is supported by the on-site measurements performed in a normally operating building. The measurement of outdoor weather data provided boundary conditions for the CFD model, while a network of wireless sensors supplied air speeds and air temperatures inside the room for the model calibration. The concepts and techniques developed here will enhance the process of achieving reliable CFD models that represent indoor spaces and provide new and valuable information for estimating the effect of the boundary conditions on the CFD model results in indoor environments. © 2012 Elsevier Ltd.

Comparison of lagrangian and eulerian simulations of slurry flows in a sudden expansion

From a review of technical literature, it was not apparent if the Lagrangian or the Eulerian dispersed phase modeling approach was more valid to simulate dilute erosive slurry flow. In this study, both modeling approaches were employed and a comparative analysis of performances and accuracy between the two models was carried out. Due to an impossibility to define, for the Eulerian model already implemented in FLUENT, a set of boundary conditions consistent with the Lagrangian impulsive equations, an Eulerian dispersed phase model was integrated in the FLUENT code using subroutines and user-defined scalar equations. Numerical predictions obtained from the two different approaches for two-phase flow in a sudden expansion were compared with the measured data. Excellent agreement was attained between the predicted and observed fluid and particle velocity in the axial direction and for the kinetic energy. Erosion profiles in a sudden expansion computed using the Lagrangian scheme yielded good qualitative agreement with measured data and predicted a maximum impact angle of 29 deg at the fluid reattachment point. The Eulerian model was adversely affected by the reattachment of the fluid phase to the wall and the simulated erosion profiles were not in agreement with the Lagrangian or measured data. Furthermore, the Eulerian model under-predicted the Lagrangian impact angle at all locations except the reattachment point. © 2010 American Society of Mechanical Engineers.

Numerical study of film cooling scheme on a blunt-nosed body in hypersonic flow

In hypersonic flights, the prediction of aerodynamic heating and the construction of a proper thermal protection system (TPS) are significantly important. In this study, the method of a film cooling technique, which is already the state of the art in cooling gas turbine engine, is proposed for a fully reusable and active TPS. Effectiveness of the film cooling scheme to reduce convective heating rates for a blunt nosed spacecraft flying at Mach number 6.56 and 40 degree angle of attack is investigated numerically. The inflow boundary conditions used the standard values at an altitude of 30 km. Computational domain consists of infinite rows of film cooling holes on the bottom of a blunt-nosed slab. Laminar and several turbulent calculations have been performed and compared each other. The influence of blowing ratios on the film cooling effectiveness is investigated. The results exhibit that the film cooling technique could be an effective method for an active cooling of blunt-nosed bodies in hypersonic flows.

When Casimir meets Kibble-Zurek

Verification of the dynamical Casimir effect (DCE) in optical systems is still elusive due to the very demanding requirements for its experimental implementation. This typically requires very fast changes in the boundary conditions of the problem. We show that an ensemble of two-level atoms collectively coupled to the electromagnetic field of a cavity, driven at low frequencies and close to a quantum phase transition, stimulates the production of photons from the vacuum. This paves the way for an effective simulation of the DCE through a mechanism that has recently found experimental demonstration. The spectral properties of the emitted radiation reflect the critical nature of the system and allow us to link the detection of DCE to the Kibble-Zurek mechanism for the production of defects when crossing a continuous phase transition.

Numerical modelling of laboratory tests for TBM cutter indentation process analysis

Considering that TBMs are nowadays used for long Trans-Alpine tunnels, the
understanding of rock breaking and chipping due to TBM cutter disks mechanism, for deep tunnelling operations, becomes very interesting. In this paper, the results from carried out laboratory tests that simulate the disk cutter action at the rock tunnel face by means of an indentation tool, acting on a rock
specimen with proper size, and the related three-dimensional and two-dimensional numerical modelling are proposed. The developed numerical models simulate the different test conditions (applied load, boundary conditions) allowing the analysis of the stresses distributions along possible breaking planes.
The influence of a confinement-free area on one side of the specimen, simulating the formation of a groove near the tool, is pointed out.
The obtained results from numerical modelling put in evidence a satisfactory agreement with the experimental observations.

Stiffness analysis of polymeric composites using the finite element method

The ability to predict the mechanical behavior of polymer composites is crucial for their design and manufacture. Extensive studies based on both macro- and micromechanical analyses are used to develop new insights into the behavior of composites. In this respect, finite element modeling has proved to be a particularly powerful tool. In this article, we present a Galerkin scheme in conjunction with the penalty method for elasticity analyses of different types of polymer composites. In this scheme, the application of Green's theorem to the model equation results in the appearance of interfacial flux terms along the boundary between the filler and polymer matrix. It is shown that for some types of composites these terms significantly affect the stress transfer between polymer and fillers. Thus, inclusion of these terms in the working equations of the scheme preserves the accuracy of the model predictions. The model is used to predict the most important bulk property of different types of composites. Composites filled with rigid or soft particles, and composites reinforced with short or continuous fibers are investigated. For each case, the results are compared with the available experimental results and data obtained from other models reported in the literature. Effects of assumptions made in the development of the model and the selection of the prescribed boundary conditions are discussed.

Preliminary analysis of organic Rankine cycles to improve vehicle efficiency

This paper presents the background rationale and key findings for a model-based study of supercritical waste heat recovery organic Rankine cycles. The paper’s objective is to cover the necessary groundwork to facilitate the future operation of a thermodynamic organic Rankine cycle model under realistic thermodynamic boundary conditions for performance optimisation of organic Rankine cycles. This involves determining the type of power cycle for organic Rankine cycles, the circuit configuration and suitable boundary conditions. The study focuses on multiple heat sources from vehicles but the findings are generally applicable, with careful consideration, to any waste heat recovery system. This paper introduces waste heat recovery and discusses the general merits of organic fluids versus water and supercritical operation versus subcritical operation from a theoretical perspective and, where possible, from a practical perspective. The benefits of regeneration are investigated from an efficiency perspective for selected subcritical and supercritical conditions. A simulation model is described with an introduction to some general Rankine cycle boundary conditions. The paper describes the analysis of real hybrid vehicle data from several driving cycles and its manipulation to represent the thermal inertia for model heat input boundary conditions. Basic theory suggests that selecting the operating pressures and temperatures to maximise the Rankine cycle performance is relatively straightforward. However, it was found that this may not be the case for an organic Rankine cycle operating in a vehicle. When operating in a driving cycle, the available heat and its quality can vary with the power output and between heat sources. For example, the available coolant heat does not vary much with the load, whereas the quantity and quality of the exhaust heat varies considerably. The key objective for operation in the vehicle is optimum utilisation of the available heat by delivering the maximum work out. The fluid selection process and the presentation and analysis of the final results of the simulation work on organic Rankine cycles are the subjects of two future publications.

The effects of global climate change on the cycling and processes of persistent organic pollutants (POPs) in the North Sea

The fate and cycling of two selected legacy persistent organic pollutants (POPs), PCB 153 and gamma-HCH, in the North Sea in the 21st century have been modelled with combined hydrodynamic and fate and transport ocean models
(HAMSOM and FANTOM, respectively). To investigate the impact of climate variability on POPs in the North Sea in the 21st century, future scenario model runs for three 10-year periods to the year 2100 using plausible levels of both in
situ concentrations and atmospheric, river and open boundary inputs are performed. This slice mode under a moderate scenario (A1B) is sufficient to provide a basis for further analysis. For the HAMSOM and atmospheric forcing, results of the IPCC A1B (SRES) 21st century scenario are utilized, where surface forcing is provided by the REMO downscaling of the ECHAM5 global atmospheric model, and open boundary conditions are provided by the MPIOM global ocean model.
Dry gas deposition and volatilization of gamma-HCH increase in the future relative to the present by up to 20% (in the spring and summer months for deposition and in summer for volatilization). In the water column, total mass of
gamma-HCH and PCB 153 remain fairly steady in all three runs. In sediment,
gamma-HCH increases in the future runs, relative to the present, while PCB 153 in sediment decreases exponentially in all three runs, but even faster in the future, due to the increased number of storms, increased duration of gale wind conditions and increased water and air temperatures, all of which are the result of climate change. Annual net sinks exceed sources at the ends of all periods.
Overall, the model results indicate that the climate change scenarios considered here generally have a negligible influence on the simulated fate and transport of the two POPs in the North Sea, although the increased number and magnitude of storms in the 21st century will result in POP resuspension and ensuing revolatilization events. Trends in emissions from primary and secondary sources will remain the key driver of levels of these contaminants over time.

On the nudging terms at open boundaries in regional ocean models

For open boundary conditions (OBCs) in regional models, a nudging term added to radiative and/or advective conditions during the wave or flow propagation outward from the model domain of interest is widely used, to prevent the predicted boundary values from evolving to become quite different from the external data, especially for a long-term integration. However, nudging time scales are basically unknown, leading to many empirical selections. In this paper, a method for objectively estimating nudging time scales during outward propagation is proposed, by using internal model dynamics near the boundary. We tested this method and other several commonly used OBCs for cases of both an idealized model domain and a realistic configuration, and model results demonstrated that the proposed method improves the model solutions. Many similarities are found between the nudging and mixing time scales, in magnitude, spatial and temporal variations, since the nudging mainly replaces the effect of the mixing terms in this study. However, the mixing time scale is not an intrinsic property of the nudging term because in other studies the nudging term might replace terms other than the mixing terms and, thus, should reflect other characteristic time scales.

Using mesh-geometry relationships to transfer analysis models between CAE tools

Integrating analysis and design models is a complex task due to differences between the models and the architectures of the toolsets used to create them. This complexity is increased with the use of many different tools for specific tasks using an analysis process. In this work various design and analysis models are linked throughout the design lifecycle, allowing them to be moved between packages in a way not currently available. Three technologies named Cellular Modeling, Virtual Topology and Equivalencing are combined to demonstrate how different finite element meshes generated on abstract analysis geometries can be linked to their original geometry. Cellular models allow interfaces between adjacent cells to be extracted and exploited to transfer analysis attributes such as mesh associativity or boundary conditions between equivalent model representations. Virtual Topology descriptions used for geometry clean-up operations are explicitly stored so they can be reused by downstream applications. Establishing the equivalence relationships between models enables analysts to utilize multiple packages for specialist tasks without worrying about compatibility issues or substantial rework.