Numerical modeling of deep soil mixing excavation support

The finite element method (FEM) is growing in popularity over the pressure diagram/hand calculation method for analysis of excavation systems in general and deep soil mixing excavations in particular. In this paper, a finite element analysis is used to study the behavior of a deep mixed excavation. Through the use of Plaxis (a FEM software program), the construction sequence is simulated by following the various construction phases allowing for deflections due to strut or anchor installation to be predicted. The numerical model used in this study simulates the soil cement columns as a continuous wall matching the bending stiffness of the actual wall. Input parameters based on laboratory tests and modeling assumptions are discussed. An example of the approach is illustrated using the Islais Creek Transport/Storage Project in San Francisco, California. Copyright ASCE 2006.

Simulations of thermal upgrading methods for oil shale reservoirs

This paper analyzes reaction and thermal front development in porous reservoirs with reacting flows, such as those encountered in shale oil extraction. A set of dimensionless parameters and a 3D code are developed in order to investigate the important physical and chemical variables of such reservoirs when heated by in situ methods. This contribution builds on a 1D model developed for the precursor study to this work. Theory necessary for this study is presented, namely shale decomposition chemical mechanisms, governing equations for multiphase flow in porous media and necessary closure models. Plotting the ratio of the thermal wave speed to the fluid speed allows one to infer that the reaction wave front ends where this ratio is at a minimum. The reaction front follows the thermal front closely, thus allowing assumptions to be made about the extent of decomposition solely by looking at thermal wave progression. Furthermore, this sensitivity analysis showed that a certain minimum permeability is required in order to ensure the formation of a traveling thermal wave. It was found that by studying the non-dimensional governing parameters of the system one can ascribe characteristic values for these parameters for given initial and boundary conditions. This allows one to roughly predict the performance of a particular method on a particular reservoir given approximate values for initial and boundary conditions. Channelling and flow blockage due to carbon residue buildup impeded each method's performance. Blockage was found to be a result of imbalanced heating. Copyright 2012, Society of Petroleum Engineers.

Matching an oscillator model to a phase response curve

The Phase Response Curve (PRC) has proven a useful tool for the reduction of complex oscillator models. It is also an information often experimentally available to the biologist. This paper introduces a numerical tool based on the sensitivity analysis of the PRC to adapt initial model parameters in order to match a particular PRC shape. We illustrate the approach on a simple biochemical model of circadian oscillator. © 2011 IEEE.

BWR fuel design options for self-sustainable Th-233U fuel cycle

In this work, we investigate a number of fuel assembly design options for a BWR core operating in a closed self-sustainable Th-233U fuel cycle. The designs rely on axially heterogeneous fuel assembly structure in order to improve fertile to fissile conversion ratio. One of the main assumptions of the current study was to restrict the fuel assembly geometry to a single axial fissile zone "sandwiched" between two fertile blanket zones. The main objective was to study the effect of the most important design parameters, such as dimensions of fissile and fertile zones and average void fraction, on the net breeding of 233U. The main design challenge in this respect is that the fuel breeding potential is at odds with axial power peaking and therefore limits the maximum achievable core power rating. The calculations were performed with BGCore system, which consists of MCNP code coupled with fuel depletion and thermo-hydraulic feedback modules. A single 3-dimensional fuel assembly with reflective radial boundaries was modeled applying simplified restrictions on maximum central line fuel temperature and Critical Power Ratio. It was found that axially heterogeneous fuel assembly design with single fissile zone can potentially achieve net breeding. In this case however, the achievable core power density is roughly one third of the reference BWR core.

Neutronic optimization in high conversion Th-233U fuel assembly with simulated annealing

