Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local time stepping”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

ARPA-E PNDIODES "PN Junctions by Ion Implantation" Project (Final Report)

Electricity accounts for nearly 40% of the energy generated in the United States, of which 65% is lost to electrical system inefficiency before reaching end-use sectors. Power electronics will play a growing role in this distribution and end-use, with up to 80% of electricity estimated to pass through power electronics from generation to consumption by 2030. This gives technical improvements in power electronics enormous potential to drive efficiency gains throughout the U.S. economy. Current Si-based power architectures including MOSFETs and IGBTs are limited by high losses, low switching frequency, and poor high-temperature performance. Given the material properties of GaN, devices can exhibit a given turn-on resistance and breakdown voltage with a minute fraction of the area of a Si device—all with higher switching speeds. While n-type doping of III-nitrides has reached a high level of maturity and integration into commercialized products, p-type doping presents distinct challenges that requires continued research. The success of the next generation of GaN-based power devices depends in large part on the ability to form selective-area p-type regions, which can be accomplished by ion implantation of the prevailing acceptor dopant, Mg. This approach is relatively inexpensive and is a fab-ready processing step for creation of current-generation power electronic devices including BJTs, IGBTs, and diode rectifiers. However, ion implantation induces lattice damage and creates point defects within the material, which negatively impacts dopant activation annealing. These forms of damage may be removed by high-temperature annealing, however short annealing time and higher pressures are required to prevent decomposition of the material. We approach this challenge through an innovative combination of approaches toward obtaining high-conductivity p-type GaN. By the novel employment of a gyrotron, a microwave source that can be shaped into a highly localized beam to dissipate very high power in a semiconductor target, a wafer of GaN can be annealed in short (<2 s) ‘pulses’ at temperatures of 1450 °C (2640 °F) to achieve p-type activation without appreciable surface degradation. The p-type activation of Mg-implanted GaN can be further aided by co-implantation with N, believed to enhance the removal of point defects that compensate the acceptors dopant. Elevation to such high temperatures for even a short time induces degradation of GaN-to overcome, we have implemented additional measures to retain high-quality crystalline GaN while removing implant damage and activating Mg dopant. A bilayer cap of AlN is deposited by a combination of chemical and physical vapor deposition to protect the GaN surface. To prevent decomposition of GaN loss of N during annealing, a custom vessel—to contain pressurized N while allowing injection of microwave energy—is employed. By detailed study of the interaction of various measures that are implemented to protect the sample surface (GaN), activate the implanted dopants, and mitigate damage of the lattice heating damage during annealing, we are able to achieve p-type GaN and observe recovery of the crystal and removal of substantial density of defects created during the implantation process.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mapping 3D grain and precipitate structure during in situ mechanical testing of open-cell metal foam using micro-computed tomography and high-energy X-ray diffraction microscopy

Open-cell metal foams are ultra-low-density cellular metals with complex hierarchical structures that span bulk, cell, ligament, and sub-ligament scales and give rise to desirable properties such as high strength-to-weight ratio and excellent energy absorption. Although literature suggests that intrinsic material structures at sub-ligament length scales (e.g., grains and precipitates) play an important role in mechanical behavior of open-cell metal foams, there are very few experimental measurements of such structures in three dimensions and for meaningful volumes of foam. This study seeks to map and track the three-dimensional (3D) grain and precipitate structures of an intact volume of open-cell aluminum foam by advancing microstructural characterization techniques that leverage X-ray micro-computed tomography (μCT) and far-field high-energy X-ray diffraction microscopy (FF-HEDM). A 6%-relative-density aluminum foam sample was mechanically tested in compression while μCT and FF-HEDM measurements were collected at interrupted loading states at beamline 1-ID of the Advanced Photon Source. Further, a new scanning strategy and reconstruction algorithm were established to enable characterization of a foam volume with diameter approximately four times wider than the nominal width of the X-ray beam. The result is a set of maps that detail both the 3D grain and precipitate structures throughout the foam volume at successive strain steps. A novel grain tracking procedure was developed to track individual grains within the foam volume by accounting for the large rigid-body motions that individual ligaments can undergo during mechanical loading. The ability to track grains and precipitate structures in three dimensions throughout large bulk deformation of ultra-low-density polycrystalline materials enables new possibilities for validating numerical models and investigating local failure mechanisms. Furthermore, the methods and procedures developed in this study could be applied to other ultra-low-density structures, such as additively manufactured lattices.

