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 mode i fracture of bitumen films

The fracture behavior of thin films of bitumen in double cantilever beam (DCB) specimens was investigated over a wide range of temperature and loading rate conditions using finite-element analysis. The model includes a phenomenological model for the mechanical behavior of bitumen, implemented into a special-purpose finite-element user material subroutine, combined with a cohesive zone model (CZM) for simulating the fracture process. The finite-element model is validated against experimental results from laboratory tests of DCB specimens by comparing measured and predicted load-line deflection histories and fracture energy release rates. Computer simulation results agreed well with experimental data of DCB joints containing bitumen films in terms of peak stress, fracture toughness, and stress-strain history response. The predicted "normalized toughness," G=2h, was found to increase in a power-law manner with effective temperaturecompensated strain rate in the ductile region as previously observed experimentally. In the brittle regime, G=2h is virtually constant. The model successfully captured the ductile and brittle failure behavior of bitumen films in opening mode (tension) for stable crack growth conditions. © 2013 American Society of Civil Engineers.

Aerodynamic design of high end wall angle turbine stages- part II: Experimental verification

An experimental investigation of a turbine stage featuring very high end wall angles is presented. The initial turbine design did not achieve a satisfactory performance and the difference between the design predictions and the test results was traced to a large separated region on the rear suction-surface. To improve the agreement between computational fluid dynamics (CFD) and experiment, it was found necessary to modify the turbulence modeling employed. The modified CFD code was then used to redesign the vane, and the changes made are described. When tested, the performance of the redesigned vane was found to have much closer agreement with the predictions than the initial vane. Finally, the flowfield and performance of the redesigned stage are compared to a similar turbine, designed to perform the same duty, which lies in an annulus of moderate end wall angles. A reduction in stage efficiency of at least 2.4% was estimated for the very high end wall angle design. © 2014 by ASME.

Corrugated all-composite sandwich structures. Part 1: Modeling

An analytical model for the compressive and shear response of monolithic and hierarchical corrugated composite cores has been developed. The stiffness model considers the contribution in stiffness from the bending- and the shear deformations of the core members in addition to the stretching deformation. The strength model is based on the normal stress and shear stress distribution over each core member when subjected to a shear or compressive load condition. The strength model also accounts for initial imperfections. In part 1 of this series, the analytical model is described and the results are compared to finite element predictions. In part 2, the analytical model is compared to experimental results and the behaviour of the corrugated structures is investigated more thoroughly using failure mechanism maps. © 2008 Elsevier Ltd. All rights reserved.

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.

Optical constants of cubic GaN/GaAs(001): Experiment and modeling

The optical constants epsilon(E)=epsilon(1)(E)+iepsilon(2)(E) of unintentionally doped cubic GaN grown on GaAs(001) have been measured at 300 K using spectral ellipsometry in the range of 1.5-5.0 eV. The epsilon(E) spectra display a structure, associated with the critical point at E-0 (direct gap) and some contribution mainly coming from the E-1 critical point. The experimental data over the entire measured spectral range (after oxide removal) has been fit using the Holden-Munoz model dielectric function [M. Munoz et al., J. Appl. Phys. 92, 5878 (2002)]. This model is based on the electronic energy-band structure near critical points plus excitonic and band-to-band Coulomb-enhancement effects at E-0, E-0 + Delta(0) and the E-1, E-1 + Delta(1), doublet. In addition to evaluating the energy of the E-0 critical point, the binding energy (R-1) of the two-dimensional exciton related to the E-1 critical point was estimated using the effective mass/k.p theory. The line, shape of the imaginary part of the cubic-GaN dielectric function shows excitonic effects at room temperature not withstanding that the exciton was not resolved. (C) 2003 American Institute of Physics.

Kinetic modeling of N incorporation in GaInNAs growth by plasma-assisted molecular-beam epitaxy