This paper reports on fuel design optimization of a PWR operating in a self sustainable Th-233U fuel cycle. Monte Carlo simulated annealing method was used in order to identify the fuel assembly configuration with the most attractive breeding performance. In previous studies, it was shown that breeding may be achieved by employing heterogeneous Seed-Blanket fuel geometry. The arrangement of seed and blanket pins within the assemblies may be determined by varying the designed parameters based on basic reactor physics phenomena which affect breeding. However, the amount of free parameters may still prove to be prohibitively large in order to systematically explore the design space for optimal solution. Therefore, the Monte Carlo annealing algorithm for neutronic optimization is applied in order to identify the most favorable design. The objective of simulated annealing optimization is to find a set of design parameters, which maximizes some given performance function (such as relative period of net breeding) under specified constraints (such as fuel cycle length). The first objective of the study was to demonstrate that the simulated annealing optimization algorithm will lead to the same fuel pins arrangement as was obtained in the previous studies which used only basic physics phenomena as guidance for optimization. In the second part of this work, the simulated annealing method was used to optimize fuel pins arrangement in much larger fuel assembly, where the basic physics intuition does not yield clearly optimal configuration. The simulated annealing method was found to be very efficient in selecting the optimal design in both cases. In the future, this method will be used for optimization of fuel assembly design with larger number of free parameters in order to determine the most favorable trade-off between the breeding performance and core average power density.

Optimization of burnable poison design for Pu incineration in fully fertile free PWR core

The design challenges of the fertile-free based fuel (FFF) can be addressed by careful and elaborate use of burnable poisons (BP). Practical fully FFF core design for PWR reactor has been reported in the past [1]. However, the burnable poison option used in the design resulted in significant end of cycle reactivity penalty due to incomplete BP depletion. Consequently, excessive Pu loading were required to maintain the target fuel cycle length, which in turn decreased the Pu burning efficiency. A systematic evaluation of commercially available BP materials in all configurations currently used in PWRs is the main objective of this work. The BP materials considered are Boron, Gd, Er, and Hf. The BP geometries were based on Wet Annular Burnable Absorber (WABA), Integral Fuel Burnable Absorber (IFBA), and Homogeneous poison/fuel mixtures. Several most promising combinations of BP designs were selected for the full core 3D simulation. All major core performance parameters for the analyzed cases are very close to those of a standard PWR with conventional UO2 fuel including possibility of reactivity control, power peaking factors, and cycle length. The MTC of all FFF cores was found at the full power conditions at all times and very close to that of the UO2 core. The Doppler coefficient of the FFF cores is also negative but somewhat lower in magnitude compared to UO2 core. The soluble boron worth of the FFF cores was calculated to be lower than that of the UO2 core by about a factor of two, which still allows the core reactivity control with acceptable soluble boron concentrations. The main conclusion of this work is that judicial application of burnable poisons for fertile free fuel has a potential to produce a core design with performance characteristics close to those of the reference PWR core with conventional UO2 fuel.

Consideration of reactivity initiated accidents for TRU recycling in LWRs

Response of a PWR core loaded with Combined Non-Fertile and UO2 (CONFU) fuel assemblies to control rod ejection accident was evaluated. A number of core arrangements and TRU fuel compositions were considered and the results were compared with the performance of a reference all-UO2 core. The comparison was based on the results of a simple point kinetics model with thermal reactivity feedbacks and temperature dependant materials properties. The reactivity coefficients and core average kinetics parameters were obtained from the full core 3D neutronic simulations. The results show that application of the CONFU assembly concept causes only minor deviation of fuel performance characteristics in reactivity initiated accidents. This is a consequence of relatively small loadings of TRU in the CONFU assembly and therefore dominating role of conventional UO2 fuel in the neutronic performance of the core.

Sensitivity analysis and optimization of the nuclear fuel cycle