36 MATERIALS SCIENCE↗

On Lateral Viscosity Contrast in the Mantle and the Rheology of Low-Frequency Geodynamics

Mantle-wide heterogeneity is largely controlled by deeply penetrating thermal convective currents. These thermal currents are likely to produce significant lateral variation in rheology, and this can profoundly influence overall material behaviour. How thermally related lateral viscosity variations impact models of glacio-isostatic and tidal deformation is largely unknown. An important step towards model improvement is to quantify, or bound, the actual viscosity variations that characterize the mantle. Simple scaling of viscosity to shear-wave velocity fluctuations yields map-views of long- wavelength viscosity variation. These give a general quantitative description and aid in estimating the depth dependence of rheological heterogeneity throughout the mantle. The upper mantle is probably characterized by two to four orders of magnitude variation (peak-to-peak). Discrepant time-scales for rebounding Holocene shorelines of Hudson Bay and southern Iceland are consistent with this characterization. Results are given in terms of a local average viscosity ratio, (Delta)eta(bar)(sub i), of volumetric concentration, phi(sub i). For the upper mantle deeper than 340 km the following reasonable limits are estimated for (delta)eta(bar) approx. equal 10(exp -2): 0.01 less than or equal to phi less than or equal to 0.15. A spectrum of ratios (Delta)eta(bar)(sub i) less than 0.1 at concentration level eta(sub i) approx. equal 10(exp -6) - 10(exp -1) in the lower mantle implies a spectrum of shorter time-scale deformational response modes for second-degree spherical harmonic deformations of the Earth. Although highly uncertain, this spectrum of spatial variation allows a purely Maxwellian viscoelastic rheology simultaneously to explain all solid tidal dispersion phenomena and long-term rebound-related mantle viscosity. Composite theory of multiphase viscoelastic media is used to demonstrate this effect.

Ivins, Erik R.↗

Historical Notes on the Expanding Universe

The article Measuring the Hubble constant by Mario Livio and Adam Riess (Physics Today, October 2013, page 41) reviewed studies of the expanding universe from the 1920s to the present. Although the history of the subject underwent considerable compression to fit the length of a magazine article, we think it may leave a misleading impression of some of the key steps to our current understanding. We therefore offer the following clarifications. Most significantly, papers by Arthur Eddington and by Willem de Sitter in 1930, who successfully promoted Georges Lematres 1927 article for the Scientific Society of Brussels, effected a paradigm shift in interpretation of extragalactic redshifts in 1930. Before then, the astronomical community was generally unaware of the existence of nonstatic cosmological solutions and did not broadly appreciate that redshifts could be thought of locally as Doppler shifts in an expanding matter distribution. Certainly, in 1929 Edwin Hubble referred only to the de Sitter solution of 1917. At the time, the relation between distance and redshift predicted in that model was generally seen purely as a manifestation of static spacetime curvature.

Hubble constant↗

Valuing Resilience for Microgrids: Challenges, Innovative Approaches, and State Needs

