3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

A Synoptic View of a Coastal Plain Estuary.

During October, 1972 the Patuxent River Estuary was monitored intensively and synoptically over two tidal cycles to determine the spatial and temporal patterns of various hydrodynamic, chemical and biological features. Forty-one depths at eleven stations along nine transects were sampled simultaneously at hourly intervals for salinity, temperature, dissolved oxygen, chlorohyll a, particulate nitrogen, nitrate, nitrite, total kjeldahl nitrogen, ammonia, particulate carbohydrate, dissolved organic carbon, total hydrolizable phosphorous, dissolved inorganic phosphorous, suspended sediment, particle size distribution, and zooplankton. Tidal velocity was continuously monitored at each depth by recording current meters. Riverine input and meteorological conditions were relatively stable for two weeks preceeding the deployment. This communication describes the calculation of the intrinsic rates of change of the observed variables from their measured distributions in the Estuary. The steady-state, one-dimensional equation of species continuity is employed to separate the advection and tidal dispersion of a hydrodynamically passive substance frbm its intrinsic rate of change at point. A new spatial transform is introduced for the purpose of interpolation and extrapolation of data.The intrinsic rate of change profiles reveal a region of heavy bloom activity in the upper estuary and a secondary bloom near the point in the River that most of the suspended material settles out. The changes in ammonia and nitrates are highly correlated to the productivity patterns. Phosphorous rates are less closely correlated to productivity. The perturbations that the Chalk Point steam electric power plant have on the heat and oxygen balances are easily discernible.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Investigations of cohesive zone properties of the interface between metal film and ceramic substrate at nano- and micro-scales

Cohesive zone characterizations of the interface between metal film and ceramic substrate at micro- and nano-scales are performed in the present research. At the nano-scale, a special potential for special material interface (Ag/MgO) is adopted to investigate the interface separation mechanism by using MD simulation, and stress-separation relationship will be obtained. At the micro-scale, peeling experiment is performed for the Al film/Al2O3 substrate system with an adhesive layer at the interface. Adhesive is a mixture of epoxy and polyimide with mass ratio 1:1, by which a brittle cohesive property is obtained. The relationships between energy release rate, the film thickness and the adhesive layer thickness are measured during the steady-state peeling. The experimental result has a similar trend as modeling result for a weak adhesion interface case.

Subsonic interface crack with crack face contact

A steady-state subsonic interface crack propagating between an elastic solid and a rigid substrate with crack face contact is studied. Two cases with respective to the contact length are considered, i.e., semi-infinite and finite crack face contact. Different from a stationary or an open subsonic interface crack, stress singularity at the crack tip in the present paper is found to be non-oscillatory. Furthermore, in the semi-infinite contact case, the singularity of the stress field near the crack tip is less than 1/2. In the finite contact case, no singularity exists near the crack tip, but less than 1/2 singularity does at the end of the contact zone. In both cases, the singularity depends on the linear contact coefficient and the crack speed. Asymptotic solutions near the crack tip are given and analyzed. In order to satisfy the contact conditions, reasonable region of the linear contact coefficient is found. In addition, the solution predicts a non-zero-energy dissipation rate due to crack face contact.

First Biennial Ocean Climate Summit: Finding Solutions for San Francisco Bay Area’s Coast and Ocean. Proceedings of the Summit: April 29, 2008, Golden Gate Club, The Presidio San Francisco, California

Executive Summary: The marine environment plays a critical role in the amount of carbon dioxide (CO2) that remains within Earth’s atmosphere, but has not received as much attention as the terrestrial environment when it comes to climate change discussions, programs, and plans for action. It is now apparent that the oceans have begun to reach a state of CO2 saturation, no longer maintaining the “steady-state” carbon cycle that existed prior to the Industrial Revolution. The increasing amount of CO2 present within the oceans and the atmosphere has an effect on climate and a cascading effect on the marine environment. Potential physical effects of climate change within the marine environment, including ocean acidification, changes in wind and upwelling regimes, increasing global sea surface temperatures, and sea level rise, can lead to dramatic, fundamental changes within marine and coastal ecosystems. Altered ecosystems can result in changing coastal economies through a reduction in marine ecosystem services such as commercial fish stocks and coastal tourism. Local impacts from climate change should be a front line issue for natural resource managers, but they often feel too overwhelmed by the magnitude of this issue to begin to take action. They may not feel they have the time, funding, or staff to take on a challenge as large as climate change and continue to not act as a result. Already, natural resource managers work to balance the needs of humans and the economy with ecosystem biodiversity and resilience. Responsible decisions are made each day that consider a wide variety of stakeholders, including community members, agencies, non-profit organizations, and business/industry. The issue of climate change must be approached as a collaborative effort, one that natural resource managers can facilitate by balancing human demands with healthy ecosystem function through research and monitoring, education and outreach, and policy reform. The Scientific Expert Group on Climate Change in their 2007 report titled, “Confronting Climate Change: Avoiding the Unmanageable and Managing the Unavoidable” charged governments around the world with developing strategies to “adapt to ongoing and future changes in climate change by integrating the implications of climate change into resource management and infrastructure development”. Resource managers must make future management decisions within an uncertain and changing climate based on both physical and biological ecosystem response to climate change and human perception of and response to the issue. Climate change is the biggest threat facing any protected area today and resource managers must lead the charge in addressing this threat. (PDF has 59 pages.)