A sensitivity study has been conducted to assess the robustness of the conclusions presented in the MIT Fuel Cycle Study. The Once Through Cycle (OTC) is considered as the base-line case, while advanced technologies with fuel recycling characterize the alternative fuel cycles. The options include limited recycling in LWRs and full recycling in fast reactors and in high conversion LWRs. Fast reactor technologies studied include both oxide and metal fueled reactors. The analysis allowed optimization of the fast reactor conversion ratio with respect to desired fuel cycle performance characteristics. The following parameters were found to significantly affect the performance of recycling technologies and their penetration over time: Capacity Factors of the fuel cycle facilities, Spent Fuel Cooling Time, Thermal Reprocessing Introduction Date, and incore and Out-of-core TRU Inventory Requirements for recycling technology. An optimization scheme of the nuclear fuel cycle is proposed. Optimization criteria and metrics of interest for different stakeholders in the fuel cycle (economics, waste management, environmental impact, etc.) are utilized for two different optimization techniques (linear and stochastic). Preliminary results covering single and multi-variable and single and multi-objective optimization demonstrate the viability of the optimization scheme.

Phase perfection in zinc Blende and Wurtzite III-V nanowires using basic growth parameters.

Controlling the crystallographic phase purity of III-V nanowires is notoriously difficult, yet this is essential for future nanowire devices. Reported methods for controlling nanowire phase require dopant addition, or a restricted choice of nanowire diameter, and only rarely yield a pure phase. Here we demonstrate that phase-perfect nanowires, of arbitrary diameter, can be achieved simply by tailoring basic growth parameters: temperature and V/III ratio. Phase purity is achieved without sacrificing important specifications of diameter and dopant levels. Pure zinc blende nanowires, free of twin defects, were achieved using a low growth temperature coupled with a high V/III ratio. Conversely, a high growth temperature coupled with a low V/III ratio produced pure wurtzite nanowires free of stacking faults. We present a comprehensive nucleation model to explain the formation of these markedly different crystal phases under these growth conditions. Critical to achieving phase purity are changes in surface energy of the nanowire side facets, which in turn are controlled by the basic growth parameters of temperature and V/III ratio. This ability to tune crystal structure between twin-free zinc blende and stacking-fault-free wurtzite not only will enhance the performance of nanowire devices but also opens new possibilities for engineering nanowire devices, without restrictions on nanowire diameters or doping.

Cement bentonite cutoff walls for polluted sites

Cement-bentonite (CB) cutoff walls have long been used to control ground water flow and contaminant migration at polluted sites. Hydraulic conductivity and unconfined compressive strength are two short-term properties often used by industry and owners in CB specification and are important parameters discussed in this paper. For polluted sites, long-term compatibility is also an important issue. These properties are coupled to a number of external factors including the mix design, construction sequence, presence/absence of contaminants at the site. Additional short-term properties for engineering assessment include the stressstrain characteristics in both drained and undrained shear in both with and without confinement as well as one-dimensional consolidation properties. Long term CB properties are affected by aging, reaction chemistry, drying, in situ stress state, and interaction with the polluted environment. © 2013 Taylor & Francis Group.

Parallel coordinates in computational engineering design

Modern Engineering Design involves the deployment of many computational tools. Re- search on challenging real-world design problems is focused on developing improvements for the engineering design process through the integration and application of advanced com- putational search/optimization and analysis tools. Successful application of these methods generates vast quantities of data on potential optimum designs. To gain maximum value from the optimization process, designers need to visualise and interpret this information leading to better understanding of the complex and multimodal relations between param- eters, objectives and decision-making of multiple and strongly conflicting criteria. Initial work by the authors has identified that the Parallel Coordinates interactive visualisation method has considerable potential in this regard. This methodology involves significant levels of user-interaction, making the engineering designer central to the process, rather than the passive recipient of a deluge of pre-formatted information. In the present work we have applied and demonstrated this methodology in two differ- ent aerodynamic turbomachinery design cases; a detailed 3D shape design for compressor blades, and a preliminary mean-line design for the whole compressor core. The first case comprises 26 design parameters for the parameterisation of the blade geometry, and we analysed the data produced from a three-objective optimization study, thus describing a design space with 29 dimensions. The latter case comprises 45 design parameters and two objective functions, hence developing a design space with 47 dimensions. In both cases the dimensionality can be managed quite easily in Parallel Coordinates space, and most importantly, we are able to identify interesting and crucial aspects of the relationships between the design parameters and optimum level of the objective functions under con- sideration. These findings guide the human designer to find answers to questions that could not even be addressed before. In this way, understanding the design leads to more intelligent decision-making and design space exploration. © 2012 AIAA.