The United States depends on the delivery of reliable, affordable, clean, and safe electricity. Electric utilities invest billions of dollars each year in generation, transmission, and distribution assets to meet this need. However, experiences with recent natural disasters of increasing frequency and duration demonstrate the shortcomings of this approach in the face of modern threats. Further, as customers rely on electricity for a broader range of important needs, such as transportation, as well as critical life-saving services and mission critical facilities such as water treatment, medical care, shelters, telecommunications, and more, the need to minimize the likelihood and impacts of outages grows. Against this backdrop, resilience has emerged as a key consideration to guide electricity spending, whether from utilities, customers, or taxpayers. Although reliability has been defined and measured for decades with broadly accepted metrics that measure how many customers lose power and at what frequency and duration, resilience considers the electricity system’s response to a disruption and its subsequent impacts on customers. Developing tools and methods to accurately assess the costs and benefits of resilience investments is a critical step toward the goal of mitigating the impacts of outages on customers and society. Today, electric system resilience is largely treated as an externality due to challenges estimating the costs of long-duration outages, impacts of outages on society, and increasing reliance on electricity for a growing set of interdependent services. These interdependencies include the water, wastewater, telecommunications, natural gas, and health sectors. Without knowing how much a given resilience investment will benefit customers or society more broadly, investors, policymakers, and regulators are less likely to make or approve such investments, and less able to prioritize those investments. State Energy Offices and public utility commissions (PUCs) lead the development of state-level energy policy and utility regulation, respectively, and each has interests in encouraging appropriate public and private investments in resilience. To this end, the National Association of State Energy Officials (NASEO) and National Association of Regulatory Utility Commissioners (NARUC), with the support of the U.S. Department of Energy (DOE) Office of Electricity (OE), formed a joint Microgrids State Working Group to explore the costs and benefits of microgrids, barriers to broader deployment of microgrids to meet resilience and other objectives, and policy and regulatory strategies to optimize investments in resilience, including but not limited to microgrids. Although no universally accepted valuation tool for resilience exists, National Laboratories, utilities, researchers, and state and federal agencies have collaborated to develop, apply, and improve a number of approaches to quantify resilience, several of which are still in progress at the time this report is published. This report seeks to share these important advances by discussing current definitions of resilience (Section 1), how microgrids are defined and used to meet resilience objectives (Section 2), new approaches to valuing resilience (Section 3), steps State Energy Offices and PUCs have taken to further resilience valuation efforts (Section 4), and finally, considerations and suggested next steps for State Energy Offices and PUCs (Section 5). Relevant examples of specific microgrid projects and resilience valuation efforts are included throughout the report. While this report is written specifically for NASEO and NARUC members, it may be useful for utilities, local governments, and individual customers interested in improving the way public and private dollars are spent to achieve resilience outcomes.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Energy considerations in computational aeroacoustics

A finite-volume multistage time-stepping Euler code is used to investigate the use of CFD algorithms for the direct calculation of acoustics. The 2D compressible inviscid flow about an accelerating or decelerating circular cylinder is used as a model problem. The time evolution of the energy transfer from the cylinder to the fluid, as the cylinder is moved from rest to some nonnegligible velocity, is clearly seen. By examining the temporal and spatial characteristics of the numerical solution, a distinction can be made between the propagating acoustic energy, the convecting energy associated with the entropy change in the fluid, and the energy contained in the local aerodynamic field. Systematic variation of the cylinder acceleration shows that the radiated acoustic energy depends strongly upon the rate of acceleration or deceleration. The computational grid has a large effect on the ratio of acoustic energy to nonphysical entropy associated energy, while the role of the explicit artificial viscosity seems to be of second order. The entropy term was nearly negligible in all cases the cylinder was started slowly.

Brentner, Kenneth S.↗

SPARCLE Optical System Design and Operational Characteristics

The SPAce Readiness Coherent Lidar Experiment (SPARCLE) is the first demonstration of a coherent Doppler wind lidar in space. SPARCLE will be flown aboard a space shuttle In the middle part of 2001 as a stepping stone towards the development and deployment of a long-life-time operational instrument in the later part of next decade. SPARCLE is an ambitious project that is intended to evaluate the suitability of coherent lidar for wind measurements, demonstrate the maturity of the technology for space application, and provide a useable data set for model development and validation. This paper describes the SPARCLE's optical system design, fabrication methods, assembly and alignment techniques, and its anticipated operational characteristics. Coherent detection is highly sensitive to aberrations in the signal phase front, and to relative alignment between the signal and the local oscillator beams. Consequently, the performance of coherent lidars is usually limited by the optical quality of the transmitter/receiver optical system. For SPARCLE having a relatively large aperture (25 cm) and a very long operating range (400 km), compared to the previously developed 2-micron coherent lidars, the optical performance requirements are even more stringent. In addition with stringent performance requirements, the physical and environment constraints associated with this instrument further challenge the limit of optical fabrication technologies.