Nutrient distributions for Georges Bank and adjacent waters in 1979

In this report we describe the temporal and spatial distributions of inorganic nutrients over Georges Bank and in adjacent waters and discuss major features with respect to tbe nutrient environments of pbytoplankton. Nitrate and orthophosphorus were rapidly depleted from the surface layer of much of the study area in spring, but major differences were found between the shallow areas on Georges Bank and the surrounding stratified waters. In the "well-mixed" area of Georges Bank, the depletion encompassed the entire water column and ammonium became the dominant form of inorganic nitrogen throughout. Dissolved silicon was depleted slowly over central Georges Bank, reaching a minimum concentration in September while orthophosphorus gradually increased during the summer. The nutrient environment of phytoplankton over central Georges Bank may be described as vertically uniform but temporally changing in the relative availability of the various nutrients. In areas that undergo stratification (e.g., the central Gulf of Maine), a quasi-steady state was established as the surface water layer formed, consisting of declining nutrient gradients from below the euphotic layer to the top of the water column. These intergrading nutrient environments are relatively stable through time. Destratification reintroduced nutrients to depleted areas beginning in October; however, dissolved silicon was again depleted over shallow Georges Bank in late autumn though nitrate remained abundant. Slope water has been found to enter the bottom layer of the Gulf of Maine via the Northeast Channel. High nutrient concentrations observed in the bottom water of the Northeast Channel are consistent with this mechanism being the nutrient source for the Gulf of Maine. (PDF file contains 40 pages.)

Wage stickiness and unemployment fluctuations: an alternative approach

28 p.

Periodic solutions of integro-differential equations which arise in population dynamics

The problem of the existence and stability of periodic solutions of infinite-lag integra-differential equations is considered. Specifically, the integrals involved are of the convolution type with the dependent variable being integrated over the range (- ∞,t), as occur in models of population growth. It is shown that Hopf bifurcation of periodic solutions from a steady state can occur, when a pair of eigenvalues crosses the imaginary axis. Also considered is the existence of traveling wave solutions of a model population equation allowing spatial diffusion in addition to the usual temporal variation. Lastly, the stability of the periodic solutions resulting from Hopf bifurcation is determined with aid of a Floquet theory.

The first chapter is devoted to linear integro-differential equations with constant coefficients utilizing the method of semi-groups of operators. The second chapter analyzes the Hopf bifurcation providing an existence theorem. Also, the two-timing perturbation procedure is applied to construct the periodic solutions. The third chapter uses two-timing to obtain traveling wave solutions of the diffusive model, as well as providing an existence theorem. The fourth chapter develops a Floquet theory for linear integro-differential equations with periodic coefficients again using the semi-group approach. The fifth chapter gives sufficient conditions for the stability or instability of a periodic solution in terms of the linearization of the equations. These results are then applied to the Hopf bifurcation problem and to a certain population equation modeling periodically fluctuating environments to deduce the stability of the corresponding periodic solutions.

Some problems in nonlinear elasticity

Two separate problems are discussed: axisymmetric equilibrium configurations of a circular membrane under pressure and subject to thrust along its edge, and the buckling of a circular cylindrical shell.

An ordinary differential equation governing the circular membrane is imbedded in a family of n-dimensional nonlinear equations. Phase plane methods are used to examine the number of solutions corresponding to a parameter which generalizes the thrust, as well as other parameters determining the shape of the nonlinearity and the undeformed shape of the membrane. It is found that in any number of dimensions there exists a value of the generalized thrust for which a countable infinity of solutions exist if some of the remaining parameters are made sufficiently large. Criteria describing the number of solutions in other cases are also given.