Comparison of Micro-Vortex generators in supersonic flows

An experimental comparison of several vortex generator geometries was conducted at Mach 1.5, 1.8, and 2.5 to better understand downstream vortex development as a function of device shape and Mach number. The devices had heights less than that of the boundary-layer ("micro"-vortex generators) and were either vane-shaped or of the alternative microramp geometry. LDV was used to measure two components of velocity at several stations downstream of the devices. The velocity data were then fitted to a vortex model so that vortex parameters such as circulation, core radius, and trajectory were estimated. Mach number dependence was seen for all parameters. Vortex height and core radius both tended to decrease slightly with increasing Mach number. A critical vane angle for maximum circulation was observed and also decreased with increasing Mach number. Circulation was seen to scale with wall-friction velocity for Mach 1.5 and 1.8 but not 2.5. © 2012 by W.R. Nolan and H. Babinsky.

Kelvin-Helmholtz-like instabilities in turbulent flows over complex surfaces

Riblets are small surface protrusions aligned with the flow direction, which confer an anisotropic roughness to the surface [6]. We have recently reported that the transitional-roughness effect in riblets, which limits their performance, is due to a Kelvin–Helmholtz-like instability of the overlying mean flow [7]. According to our DNSs, the instability sets on as the Reynolds number based on the roughness size of the riblets increases, and coherent, elongated spanwise vortices begin to develop immediately above the riblet tips, causing the degradation of the drag-reduction effect. This is a very novel concept, since prior studies had proposed that the degradation was due to the interaction of riblets with the flow as independent units, either to the lodging of quasi-streamwise vortices in the surface grooves [2] or to the shedding of secondary streamwise vorticity at the riblet peaks [9]. We have proposed an approximate inviscid analysis for the instability, in which the presence of riblets is modelled through an average boundary condition for an overlying, spanwise-independent mean flow. This simplification lacks the accuracy of an exact analysis [4], but in turn applies to riblet surfaces in general. Our analysis succeeds in predicting the riblet size for the onset of the instability, while qualitatively reproducing the wavelengths and shapes of the spanwise structures observed in the DNSs. The analysis also connects the observations with the Kelvin–Helmholtz instability of mixing layers. The fundamental riblet length scale for the onset of the instability is a ‘penetration length,’ which reflects how easily the perturbation flow moves through the riblet grooves. This result is in excellent agreement with the available experimental evidence, and has enabled the identification of the key geometric parameters to delay the breakdown. Although the appearance of elongated spanwise vortices was unexpected in the case of riblets, similar phenomena had already been observed over other rough [3], porous [1] and permeable [11] surfaces, as well as over plant [5,14] and urban [12] canopies, both in the transitional and in the fully-rough regimes. However, the theoretical analyses that support the connection of these observations with the Kelvin–Helmholtz instability are somewhat scarce [7, 11, 13]. It has been recently proposed that Kelvin–Helmholtz-like instabilities are a dominant feature common to “obstructed” shear flows [8]. It is interesting that the instability does not require an inflection point to develop, as is often claimed in the literature. The Kelvin-Helmholtz rollers are rather triggered by the apparent wall-normal-transpiration ability of the flow at the plane immediately above the obstructing elements [7,11]. Although both conditions are generally complementary, if wall-normal transpiration is not present the spanwise vortices may not develop, even if an inflection point exists within the roughness [10]. REFERENCES [1] Breugem, W. P., Boersma, B. J. & Uittenbogaard, R. E. 2006 J. Fluid Mech. 562, 35–72. [2] Choi, H., Moin, P. & Kim, J. 1993 J. Fluid Mech. 255, 503–539. [3] Coceal, O., Dobre, A., Thomas, T. G. & Belcher, S. E. 2007 J. Fluid Mech. 589, 375–409. [4] Ehrenstein, U. 2009 Phys. Fluids 8, 3194–3196. [5] Finnigan, J. 2000 Ann. Rev. Fluid Mech. 32, 519–571. [6] Garcia-Mayoral, R. & Jimenez, J. 2011 Phil. Trans. R. Soc. A 369, 1412–1427. [7] Garcia-Mayoral, R. & Jimenez, J. 2011 J. Fluid Mech. doi: 10.1017/jfm.2011.114. [8] Ghisalberti, M. 2009 J. Fluid Mech. 641, 51–61. [9] Goldstein, D. B. & Tuan, T. C. 1998 J. Fluid Mech. 363, 115–151. [10] Hahn, S., Je, J. & Choi, H. 2002 J. Fluid Mech. 450, 259–285. [11] Jimenez, J., Uhlman, M., Pinelli, A. & G., K. 2001 J. Fluid Mech. 442, 89–117. [12] Letzel, M. O., Krane, M. & Raasch, S. 2008 Atmos. Environ. 42, 8770–8784. [13] Py, C., de Langre, E. & Moulia, B. 2006 J. Fluid Mech. 568, 425–449. [14] Raupach, M. R., Finnigan, J. & Brunet, Y. 1996 Boundary-Layer Meteorol. 78, 351–382.

