Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oscillatory response”

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.

266 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Characteristic terahertz emissions induced by optically excited collective orbital modes

Here, we study the generation of collective orbital modes, their evolution, and the characteristic nonlinear optical response induced by them in a photoinduced orbital-ordered correlated oxide using real-time simulations based on an interacting multiband tight-binding (TB) model. The d-d optical transitions under femtoseconds light-pulse in an orbital-ordered state excite collective orbital modes, also known as “orbitons”. Consistently incorporating electronic interactions and the interplay between charge, spin, and lattice degrees of freedom in the TB-model provides a clearer understanding of how these factors influence the generation and evolution of collective orbital modes. The dynamics of Jahn-Teller vibrational modes in the photoinduced state modify the intersite orbital interaction, which further amplifies these orbital modes. In the presence of weak ferroelectricity, the excitation of collective orbital modes induces a strong THz oscillatory photocurrent, which is long-lived. This suggests an alternative way to experimentally detect low-energy collective modes through THz-emission studies in the photoinduced state. Our study also elucidates that quasiparticle dynamics in improper ferroelectric oxides can be exploited to achieve highly interesting and nontrivial optoelectronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗
NASA NTRS2002

Characterization of Truncated Tumor-Associated NADH Oxidase (ttNOX)

Bacterial, plant and animal cells possess novel surface proteins that exhibit both NADH oxidation (NOX) or hydroquinone and protein disulfide-thiol interchange. These enzymatic activities alternate to yield oscillating patterns wjth period lengths of approximately 24 minutes. The catalytic period of NOX proteins are temperature compensated and gravity responsive. We report the cloning, expression and characterization of truncated tumor-associated NADH oxidase (ttNOX), in which the membrane spanning region has been deleted. The cDNA (originated from HeLa cells) was cloned into pET-34b and pET-14b (Novagen) vectors for E. coli expression. Optimized expression and purification protocols yielded greater than 300mg per liter of culture with greater than 95% purity. Circular dichroism data was collected from a 2.7mg/ml solution in a 0.1mm cuvette with variable scanning using an Olis RSM CD spectrophotometer. The ellipticity values were scanned from 190 to 260nm. The spectra recorded have characteristics for alpha proteins with band maxima at 216nm and a possible shoulder at 212nm at 12OC and 250 C. Protein crystal screens are in progress and, to date, only small crystals have been observed. The regular periodic oscillatory change in the ttNOX protein is indicative of a possible time-keeping functional role. A single protein possessing alternating catalytic activities, with a potential biological clock function, is unprecedented and structural determination is paramount to understanding this role.

Karr, Laurel J.↗
NASA NTRSDate not supplied

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗
DOE OSTI2023

Neutron calorimeter development at fission and fusion pulsed neutron sources

A neutron calorimeter fielded in pulsed neutron irradiation experiments at the National Ignition Facility (NIF) and White Sands Missile Range Fast Burst Reactor (WSMR FBR) is described, and results are provided and discussed. The calorimeter is modular in design to allow fielding of various neutron absorber materials and thicknesses, as well as use at both fission and fusion sources, all of which has been done. Neutron calorimeters designed for fission sources have been fielded at fusion sources as well, but these instances are not documented in available literature—making this work novel. The functional calorimeter consists of a neutron absorber and two sheathed, grounded thermocouples secured within a subassembly made of thermally insulating materials. This subassembly is placed into a sealed vessel for material and facility safety considerations. The device was fielded at NIF in four polar-direct-drive exploding pusher shots and eleven WSMR FBR pulses. Calorimeter data from NIF experiments show a discernable thermal response accredited to neutron heating within the first 10 s, but the data are dominated by noise from ambient X-ray heating after time = 10 s post-pulse. This response is highly consistent between shots and linearly related to incident neutron fluence and neutron absorber thickness. A temperature rise of 1.1E-2 K/(10 13 neutrons/cm 2 ) per millimeter of absorber thickness was determined for both 0.4 mm and 1 mm thick depleted uranium absorbers with linear best fit R 2 of 0.97 and 0.98 for each thickness and y-axis offsets of 8.8E-3 and 1.3E-2, respectively. WSMR FBR experiment data show a temperature rise of 6.2E-2 K/(10 13 neutrons/cm 2 ) with linear best fit R 2 of 0.14 and y-axis offset of 7.6E-3 for a calorimeter with a 1 mm thick enriched boron carbide ( 10 B 4 C) absorber, though it must be noted that device-to-device variability is significant at both facilities. So these data are also affected by noise on a similar timescale as at NIF, though the noise is significantly lesser in magnitude but oscillatory.

