Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave generator”

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 559 records · Page 31

Effects of Forward- and Backward-Facing Steps on the Crossflow Receptivity and Stability in Supersonic Boundary Layers

The effects of forward- and backward-facing steps on the receptivity and stability of three-dimensional supersonic boundary layers over a swept wing with a blunt leading edge are numerically investigated for a freestream Mach number of 3 and a sweep angle of 30 degrees. The flow fields are obtained by solving the full Navier-Stokes equations. The evolution of instability waves generated by surface roughness is simulated with and without the forward- and backward-facing steps. The separation bubble lengths are about 5-10 step heights for the forward-facing step and are about 10 for the backward-facing step. The linear stability calculations show very strong instability in the separated region with a large frequency domain. The simulation results show that the presence of backward-facing steps decreases the amplitude of the stationary crossflow vortices with longer spanwise wavelengths by about fifty percent and the presence of forward-facing steps does not modify the amplitudes noticeably across the steps. The waves with the shorter wavelengths grow substantially downstream of the step in agreement with the linear stability prediction.

Balakumar, P.↗

Explosions of infalling comets in Jupiter's atmosphere

In view of the expected collision of comet Shoemaker-Levy 9 (1993e) with Jupiter in 1994 July, we calculate basic properties of the initial interaction for a simplified Jovian atmosphere. The comet is expected to impact Jupiter at 60 km/sec and at an angle of 45 deg to the zenith. The shock wave generated by the bolide should be optically thick once it has penetrated to an atmospheric density approximately 10(exp -6) gr/cm(exp 3), and we calculate the post-shock conditions assuming local thermodynamic equilibrium (LTE) for shock velocities v(sub sh) in the range 10 to 60 km/sec and preshock densities rho(sub a) = 10(exp -6) to 10(exp -2) gr/cm(exp 3). Our shock calculations include molecular hydrogen, atomic hydrogen, ionized hydrogen, neutral helium, and singly ionized helium. Even at the highest shock velocity, the gas is only partially ionized and the postshock temperature rises with preshock density in order to maintain the ionization. The value of the effective shock adiabatic index gamma(sub sh) varies from 1.17 (at low v(sub sh) and rho(sub a) to 1.40 (at high v(sub sh) and rho(sub a). The ablation rate is limited by the radiative flux that reaches the bolide surface. We argue that the ablated gas does not efficiently transfer its kinetic energy to the atmosphere, and it ultimately slows in a similar fashion to the comet material. As the bolide initially falls through the atmosphere, the character of the shock emission changes. At rho(sub a) approximately 10(exp -8) gr/cm (exp 3), the gas is optically thin and we expect line emission; in the optical spectrum, Balmer emission is expected from the shocked atmosphere and low-ionization metal lines from ablated cometary material. At rho(sub a) approximately 10(exp -6) gr/m(exp 3), the shocked gas is optically thick and the shock front near the bolide produces a blackbody spectrum. The temperature is favorable for ultraviolet (1000 to 3000 A) emission and the luminosity may be approximately 5 x 10(exp 23) ergs/sec for approximately 0.6 sec for a bolide 1 km in radius. At rho(sub a) approximately 10(exp -4) gr/cm(exp 3), the bolide has passed below the ultraviolet photosphere. The shock front emits considerable ionizing radiation, but it is absorbed in a narrow preshock region. The bolometric correction for the optical luminosity is large and we expect a 3000 to 8000 A luminosity of approximately 3 x 10(exp 23) ergs/sec for approximately 1 sec. The optical emission is strongly peaked in the vicinity of the bolide. The bolide does have a somewhat less luminous, optically thick trail extending greater than or equal to 10 km, but the radiation is characterized by a temperature of 4000 to 5000 K. From the fragmentation model of Chyba, Thomas, & Zahnle (1993), the bolide deposits most of its kinetic energy at rho(sub a) approximately 10(exp -3) gr/cm(exp 3) and this is the effective explosion site. The shock wave from such an explosion can move up about one density scale height. We examine the breakout of the shock front from the Jovian atmosphere and find that the shock acceleration in the decreasing density region is slow, so that the energy flux in the shock front is small. Higher velocities might be generated by shock acceleration along the channel left by the bolide if the shock motion can occur before the channel closes off as a result of radiative cooling. Hot gas created by the explosion ultimately rises due to buoyancy on a timescale of a minute. The luminosity is highest when the bubble first rises into the optically thin part of the atmosphere and may be approximately 1 x 10(exp 25) ergs/sec in the near-infrared. Roughly 1% of the initial bolide energy may be radiated in this way; the rest of the energy is lost to sound waves from the initial explosion and to work done by the bubble on the surrounding atmosphere.

Chevalier, Roger A.↗

A nonlinear steady model for moist hydrostatic mountain waves

The dynamics of hydrostatic gravity waves generated by the passage of a steady, stably stratified, moist flow over a two-dimensional topography is considered. Coriolis effects are neglected. The cloud region is determined by the dynamics, and within that region the Brunt-Vaisala frequency takes on a value smaller than the outside value. In both the dry and cloudy regions the Brunt-Vaisala frequency is constant with height. The moist layer is considered to be either next to the mountain or at midlevels and to be deep enough so that an entire cloud forms in that layer. The nonlinearity in the flow and lower boundary affects the dynamics of these waves and wave drag. The latter is found to depend upon: (1) the location of the moist layer with respect to the ground, (2) the amount of moisture, (3) the degree of nonlinearity and (4) the departure from symmetry in the bottom topography.

Barcilon, A.↗

Investigation of a hypersonic crossing shock wave/turbulent boundary layer interaction

A combined theoretical and experimental study is presented for the interaction between crossing shock waves generated by (10 deg, 10 deg) sharp fins and a flat plate turbulent boundary layer at Mach 8.3. The theoretical model is the full 3D mean compressible Reynolds-averaged Navier-Stokes (RANS) equations incorporating the algebraic turbulent eddy viscosity model of Baldwin and Lomax (1978). A grid refinement study indicated that adequate resolution of the flow field has been achieved. Computed results agree well with experiment for surface pressure and surface flow patterns and for pitot pressure and yaw angle profiles in the flow field. The computations, however, significantly overpredict surface heat transfer. Analysis of the computed flow field results indicates the formation of complex streamline and wave structures within the interaction region.

Narayanswami, N↗

The Stability of Radiatively Cooling Jets. 2: Nonlinear Evolution

We use two-dimensional time-dependent hydrodynamical simulations to follow the growth of the Kelvin-Helmholtz (K-H) instability in cooling jets into the nonlinear regime. We focus primarily on asymmetric modes that give rise to transverse displacements of the jet beam. A variety of Mach numbers and two different cooling curves are studied. The growth rates of waves in the linear regime measured from the numerical simulations are in excellent agreement with the predictions of the linear stability analysis presented in the first paper in this series. In the nonlinear regime, the simulations show that asymmetric modes of the K-H instability can affect the structure and evolution of cooling jets in a number of ways. We find that jets in which the growth rate of the sinusoidal surface wave has a maximum at a so-called resonant frequency can be dominated by large-amplitude sinusoidal oscillations near this frequency. Eventually, growth of this wave can disrupt the jet. On the other hand, nonlinear body waves tend to produce low-amplitude wiggles in the shape of the jet but can result in strong shocks in the jet beam. In cooling jets, these shocks can produce dense knots and filaments of cooling gas within the jet. Ripples in the surface of the jet beam caused by both surface and body waves generate oblique shock "spurs" driven into the ambient gas. Our simulations show these shock "spurs" can accelerate ambient gas at large distances from the jet beam to low velocities, which represents a new mechanism by which low-velocity bipolar outflows may be driven by high-velocity jets. Rapid entrainment and acceleration of ambient gas may also occur if the jet is disrupted. For parameters typical of protostellar jets, the frequency at which K-H growth is a maximum (or highest frequency to which the entire jet can respond dynamically) will be associated with perturbations with a period of - 200 yr. Higher frequency (shorter period) perturbations excite waves associated with body modes that produce internal shocks and only small-amplitude wiggles within the jet. The fact that most observed systems show no evidence for large-amplitude sinusoidal oscillation leading to disruption is indicative that the perturbation frequencies are generally large, consistent with the suggestion that pro- tostellar jets arise from the inner regions (r less than 1 AU) of accretion disks.

Stone, James M.↗

Infrasound Detections of Low-Magnitude Earthquakes: Preliminary Results of the West Texas Acoustic Experiment

Infrasound observations have grown increasingly important for the monitoring of earthquakes. While large earthquakes generate infrasound that can be detected thousands of kilometers away, there are few near-field observations of infrasound generated by low-magnitude events. We describe preliminary results of the West Texas Acoustic Experiment, during which infrasound sensors collected continuous data in the Permian Basin for a six-month period spanning January—June 2023. During this time, more than 1000 earthquakes with magnitudes between 1.2 and 4.2 occurred within 50 km of the network. We used spectral analysis, array processing, and manual inspection of waveforms to evaluate arrivals of infrasound signals following 84 events with magnitudes between 2.5 and 4.2. Here, we describe eight such events and the infrasound signals associated with each. We find detections of seismic-to-acoustic infrasound signals associated with seven events. We also find strong evidence of a laterally-propagating, purely acoustic wave generated by an M2.9 earthquake.

58 GEOSCIENCES↗

Proximal impact deposits at the Cretaceous-Tertiary boundary in the Gulf of Mexico: a restudy of DSDP Leg 77 Sites 536 and 540

Restudy of Deep Sea Drilling Project Sites 536 and 540 in the southeast Gulf of Mexico gives evidence for a giant wave at Cretaceous-Tertiary boundary time. Five units are recognized: (1) Cenomanian limestone underlies a hiatus in which the five highest Cretaceous stages are missing, possibly because of catastrophic K-T erosion. (2) Pebbly mudstone, 45 m thick, represents a submarine landslide possibly of K-T age. (3) Current-bedded sandstone, more than 2.5 m thick, contains anomalous iridium, tektite glass, and shocked quartz; it is interpreted as ejecta from a nearby impact crater, reworked on the deep-sea floor by the resulting tsunami. (4) A 50-cm interval of calcareous mudstone containing small Cretaceous planktic foraminifera and the Ir peak is interpreted as the silt-size fraction of the Cretaceous material suspended by the impact-generated wave. (5) Calcareous mudstone with basal Tertiary forams and the uppermost tail of the Ir anomaly overlies the disturbed interval, dating the impact and wave event as K-T boundary age. Like Beloc in Haiti and Mimbral in Mexico, Sites 536 and 540 are consistent with a large K-T age impact at the nearby Chicxulub crater.

NASA Discipline Exobiology↗

The Impact of the Lightning Surge on SiC-based Medium-voltage Three-phase Four-wire Grid-connected Converters

The surge voltage induced by lightning will generate an inrush current in a grid-connected converter. It may not be large in Si-based converters because of the large grid side inductive filter. However, for SiC-based multi-level converters, the inrush current could be very large due to the small filter inductance, which is selected based on the current harmonic consideration. Most previous work focuses on system-level lightning protection while the converter-level impact, especially the internal components of the converter, is rarely discussed. This paper presents the analysis of the impact of the lightning surge on a 10 kV SiC MOSFET-based 13.8 kV three-phase four-wire grid-connected converter. The IEEE arrester model and the combined wave generator model in IEC 61000-4-5 standard are used to evaluate the inrush current, DC-link overvoltage, and component-to-ground potential of the converter. Analysis and simulation results have been provided to demonstrate the impact of lightning on the grid-connected converter. It is found that large inrush current occurs due to the small filter inductance, and the neutral-to-ground impedance impacts both the inrush current and the neutral point-to-ground potential.

Li, Haiguo↗

Electrodynamic interactions

The electrodynamic interaction of long metallic tethers with the ionosphere is introduced. Due to the tether's motion across the Earth's field, a polarization electric field is generated along the tether. In a system like the Tethered Satellite System (TSS) where the tether is coated with dielectric and the electrical contact is between its two terminations and the ionosphere, the upper end termination will collect electrons and, in a passive system, the lower end termination will collect ions. Space charge regions will develop around the two terminations of the system, which can be at considerable potential with respect to the unperturbed ionospheric plasma and the current in the tether will be most significantly determined by local processes in such space charge regions. A qualitative view of the global perturbation induced by TSS in the ionosphere is outlined. One of the applications of electrodynamic tethers is that of using them as low frequency wave generators to communicate to the Earth. This application is briefly discussed. An equivalent circuit of the tether system is described. The basic equation of the circuit is given.

Dobrowolny, M.↗

Discovery of a nearby radio relic in the low-mass, merging cluster Abell 4067

Shock waves generated during cluster mergers offer a powerful probe of how large-scale structure grows and evolves in the Universe. As part of the MeerKAT-South Pole Telescope (SPT) survey, we report the discovery of a single arc-like radio relic in the galaxy cluster Abell 4067 ($z=0.099$), one of the lowest-mass clusters known to host such a structure. MeerKAT UHF-band (0.58–1.09 GHz) observations reveal a relic with a largest linear size of $\sim 1.48\pm 0.02$ Mpc, located at a projected distance of 0.95 Mpc from the cluster centre. XMM–Newton X-ray data show that the relic’s position and orientation relative to the intracluster medium (ICM) elongation are consistent with a merger-driven shock-wave scenario. The relic has an estimated radio power of $3.10\pm 0.03\times 10^{24}$ W Hz$^{-1}$ at 150 MHz. When placed in the $P_{150\, \mathrm{MHz}}$–$M_{500}$ scaling relation, the Abell 4067 relic appears less luminous compared to relics in more massive clusters, suggesting an association with weak merger shocks. This finding supports the idea that relics in low-mass clusters may form through less energetic merger events, leading to weak merger shocks. The latter is supported by the absence of a detectable central radio halo in Abell 4067, which reinforces the idea that luminous radio haloes are not a universal outcome of cluster mergers and highlights the role of cluster mass, merger energetics, and evolutionary stage in shaping diffuse radio emission in the ICM.

79 ASTRONOMY AND ASTROPHYSICS↗

Dependence of wind-farm-induced gravity waves and wind farm performance on non-dimensional atmospheric parameters and simulation configuration

This large-eddy simulation (LES) study examines how wind-farm-induced atmospheric gravity waves (AGWs) and wind farm performance depend on non-dimensional atmospheric parameters and simulation configuration. A hypothetical aligned wind farm of actuator disks is simulated under neutral surface conditions, with a stable capping inversion and a mildly stable free atmosphere, to assess the effects of stratification beyond the atmospheric boundary layer (ABL) on ABL flow. Simulation set-ups fully resolving AGWs are validated to minimize spurious wave generation and reflection from the domain boundaries. The validated set-up is then used to analyze AGW types and characteristics, as well as stratification impacts under conventionally neutral boundary layer (CNBL) conditions. These conditions are governed by four non-dimensional parameters: the Froude numbers of the free atmosphere and capping inversion (Fr, Fr i ), and the aspect ratios of the ABL and wind farm (H̃ i , S h ). Simulation configurations that fully resolve AGWs – capturing at least one wavelength both horizontally and vertically – yield the most realistic stratification effects on ABL flow, whereas partial or unresolved configurations produce non-physical, channel-like behavior. A coherent description of the AGW phenomena is provided, highlighting the central role of capping inversion displacement in linking ABL fluctuations with AGWs. Trapped waves are confined within the capping inversion, while interfacial and internal waves aloft are identified as the AGW types most relevant to wind farm performance. The wavy inversion, analogous to an interfacial wave, forms converging and diverging zones that drive power fluctuations across the farm. The interfacial wavelength, measured over the wind farm, corresponds to one diverging, one converging, and one mildly diverging zone. As the interfacial wavelength decreases with Fr i , multiple convergence–divergence zones develop under sub-critical conditions (Fr i <1.0), while for super-critical conditions (Fr i > 1.0), the wavelength approaches the farm length. Wave amplitude increases with decreasing H̃ i (i.e., shallower capping inversions). Wind farm performance is most sensitive to H̃ i : shallow boundary layers increase blockage and reduce efficiency, while deeper layers enhance efficiency. Increasing Fr and Fr i mitigates blockage, and increasing S h mainly improves wake recovery. Although local power fluctuations arise from AGWs, overall wind farm efficiency remains nearly constant with Fr and Fr i , improving primarily with larger H̃ i and S h .

17 WIND ENERGY↗

Stochastic particle acceleration in solar flares

It is proposed that particles during the second phase of solar flares are accelerated by stochastic resonant scattering off hydromagnetic waves and first order Fermi acceleration in shock waves generated in the impulsive phase of the flare. Solutions allow arbitrary power law momentum dependences of the momentum diffusion coefficient as well as the momentum diffusion coefficient as well as the momentum loss time. The acceleration time scale to a characteristic energy approximately 100 keV for protons can be as short as 5s. The resulting electron spectra show a characteristic double power law with a transition around 200 keV and are correlated to the proton spectra evaluated under equal boundary conditions, indicating that electrons and protons are accelerated by the same mechanism. The correlation between the different spectral indices in the electron double power law and between electron and proton spectra are governed by the ratio of first to second order acceleration and therefore allow a determination of the Alfven Mach number of the shock wave.

Droege, W.↗

The Space Technology 7 disturbance reduction system

The NASA New Millennium Program Space Technology 7 (ST7) project was selected for the purpose of flight validation of technology for minimizing disturbances on spacecraft and freely-floating test masses. Measurements of the motion of test masses are used for determining the gravitational field in its location, providing information on the nearby mass distribution or looking gravitational waves generated by distant star systems. Any unmeasured forces acting on the test mass cause a deviation from a purely gravitational trajectory and so limit the accuracy of the gravity field estimation. The ST7 project will validate the capability to eliminate unwanted forces to a level 3000 times lower than the current state of the art. The ST7 design is based on the concept of a freely-floating test mass surrounded by a spacecraft which shields the test mass from unwanted disturbances. since some disturbance of the test mass can be caused by motion of the spacecraft itself with respect to the test mass, the disturbance reduction system includes the capability of precisely measuring the position of the spacecraft with respect to the test mass and control of the spacecraft position to within 1/100 of a wavelength of light.. The overall system performance of the ST7 project will be infused into space observatories for gravitational waves such as the Laser Interferometer Space Antenna project, and into future mission for mapping the time-variable Earth gravity field. The ST7 precision control capability also has applications for separated-spacecraft interferometric imaging projects. ST7 will consist of an instrument package and a set of microthrusters, which will be attached to the European Space Agency’s LISA Pathfinder spacecraft with launch scheduled for 2008.

Folkner, W. M.↗

A Curated Dataset of Regional Meteor Events with Simultaneous Optical and Infrasound Observations (2006–2011)

We present a curated, openly accessible dataset of 71 regional meteor events simultaneously recorded by optical and infrasound instrumentation between 2006 and 2011. These events were captured during an observational campaign using the all-sky cameras of the Southern Ontario Meteor Network and the co-located Elginfield Infrasound Array. Each entry provides optical trajectory measurements, infrasound waveforms, and atmospheric specification profiles. The integration of optical and acoustic data enables robust linkage between observed acoustic signals and specific points along meteor trajectories, offering new opportunities to examine shock wave generation, propagation, and energy deposition processes. This release fills a critical observational gap by providing the first validated, openly accessible archive of simultaneous optical–infrasound meteor observations that supports trajectory reconstruction, acoustic propagation modeling, and energy deposition analyses. By making these data openly available in a structured format, this work establishes a durable reference resource that advances reproducibility, fosters cross-disciplinary research, and underpins future developments in meteor physics, atmospheric acoustics, and planetary defense.

astrometry↗

In-Situ Detection of Process-Induced Porosity During Cure of Out-of-Autoclave Composites

Composite materials offer unique benefits in aerospace applications, including high strength-to-weight ratio, and are becoming increasingly used by industry manufacturers. Current manufacturing and processing methods can lead to defects in the composites, which are identified after fabrication using inspection methods. This study utilized a high-temperature ultrasonic inspection system to detect process-induced porosity in an out-of-autoclave (OOA) composite panel during cure. In order to perform ultrasonic scans in-situ during cure, the system needed to be able to operate at oven temperatures up to 200 °C. There are no commercially available ultrasonic scanning systems that can operate continuously at these temperatures, so an enclosure was fabricated to insulate the ultrasonic system in the curing oven. The enclosure housed a MISTRAS® motorized X-Y raster scanner and contact transducer (Olympus® X2002). The Olympus® X2002 contact transducer is a continuous high-temperature delay line transducer with a center frequency of 2.25 MHz. In order to prevent the inside of the enclosure from reaching high temperatures, a sparge pipe inside the cooling enclosure was connected to a liquid nitrogen (LN2) tank outside of the oven via insulated hoses. The cooling system was automatically controlled by a temperature controller that was programmed to open a solenoid valve along the hose at 31 ℃ and close at 28 ℃ inside the cooling enclosure. When the solenoid valve was open, LN2 would flow from the LN2 tank and vaporize prior to reaching cooling enclosure. A U-shaped sparge pipe had holes drilled along its length to uniformly introduce cold nitrogen gas into the cooling enclosure. This prevented the motors from ever reaching the maximum desired operating temperature of 38 ℃. In the bottom of the enclosure, the tool plate was placed with the composite facing the outside of the enclosure, exposing it to the temperatures of the oven. The tool plate was 63.5 mm thick borosilicate glass, which was chosen for its low ultrasonic attenuation and moderate thermal conductivity. A thin layer of ultrasonic couplant was placed on the tool plate. The transducer sent and received ultrasonic waves through the tool plate into the composite. A MISTRAS® Remote UT System performed all wave generation, analog to digital conversion, and data acquisition and processing. The transducer was connected directly to the Remote UT System via a high-temperature BNC cable (CD International, RG-316). Polyimide insulated copper wire was used to connect the scanner to a 2-axis motor drive, which was connected to the Remote UT System. The motor drive and Remote UT System were located outside the industrial oven. A diagram of the ultrasonic system is shown in Figure 1. Figure 1: Diagram of ultrasonic in-situ defect detection system. In this study, a porosity gradient was introduced into the carbon fiber reinforced polymer (CFRP) composite through a misfit of the part and caul plate. A 24-ply quasi-isotropic ([0/90, ±45]12S) composite panel was laid up using Newport® AS4C-M/NB321 out-of-autoclave (OOA) plain weave (PW) prepreg (42% resin content, 195 g/m2). The first 16 plies were 203 mm × 203 mm. Ply drops were incorporated into the panel by decreasing the length of plies 17-24. The sizes of each ply are included in Table 1. All plies were centered along their length (x-axis) creating a trapezoidal type cross-section. Table 1: Ply sizes Plies Length (mm) Width (mm) 1-16 203 203 17-18 102 203 19-20 85 203 21-22 68 203 23-24 51 203 A 15-5 stainless steel caul plate (Size: 203 × 203 × 1.2 mm) was placed on the laminate. The composite and caul plate were vacuum bagged to the tool plate. Per the product datasheet, Newport 321 can be cured between 121 °C and 149 °C, depending on service temperature, with a hold for 90-120 minutes. Based on prior thermal testing of the cooling enclosure, the desired composite temperature cycle could be attained when the air temperature of the oven was ramped to 160 °C at a rate of 1.7 °C/min with a 3 hour hold. The extended hold was not necessary to cure the composite, but did allow for additional scans at the cure temperature during preliminary testing. The composite panel was under vacuum pressure for the entire cure cycle. Under vacuum, the caul plate deformed causing high and low pressure regions. The low-pressure regions formed high porosity in the composite part. The ultrasonic data collected showed that the inspection system was able to monitor the evolution of porosity within the composite throughout the cure cycle. Fifty scans of a 203 mm × 51 mm (8 in. × 2 in.) area of the composite were completed throughout the cure cycle. The scan speed was set to 43 mm/s in the x-direction and 36 mm/s in the y-direction, with a 1.0 mm/pixel resolution in both directions. Each scan took approximately five minutes to complete. Figure 2 shows C-scan amplitude data (top view) from one scan during the cure cycle. The color map indicates the maximum amplitude (%) measured with a time gate that contains reflections from the back surface (furthest from the transducer in contact with the caul plate) of the composite. Because porosity in the composite increases the ultrasonic attenuation in the composite part, a higher amplitude (red) indicates low porosity and a lower amplitude (blue) indicates high porosity. The regions where the caul plate was applying increased pressure (center and ends of the panel) resulted in reduced porosity, whereas the ply drop regions had reduced pressure and thus had higher porosity. Figure 2: Amplitude C-scan (top view) of OOA composite from one scan during the 160 °C temperature hold. Higher amplitude (red) indicates low porosity. Lower amplitude (blue) indicates high porosity. A B-scan is a two-dimensional image plotting the ultrasonic signals from a single y-position at every x-location of the composite panel. The travel time (μs) of the ultrasonic wave (representing the through-thickness location) is plotted along the y-axis, and the ultrasonic wave signal is represented as a contour map. Based on the B-scan data, the location of the defects with respect to the composite thickness was able to be determined. The ply drop regions had high porosity near the front (tool plate) surface of the composite (Figure 3). Figure 3: B-scan (cross-section view) of OOA composite from one scan during the 160 °C temperature hold. Using the data obtained from amplitude C-scans and B-scans throughout the cure cycle, the high and low porosity regions within the composite laminate were detected and localized with high spatial resolution. This paper will discuss the evolution of porosity content in the OOA composite during cure. OOA tests aided the transition of this system to an autoclave, which is the primary method of curing aerospace-grade thermoset composites.

Tyler B. Hudson↗

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗

The role of cool-flame fluctuations in high-pressure spray flames, studied using high-speed optical diagnostics and Large-Eddy Simulations

Here this work investigates the low- and high-temperature ignition and combustion processes, applied to the Engine Combustion Network Spray A flame, combining advanced optical diagnostics and large-eddy simulations (LES). Simultaneous high-speed (50 kHz) formaldehyde (CH 2 0) planar laser-induced fluorescence (PLIF) and line-of-sight OH* chemiluminescence imaging were used to measure the low- and high-temperature flame, during ignition as well as during quasi-steady combustion. While tracking the cool flame at the laser sheet plane, the present experimental setup allows detection of distinct ignition spots and dynamic fluctuations of the lift-off length over time, which overcomes limitations for flame tracking when using schlieren imaging. After significant development to improve LES prediction of the low-and high-temperature flame position, both during the ignition processes and quasi-steady combustion, the simulations were analyzed to gain understanding of the mixture variance and how this variance affects formation/consumption of CH 2 0. Analysis of the high-temperature ignition period shows that a key improvement in the LES is the ability to predict heterogeneous ignition sites, not only in the head of the jet, but in shear layers at the jet edge close to the position where flame lift-off eventually stabilizes. The LES analysis also shows concentrated pockets of CH 2 0, in the center of jet and at 20 mm downstream of the injector (in regions where the equivalence ratio is greater than 6), that are of similar length scale and frequency as the experiment (approximately 5–6 kHz). The periodic oscillation of CH 2 0 match the frequency of pressure waves generated during auto-ignition and reflected within the constant-volume vessel throughout injection. The ability of LES to capture the periodic appearance and destruction of CH 2 0 is particularly important because these structures travel downstream and become rich premixed flames that affect soot production.

42 ENGINEERING↗