We have studied the growth of GaInNAs by a plasma-assisted molecular-beam epitaxy (MBE). It was found that the N-radicals were incorporated into the epitaxial layer like dopant atoms. In the range of 400-500 degrees C, the growth temperature (T-g) mainly affected the crystal quality of GaInNAs rather than the N concentration. The N concentration dropped rapidly when T-g exceeded 500 degrees C. Considering N desorption alone is insufficient to account for the strong falloff of the N concentration with T-g over 500 degrees C, the effect of thermally-activated N surface segregation must be taken into account. The N concentration was independent of the arsenic pressure and the In concentration in GaInNAs layers, but inversely proportional to the growth rate. Based on the experimental results, a kinetic model including N desorption and surface segregation was developed to analyze quantitatively the N incorporation in MBE growth. (C) 2000 American Institute of Physics. [S0003-6951(00)00928-1].

Modeling and optimization of semiconductor manufacturing process with neural networks

A neural network-based process model is proposed to optimize the semiconductor manufacturing process. Being different from some works in several research groups which developed neural network-based models to predict process quality with a set of process variables of only single manufacturing step, we applied this model to wafer fabrication parameters control and wafer lot yield optimization. The original data are collected from a wafer fabrication line, including technological parameters and wafer test results. The wafer lot yield is taken as the optimization target. Learning from historical technological records and wafer test results, the model can predict the wafer yield. To eliminate the "bad" or noisy samples from the sample set, an experimental method was used to determine the number of hidden units so that both good learning ability and prediction capability can be obtained.

Modeling of a coal-fired slagging combustor: Development of a slag submodel

In a slagging combustor or furnace, the high combustion temperature makes the molten slag layer cover the wall and capture the particles. If these particles contain combustible matter, they will continue to burn on the running slag. As a result, the total amount of ash deposition will be much greater than that in dry-wall combustors and the total heat flux through the deposition surface will change greatly. Considering the limitations of existing simulation methods for slagging combustion, this paper introduces a new wall burning model and slag flow model from the analysis; of particle deposition phenomena. Combined with a conventional combustion simulation program, the total computational frame is introduced. From comparisons of simulation results from several kinds of methods with experimental data, the conclusion is drawn that the conventional simulation methods are not very suitable for slagging combustion and the wall burning mechanism should be considered more thoroughly.

MODELING OF ELECTRON-TUNNELING TIMES IN ASYMMETRIC DOUBLE QUANTUM-WELLS

A scattering process modeled by an imaginary potential V(I) in the wide well of an asymmetric double quantum well structure (DQWS) is used to model the electron tunneling from the narrow well. Taking V(I) approximately -5 meV, the ground resonant level lifetimes of the narrow well in the DQWS are in quantitative agreement with the experimental resonance and non-resonance tunneling times. The corresponding scattering time 66 fs is much faster than the intersubband scattering time of LO-photon emission.

A finite-element algorithm for modeling variably saturated flows

A general numerical algorithm in the context of finite element scheme is developed to solve Richards’ equation, in which a mass-conservative, modified head based scheme (MHB) is proposed to approximate the governing equation, and mass-lumping techniques are used to keep the numerical simulation stable. The MHB scheme is compared with the modified Picard iteration scheme (MPI) in a ponding infiltration example. Although the MHB scheme is a little inferior to the MPI scheme in respect of mass balance, it is superior in convergence character and simplicity. Fully implicit, explicit and geometric average conductivity methods are performed and compared, the first one is superior in simulation accuracy and can use large time-step size, but the others are superior in iteration efficiency. The algorithm works well over a wide variety of problems, such as infiltration fronts, steady-state and transient water tables, and transient seepage faces, as demonstrated by its performance against published experimental data. The algorithm is presented in sufficient detail to facilitate its implementation.

Modeling study on the flow, heat transfer and energy conversion characteristics of low-power arc-heated hydrogen/nitrogen thrusters