42 ENGINEERING↗
NASA NTRS2004

Oscillatory Hierarchy Controlling Cortical Excitability and Stimulus Integration

Cortical gamma band oscillations have been recorded in sensory cortices of cats and monkeys, and are thought to aid in perceptual binding. Gamma activity has also been recorded in the rat hippocampus and entorhinal cortex, where it has been shown, that field gamma power is modulated at theta frequency. Since the power of gamma activity in the sensory cortices is not constant (gamma-bursts). we decided to examine the relationship between gamma power and the phase of low frequency oscillation in the auditory cortex of the awake macaque. Macaque monkeys were surgically prepared for chronic awake electrophysiological recording. During the time of the experiments. linear array multielectrodes were inserted in area AI to obtain laminar current source density (CSD) and multiunit activity profiles. Instantaneous theta and gamma power and phase was extracted by applying the Morlet wavelet transformation to the CSD. Gamma power was averaged for every 1 degree of low frequency oscillations to calculate power-phase relation. Both gamma and theta-delta power are largest in the supragranular layers. Power modulation of gamma activity is phase locked to spontaneous, as well as stimulus-related local theta and delta field oscillations. Our analysis also revealed that the power of theta oscillations is always largest at a certain phase of delta oscillation. Auditory stimuli produce evoked responses in the theta band (Le., there is pre- to post-stimulus addition of theta power), but there is also indication that stimuli may cause partial phase re-setting of spontaneous delta (and thus also theta and gamma) oscillations. We also show that spontaneous oscillations might play a role in the processing of incoming sensory signals by 'preparing' the cortex.

Shah, A. S.↗
DOE OSTI2020

Optimization-Based Fast-Frequency Estimation and Control of Low-Inertia Microgrids

The lack of inertial response from non-synchronous, inverter-based generation in microgrids makes the power system vulnerable to a large rate of change of frequency (ROCOF) and frequency excursions. Energy storage systems (ESSs) can be utilized to provide fast-frequency support to prevent such large excursions in the system. However, fast-frequency support is a power-intensive application that has a significant impact on the ESS lifetime. In this paper, a framework that allows the ESS operator to provide fast-frequency support as a service is proposed. The framework maintains the desired quality-of-service (limiting the ROCOF and frequency) while taking into account the ESS lifetime and physical limits. The framework utilizes moving horizon estimation (MHE) to estimate the frequency deviation and ROCOF from noisy phase-locked loop (PLL) measurements. These estimates are employed by a model predictive control (MPC) algorithm that computes control actions by solving a finite-horizon, online optimization problem. Additionally, this approach avoids oscillatory behavior induced by delays that are common when using low-pass filters as with traditional derivative-based (virtual inertia) controllers. MATLAB/Simulink simulations on a test system from Cordova, Alaska, show the effectiveness of the MHE-MPC approach to reduce frequency deviations and ROCOF of a low-inertia microgrid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2005

Dynamic Weakening (Extinction) of Simple Hydrocarbon-air Counterflow Diffusion Flames by Oscillatory Inflows