Amzajerdian, Farzin↗

Forecasting Global Fire Emissions on Subseasonal to Seasonal (S2S) Time Scales

Fire emissions of gases and aerosols alter atmospheric composition and have substantial impacts on climate, ecosystem function, and human health. Warming climate and human expansion in fire‐prone landscapes exacerbate fire impacts and call for more effective management tools. Here we developed a global fire forecasting system that predicts monthly emissions using past fire data and climate variables for lead times of 1 to 6 months. Using monthly fire emissions from the Global Fire Emissions Database (GFED) as the prediction target, we fit a statistical time series model, the Autoregressive Integrated Moving Average model with eXogenous variables (ARIMAX), in over 1,300 different fire regions. Optimized parameters were then used to forecast future emissions. The forecast system took into account information about region‐specific seasonality, long‐term trends, recent fire observations, and climate drivers representing both large‐scale climate variability and local fire weather. We cross‐validated the forecast skill of the system with different combinations of predictors and forecast lead times. The reference model, which combined endogenous and exogenous predictors with a 1 month forecast lead time, explained 52% of the variability in the global fire emissions anomaly, considerably exceeding the performance of a reference model that assumed persistent emissions during the forecast period. The system also successfully resolved detailed spatial patterns of fire emissions anomalies in regions with significant fire activity. This study bridges the gap between the efforts of near‐real‐time fire forecasts and seasonal fire outlooks and represents a step toward establishing an operational global fire, smoke, and carbon cycle forecasting system.

vapor pressure↗

Improvements to Stitching Controls for Manufacturing Advanced Stitched Composites

Stitched composites have been extensively studied and have been shown to have benefits over unstitched composites for stiffened structures, including improved damage tolerance, reduced weight, and fewer fasteners. These benefits make stitched composites attractive from a performance perspective, however, the nature of stitching the preforms prior to infusion includes additional manufacturing steps and equipment compared to conventional resin infusion. Stitching of composite material structures at NASA Langley Research Center (LaRC) is performed at the Integrated Structural Assembly of Advanced Composites (ISAAC) facility using control algorithms developed by the lead author. Two deficiencies have been identified in the current control algorithms, and mitigations for these deficiencies were developed within the Hi-Rate Composite Aircraft Manufacturing (HiCAM) project and are presented herein. First, decoupling stitch translation and rotation takes more time for stitching since two robot operation control steps are executed for each stitch. Therefore, a method was developed to determine when to decouple the translation and rotation for stitching, which reduces stitch time by decoupling the operation control steps only when necessary to execute the stitch properly. Second, a correction to the sideslip stitching approach that eliminates a significant portion of the pre- and post-stitching activities was developed. The developed method adjusts the sideslip stitch angle based on the rotations of previous stitches to ensure that the width of the stitching seam remains constant along its length. Constant seam width eliminates observed overlap of the insertion and catcher thread portions for highly curved paths, which maintains the desired uniformity of through-the-thickness architecture along the seam and eliminates the local degradation of properties where the insertion needle is too close to, or crosses, the catcher needle portion of the seam.

stitching↗

Why Seeding Works When Nucleation Barriers Vanish

Crystallization is a process governed by the interplay between nucleation and growth. While crystalline seeds are known to reduce nucleation barriers and accelerate crystallization under nucleation-limited conditions, their influence when nucleation is not the limiting step remains poorly understood. This creates a mechanistic puzzle in systems where nucleation barriers are already negligible, yet seeding still accelerates crystallization. Here, the synthesis of zeolites is a quintessential example of growth-limited crystallization in which seeds accelerate the process despite negligible homogeneous nucleation barriers. Using coarse-grained molecular dynamics simulations─validated across two zeolites and the unrelated case of ice crystallization─we establish that growth-limited crystallization produces many small, misoriented crystallites whose slow coarsening into larger domains controls the emergence of X-ray-detectable crystallinity. Local crystalline order, structural coherence, and X-ray detectability are therefore kinetically decoupled milestones: the first can be reached rapidly while the latter two lag significantly. Seeds resolve this lag by imposing a common orientational registry on nascent crystallites, enabling their coherent coalescence into large seed-bound domains. This early coherence-building step produces a crystallite-size asymmetry that accelerates subsequent coarsening, advancing the onset of X-ray-detectable crystallinity without necessarily increasing the nucleated fraction. We conclude that under growth-limited conditions, the apparent induction period observed in powder X-ray diffraction reflects the time required to build long-range coherence, not the time to form crystalline material. Accordingly, seeds function not primarily by reducing nucleation barriers but by enforcing spatial coherence, thereby shortening the time required to develop long-range order detectable by X-ray diffraction.