Donnell-type equations are used to model a circular cylindrical shell. The static problem of bifurcation of buckled modes from Poisson expansion is analyzed using an iteration scheme and pertubation methods. Analysis shows that although buckling loads are usually simple eigenvalues, they may have arbitrarily large but finite multiplicity when the ratio of the shell's length and circumference is rational. A numerical study of the critical buckling load for simple eigenvalues indicates that the number of waves along the axis of the deformed shell is roughly proportional to the length of the shell, suggesting the possibility of a "characteristic length." Further numerical work indicates that initial post-buckling curves are typically steep, although the load may increase or decrease. It is shown that either a sheet of solutions or two distinct branches bifurcate from a double eigenvalue. Furthermore, a shell may be subject to a uniform torque, even though one is not prescribed at the ends of the shell, through the interaction of two modes with the same number of circumferential waves. Finally, multiple time scale techniques are used to study the dynamic buckling of a rectangular plate as well as a circular cylindrical shell; transition to a new steady state amplitude determined by the nonlinearity is shown. The importance of damping in determining equilibrium configurations independent of initial conditions is illustrated.

Anatomical and developmental patterns of interleukin-2 gene expression in the mouse: analysis of IL-2 expressing cells and partial characterization of interactions involved in mediating IL-2 gene expression in vivo

Interleukin-2 (IL-2) is an important mediator in the vertebrate immune system. IL-2 is a potent growth factor that mature T lymphocytes use as a proliferation signal and the production of IL-2 is crucial for the clonal expansion of antigen-specific T cells in the primary immune response. IL-2 driven proliferation is dependent on the interaction of the lymphokine with its cognate multichain receptor. IL-2 expression is induced only upon stimulation and transcriptional activation of the IL-2 gene relies extensively on the coordinate interaction of numerous inducible and constitutive trans-acting factors. Over the past several years, thousands of papers have been published regarding molecular and cellular aspects of IL-2 gene expression and IL-2 function. The vast majority of these reports describe work that has been carried out in vitro. However, considerably less is known about control of IL-2 gene expression and IL-2 function in vivo.

To gain new insight into the regulation of IL-2 gene expression in vivo, anatomical and developmental patterns of IL-2 gene expression in the mouse were established by employing in situ hybridization and immunohistochemical staining methodologies to tissue sections generated from normal mice and mutant animals in which T -cell development was perturbed. Results from these studies revealed several interesting aspects of IL-2 gene expression, such as (1) induction of IL-2 gene expression and protein synthesis in the thymus, the primary site of T-cell development in the body, (2) cell-type specificity of IL-2 gene expression in vivo, (3) participation of IL-2 in the extrathymic expansion of mature T cells in particular tissues, independent of an acute immune response to foreign antigen, (4) involvement of IL-2 in maintaining immunologic balance in the mucosal immune system, and (5) potential function of IL-2 in early events associated with hematopoiesis.

Extensive analysis of IL-2 mRNA accumulation and protein production in the murine thymus at various stages of development established the existence of two classes of intrathymic IL-2 producing cells. One class of intrathymic IL-2 producers was found exclusively in the fetal thymus. Cells belonging to this subset were restricted to the outermost region of the thymus. IL-2 expression in the fetal thymus was highly transient; a dramatic peak ofiL-2 mRNA accumulation was identified at day 14.5 of gestation and maximal IL-2 protein production was observed 12 hours later, after which both IL-2 mRNA and protein levels rapidly decreased. Significantly, the presence of IL-2 expressing cells in the day 14-15 fetal thymus was not contingent on the generation of T-cell receptor (TcR) positive cells. The second class of IL-2 producing cells was also detectable in the fetal thymus (cells found in this class represented a minority subset of IL-2 producers in the fetal thymus) but persist in the thymus during later stages of development and after birth. Intrathymic IL-2 producers in postnatal animals were located in the subcapsular region and cortex, indicating that these cells reside in the same areas where immature T cells are consigned. The frequency of IL-2 expressing cells in the postnatal thymus was extremely low, indicating that induction of IL-2 expression and protein synthesis are indicative of a rare activation event. Unlike the fetal class of intrathymic IL-2 producers, the presence of IL-2 producing cells in the postnatal thymus was dependent on to the generation of TcR+ cells. Subsequent examination of intrathymic IL-2 production in mutant postnatal mice unable to produce either αβ or γδ T cells showed that postnatal IL-2 producers in the thymus belong to both αβ and γδ lineages. Additionally, further studies indicated that IL-2 synthesis by immature αβ -T cells depends on the expression of bonafide TcR αβ-heterodimers. Taken altogether, IL-2 production in the postnatal thymus relies on the generation of αβ or γδ-TcR^+ cells and induction of IL-2 protein synthesis can be linked to an activation event mediated via the TcR.

