Non-reflective boundary conditions for non-uniform mean velocity inlet and outlet flows

The numerical solution of problems in unbounded physical space requires a truncation of the computational domain to a reasonable size. As a result, the conditions on the artificial boundaries are generally unknown. Assumptions like constant pressure or velocities are only valid in the far field and lead to spurious reflections if applied on the boundaries of the truncated domain. A number of attempts have been made over the past decades to design conditions that prevent such reflections. One approach is based on characteristics. The standard analysis assumes a spatially uniform mean flow field but this is often impractical. In the present paper we show how to extend the formulation to the more general case of a non-uniform mean velocity field. A number of test cases are provided and our results compare favourably with other boundary conditions. In principle the present approach can be extended to include non-uniformities in all variables.

Sufficient conditions for zeno behavior in lagrangian hybrid systems

This paper presents easily verifiable sufficient conditions for the existence of Zenobehavior in Lagrangian hybrid systems, i.e., hybrid systems modeling mechanical systemsundergoing impacts. © 2008 Springer-Verlag Berlin Heidelberg.

Lyapunov-like conditions for the existence of zeno behavior in hybrid and lagrangian hybrid systems

Lyapunov-like conditions that utilize generalizations of energy and barrier functions certifying Zeno behavior near Zeno equilibria are presented. To better illustrate these conditions, we will study them in the context of Lagrangian hybrid systems. Through the observation that Lagrangian hybrid systems with isolated Zeno equilibria must have a onedimensional configuration space, we utilize our Lyapunov-like conditions to obtain easily verifiable necessary and sufficient conditions for the existence of Zeno behavior in systems of this form. © 2007 IEEE.

Effect of soil conditions on uplift of underground structures in liquefied soil

In an earthquake, underground structures located in liquefiable soil deposits are susceptible to floatation following an earthquake event due to their lower unit weight relative to the surrounding saturated soil. Such uplift response of the buoyant structure is influenced by the soil it is buried in. In the case of a liquefiable soil deposit, the soil can lose its shear strength significantly in the event of an earthquake. If the soil liquefies fully, the buoyant structure can float towards the soil surface. However, a partly liquefied soil deposit retains some of its initial shear strength and resists the uplift. This paper discusses the different soil conditions and their influence on the uplift response of buoyant structures. © 2012 World Scientific Publishing Company.

Dynamical structure function identifiability conditions enabling signal structure reconstruction

Networks of controlled dynamical systems exhibit a variety of interconnection patterns that could be interpreted as the structure of the system. One such interpretation of system structure is a system's signal structure, characterized as the open-loop causal dependencies among manifest variables and represented by its dynamical structure function. Although this notion of structure is among the weakest available, previous work has shown that if no a priori structural information is known about the system, not even the Boolean structure of the dynamical structure function is identifiable. Consequently, one method previously suggested for obtaining the necessary a priori structural information is to leverage knowledge about target specificity of the controlled inputs. This work extends these results to demonstrate precisely the a priori structural information that is both necessary and sufficient to reconstruct the network from input-output data. This extension is important because it significantly broadens the applicability of the identifiability conditions, enabling the design of network reconstruction experiments that were previously impossible due to practical constraints on the types of actuation mechanisms available to the engineer or scientist. The work is motivated by the proteomics problem of reconstructing the Per-Arnt-Sim Kinase pathway used in the metabolism of sugars. © 2012 IEEE.

Complex chemistry DNS of n-heptane spray autoignition at high pressure and intermediate temperature conditions

Direct Numerical Simulations (DNS) of turbulent n-heptane sprays autoigniting at high pressure (P=24bar) and intermediate air temperature (Tair=1000K) have been performed to investigate the physical mechanisms present under conditions where low-temperature chemistry is expected to be important. The initial turbulence in the carrier gas, the global equivalence ratio in the spray region, and the initial droplet size distribution of the spray were varied. Results show that spray ignition exhibits a spotty nature, with several kernels developing independently in those regions where the mixture fraction is close to its most reactive value ξMR (as determined from homogeneous reactor calculations) and the scalar dissipation rate is low. Turbulence reduces the ignition delay time as it promotes mixing between air and the fuel vapor, eventually resulting in lower values of scalar dissipation. High values of the global equivalence ratio are responsible for a larger number of ignition kernels, due to the higher probability of finding regions where ξ=ξMR. Spray polydispersity results in the occurrence of ignition over a wider range of mixture fraction values. This is a consequence of the inhomogeneities in the mixing field that characterize these sprays, where poorly mixed rich spots are seen to alternate with leaner ones which are well-mixed. The DNS simulations presented in this work have also been used to assess the applicability of the Conditional Moment Closure (CMC) method to the simulation of spray combustion. CMC is found to be a valid method for capturing spray autoignition, although care should be taken in the modelling of the unclosed terms appearing in the CMC equations. © 2013 The Combustion Institute.

Effect of operating conditions and fuel type on crevice hc emissions: Model results and comparison with experiments

A one-dimensional model for crevice HC post-flame oxidation is used to calculate and understand the effect of operating parameters and fuel type (propane and isooctane) on the extent of crevice hydrocarbon and the product distribution in the post flame environment. The calculations show that the main parameters controlling oxidation are: bulk burned gas temperatures, wall temperatures, turbulent diffusivity, and fuel oxidation rates. Calculated extents of oxidation agree well with experimental values, and the sensitivities to operating conditions (wall temperatures, equivalence ratio, fuel type) are reasonably well captured. Whereas the bulk gas temperatures largely determine the extent of oxidation, the hydrocarbon product distribution is not very much affected by the burned gas temperatures, but mostly by diffusion rates. Uncertainties in both turbulent diffusion rates as well as in mechanisms are an important factor limiting the predictive capabilities of the model. However, it seems well suited to sensitivity calculations about a baseline. Copyright © 1999 Society of Automotive Engineers, Inc.

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.

Effect of high temperature post-annealing on sidewalls of GaAs NWs grown by MOCVD

The sidewall facets of GaAs nanowires (NWs) were studied. It has been found that the sidewalls of GaAs NWs grown at 450 °C are {112} facets. However, the sidewalls of GaAs NWs start to become {110} during the postannealing at 650 °C for 30 min. © 2010 IEEE.

Off-wall boundary conditions for turbulent flows obtained from buffer-layer minimal flow units

There is strong evidence that the transport processes in the buffer region of wall-bounded turbulence are common across various flow configurations, even in the embryonic turbulence in transition (Park et al., Phys. Fl. 24). We use this premise to develop off-wall boundary conditions for turbulent simulations. Boundary conditions are constructed from DNS databases using periodic minimal flow units and reduced order modeling. The DNS data was taken from a channel at Reτ=400 and a zero-pressure gradient transitional boundary layer (Sayadi et al., submitted to J. Fluid Mech.). Both types of boundary conditions were first tested on a DNS of the core of the channel flow with the aim of extending their application to LES and to spatially evolving flows.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.