Growth-Limited Crystallization↗

Knowledge-driven design of solid-electrolyte interphases on lithium metal via multiscale modelling

Abstract Due to its high energy density, lithium metal is a promising electrode for future energy storage. However, its practical capacity, cyclability and safety heavily depend on controlling its reactivity in contact with liquid electrolytes, which leads to the formation of a solid electrolyte interphase (SEI). In particular, there is a lack of fundamental mechanistic understanding of how the electrolyte composition impacts the SEI formation and its governing processes. Here, we present an in-depth model-based analysis of the initial SEI formation on lithium metal in a carbonate-based electrolyte. Thereby we reach for significantly larger length and time scales than comparable molecular dynamic studies. Our multiscale kinetic Monte Carlo/continuum model shows a layered, mostly inorganic SEI consisting of LiF on top of Li 2 CO 3 and Li after 1 µs. Its formation is traced back to a complex interplay of various electrolyte and salt decomposition processes. We further reveal that low local Li + concentrations result in a more mosaic-like, partly organic SEI and that a faster passivation of the lithium metal surface can be achieved by increasing the salt concentration. Based on this we suggest design strategies for SEI on lithium metal and make an important step towards knowledge-driven SEI engineering.

25 ENERGY STORAGE↗

Microstrain screening towards defect-less layered transition metal oxide cathodes

Microstrain and the associated surface-to-bulk propagation of structural defects are known to be major roadblocks to developing high-energy and long-life batteries. However, the origin and effects of microstrain during the synthesis of battery materials remain largely unknown. Here, in this study, we perform microstrain screening during real-time and realistic synthesis of sodium layered oxide cathodes. Evidence gathered from multiscale in situ synchrotron X-ray diffraction and microscopy characterization collectively reveals that the spatial distribution of transition metals within individual precursor particles strongly governs the nanoscale phase transformation, local charge heterogeneity and accumulation of microstrain during synthesis. This unexpected dominance of transition metals results in a counterintuitive outward propagation of defect nucleation and growth. These insights direct a more rational synthesis route to reduce the microstrain and crystallographic defects within the bulk lattice, leading to significantly improved structural stability. The present work on microstrain screening represents a critical step towards synthesis-by-design of defect-less battery materials.

25 ENERGY STORAGE↗

Pursuing graphite-based K-ion O 2 batteries: a lesson from Li-ion batteries

The replacement of lithium metal by lithium intercalated graphite was crucial for developing safe lithium-ion batteries. Superoxide-based potassium–oxygen batteries represent an exciting metal–oxygen battery with the highest energy efficiencies. However, using potassium metal is a more serious safety threat than lithium. Herein, we explore the possibility of graphite-intercalation anodes for potassium-ion oxygen batteries (PIOBs) with enhanced safety. This work demonstrates for the first time that establishing an artificial potassium salt-rich solid electrolyte interphase (SEI) enables a reversible graphite-intercalation anode (249.6 mA h g -1 after 600 cycles) in a potassium bis(trifluoromethanesulfonyl)imide (KTFSI)-based localized high-concentration electrolyte. Such an electrolyte is stable with the superoxide cathode. The PIOB delivers energy efficiencies above 90% at a depth of discharge (DOD) of 25% for 80 cycles. A three-electrode measurement shows that its overpotential mainly comes from the anode. The lifespan is limited to the gradual degradation of the artificial SEI caused by oxygen crossover. This work represents a step towards achieving holistic anode–electrolyte–cathode compatibility in a PIOB and its realistic evaluation under a controlled DOD.