A canonical normal SBLI flow relevant to external compression inlets

The normal shock wave / boundary layer interaction (normal SBLI) is important to the operation and performance of a supersonic inlet, and the normal SBLI is particularly prominent in external compression inlets. To improve our understanding of such interactions, it is helpful to make use of fundamental flows which capture the main elements of inlets, without resorting to the level of complexity and system integration associated with full-geometry inlets. In this paper, several fundamental fiow-fleld configurations have been considered as possible test cases to represent the normal SBLI aspects found in typical external compression inlets, and it was found that the spillage-diffuser more closely retains the basic flow features of an external compression inlet than the other configurations. In particular, this flow-fleld allows the normal shock Mach number as well as the amount and rate of subsonic diffusion to be all held approximately constant mid independent of the application of flow control. In addition, a survey of several external compression inlets was conducted to quantify the flow and geometric parameters of the spillage-diffuser relevant to actual inlets. The results indicated that such a flow may be especially relevant if the terminal Mach number is about 1.3 to 1.4, the confinement parameter is around 10%, the width around twice or three times the height, and with the area expansion just downstream of the shock on the conservative side of the stall limit for incompressible diffusers. © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

An investigation of higher-order multi-objective optimisation for 3D aerodynamic shape design

We investigate the performance of different variants of a suitably tailored Tabu Search optimisation algorithm on a higher-order design problem. We consider four objective func- tions to describe the performance of a compressor stator row, subject to a number of equality and inequality constraints. The same design problem has been previously in- vestigated through single-, bi- and three-objective optimisation studies. However, in this study we explore the capabilities of enhanced variants of our Multi-objective Tabu Search (MOTS) optimisation algorithm in the context of detailed 3D aerodynamic shape design. It is shown that with these enhancements to the local search of the MOTS algorithm we can achieve a rapid exploration of complicated design spaces, but there is a trade-off be- tween speed and the quality of the trade-off surface found. Rapidly explored design spaces reveal the extremes of the objective functions, but the compromise optimum areas are not very well explored. However, there are ways to adapt the behaviour of the optimiser and maintain both a very efficient rate of progress towards the global optimum Pareto front and a healthy number of design configurations lying on the trade-off surface and exploring the compromise optimum regions. These compromise solutions almost always represent the best qualitative balance between the objectives under consideration. Such enhancements to the effectiveness of design space exploration make engineering design optimisation with multiple objectives and robustness criteria ever more practicable and attractive for modern advanced engineering design. Finally, new research questions are addressed that highlight the trade-offs between intelligence in optimisation algorithms and acquisition of qualita- tive information through computational engineering design processes that reveal patterns and relations between design parameters and objective functions, but also speed versus optimum quality. © 2012 AIAA.