Using Level Sets as the basis for a scalable, parallel geometry engine and mesh generation system

The background to this review paper is research we have performed over recent years aimed at developing a simulation system capable of handling large scale, real world applications implemented in an end-to-end parallel, scalable manner. The particular focus of this paper is the use of a Level Set solid modeling geometry kernel within this parallel framework to enable automated design optimization without topological restrictions and on geometries of arbitrary complexity. Also described is another interesting application of Level Sets: their use in guiding the export of a body-conformal mesh from our basic cut-Cartesian background octree - mesh - this permits third party flow solvers to be deployed. As a practical demonstrations meshes of guaranteed quality are generated and flow-solved for a B747 in full landing configuration and an automated optimization is performed on a cooled turbine tip geometry. Copyright © 2009 by W.N.Dawes.

Evaluating the impacts of iteration on PD processes by transforming task network models into System Dynamics models

Iteration is unavoidable in the design process and should be incorporated when planning and managing projects in order to minimize surprises and reduce schedule distortions. However, planning and managing iteration is challenging because the relationships between its causes and effects are complex. Most approaches which use mathematical models to analyze the impact of iteration on the design process focus on a relatively small number of its causes and effects. Therefore, insights derived from these analytical models may not be robust under a broader consideration of potential influencing factors. In this article, we synthesize an explanatory framework which describes the network of causes and effects of iteration identified from the literature, and introduce an analytic approach which combines a task network modeling approach with System Dynamics simulation. Our approach models the network of causes and effects of iteration alongside the process architecture which is required to analyze the impact of iteration on design process performance. We show how this allows managers to assess the impact of changes to process architecture and to management levers which influence iterative behavior, accounting for the fact that these changes can occur simultaneously and can accumulate in non-linear ways. We also discuss how the insights resulting from this analysis can be visualized for easier consumption by project participants not familiar with simulation methods. Copyright © 2010 by ASME.

Efficient truncation scheme for modeling deepwater mooring lines

This paper is aimed at enabling the confident use of existing model test facilities for ultra deepwater application without having to compromise on the widely accepted range of scales currently used by the floating production industry. Passive line truncation has traditionally been the preferred method of creating an equivalent numerical model at reduced depth; however, these techniques tend to suffer in capturing accurately line dynamic response and so reproducing peak tensions. In an attempt to improve credibility of model test data the proposed truncation procedure sets up the truncated model, based on line dynamic response rather than quasi-static system stiffness. The upper sections of each line are modeled in detail, capturing the wave action zone and all coupling effects with the vessel. These terminate to an approximate analytical model that aims to simulate the remainder of the line. Stages 1 & 2 are used to derive a water depth truncation ratio. Here vibration decay of transverse elastic waves is assessed and it is found that below a certain length criterion, the transverse vibrational characteristics for each line are inertia driven, hence with respect to these motions the truncated model can assume a linear damper whose coefficient depends on the local line properties and vibration frequency. Stage 3 endeavors to match the individual line stiffness between the full depth and truncated models. In deepwater it is likely that taut polyester moorings will be used which are predominantly straight and have high axial stiffness that provides the principal restoring force to static and low frequency vessel motions. Consequently, it means that the natural frequencies of axial vibrations are above the typical wave frequency range allowing for a quasi-static solution. In cases of exceptionally large wave frequency vessel motions, localized curvature at the chain seabed segment and tangential skin drag on the polyester rope can increase dynamic peak tensions considerably. The focus of this paper is to develop an efficient scheme based on analytic formulation, for replicating these forces at the truncation. The paper will close with an example case study of a single mooring under extreme conditions that replicates exactly the static and dynamic characteristics of the full depth line. Copyright © 2012 by the International Society of Offshore and Polar Engineers (ISOPE).

3D modeling of high-T c superconductors by finite element software

A three-dimensional (3D) numerical model is proposed to solve the electromagnetic problems involving transport current and background field of a high-T c superconducting (HTS) system. The model is characterized by the E-J power law and H-formulation, and is successfully implemented using finite element software. We first discuss the model in detail, including the mesh methods, boundary conditions and computing time. To validate the 3D model, we calculate the ac loss and trapped field solution for a bulk material and compare the results with the previously verified 2D solutions and an analytical solution. We then apply our model to test some typical problems such as superconducting bulk array and twisted conductors, which cannot be tackled by the 2D models. The new 3D model could be a powerful tool for researchers and engineers to investigate problems with a greater level of complicity.