25 ENERGY STORAGE↗

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extreme Scale Infrasound Inversion and Prediction for Weather Characterization and Acute Event Detection

Accurate and timely weather predictions are critical to many aspects of society with a profound impact on our economy, general well-being, and national security. In particular, our ability to forecast severe weather systems is necessary to avoid injuries and fatalities, but also important to minimize infrastructure damage and maximize mitigation strategies. The weather community has developed a range of sophisticated numerical models that are executed at various spatial and temporal scales in an attempt to issue global, regional, and local forecasts in pseudo real time. The accuracy however depends on the time period of the forecast, the nonlinearities of the dynamics, and the target spatial resolution. Significant uncertainties plague these predictions including errors in initial conditions, material properties, data, and model approximations. To address these shortcomings, a continuous data collection occurs at an effort level that is even larger than the modeling process. It has been demonstrated that the accuracy of the predictions depends on the quality of the data and is independent to a certain extent on the sophistication of the numerical models. Data assimilation has become one of the more critical steps in the overall weather prediction business and consequently substantial improvements in the quality of the data would have transformational benefits. This paper describes the use of infrasound inversion technology, enabled through exascale computing, that could potentially achieve orders of magnitude improvement in data quality and therefore transform weather predictions with significant impact on many aspects of our society.

58 GEOSCIENCES↗

AMPPED LDRD: Advanced Models of the Physics & Phenomena of Electrical Discharge.

This Laboratory Directed Research and Development (LDRD) project aimed at addressing an unmet need for physics-based models that enable prediction of electrical breakdown in the presence of high voltage gradients and real-world environments. This project focused on understanding coupled plasma-surface feedback phenomena critical to electrical discharge processes and design of breakdown voltages, breakdown times, and breakdown paths. New secondary electron emission and coupled plasma-surface measurement methods and models were developed. Highly reproducible discharge experiments were developed to provide current-voltage data as a function of time, and methods were developed to condition electrodes for predictable discharge. To assess plasma-surface feedback mechanisms, new surface physics capabilities, including measurement systems for local and global photoemission yield and ion-induced secondary electron emission yield were developed. Discharge models were developed to incorporate these yields for predicting discharge behavior to understand influence of surface physics on discharge current and time evolution, and methods to accelerate these simulations were demonstrated. A key finding was that the initiation process of a non-vacuum discharge appears highly dependent on photoemission at the cathode. Advanced diagnostics were developed, including spatial electrical field imaging, spectroscopy of plasmas, and high-speed imaging of discharge event s in metal and dielectric environments. These diagnostics and improved electrical models are a step toward enabling quantitative design margins, assurance of safety architectures, and implementation of novel materials and technologies for wide ranging energy and electronics applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗

Laser-induced digital oxidation for copper-based flexible photodetectors

Copper oxide compounds (Cu x O) with bandgaps of 1.3–2.1 eV (CuO) and 2.1–2.6 eV (Cu 2 O) have been investigated as promising p-type semiconducting materials. Cu x O is generally obtained by deposition or thermal oxidation, but those methods are not optimal for flexible substrates. Furthermore, additional patterning steps are required to fabricate devices. In this work, we present an easy, controllable method to fabricate a metal-semiconductor-metal (MSM) photodetector using laser-induced oxidation of a thin Cu film. After laser irradiation, the Cu film is heated under ambient conditions, and this leads to a thermal oxidation reaction, in which Cu oxide is monolithically formed in the Cu film and a Cu-Cu x O-Cu MSM structure is produced. Since the laser offers localized heating, an arbitrary CuxO pattern can be written in the Cu film by spatially controlled heating. In addition, by optimizing the heating time, the laser-induced oxidation can be successfully performed even on a flexible substrate. To study the laser-induced oxidation, we examined the correlation between laser parameters and the oxidation pattern and analyzed the composition using scanning electron microscopy, Raman spectroscopy, and X-ray photoemission spectroscopy. Furthermore, we measured the transient photoresponse and employed scanning photocurrent microscopy to investigate the mechanism of carrier transport behavior.

36 MATERIALS SCIENCE↗