The CTH Eulerian hydrocode, and the SPHINX smooth particle hydrodynamics (SPH) code were used to model a shock tube, two long rod penetrations into semi-infinite steel targets, and a long rod penetration into a spaced plate array. The results were then compared to experimental data. Both SPHINX and CTH modeled the one-dimensional shock tube problem well. Both codes did a reasonable job in modeling the outcome of the axisymmetric rod impact problem. Neither code correctly reproduced the depth of penetration in both experiments. In the 3-D problem, both codes reasonably replicated the penetration of the rod through the first plate. After this, however, the predictions of both codes began to diverge from the results seen in the experiment. In terms of computer resources, the run times are problem dependent, and are discussed in the text.
As part of the FAA`s National Aging Aircraft Research Program to foster new technologies for civil aircraft maintenance and repair, use of bonded composite doublers on metal aircraft structures has been advanced. Research and validation of such doubler applications on US certified commercial aircraft has begun. A specific composite application to assess the capabilities of composite doublers was chosen on a L-1011 aircraft for reinforcement of the comer of a cargo door frame where a boron-epoxy repair patch of up to 72 plies was installed. A primary inspection requirement for these doublers is the identification of disbonds between the composite laminate and the aluminum parent material. This paper describes the development of an ultrasonic pulse echo technique using a modified immersion focus transducer where a robust signal amplitude signature of the composite aluminum interface is obtained to characterize the condition of the bond. Example waveforms and C-scan images are shown to illustrate the ultrasonic response for various transducer configurations using a boron-epoxy aluminum skin calibration test sample where disbonds and delaminations were built-in. The modified focus transducer is compatible with portable ultrasonic scanning systems that utilize the weeper or dripless bubbler technologies when an ultrasonic inspection of the boron-epoxy composite doublers installed on aircraft is implemented.
An eddy current inspection method was developed at the Federal Aviation Administration`s Airworthiness Assurance NDI Validation Center (AANC) to easily and rapidly detect subsurface fatigue cracks in the wheel well fairing on the US Coast Guard (USCG) HC-130H aircraft caused by fatigue. The inspection procedure locates cracks as small as 10.2 millimeters in length at 2.54 mm below the skin surface at raised fastener sites. The test procedure developed baseline three USCG aircraft. Inspection results on the three aircraft reveals good correlation with results made during subsequent structural disassembly.
GaN MIS diodes were demonstrated utilizing AlN and Ga{sub 2}O{sub 3}(Gd{sub 2}O{sub 3}) as insulators. A 345 {angstrom} of AlN was grown on the MOCVD grown n-GaN in a MOMBE system using trimethylamine alane as Al precursor and nitrogen generated from a wavemat ECR N2 plasma. For the Ga{sub 2}O{sub 3}(Gd{sub 2}O{sub 3}) growth, a multi MBE chamber was used and a 195 {angstrom} oxide is E-beam evaporated from a single crystal source of Ga{sub 5}Gd{sub 3}O{sub 12}. The forward breakdown voltage of AlN and Ga{sub 2}O{sub 3}(Gd{sub 2}O{sub 3}) diodes are 5V and 6V, respectively, which are significantly improved from {approximately} 1.2 V of schottky contact. From the C-V measurements, both kinds of diodes showed good charge modulation from accumulation to depletion at different frequencies. The insulator GaN interface roughness and the thickness of the insulator were measured with x-ray reflectivity.
In the past thirty-six months, tremendous strides have been made in x-ray production using high-current z-pinches. Today, the x-ray energy (1.9 MJ) and power (200 TW) output of the Z accelerator (formerly PBFA-II) is the largest available in the laboratory. These z-pinch x-ray sources are being developed for research into the physics of high energy density plasmas of interest in weapon behavior and in inertial confinement fusion. Beyond the Z accelerator current of 20 MA, an extrapolation to the X-1 accelerator level of 60 MA may have the potential to drive high-yield ICF reactions at affordable cost if several challenging technical problems can be overcome. New developments have also taken place at Sandia in the area of high current, mm-diameter electron beams for advanced hydrodynamic radiography. On SABRE, x-ray spot diameters were less than 2 mm with a dose of 100 R at 1 meter in a 40 ns pulse.
Prosperity Games TM are an outgrowth and adaptation move/countermove and seminar War Games. Prosperity Games TM are simulations that explore complex issues in a variety of areas including economics, politics, sociology, environment, education, and research. These issues can be examined from a variety of perspectives ranging from a global, macroeconomic and geopolitical viewpoint down to the details of customer/supplier/market interactions in specific industries. All Prosperity Games TM are unique in that both the game format and the player contributions vary from game to game. This report documents the Industry Partnership Prosperity Game sponsored by the Technology Partnerships and Commercialization Center at Sandia National Laboratories. Players came from the Sandia line organizations, the Sandia business development and technology partnerships organizations, the US Department of Energy, academia, and industry The primary objectives of this game were to: explore ways to increase industry partnerships to meet long-term Sandia goals; improve Sandia business development and marketing strategies and tactics; improve the process by which Sandia develops long-term strategic alliances. The game actions and recommendations of these players provided valuable insights as to what Sandia can do to meet these objectives.
LDRD research activities have focused on increasing the robustness and efficiency of optimization studies for computationally complex engineering problems. Engineering applications can be characterized by extreme computational expense, lack of gradient information, discrete parameters, non-converging simulations, and nonsmooth, multimodal, and discontinuous response variations. Guided by these challenges, the LDRD research activities have developed application-specific techniques, fundamental optimization algorithms, multilevel hybrid and sequential approximate optimization strategies, parallel processing approaches, and automatic differentiation and adjoint augmentation methods. This report surveys these activities and summarizes the key findings and recommendations.
During the 1980s Sandia designed, developed, fabricated, tested, and delivered hundreds of thousands of radiation hardened Integrated Circuits (IC) for use in weapons and satellites. Initially, Sandia carried out all phases, design through delivery, so that development of next generation ICs and production of current generation circuits were carried out simultaneously. All this changed in the mid-eighties when an outside contractor was brought in to produce ICs that Sandia developed, in effect creating a crisp separation between development and production. This partnership had a severe impact on operations, but its more damaging effect was the degradation of Sandia`s microelectronics capabilities. This report outlines microelectronics development and production in the early eighties and summarizes the impact of changing to a separate contractor for production. This record suggests that low volume production be best accomplished within the development organization.
This report describes the results of composites fabrication research sponsored by the Laboratory Directed Research and Development (LDRD) program at Sandia National Laboratories. They have developed, prototyped, and demonstrated the feasibility of a novel robotic technique for rapid fabrication of composite structures. Its chief innovation is that, unlike all other available fabrication methods, it does not require a mold. Instead, the structure is built patch by patch, using a rapidly reconfigurable forming surface, and a robot to position the evolving part. Both of these components are programmable, so only the control software needs to be changed to produce a new shape. Hence it should be possible to automatically program the system to produce a shape directly from an electronic model of it. It is therefore likely that the method will enable faster and less expensive fabrication of composites.
In the assessment of a system, understanding the system is central. Even so, most of the current literature takes a narrow view of understanding, making only the catalog of system ``assets`` explicit, while maintaining the balance of the analyst`s understanding inside the analyst`s head. This can lead to problems with non-repeatability and incompleteness of assessment results. This paper introduces the notion of using explicit system models to document the analyst`s understanding of the system and shows that, from these models, standard assessment products, such as fault trees and event trees, can be automatically derived. This paper also presents five ``views`` of a system that can be used to document the analyst`s understanding of the system. These views go well beyond the standard instruction to identify the system`s assets to show that a much richer understanding of the system can be required for effective assessment.
In parameter estimation considerable insight is provided by examining sensitivity coefficients. This paper focuses on the use of sensitivity coefficients in connection with estimating thermal properties in the heat conduction equation. A general methodology for computing sensitivity coefficients can be an important design tool. The use of such a tool is demonstrated in this paper. A control volume, finite element program is used, and briefly described, to implement numerical sensitivity coefficient calculations. In this approach general problems can be studied. Several example problems are presented to demonstrate the insight gained from sensitivity coefficients. The problems are selected from experimental studies to characterize the thermal properties of carbon-carbon composite. Sensitivity coefficients show that in an experiment that is not well designed, additional materials in the experimental configuration can have a larger impact on the temperature than the material of interest. Two-dimensional configurations demonstrate that there can be isolated areas of insensitivity and the difficulty of estimating multiple parameters.
Effective use of resources that are shared among multiple products or processes is critical for agile manufacturing. This paper describes the development and implementation of a computerized model to support production planning in a complex manufacturing system at the Pantex Plant, a US Department of Energy facility. The model integrates two different production processes (nuclear weapon disposal and stockpile evaluation) that use common facilities and personnel at the plant. The two production processes are characteristic of flow-shop and job shop operations. The model reflects the interactions of scheduling constraints, material flow constraints, and the availability of required technicians and facilities. Operational results show significant productivity increases from use of the model.
Capacitance-voltage and thermally-stimulated-current techniques are used to estimate trapped hole and electron densities in MOS oxides as functions of irradiation and isochronal anneal temperature. Trapped-charge annealing and compensation effects are discussed.
Highly porous sol-gel films have potential applications as electrical and thermal insulators, catalyst supports, sensors, and membranes for gas separations. Pore dimensions in these sol-gel films are usually small e.g., on the order of tens of nanometers or less. Their successful fabrications, however, greatly depend on the fundamental understanding of mechanisms that underlie the phenomena of pore evolution, network shrinkage, and stress development since the final microstructure of a solid gel film is strongly affected by composition of its starting sol and its processing conditions. This report documents a simplified one-dimensional analysis of drying a solidifying sol-gel thin film coating supported by an impermeable solid substrate. Portions of this work were presented at the 1994 Annual Joint Meeting of the New Mexico Section of the American Ceramic Society and Materials Research Society in Albuquerque. The authors considered the solid/liquid two phase coexistent regime during the drying solidifying process in which solvent is removed continuously via evaporation, the solid phase grows significantly in mechanical strength, and pore space shrinks appreciably. From overall and differential mass balances and a force balance at equilibrium, coupled with empirical correlations of solid phase modulus and permeability to strain or deformation, the authors followed the evolution of pore space, solid phase elastic stress, and liquid phase hydrodynamic pressure; they also determined their respective values at equilibrium. By assuming microscopic pore shape models, they estimated and compared the predicted mean pore radii. Their simplified one-dimensional analysis shows that the final mean pore radius is controlled by four parameters: pore-liquid surface tension, solid phase modulus, mean pore radius, and porosity at the initial stress-free state. The one-dimensional model can be employed to guide process design and optimization in sol-gel film fabrications.
Practical work of adhesion measurements are being studied for several types of polymer/metal combinations in order to obtain a better understanding of the adhesive failure mechanisms for systems containing encapsulated and bonded components. The primary question is whether studies of model systems can be extended to systems of technological interest. The authors report on their first attempts to obtain the work of adhesion between a PDMS polymer and stainless steel. The work of adhesion measurements were made using three techniques -- contact angle, adhesive fracture energy at low deformation rates and JKR. Previous work by Whitesides` group show a good correlation between JKR and contact angle measurements for PDMS. Their initial work focused on duplicating the PDMS measurements of Chaudury. In addition, in this paper the authors extend the work of adhesion measurement to third technique -- interfacial failure energy. The ability to determine the reversible work of adhesion for practical adhesive joints allows understanding of several issues that control adhesion: surface preparation, nature of the interphase region, and bond durability.
Gate oxide electric fields are expected to increase to greater than 5 MV/cm as feature size approaches 0.1 micrometers in advanced integrated circuit (IC) technologies. Work by Johnston, et al. raised the concern that single event gate rupture (SEGR) may limit the scaling of advanced ICs for space applications. SEGR has also been observed in field programmable gate arrays, which rely on thin dielectrics for electrical programming at very high electric fields. The focus of this effort is to further explore the mechanisms for SEGR in thin gate oxides. The authors examine the characteristics of heavy ion induced breakdown and compare them to ion induced damage in thin gate oxides. Further, the authors study the impact of precursor damage in oxides on SEGR threshold. Finally, they compare thermal and nitrided oxides to see if SEGR is improved by incorporating nitrogen in the oxide.
Sandia National Laboratories has initiated many joint research and development projects with the two premier Russian nuclear laboratories, VNIIEF and VNIITF, (historically known as Arzamas-16 and Chelyabinsk-70) in a wide spectrum of areas. One of the areas in which critical dialogue and technical exchange is continuing to take place is in the realm of system surety. Activities primarily include either safety or security methodology development, processes, accident environment analyses and testing, accident data-bases, assessments, and product design. Furthermore, a continuing dialog has been established between the organizations with regard to developing a better understanding of how risk is perceived and analyzed in Russia versus that in the US. The result of such efforts could reduce the risk of systems to incur accidents or incidents resulting in high consequences to the public. The purpose of this paper is to provide a current overview of the Sandia surety program and its various initiatives with the Russian institutes, with an emphasis on the program scope and rationale. The historical scope of projects will be indicated. A few specific projects will be discussed, along with results to date. The extension of the joint surety initiatives to other government and industry organizations will be described. This will include the current status of a joint Sandia/VNIIEF initiative to establish an International Surety Center for Energy Intensive and High Consequence Systems and Infrastructures.
Preliminary work is presented on an effort to generate synthetic constitutive data for random composite materials. The long-ranged goal is to use the overall response determined from finite element simulations of representative volumes (RV) of the heterogeneous material to construct a homogenized constitutive model. A simple composite of a matrix containing polydispersed spheres was chosen as the first configuration to simulate. Here the accuracy of the numerical simulation tools is tested by determining effective elastic constants of the ordered elastic composite in which equal-sized spheres are arranged in each of three cubic lattice configurations. The resulting anisotropic effective elastic constant values agree with theoretical results to better than 10%, with typical agreement being better than 4%.
The thermal conduction of a portion of an enhanced surface heat exchanger for a gas fired heat pipe solar receiver was modeled using the boundary element and finite element methods (BEM and FEM) to determine the effect of weld fillet size on performance of a stud welded pin fin. A process that could be utilized by others for designing the surface mesh on an object of interest, performing a conversion from the mesh into the input format utilized by the BEM code, obtaining output on the surface of the object, and displaying visual results was developed. It was determined that the weld fillet on the pin fin significantly enhanced the heat performance, improving the operating margin of the heat exchanger. The performance of the BEM program on the pin fin was measured (as computational time) and used as a performance comparison with the FEM model. Given similar surface element densities, the BEM method took longer to get a solution than the FEM method. The FEM method creates a sparse matrix that scales in storage and computation as the number of nodes (N), whereas the BEM method scales as N{sup 2} in storage and N{sup 3} in computation.
This paper explores one woman`s journey through her recent promotion into management, and will identify key factors that helped prepare and position her to be ready to exercise leadership through a formal management role. It discusses assessment of qualifications and skills, acquisition of needed skills, the influence of luck and timing, and the use of mentors and delegation as survival skills to get through the transition period and become fully functional as a manager. It also includes insights into sensitive issues such as how to relate to former peers, how to gain credibility as the junior member of the management team, and how to juggle family responsibilities with increased time commitments at work. It emphasizes the importance of acknowledging the help the authors receive in reaching their own career goals and offering the same kind of help and support to those in the early stages of their careers.
The traditional method for determination of alpha-emitting isotopes on air filters has been to process the samples by radiochemical methods. However, this method is too slow for incidents involving radioactive materials where the determination of personnel dose is urgent. A method is developed to directly analyze the air filters taken from personal and area air monitors. The site knowledge is used in combination with alpha-spectral information to identify isotopes. A mathematical function is developed to estimate the activity for each isotope. The strengths and weaknesses of the method are discussed.
The reliability of microengines is a function of the design of the mechanical linkage used to connect the electrostatic actuator to the drive. We have completed a series of reliability stress tests on surface micromachined microengines driving an inertial load. In these experiments, we used microengines that had pin mechanisms with guides connecting the drive arms to the electrostatic actuators. Comparing this data to previous results using flexure linkages revealed that the pin linkage design was less reliable. The devices were stressed to failure at eight frequencies, both above and below the measured resonance frequency of the microengine. Significant amounts of wear debris were observed both around the hub and pin joint of the drive gear. Additionally, wear tracks were observed in the area where the moving shuttle rubbed against the guides of the pin linkage. At each frequency, we analyzed the statistical data yielding a lifetime (t50) for median cycles to failure and σ, the shape parameter of the distribution. A model was developed to describe the failure data based on fundamental wear mechanisms and forces exhibited in mechanical resonant systems. The comparison to the model will be discussed.
The monolithic integration of MicroElectroMechanical Systems (MEMS) with the driving, controlling, and signal processing electronics promises to improve the performance of micromechanical devices as well as lower their manufacturing, packaging, and instrumentation costs. Key to this integration is the proper interleaving, combining, and customizing of the manufacturing processes to produce functional integrated micromechanical devices with electronics. We have developed a MEMS-first monolothic integrated process that first seals the micromechanical devices in a planarized trench and then builds the electronics in a conventional CMOS process. To date, most of the research published on this technology has focused on the performance characteristics of the mechanical portion of the devices, with little information on the attributes of the accompanying electronics. This work attempts to reduce this information void by presenting the results of SPICE Level 3 and BSIM3v3.1 model parameters extracted for the CMOS portion of the MEMS-first process. Transistor-level simulations of MOSFET current, capacitance, output resistance, and transconductance versus voltage using the extracted model parameters closely match the measured data. Moreover, in model validation efforts, circuit-level simulation values for the average gate propagation delay in a 101-stage ring oscillator are within 13-18% of the measured data. These results establish the following: (1) the MEMS-first approach produces functional CMOS devices integrated on a single chip with MEMS devices and (2) the devices manufactured in the approach have excellent transistor characteristics. Thus, the MEMS-first approach renders a solid technology foundation for customers designing in the technology.
Proceedings International Carnahan Conference on Security Technology
Garcia, M.L.
The Southwest Surety Institute was formed in June, 1996 by Arizona State University (ASU), New Mexico Institute of Mining and Technology (NM Tech), New Mexico State University (NMSU), and Sandia National Laboratories (SNL) to provide new educational programs in Security Engineering. This is the first science-based program of its kind in the United States, directed at educating Security Engineers to help government and industry address their security needs. Current courses include security system design, evaluation, principles, and technology, the criminal justice system, and each member brings a unique educational capability to the institute. NMSU provides a security technology minor, merging programs in Criminal Justice and Electronics Technology. NM Tech has a formidable explosives testing and evaluation facility. ASU is developing a masters program in Security Engineering at their School of Technology located on a new campus in Mesa, Arizona. The Sandia National Laboratories security system design and evaluation process forms the basis for the security engineering curricula. In an effort to leverage the special capabilities of each university, distance education will be used to share courses among institute members and eventually with other sites across the country. The Institute will also pursue research and development funding in the areas of physical security information security, computer modeling and analysis, and counter-terrorist technology. Individual Institute members are currently working with sponsors from government and industry in areas such as counter-terrorism, microelectronics, banking, aviation, and sensor development.
Design and analysis of physical protection systems requires (1) identification of mission critical assets, (2) identification of potential threats that might undermine mission capability; (3) identification of the consequences of loss of mission-critical assets (e.g., time and cost to recover required capability and impact on operational readiness), and (4) analysis of the effectiveness of physical protection elements. CPA (cost and performance analysis) addresses the fourth of these four issues. CPA is a methodology that joins activity based cost estimation with performance-based analysis of physical protection systems. CPA offers system managers an approach that supports both tactical decision making and strategic planning. Current exploratory applications of the CPA methodology address analysis of alternative conceptual designs. Hypothetical data is used to illustrate this process.
We discuss the development, design and operation of a walk-through trace detection portal designed to screen personnel for explosives. Developed at Sandia National Laboratories (SNL) with primary funding from the Federal Aviation Administration (FAA) and additional support from the Department of Energy Office of Safeguards and Security, this portal is intended primarily for use in airport terminals and in other localities where a very high throughput of pedestrian traffic is combined with stringent security requirements. The portal is capable of detecting both vapor and particulate contamination, with the collection of explosive material being based upon the entrainment of that material in air flows over the body of the person being screened. This portal is capable of detecting high explosives of interest to the FAA. We discuss the results of field testing of the portal in the Albuquerque International Airport in September, 1997 and more recent steps towards commercialization of the portal.
ASME International Mechanical Engineering Congress and Exposition, Proceedings (IMECE)
Wong, C.C.; Hudson, M.L.; Potter, D.L.; Bartel, T.J.
A reliable micro gas pump is an essential element to the development of many micro-systems for chemical gas analyses. At Sandia, we are exploring a different pumping mechanism, gas transport by thermal transpiration. Thermal transpiration refers to the rarefied gas dynamics developed in a micro-channel with a longitudinal temperature gradient. To investigate the potential of thermal transpiration for gas pumping in micro-systems, we have performed simulations and model analysis to design micro-devices and to assess their design performance before the fabrication process. Our effort is to apply ICARUS (a Direct Simulation Monte Carlo code developed at Sandia) to characterize the fluid transport and evaluate the design performance. The design being considered has two plenums at different temperatures (hot and cold) separated by a micro-channel of 0.1 micron wide and 1 micron long. The temperature difference between the two plenums is 30 Kelvin. ICARUS results, a quasi-steady analysis, predicts a net flow through the micro-channel with a velocity magnitude of about 0.4 m/s due to temperature gradient at the wall (thermal creep flow) at the early time. Later as the pressure builds up in the hot plenum, flow is reversed. Eventually when the system reaches steady state equilibrium, the net flow becomes zero. The thermal creep effect is compensated by the thermo-molecular pressure effect. This result demonstrates that it is important to include the thermo-molecular pressure effect when designing a pumping mechanism based on thermal transpiration. The DSMC technique can model this complex thermal transpiration problem.
This paper summarizes results of metal cutting tests using an actively damped boring bar to suppress regenerative chatter. PZT stack actuators were integrated into a commercially available two-inch diameter boring bar to suppress bending vibrations. Since the modified tool requires no specialized mounting hardware, it can be readily mounted on a variety of machines. A cutting test using the prototype bar to remove metal from a hardened steel workpiece verifies that the actively damped tool yields significant vibration reduction and improved surface finish as compared to the open-loop case. In addition, the overall performance of the prototype bar is compared to that of an unmodified bar of pristine geometry, revealing that a significant enlargement of the stable machining envelope is obtained through application of feedback control.
This paper summarizes the design, modeling, and initial evaluation of a hinged structure for friction measurement in surface micromachining technology. While the area requirements are small, the present structure allows a much larger velocity and pressure range to be evaluated as compared to comb drive structures. The device consists of a cantilevered driver beam connected to a friction pad through a strategically located hinge. AC excitation of the beam flexure forces axial sliding of the friction pad due to beam foreshortening. Normal force is controlled by DC voltage on wings adjacent to the friction pad. While the achievable slip is small (10-30 nm), it is sufficient to disengage contacting asperities and engage new points of contact, and thus should be representative of frictional processes. Furthermore, the design enables the friction pad contact area to remain relatively constant over the excitation cycle. Computer simulation results are provided to mimic on-going experimental work. Increased friction forces are shown to enhance the size of hysteresis loops relating beam deflection to driver voltage.
A variety of tomographic techniques that have been applied to multiphase flows are described. The methods discussed include electrical impedance tomography (EIT), magnetic resonance imaging (MRI), positron emission tomography (PET), gamma-densitometry tomography (GDT), radiative particle tracking (RDT), X-ray imaging, and acoustic tomography. Also presented is a case study in which measurements were made with EIT and GDT in two-phase flows. Both solid-liquid and gas-liquid flows were examined. EIT and GDT were applied independently to predict mean and spatially resolved phase volume fractions. The results from the two systems compared well.
ASME International Mechanical Engineering Congress and Exposition, Proceedings (IMECE)
Hudson, Mary L.
A microscale gas chromatography column is one component in a microscale chemistry laboratory for detecting chemical agents. Several columns were fabricated using the Bosch etch process which allows deep, high aspect ratio channels of rectangular cross-section. A design tool, based on analytical models, was developed to evaluate the effects of operating conditions and column specifications on separation resolution and time. The effects of slip flow, channel configuration, and cross-sectional shape were included to evaluate the differences between conventional round, straight columns and the microscale rectangular, spiral columns. Experimental data were obtained and compared with the predicted flowrates and theoretical number of plates. The design tool was then employed to select more optimum channel dimensions and operating conditions for high resolution separations.
The response and operation of a heat flux gage is studied using sensitivity analysis. Sensitivity analysis is the process by which one determines the sensitivity of a model output to changes in the model parameters. This process uses sensitivity coefficients, which are defined as partial derivatives of field variables (e.g. temperature) with respect to model parameters (e.g. thermal properties and boundary conditions). Computing sensitivity coefficients, in addition to the response of a heat flux gage, AIDS in identifying model parameters that significantly impact the temperature response. A control volume finite element based code is used to implement numerical sensitivity coefficient calculations, allowing general problems to be studied. Sensitivity coefficients are discussed for the well known Gardon gage.
Equations are presented to describe the sensitivity of the temperature field in a heat conducting body to changes in the volumetric heat source and volumetric heat capacity. These sensitivity equations, along with others not presented, are applied to a thermal battery problem to compute the sensitivity of the temperature field to 19 model input parameters. Sensitivity coefficients, along with assumed standard deviation in these parameters, are used to estimate the uncertainty in the temperature prediction. From the 19 parameters investigated, the battery cell heat source and volumetric heat capacity were clearly identified as being the major contributors to the overall uncertainty in the temperature predictions. The predicted operational life of the thermal battery was shown to be very sensitive to uncertainty in these parameters.
Research is underway to develop a 75-kW heat pipe to transfer solar energy from the focus of a parabolic dish concentrator to the heater tubes of a Stirling engine. The high flux levels and high total power level encountered in this application have made it necessary to use a high-performance wick structure with fibers on the order of 4 to 8 microns in diameter. This fine wick structure is highly susceptible to corrosion damage and plugging, as dissolved contaminants plate out on the evaporator surface. Normal operation of the heat pipe also tends to concentrate contaminants in localized areas of the evaporator surface where heat fluxes are the highest. Sandia National Laboratories is conducting a systematic study to identify procedures that reduce corrosion and contamination problems in liquid-metal heat pipes. A series of heat pipes are being tested to explore different options for cleaning heat-pipe systems. Models are being developed to help understand the overall importance of operating parameters on the life of heat-pipe systems. In this paper, we present our efforts to reduce corrosion damage.
This study compares the moduli and stresses obtained from dynamic measurements (e.g., logs or tomograms) and static tests (microfracture stress tests, core tri-axial compression tests) at M-Site, where there is a full suite of both types of data as well as other supporting information. The study shows that the dynamic moduli and log-derived stresses are considerably different from the measured in situ values as determined from microfracture stress tests. 2-D images of moduli and stress were also calculated from p-wave and s-wave tomograms, but the primary value of these results is in the qualitative description of the reservoir. The choice of modulus and stress values has a significant effect on processes such as hydraulic fracturing.
We are interested in the stability of three-dimensional fluid flows to small disturbances. One computational approach is to solve a sequence of large sparse generalized eigenvalue problems for the leading modes that arise from discretizating the differential equations modeling the flow. The modes of interest are the eigenvalues of largest real part and their associated eigenvectors. We discuss our work to develop an efficient and reliable eigensolver for use by the massively parallel simulation code MPSalsa. MPSalsa allows simulation of complex 3D fluid flow, heat transfer, and mass transfer with detailed bulk fluid and surface chemical reaction kinetics.
A linear least-squares procedure for the determination of modal residues using time-domain system realization theory is presented. The present procedure is intended to complement existing techniques for time-domain system identification and is shown to be theoretically equivalent to residue determination in realization algorithms such as the eigensystem realization algorithm and Q-Markov covariance equivalent realization method. However, isolating the optimal residue estimation problem from the general realization problem affords several alternative strategies as compared to standard realization algorithms for structural dynamics identification. Primary among these are alternative techniques for handling data sets with large numbers of sensors using small numbers of reference point responses and the inclusion of terms that accurately model the effects of residual flexibility. The accuracy and efficiency of the present realization theory-based procedure is demonstrated for both simulated and experimental data.
In order for the rapidly emerging field of MicroElectroMechanical Systems (MEMS) to meet its extraordinary expectations regarding commercial impact, issues pertaining to how they fail must be understood. We identify failure modes common to a broad range of MEMS actuators, including adhesion (stiction) and friction-induced failures caused by improper operational methods, mechanical instabilities, and electrical instabilities. Demonstrated methods to mitigate these failure modes include implementing optimized designs, model-based operational methods, and chemical surface treatments.
The fluid and particle dynamics of a high-velocity oxygen fuel (HVOF) thermal spray torch are analyzed using computational and experimental techniques. Three-dimensional computational fluid dynamics (CFD) results are presented for a curved aircap used for coating interior surfaces such as engine cylinder bores. The device analyzed is similar to the Metco diamond jet rotating wire (DJRW) torch. The feed gases are injected through an axisymmetric nozzle into the curved aircap. Premixed propylene and oxygen are introduced from an annulus in the nozzle, while cooling air is injected between the nozzle and the interior wall of the aircap. The combustion process is modeled using a single-step, finite-rate chemistry model with a total of nine gas species which includes dissociation of combustion products. A continually fed steel wire passes through the center of the nozzle, and melting occurs at a conical tip near the exit of the aircap. Wire melting is simulated computationally by injecting liquid steel particles into the flow field near the tip of the wire. Experimental particle velocity measurements during wire feed were also taken using a laser two-focus (L2F) velocimeter system. Flow fields inside and outside the aircap are presented, and particle velocity predictions are compared with experimental measurements outside of the aircap.
A methodology is presented that allows the derivation of low-truncation-error finite difference equations for photonics simulation. This methodology is applied to the case of wide-angle beam propagation in two dimensions, resulting in finite difference equations for both TE and TM polarization that are quasi-fourth-order accurate even in the presence of interfaces between dissimilar dielectrics. This accuracy is accomplished without an appreciable increase in numerical overhead and is concretely demonstrated for two test problems having known solutions. These finite difference equations facilitate an approach to the ideal of grid-independent computing and should allow the simulation of relevant photonics devices on personal computers.
Although they appear deceptively simple, batteries embody a complex set of interacting physical and chemical processes. While the discrete engineering characteristics of a battery, such as the physical dimensions of the individual components, are relatively straightforward to define explicitly, their myriad chemical and physical processes, including interactions, are much more difficult to accurately represent. For this reason, development of analytical models that can consistently predict the performance of a battery has only been partially successful, even though significant resources have been applied to this problem. As an alternative approach, we have begun development of non-phenomenological models for battery systems based on artificial neural networks. This paper describes initial feasibility studies as well as current models and makes comparisons between predicted and actual performance.
In this paper, a support and preload system is presented in which the frequencies and damping of the test article are affected by the stiffness and damping of the supporting structure. A dynamic model is derived for the support system that includes the damping as well as the mass and stiffness of the supports. The frequencies, damping, and mode shapes are compared with the experimentally determined parameters. It is shown that for a seemingly simple support system, deriving a predictive model is not a trivial task.
A method is presented for estimating uncertain or unknown parameters in a mathematical model using measurements of transient response. The method is based on a least squares formulation in which the differences between the model and test-based responses are minimized. An application of the method is presented for a nonlinear structural dynamic system. The method is also applied to a model of the Department of Energy armored tractor trailer. For the subject problem, the transient response was generated by driving the vehicle over a bump of prescribed shape and size. Results from the analysis and inspection of the test data revealed that a linear model of the vehicle's suspension is not adequate to accurately predict the response caused by the bump.
Recent advances in Vertical-Cavity Surface-Emitting Laser (VCSEL) technology that have led to higher efficiencies and lower thresholds have opened up a new realm of applications for these devices. In particular, phase-locked arrays of VCSELs1, previously thought to be impractical due to thermal considerations, now look extremely attractive as high-power and highbrightness sources. In addition, a new understanding of waveguiding in VCSELs2 has led to practical methods for designing phase-locked arrays employing either evansecent or leaky-mode (antiguided) coupling. The latter type of coupling is particularly attractive in light of previous calculations1 that predict especially strong mode discrimination against higher-order lateral modes. In this paper we report the first detailed simulation of leaky-mode coupling between two VCSEL pixels performed without the use of simplifying assumptions such as the effective index model. The results of this simulation are, however, found to be in good agreement with previously-developed simple theories3 of leaky-mode coupling.