This study of laminar non-premixed HC-air flames used an Oscillatory-input Opposed Jet Burner (OOJB) system developed from a previously well-characterized 7.2-mm Pyrex-nozzle OJB system. Over 600 dynamic Flame Strength (FS) measurements were obtained on unanchored (free-floating) laminar Counterflow Diffusion Flames (CFDF's). Flames were stabilized using plug inflows having steady-plus-sinusoidal axial velocities of varied magnitude, frequency, f, up to 1600 Hz, and phase angle from 0 (most data) to 360 degrees. Dynamic FS is defined as the maximum average air input velocity (U(sub air), at nozzle exit) a CFDF can sustain before strain-induced extinction occurs due to prescribed oscillatory peak-to-peak velocity inputs superimposed on steady inputs. Initially, dynamic flame extinction data were obtained at low f, and were supported by 25-120 Hz Hot-Wire cold-flow velocity data at nozzle exits. Later, expanded extinction data were supported by 4-1600 Hz Probe Microphone (PM) pk/pk P data at nozzle exits. The PM data were first obtained without flows, and later with cold stagnating flows, which better represent speaker-diaphragm dynamics during runs. The PM approach enabled characterizations of Dynamic Flame Weakening (DFW) of CFDF's from 8 to 1600 Hz. DFW was defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = - 100 d(U(sub air) / U(sub air),0Hz) / d(pkpk P). The linear normalization with respect to acoustic pressure magnitude (and steady state (SS) FS) led to a DFW unaffected by strong internal resonances. For the C2H4/N2-air system, from 8 to 20 Hz, DFW is constant at 8.52 plus or minus 0.20 (% weakening)/Pa. This reflects a quasi-steady flame response to an acoustically induced dU(sub air)/dP. Also, it is surprisingly independent of C2H4/N2 mole fraction due to normalization by SS FS. From 20 to approximately 150 Hz, the C2H4/N2 air-flames weakened progressively less, with an inflection at approximately 70 Hz, and became asymptotically insensitive (DFW approximately 0) at approximately 300 Hz, which continued to 1600 Hz. The DFW of CH4-air flames followed a similar pattern, but showed much greater weakening than C2H4/N2-air flames; i.e., the quasi-steady DFW (8 to approximately 15 Hz) was 44.3 %/Pa, or approximately 5x larger, even though the 0 Hz (SS) FS was only 3.0 x smaller. The quasi-steady DFW's of C3H8-air and C2H6-air were intermediate at 34.8 and 20.9 %Pa, respectively. The DFW profiles of all four fuels, at various frequencies, correlated well but non-linearly with respective SS FS's. Notably, the DFW profile for C3H8 air fell more rapidly in the range greater than 15 to 60 Hz, compared with the 1- and 2-carbon fuels. This may indicate a shift in chemical kinetics, and/or O2 transport to a flame that moved closer to the fuel-side. In conclusion, Dynamic Flame Weakening limits appear significant and unique for each fuel, and correlate closely, but non-linearly, with Steady-State Flame Strengths at any given frequency. For reasons unknown, the dynamic flames didn't weaken more at intermediate frequencies (e.g., at 20-50 Hz) than they did at low frequencies (less than 15 Hz), where quasi-steady weakening appears to dominate. Quasi-steady flame weakening ostensibly represents a transient input strain rate maximum that just exceeds the steady-state strain-rate-limited extinction limit for a few cycles. Clearly, further detailed mechanistic understanding is needed in the fall-off region.

Pellett, G.↗
NASA NTRS1991

The sun and interplanetary magnetic field

The interplanetary magnetic field (IMF) serves as a link between the sun, the response of the earth to solar activity and variations in galactic cosmic radiation. The IMF originates as a solar-coronal magnetic field that is transported into space by the solar wind. The close connection between solar magnetic fields and the origin and structure of the solar wind is described. The solar wind forms the heliosphere, a cavity containing the magnetized solar plasma from which the interstellar plasma and field are excluded. The entry of galactic cosmic rays into the heliosphere and their strong interaction with the IMF are discussed, this topic being of primary importance to the production and temporal variations of radiogenic elements. The profound influence of the IMF on geomagnetic activity and the aurora is discussed within the context of merging or reconnection with the planetary field. The physical connection is thus established between solar magnetic fields, magnetic storms and aurora. The state of the solar wind and IMF during the Maunder minimum is considered and an explanation for the (relative) absence of sunspots and aurora is proposed. The mechanism is an interruption of the oscillatory solar dynamo, a consequent reduction in the heating of the corona, a cessation of the supersonic solar wind and a weakening or absence of southward-directed magnetic fields in the vicinity of the earth.

Smith, Edward J.↗
NASA NTRSDate not supplied

Dynamics of A Vibration Isolation System Including Inertia of the Human Body

Using an exercise device in a spacecraft is liable to transmit an unacceptable amount of vibration to that vehicle. This is commonly mitigated by a Vibration Isolation System (VIS), whose dynamics must be analyzed to confirm that the oscillatory forces on the spacecraft remain within allowed range, both from a structural and microgravity perspective (see, e.g., [1]). When modeling a VIS for countermeasures devices, one common approach is to record forces and moments applied on the floor while exercising, and then drive the VIS simulation by applying these recorded loads to the exercise platform part of the VIS mechanical model. This approach misses the fact that when exercising on a moving platform, the force and moment on it will differ from that on the stationary floor due to inertial effects involving the human body. For example, standing up on a platform as it gives under the subject’s feet reduces the foot force on it, and such inertial effects are especially complex for rotational motion. In principle, one could model both the motion of the human body and dynamics of the VIS mechanism in a single combined simulation, e.g., employing a tool such as the commonly used biomechanical simulation OpenSim. Here, the joints of the human body would be driven kinematically along prescribed exercise trajectories while the dynamics engine computed the response of the VIS degrees of freedom. However, mechanism designers and biomechanics experts have their own established tools, making it very desirable to have a way of decoupling the biomechanics from the VIS modeling, simulation, and analyses. We have derived a set of equations that rigorously accomplishes this goal, and have implemented them as an interface function that provides an alternative driving mechanism for an existing force-based VIS analysis simulation. When enabled, the simulated human/VIS system dynamics is now driven by this function, instead of the recorded force methodology described above. The function is designed to accept input from a data file containing the required time-stamped human motion and inertia terms corresponding to the specific exercise in question. This data file is generated by an OpenSim plugin written for that purpose. The existing VIS analytical simulation is developed using NASA’s Trick Simulation Environment [2], as well as its MBDyn multibody dynamics [3] package. The presentation will provide a detailed overview of the mathematical formulation, assumptions, plugin implementation, software interfaces, and results for a sample set of representative exercises. The results from this work aim to better inform VIS design efforts, as well as countermeasure device/protocol designs with respect to exercise type and frequency effects on vehicle structural and microgravity restrictions.

Countermeasures↗
NASA NTRS1992

A first-order global model of Late Cenozoic climatic change: Orbital forcing as a pacemaker of the ice ages

The development of a theory of the evolution of the climate of the earth over millions of years can be subdivided into three fundamental, nested, problems: (1) to establish by equilibrium climate models (e.g., general circulation models) the diagnostic relations, valid at any time, between the fast-response climate variables (i.e., the 'weather statistics') and both the prescribed external radiative forcing and the prescribed distribution of the slow response variables (e.g., the ice sheets and shelves, the deep ocean state, and the atmospheric CO2 concentration); (2) to construct, by an essentially inductive process, a model of the time-dependent evolution of the slow-response climatic variables over time scales longer than the damping times of these variables but shorter than the time scale of tectonic changes in the boundary conditions (e.g., altered geography and elevation of the continents, slow outgassing, and weathering) and ultra-slow astronomical changes such as in the solar radiative output; and (3) to determine the nature of these ultra-slow processes and their effects on the evolution of the equilibrium state of the climatic system about which the above time-dependent variations occur. All three problems are discussed in the context of the theory of the Quaternary climate, which will be incomplete unless it is embedded in a more general theory for the fuller Cenozoic that can accommodate the onset of the ice-age fluctuations. We construct a simple mathematical model for the Late Cenozoic climatic changes based on the hypothesis that forced and free variations of the concentration of atmospheric greenhouse gases (notably CO2), coupled with changes in the deep ocean state and ice mass, under the additional 'pacemaking' influence of earth-orbital forcing, are primary determinants of the climate state over this period. Our goal is to illustrate how a single model governing both very long term variations and higher frequency oscillatory variations in the Pleistocene can be formulated with relatively few adjustable parameters.

Saltzman, Barry↗
NASA NTRS1995

Decadal variability of the Tropical Atlantic Ocean Surface Temperature in shipboard measurements and in a Global Ocean-Atmosphere model

Sea surface temperature (SST) variability was investigated in a 200-yr integration of a global model of the coupled oceanic and atmospheric general circulations developed at the Geophysical Fluid Dynamics Laboratory (GFDL). The second 100 yr of SST in the coupled model's tropical Atlantic region were analyzed with a variety of techniques. Analyses of SST time series, averaged over approximately the same subregions as the Global Ocean Surface Temperature Atlas (GOSTA) time series, showed that the GFDL SST anomalies also undergo pronounced quasi-oscillatory decadal and multidecadal variability but at somewhat shorter timescales than the GOSTA SST anomalies. Further analyses of the horizontal structures of the decadal timescale variability in the GFDL coupled model showed the existence of two types of variability in general agreement with results of the GOSTA SST time series analyses. One type, characterized by timescales between 8 and 11 yr, has high spatial coherence within each hemisphere but not between the two hemispheres of the tropical Atlantic. A second type, characterized by timescales between 12 and 20 yr, has high spatial coherence between the two hemispheres. The second type of variability is considerably weaker than the first. As in the GOSTA time series, the multidecadal variability in the GFDL SST time series has approximately opposite phases between the tropical North and South Atlantic Oceans. Empirical orthogonal function analyses of the tropical Atlantic SST anomalies revealed a north-south bipolar pattern as the dominant pattern of decadal variability. It is suggested that the bipolar pattern can be interpreted as decadal variability of the interhemispheric gradient of SST anomalies. The decadal and multidecadal timescale variability of the tropical Atlantic SST, both in the actual and in the GFDL model, stands out significantly above the background 'red noise' and is coherent within each of the time series, suggesting that specific sets of processes may be responsible for the choice of the decadal and multidecadal timescales. Finally, it must be emphasized that the GFDL coupled ocean-atmosphere model generates the decadal and multidecadal timescale variability without any externally applied force, solar or lunar, at those timescales.

Mehta, Vikram M.↗
DOE OSTI2022

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Unsteady Extinction of Opposed Jet Ethylene/Methane HIFiRE Surrogate Fuel Mixtures vs Air

A unique idealized study of the subject fuel vs. air systems was conducted using an Oscillatory-input Opposed Jet Burner (OOJB) system and a newly refined analysis. Extensive dynamic-extinction measurements were obtained on unanchored (free-floating) laminar Counter Flow Diffusion Flames (CFDFs) at 1-atm, stabilized by steady input velocities (e.g., U(sub air)) and perturbed by superimposed in-phase sinusoidal velocity inputs at fuel and air nozzle exits. Ethylene (C2H4) and methane (CH4), and intermediate 64/36 and 15/85 molar percent mixtures were studied. The latter gaseous surrogates were chosen earlier to mimic ignition and respective steady Flame Strengths (FS = U(sub air)) of vaporized and cracked, and un-cracked, JP-7 "like" kerosene for a Hypersonic International Flight Research Experimentation (HIFiRE) scramjet. For steady idealized flameholding, the 100% C2H4 flame is respectively approx. 1.3 and approx.2.7 times stronger than a 64/36 mix and CH4; but is still 12.0 times weaker than a 100% H2-air flame. Limited Hot-Wire (HW) measurements of velocity oscillations at convergent-nozzle exits, and more extensive Probe Microphone (PM) measurements of acoustic pressures, were used to normalize Dynamic FSs, which decayed linearly with pk/pk U(sub air) (velocity magnitude, HW), and also pk/pk P (pressure magnitude, PM). Thus Dynamic Flame Weakening (DFW) is defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = -100 d(U(sub air)/U(sub air),0Hz)/d(pkpk P). Key findings are: (1) Ethylene flames are uniquely strong and resilient to extinction by oscillating inflows below 150 Hz; (2) Methane flames are uniquely weak; (3) Ethylene / methane surrogate flames are disproportionately strong with respect to ethylene content; and (4) Flame weakening is consistent with limited published results on forced unsteady CFDFs. Thus from 0 to approx. 10 Hz and slightly higher, lagging diffusive responses of key species led to progressive phase lags (relative to inputs) in the oscillating flames, and caused maximum weakening. At 20 to 150 Hz, diffusion-rate-limited effects diminished, causing flames to "regain strengnth," and eventually become completely insensitive beyond 300 Hz. Detailed mechanistic understanding is needed. Overall, ethylene flames are remarkably resilient to dynamic extinction by oscillating inflows. They are the strongest, with the notable exception of H2. For HIFiRE tests, the 64%/36% surrogate disproportionally retains the high dynamic FS of ethylene, so the potential for loss of scramjet flameholding (flameout) due to low frequency oscillations is significantly mitigated.

Vaden, Sarah N.↗