FPGA, physics-based modeling of IGBT and pin diode for hardware co-simulation of complex power electronic converters and systems

This paper presents the steps and the challenges for implementing analytical, physics-based models for the insulated gate bipolar transistor (IGBT) and the PIN diode in hardware and more specifically in field programmable gate arrays (FPGAs). The models can be utilised in hardware co-simulation of complex power electronic converters and entire power systems in order to reduce the simulation time without compromising the accuracy of results. Such a co-simulation allows reliable prediction of the system's performance as well as accurate investigation of the power devices' behaviour during operation. Ultimately, this will allow application-specific optimisation of the devices' structure, circuit topologies as well as enhancement of the control and/or protection schemes.

Development of a new apparatus for modeling deep excavation related problems in geotechnical centrifuge

Underground space is commonly exploited both to maximise the utility of costly land in urban development and to reduce the vertical load acting on the ground. Deep excavations are carried out to construct various types of underground infrastructure such as deep basements, subways and service tunnels. Although the soil response to excavation is known in principle, designers lack practical calculation methods for predicting both short- and long-term ground movements. As the understanding of how soil behaves around an excavation in both the short and long term is insufficient and usually empirical, the judgements used in design are also empirical and serious accidents are common. To gain a better understanding of the mechanisms involved in soil excavation, a new apparatus for the centrifuge model testing of deep excavations in soft clay has been developed. This apparatus simulates the field construction sequence of a multi-propped retaining wall during centrifuge flight. A comparison is given between the new technique and the previously used method of draining heavy fluid to simulate excavation in a centrifuge model. The new system has the benefit of giving the correct initial ground conditions before excavation and the proper earth pressure distribution on the retaining structures during excavation, whereas heavy fluid only gives an earth pressure coefficient of unity and is unable to capture any changes in the earth pressure coefficient of soil inside the zone of excavation, for example owing to wall movements. Settlements of the ground surface, changes in pore water pressure, variations in earth pressure, prop forces and bending moments in the retaining wall are all monitored during excavation. Furthermore, digital images taken of a cross-section during the test are analysed using particle image velocimetry to illustrate ground deformation and soil–structure interaction mechanisms. The significance of these observations is discussed.

Modeling natural sounds with modulation cascade processes

Natural sounds are structured on many time-scales. A typical segment of speech, for example, contains features that span four orders of magnitude: Sentences ($\sim1$s); phonemes ($\sim10$−$1$ s); glottal pulses ($\sim 10$−$2$s); and formants ($\sim 10$−$3$s). The auditory system uses information from each of these time-scales to solve complicated tasks such as auditory scene analysis [1]. One route toward understanding how auditory processing accomplishes this analysis is to build neuroscience-inspired algorithms which solve similar tasks and to compare the properties of these algorithms with properties of auditory processing. There is however a discord: Current machine-audition algorithms largely concentrate on the shorter time-scale structures in sounds, and the longer structures are ignored. The reason for this is two-fold. Firstly, it is a difficult technical problem to construct an algorithm that utilises both sorts of information. Secondly, it is computationally demanding to simultaneously process data both at high resolution (to extract short temporal information) and for long duration (to extract long temporal information). The contribution of this work is to develop a new statistical model for natural sounds that captures structure across a wide range of time-scales, and to provide efficient learning and inference algorithms. We demonstrate the success of this approach on a missing data task.

Modeling nonlinearities in MEMS oscillators

We present a mathematical model of a microelectromechanical system (MEMS) oscillator that integrates the nonlinearities of the MEMS resonator and the oscillator circuitry in a single numerical modeling environment. This is achieved by transforming the conventional nonlinear mechanical model into the electrical domain while simultaneously considering the prominent nonlinearities of the resonator. The proposed nonlinear electrical model is validated by comparing the simulated amplitude¿frequency response with measurements on an open-loop electrically addressed flexural silicon MEMS resonator driven to large motional amplitudes. Next, the essential nonlinearities in the oscillator circuit are investigated and a mathematical model of a MEMS oscillator is proposed that integrates the nonlinearities of the resonator. The concept is illustrated for MEMS transimpedance-amplifier-based square-wave and sine-wave oscillators. Closed-form expressions of steady-state output power and output frequency are derived for both oscillator models and compared with experimental and simulation results, with a good match in the predicted trends in all three cases. © 1986-2012 IEEE.