A modeling study is conducted to investigate the effect of hydrogen content in propellants on the plasma flow, heat transfer and energy conversion characteristics of low-power (kW class) arc-heated hydrogen/nitrogen thrusters (arcjets). 1:0 (pure hydrogen), 3:1 (to simulate decomposed ammonia), 2:1 (to simulate decomposed hydrazine) and 0:1 (pure nitrogen) hydrogen/nitrogen mixtures are chosen as the propellants. Both the gas flow region inside the thruster nozzle and the anode-nozzle wall are included in the computational domain in order to better treat the conjugate heat transfer between the gas flow region and the solid wall region. The axial variations of the enthalpy flux, kinetic energy flux, directed kinetic-energy flux, and momentum flux, all normalized to the mass flow rate of the propellant, are used to investigate the energy conversion process inside the thruster nozzle. The modeling results show that the values of the arc voltage, the gas axial-velocity at the thruster exit, and the specific impulse of the arcjet thruster all increase with increasing hydrogen content in the propellant, but the gas temperature at the nitrogen thruster exit is significantly higher than that for other three propellants. The flow, heat transfer, and energy conversion processes taking place in the thruster nozzle have some common features for all the four propellants. The propellant is heated mainly in the near-cathode and constrictor region, accompanied with a rapid increase of the enthalpy flux, and after achieving its maximum value, the enthalpy flux decreases appreciably due to the conversion of gas internal energy into its kinetic energy in the divergent segment of the thruster nozzle. The kinetic energy flux, directed kinetic energy flux and momentum flux also increase at first due to the arc heating and the thermodynamic expansion, assume their maximum inside the nozzle and then decrease gradually as the propellant flows toward the thruster exit. It is found that a large energy loss (31-52%) occurs in the thruster nozzle due to the heat transfer to the nozzle wall and too long nozzle is not necessary. Modeling results for the NASA 1-kW class arcjet thruster with hydrogen or decomposed hydrazine as the propellant are found to compare favorably with available experimental data.

Modeling study to compare the flow and heat transfer characteristics of low-power hydrogen, nitrogen and argon arc-heated thrusters

A modelling study is performed to compare the plasma °ow and heat transfer char- acteristics of low-power arc-heated thrusters (arcjets) for three di®erent propellants: hydrogen, nitrogen and argon. The all-speed SIMPLE algorithm is employed to solve the governing equa- tions, which take into account the e®ects of compressibility, Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. The temperature, veloc- ity and Mach number distributions calculated within the thruster nozzle obtained with di®erent propellant gases are compared for the same thruster structure, dimensions, inlet-gas stagnant pressure and arc currents. The temperature distributions in the solid region of the anode-nozzle wall are also given. It is found that the °ow and energy conversion processes in the thruster nozzle show many similar features for all three propellants. For example, the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appear- ing near the cathode tip; the °ow transition from the subsonic to supersonic regime occurs within the constrictor region; the highest axial velocity appears inside the nozzle; and most of the input propellant °ows towards the thruster exit through the cooler gas region near the anode-nozzle wall. However, since the properties of hydrogen, nitrogen and argon, especially their molecular weights, speci¯c enthalpies and thermal conductivities, are di®erent, there are appreciable di®er- ences in arcjet performance. For example, compared to the other two propellants, the hydrogen arcjet thruster shows a higher plasma temperature in the arc region, and higher axial velocity but lower temperature at the thruster exit. Correspondingly, the hydrogen arcjet thruster has the highest speci¯c impulse and arc voltage for the same inlet stagnant pressure and arc current. The predictions of the modelling are compared favourably with available experimental results.

Numerical modeling of heat transfer and plasma flow in a low power arcjet thruster

Modeling studies are preformed to investigate the plasma and heat transfer characteristics of a low power argon arcjet thruster. Computed temperature, velocity, static pressure, and Mach number distribution in arcjet thruster under typical operating condition are presented in this paper. It shows that the performance data from numerical modeling results are basically consistent with the experimental measured values.

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.