With regard to tissue specificity of IL-2 gene expression in vivo, analysis of whole body sections obtained from normal neonatal mouse pups by in situ hybridization demonstrated that IL-2 mRNA^+ cells were found in both lymphoid and nonlymphoid tissues with which T cells are associated, such as the thymus (as described above), dermis and gut. Tissues devoid of IL-2 mRNA^+ cells included brain, heart, lung, liver, stomach, spine, spinal cord, kidney, and bladder. Additional analysis of isolated tissues taken from older animals revealed that IL-2 expression was undetectable in bone marrow and in nonactivated spleen and lymph nodes. Thus, it appears that extrathymic IL-2 expressing cells in nonimmunologically challenged animals are relegated to particular epidermal and epithelial tissues in which characterized subsets of T cells reside and thatinduction of IL-2 gene expression associated with these tissues may be a result of T-cell activation therein.

Based on the neonatal in situ hybridization results, a detailed investigation into possible induction of IL-2 expression resulting in IL-2 protein synthesis in the skin and gut revealed that IL-2 expression is induced in the epidermis and intestine and IL-2 protein is available to drive cell proliferation of resident cells and/or participate in immune function in these tissues. Pertaining to IL-2 expression in the skin, maximal IL-2 mRNA accumulation and protein production were observed when resident Vγ_3^+ T-cell populations were expanding. At this age, both IL-2 mRNA^+ cells and IL-2 protein production were intimately associated with hair follicles. Likewise, at this age a significant number of CD3ε^+ cells were also found in association with follicles. The colocalization of IL-2 expression and CD3ε^+ cells suggests that IL-2 expression is induced when T cells are in contact with hair follicles. In contrast, neither IL-2 mRNA nor IL-2 protein were readily detected once T-cell density in the skin reached steady-state proportions. At this point, T cells were no longer found associated with hair follicles but were evenly distributed throughout the epidermis. In addition, IL-2 expression in the skin was contingent upon the presence of mature T cells therein and induction of IL-2 protein synthesis in the skin did not depend on the expression of a specific TcR on resident T cells. These newly disclosed properties of IL-2 expression in the skin indicate that IL-2 may play an additional role in controlling mature T-cell proliferation by participating in the extrathymic expansion of T cells, particularly those associated with the epidermis.

Finally, regarding IL-2 expression and protein synthesis in the gut, IL-2 producing cells were found associated with the lamina propria of neonatal animals and gut-associated IL-2 production persisted throughout life. In older animals, the frequency of IL-2 producing cells in the small intestine was not identical to that in the large intestine and this difference may reflect regional specialization of the mucosal immune system in response to enteric antigen. Similar to other instances of IL-2 gene expression in vivo, a failure to generate mature T cells also led to an abrogation of IL-2 protein production in the gut. The presence of IL-2 producing cells in the neonatal gut suggested that these cells may be generated during fetal development. Examination of the fetal gut to determine the distribution of IL-2 producing cells therein indicated that there was a tenfold increase in the number of gut-associated IL-2 producers at day 20 of gestation compared to that observed four days earlier and there was little difference between the frequency of IL-2 producing cells in prenatal versus neonatal gut. The origin of these fetally-derived IL-2 producing cells is unclear. Prior to the immigration of IL-2 inducible cells to the fetal gut and/or induction of IL-2 expression therein, IL-2 protein was observed in the fetal liver and fetal omentum, as well as the fetal thymus. Considering that induction of IL-2 protein synthesis may be an indication of future functional capability, detection of IL-2 producing cells in the fetal liver and fetal omentum raises the possibility that IL-2 producing cells in the fetal gut may be extrathymic in origin and IL-2 producing cells in these fetal tissues may not belong solely to the T lineage. Overall, these results provide increased understanding of the nature of IL-2 producing cells in the gut and how the absence of IL-2 production therein and in fetal hematopoietic tissues can result in the acute pathology observed in IL-2 deficient animals.

Energy transfer between silicon-oxygen-related defects and Yb3+ in transparent glass ceramics containing Ba2TiSi2O 8 nanocrystals