Modeling the carbon nanofiber addressed liquid crystal microlens array from experimentally observed optical phenomena

This paper presents a novel method of using experimentally observed optical phenomena to reverse-engineer a model of the carbon nanofiber-addressed liquid crystal microlens array (C-MLA) using Zemax. It presents the first images of the optical profile for the C-MLA along the optic axis. The first working optical models of the C-MLA have been developed by matching the simulation results to the experimental results. This approach bypasses the need to know the exact carbon nanofiber-liquid crystal interaction and can be easily adapted to other systems where the nature of an optical device is unknown. Results show that the C-MLA behaves like a simple lensing system at 0.060-0.276 V/μm. In this lensing mode the C-MLA is successfully modeled as a reflective convex lens array intersecting with a flat reflective plane. The C-MLA at these field strengths exhibits characteristics of mostly spherical or low order aspheric arrays, with some aspects of high power aspherics. It also exhibits properties associated with varying lens apertures and strengths, which concur with previously theorized models based on E-field patterns. This work uniquely provides evidence demonstrating an apparent "rippling" of the liquid crystal texture at low field strengths, which were successfully reproduced using rippled Gaussian-like lens profiles. © 2014 Published by Elsevier B.V.

Design and modeling of an integrated device for acoustic resonance spectroscopy

This paper investigates the design and modelling of an integrated device for acoustic resonance spectroscopy (ARS). Miniaturisation of such platforms can be achieved using MEMS technology thereby enabling scaling of device dimensions to investigate smaller specimens while simultaneously operating at higher frequencies. We propose an integrated device where the transducers are mounted in close proximity with the specimen to be analysed (e.g. by integrating ultrasound transducers within a microfluidic channel). A finite element (FE) model and a simplified analytical model have been constructed to predict the acoustic response of a sample embedded in such a device configuration. A FE simulation is performed in COMSOL by embedding the piezoelectric transducers in representative fluid media. Resonant frequencies associated with the measurement can be extracted from this data. The response of various media modelled through FEA matches with analytical predictions for a range of biological media. A variety of biological media may be identified by using the measured resonant frequencies as a signature of relevant physical characteristics. The paper establishes the modelling basis of an integrated acoustic resonant spectrometer that is then applied to examine the impact of geometrical scaling on system resolution. © 2013 IEEE.

Kinetic modeling of growth and biodegradation of phenol and m-cresol using Alcaligenes faecalis

A phenol-degrading. microorganism, Alcaligenes faecalis, was used to study the substrate interactions during cell growth on phenol and m-cresol dual substrates. Both phenol and m-cresol could be utilized by the bacteria as,the sole carbon and energy sources. When cells grew on the mixture of phenol and m-cresol, strong substrate interactions were observed. m-Cresol inhibited the degradation of phenol, on the other hand, phenol also inhibited the utilization of m-cresol, the overall cell growth rate was the co-action of phenol and m-cresol. In addition, the cell growth and substrate degradation kinetics of phenol, m-cresol as single and mixed substrates for A. faecalis in batch cultures were also investigated over a wide range of initial phenol concentrations (10-1400 mg L-1) and initial m-cresol concentrations (5-200 mg L-1). The single-substrate kinetics was described well using the Haldane-type kinetic models, with model constants of it mu(m1) = 0.15 h(-1), K-S1 = 2.22 mg L-1 and K-i1 = 245.37 mg L-1 for cell growth on phenol and mu(m2) = 0.0782 h(-1), K-S2 = 1.30 mg L-1 and K-i2 = 71.77 mgL(-1), K-i2' = 5480 (mg L-1)(2) for cell growth on m-cresol. Proposed cell growth kinetic model was used to characterize the substrates interactions in the dual substrates system, the obtained parameters representing interactions between phenol and m-cresol were, K = 1.8 x 10(-6), M = 5.5 x 10(-5), Q = 6.7 x 10(-4). The results received in the experiments demonstrated that these models adequately described the dynamic behaviors of phenol and m-cresol as single and mixed substrates by the strain of A. faecalis.

Diffusion-limited kinetics modeling of one-step polyimide formation

A diffusion-limited kinetic model was developed to describe the imidization of one-step polythioetherimide formation based on an endgroup diffusion model. The changes of conversion and viscosity during the imidization were monitored with thermogravimetric analysis and dynamic stress rheometry, respectively. It was observed that the imidization rate began to decelerate after a fast early stage, whereas the viscosity in the system increased dramatically after a period of low value. Amic acid and imide formations concurrently take place in the one-step polyimide formation, but the formation of amic acid is much slower than that of imide and is the rate-limiting step of imidization. When a second-order kinetic model was used to describe the imidization, the effect of viscosity on the diffusion resistance of reactive groups needed to be included. In order to predict the change of viscosity during the imidization, the Lipshitz-Macosko model was modified and introduced into the diffusion-limited kinetic model by the Stokes-Einstein equation. The comparison of the modeled results with experimental data indicated that the diffusion-limited kinetic model and the modified Lipshitz-Macosko model were able to efficiently predict the changes of conversion and viscosity with temperature and time during the one-step polythioetherimide formation. (C) 2001 John Wiley & Sons, Inc.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone) .3. Determination of the time of pressure application

The curing temperature, pressure, and curing time have significant influence on finished thermosetting composite products. The time of pressure application is one of the most important processing parameters in the manufacture of a thermosetting composite. The determination of the time of pressure application relies on analysis of the viscosity variation of the polymer, associated with curing temperature and curing time. To determine it, the influence of the time of pressure application on the physical properties of epoxy-terminated poly(phenylene ether ketone) (E-PEK)-based continuous carbon fiber composite was studied. It was found that a stepwise temperature cure cycle is more suitable for manufacture of this composite. There are two viscosity valleys, in the case of the E-PEK system, associated with temperature during a stepwise cure cycle. The analysis on the effects of reinforcement fraction and defect content on the composite sheet quality indicates that the width-adjustable second viscosity valley provides a suitable pressing window. The viscosity, ranging from 400 to 1200 Pa . s at the second viscosity valley, is the optimal viscosity range for applying pressure to ensure appropriate resin flow during curing process, which enables one to get a finished composite with optimal fiber volume fraction and low void content. (C) 1997 John Wiley & Sons, Inc.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone) .4. Cure cycle simulation

Epoxy-terminated poly(phenylene ether ketone) (E-PEK) developed in this Institute is a candidate matrix resin for polymer composites as structural materials. Cure cycles for this reaction system were simulated according to the previously established processing model. It is found that for the E-PEK system, the curing process is best completed by a stepwise cure cycle comprising two isothermal processes at different temperatures, T-1 and T-2. The cure cycles over a wide range of processing parameters simulated, based on the established processing model, indicate that the processing window is width-adjustable. Analysis of the mechanical properties of the composite sheets showed that the simulated cure cycles are acceptable and reliable. (C) 1997 John Wiley & Sons, Inc.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone) .1. DSC characterization of curing reaction

The curing reaction process of epoxy-terminated poly(phenylene ether ketone) (E-PEK) with 4,4'-diaminodiphenyl sulfone (DDS) and hexahydrophthalic acid anhydride (Nadic) as curing agents was investigated using isothermal differential scanning calorimetry (IDSC) and nonisothermal differential scanning calorimetry (DDSC) techniques. It was found that the curing reactions of E-PEK/DDS and E-PEK/Nadic are nth-order reactions but not autoaccelerating. The experimental results revealed that the curing reaction kinetics parameters measured from IDSC and DDSC are not equivalent. This means that, in the curing reaction kinetics model for our E-PEK system, both isothermal and nonisothermal reaction kinetics parameters are needed to describe isothermal and nonisothermal curing processes, The isothermal and nonisothermal curing processes were successfully simulated using this model. A new extrapolation method was suggested. On the basis of this method the maximum extent of the curing reaction (A(ult)) that is able to reach a certain temperature can be predicted. The A(ult) for the E-PEK system estimated by the new method agrees well with the results obtained from another procedure reported in the literature. (C) 1997 John Wiley & Sons, Inc.