We report on the conversion of near-ultraviolet radiation of 250-350 nm into near-infrared emission of 970-1100 nm in Yb3+-doped transparent glass ceramics containing Ba2TiSi2O8 nanocrystals due to the energy transfer from the silicon-oxygen-related defects to Yb3+ ions. Efficient Yb3+ emission (F-2(5/2)-> F-2(7/2)) was detected under the excitation of defects absorption at 314 nm. The occurrence of energy transfer is proven by both steady state and time-resolved emission spectra, respectively, at 15 K. The Yb2O3 concentration dependent energy transfer efficiency has also been evaluated, and the maximum value is 65% for 8 mol % Yb2O3 doped glass ceramic. These materials are promising for the enhancement of photovoltaic conversion efficiency of silicon solar cells via spectra modification.

Flight dynamics in Drosophila through a dynamically-scaled robotic approach

Flies are particularly adept at balancing the competing demands of delay tolerance, performance, and robustness during flight, which invites thoughtful examination of their multimodal feedback architecture. This dissertation examines stabilization requirements for inner-loop feedback strategies in the flapping flight of Drosophila, the fruit fly, against the backdrop of sensorimotor transformations present in the animal. Flies have evolved multiple specializations to reduce sensorimotor latency, but sensory delay during flight is still significant on the timescale of body dynamics. I explored the effect of sensor delay on flight stability and performance for yaw turns using a dynamically-scaled robot equipped with a real-time feedback system that performed active turns in response to measured yaw torque. The results show a fundamental tradeoff between sensor delay and permissible feedback gain, and suggest that fast mechanosensory feedback provides a source of active damping that compliments that contributed by passive effects. Presented in the context of these findings, a control architecture whereby a haltere-mediated inner-loop proportional controller provides damping for slower visually-mediated feedback is consistent with tethered-flight measurements, free-flight observations, and engineering design principles. Additionally, I investigated how flies adjust stroke features to regulate and stabilize level forward flight. The results suggest that few changes to hovering kinematics are actually required to meet steady-state lift and thrust requirements at different flight speeds, and the primary driver of equilibrium velocity is the aerodynamic pitch moment. This finding is consistent with prior hypotheses and observations regarding the relationship between body pitch and flight speed in fruit flies. The results also show that the dynamics may be stabilized with additional pitch damping, but the magnitude of required damping increases with flight speed. I posit that differences in stroke deviation between the upstroke and downstroke might play a critical role in this stabilization. Fast mechanosensory feedback of the pitch rate could enable active damping, which would inherently exhibit gain scheduling with flight speed if pitch torque is regulated by adjusting stroke deviation. Such a control scheme would provide an elegant solution for flight stabilization across a wide range of flight speeds.

Optical bistability via amplitude and phase control of a microwave field

We investigate the steady-state optical bistability behavior in a three-level A-type atomic system closed by a microwave field under the condition that the applied fields are in resonance with corresponding atomic transitions. It is shown that the bistable hysteresis cycles can be controlled by both the amplitude and the phase of the microwave field. (c) 2006 Elsevier B.V. All rights reserved.

Modeling and analysis of hysteretic structural behavior

For damaging response, the force-displacement relationship of a structure is highly nonlinear and history-dependent. For satisfactory analysis of such behavior, it is important to be able to characterize and to model the phenomenon of hysteresis accurately. A number of models have been proposed for response studies of hysteretic structures, some of which are examined in detail in this thesis. There are two popular classes of models used in the analysis of curvilinear hysteretic systems. The first is of the distributed element or assemblage type, which models the physical behavior of the system by using well-known building blocks. The second class of models is of the differential equation type, which is based on the introduction of an extra variable to describe the history dependence of the system.

Owing to their mathematical simplicity, the latter models have been used extensively for various applications in structural dynamics, most notably in the estimation of the response statistics of hysteretic systems subjected to stochastic excitation. But the fundamental characteristics of these models are still not clearly understood. A response analysis of systems using both the Distributed Element model and the differential equation model when subjected to a variety of quasi-static and dynamic loading conditions leads to the following conclusion: Caution must be exercised when employing the models belonging to the second class in structural response studies as they can produce misleading results.

The Massing's hypothesis, originally proposed for steady-state loading, can be extended to general transient loading as well, leading to considerable simplification in the analysis of the Distributed Element models. A simple, nonparametric identification technique is also outlined, by means of which an optimal model representation involving one additional state variable is determined for hysteretic systems.