Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing front”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A new soil mixing layer model for simulating conservative solute loss from initially saturated soil to surface runoff

Surface water pollution due to solute loss from soil to surface runoff is a serious threat to the environment, and it is of great significance to predict the solute loss. While soil mixing layer theory has been widely used to predict the solute loss, it has been observed that the conventional soil mixing layer models with constant mixing layer depth underestimate solute concentration in surface runoff, and thus underestimate the solute loss after surface runoff starts. This study introduced a new soil mixing layer model with a new concept that the mixing layer depth increases downward with time after surface runoff occurs. This new concept was implemented by adding a mixing front between the soil mixing layer and its underlying soil layer, which is conceptually similar to the wetting front of soil infiltration given by the Green-Ampt equation. This is equivalent to treating the soil mixing layer as a deforming control volume. A mechanistic explanation for the new model was presented. This study also developed a closed form solution of the new model. The new model was evaluated by simulating five laboratory experiments conducted by two groups of researchers. The first two experiments have ponding water, and consider different drainage and runoff conditions at both early and late runoff times. The other three experiments do not have ponding water, and consider different infiltration conditions at the bottom of soil boxes. Simulation results of the five experiments show that the new model resolves the underestimation problem of conventional soil mixing layer models. The new model is subject to a total of six limits such as that the new model only considers conservative solute. Three limits are related to the assumptions used to derive the closed form solution, and they may be alleviated in a future research as briefly discussed at the end of the paper.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

Thermohaline circulation below the Ross Ice Shelf - A consequence of tidally induced vertical mixing and basal melting

The warmest water below parts of the Ross Ice Shelf resides in the lowest portion of the water column because of its high salinity. Vertical mixing caused by tidal stirring can thus induce ablation by lifting the warm but dense water into contact with the ice shelf. A numerical tidal simulation indicates that vertically well-mixed conditions predominate in the southeastern part of the sub-ice shelf cavity, where the water column thickness is small. Basal melting in this region is expected to be between 0.05 and 0.5 m/yr and will drive a thermohaline circulation having the following characteristics: high salinity shelf water (at - 1.8 C), formed by winter sea ice production in the open Ross Sea, flows along the seabed toward the tidal mixing fronts below the ice shelf; and meltwater (at -2.2 C), produced in the well-mixed region, flows out of the sub-ice shelf cavity along the ice shelf bottom. Sensitivity of this ablation process to climatic change is expected to be small because high salinity shelf water is constrained to have the sea surface freezing temperature.

Macayeal, D. R.↗
NASA NTRS2023

Anomalous DOC Signatures Reveal Iron Control on Export Dynamics in the Pacific Southern Ocean

Here we shed light on two mechanisms that stimulate deep particle export via upper-ocean iron fertilization in the Southern Ocean: deep frontal mixing and melting of sea ice. We present data collected a decade apart in the Pacific sector of the Southern Ocean when, serendipitously, seasonal Antarctic ice melt was anomalously low (2008) and anomalously high (2017). In 2008, the low ice melt year, we concluded that vertical mixing of iron into the euphotic zone via deep-mixing fronts was the primary stimulant of export that reached depths of ~1500 meters. This process was evidenced by localized enhancements of dissolved organic carbon (DOC) concentrations up to 4 μmol C kg-1 beneath seven branches of fronts embedded within the Antarctic Circumpolar Current (ACC). We used these enhanced DOC concentrations in the bathypelagic as primary indications of the depths and locations of recent export, as it is a logical residue of such. In 2017, the year in which sea ice melt was anomalously high, we concluded that the main driver of a widespread export event to the seafloor was the lateral influx of iron within the melt. Indications of this event included substantial enhancements of DOC concentrations (2 - 31 6 μmol C kg-1), elevated beam attenuation, and enhanced surface iron concentrations associated with a layer of low salinity water at a nearby station. Further, significant deficits of upper ocean silicic acid during the 2017 occupation indicated that deep export was likely stimulated by an iron-fueled diatom bloom. This analysis highlights the impact of iron supplied from frontal vertical mixing and sea ice melt on export and ultimately for long-term carbon sequestration in the Southern Ocean, as well as the utility of deep DOC enrichments as signatures of particle export. Understanding the impact that ice melt events have on carbon export is crucial given that anomalous events are occurring more often as our climate changes.

Dissolved Organic Carbon↗
DOE OSTI2022

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗
NASA NTRS1987

Mixing levels, the Apennine Front soil component, and compositional trends in the Apollo 15 soils

New compositional data (by INAA) for 29 samples of Apollo 15 soil fractions are presented along with data for six depth intervals of the bottom half of the double drive tube at station 2 on the Apennine Front. The compositional trends in the Apollo 15 soils are discussed in terms of four techniques used to determine the important components of the regolith: graphical techniques, factor analysis, chemical analysis of individual soil particles, and mixing models. The concept of 'mixing levels' is introduced, which states that a sample of lunar regolith can be considered as a mixture on several different compositional levels (I, igneous rock types; II, local rock types; III, soils), and it is shown that most attempts at modeling lunar soils as mixtures have been directed at level II. It is argued that only four components are required to account for the variation in concentrations of lithophile elements in the Apollo 15 samples: mare basalt, KREEP basalt, green glass, and Apennine front soil (best represented by the station 2 samples).

Korotev, Randy L.↗
DOE OSTI2021

Advanced Energy Approaches for Mixed-Use Developments in the Front Range Urban Corridor

The primary objective of the document is to present example technical approaches and strategies that planners and developers can reference for all-electric mixed-use developments in the Rocky Mountain Front Range urban corridor that are considering advanced energy performance goals such as high levels of energy efficiency and renewable integration. The strategies are tailored toward the vision for a zero energy, all-electric, district thermal system-connected starter district. In this example, on-site solar generation is emphasized as a strategy that can be combined with electrical and thermal energy storage for enhanced community resilience and demand flexibility. Opportunities for utilizing or adapting these strategies include: sharing with vertical developers for potential integration of strategies into vertical design guidelines; incorporation into Owner Project Requirements (OPRs); providing example strategies to design teams as a starting point for design; informing structure of agenda for a building design charrette; comparing building designs to example strategies at various points in the design process to evaluate design progression. While the technical strategies outlined in this document can help inform the development of guidelines for projects in the Front Range urban corridor, this document does not set requirements for any rating, certification, program, or policy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1986

The interaction of a surface cold front with a prefrontal thermodynamically well-mixed boundary layer

The dynamics of a surface cold front as it proceeds from sea to shore and interacts with a thermodynamically well-mixed boundary layer are investigated by means of numerical simulations using a modified version of the two-dimensional model of Reeder and Smith (1987). The theoretical basis of the model is explained, and the results are presented graphically and characterized in detail. Good agreement is found with observational data on Australian summer cold fronts, and the transverse frontal circulation is shown to be significantly affected by diabatic heating.

Reeder, Michael J.↗
NASA NTRS2000

Stereo Imaging Velocimetry of Mixing Driven by Buoyancy Induced Flow Fields

Mixing of two fluids generated by steady and particularly g-jitter acceleration is fundamental towards the understanding of transport phenomena in a microgravity environment. We propose to carry out flight and ground-based experiments to quantify flow fields due to g-jitter type of accelerations using Stereo Imaging Velocimetry (SIV), and measure the concentration field using laser fluorescence. The understanding of the effects of g-jitter on transport phenomena is of great practical interest to the microgravity community and impacts the design of experiments for the Space Shuttle as well as the International Space Station. The aim of our proposed research is to provide quantitative data to the community on the effects of g-jitter on flow fields due to mixing induced by buoyancy forces. The fundamental phenomenon of mixing occurs in a broad range of materials processing encompassing the growth of opto-electronic materials and semiconductors, (by directional freezing and physical vapor transport), to solution and protein crystal growth. In materials processing of these systems, crystal homogeneity, which is affected by the solutal field distribution, is one of the major issues. The understanding of fluid mixing driven by buoyancy forces, besides its importance as a topic in fundamental science, can contribute towards the understanding of how solutal fields behave under various body forces. The body forces of interest are steady acceleration and g-jitter acceleration as in a Space Shuttle environment or the International Space Station. Since control of the body force is important, the flight experiment will be carried out on a tunable microgravity vibration isolation mount, which will permit us to precisely input the desired forcing function to simulate a range of body forces. To that end, we propose to design a flight experiment that can only be carried out under microgravity conditions to fully exploit the effects of various body forces on fluid mixing. Recent flight experiments, by the P.I. through collaboration with the Canadian Space Agency (STS-85, August 1997), aimed at determining the stability of the interface between two miscible liquids inside an enclosure show that a long liquid column (5 cm) under microgravity isolation conditions can be stable, i.e. the interface remains sharp and vertical over a short time scale; thus transport occurs by molecular mass diffusion. On the other hand, when the two liquids were excited from a controlled vibration source (Microgravity Vibration Isolation Mount) two to four mode large amplitude quasi-stationary waves were observed. The data was limited to CCD recording of the dynamics of the interface between the two fluids. We propose to carry out flight experiments to quantify the dynamics of the flow field using Stereo Imaging Velocimetry and measure the concentration field using laser fluorescence. The results will serve as a basis to understand effects of g-jitter on transport phenomena, in this case mass diffusion. As the measurement of the kinematics of the flow field will shed light on the instability mechanism. The research will allow measurement of the flow field in microgravity environment to prove two hypotheses: (1) Maxwell's hypothesis: finite convection always exists in diffusing systems, and (2) Quasi-stationary waves inside a bounded enclosure in a microgravity environment is generated by Kelvin-Helmholtz instability; resonance of the interface which produces incipient mixing is due to Rayleigh-Taylor instability. The first hypothesis can be used as a benchmark experiment to illustrate diffusive mixing. The second hypothesis will lead to the understanding of g-jitter effects on buoyancy driven flow fields which occur in many situations involving materials processing, and other basic fluid physics phenomena. In addition, the second hypothesis will also provide insight in how Rayleigh-Taylor and Kelvin-Helmholtz instabilities propagate concentration fronts during mixing. Measurement of the flow field using SIV is important because it is the flow field which causes instability at the interface between the two fluids. Mixing driven by buoyancy induced flow fields will be addressed both experimentally and computationally. The experimental effort will address the kinematics of mixing: stretching, transport and chaos. Quantification of the mechanisms of mixing will consists of measuring the flow field using the SIV system at Glenn and capturing the dynamics of the interface, to measure mass transport, using a CCD camera. These experiments will be carried out within the framework of Earth's gravity and g-jitter microgravity acceleration as in a Space Shuttle environment or the International Space Station. The g-jitter will be induced and controlled using a tunable vibration isolation platform to isolate against vibration as well as input periodic and random vibration to the system. The parametric range of the microgravity experiment will be extended from the experiments on STS-85 to investigate higher mode quasi-stationary waves (8 to 12), as well as resonance regions which leads to chaos and turbulence. Ground-based experiments will focus on effects of vibration on stably stratified fluid layers in order to scale for possible scenarios in a microgravity environment. These vibrations will be subjected perpendicular to the concentration field on the ground since the parallel case can only be carried out in a microgravity environment. The concept of dynamical similarity will be applied to tune the experiments as closely as possible to a Space Shuttle environment or the International Space Station. The computational effort will take advantage of the Computational Laboratory at Glenn to corroborate the experimental findings with predictions of the dynamics of the flow field using the codes FLUENT (finite difference based) and FIDAP (finite element based). We will investigate two important cases, single-fluid model to address dilute systems with negligible jump in viscosity and the more general two-fluid model which accounts for finite jump in viscosity. Apart from its microgravity relevance, this experiment is well suited to study dynamics in nonlinear systems.

Duval, W. M. B.↗
DOE OSTI2020

Experimental study of vorticity-strain interactions in turbulent premixed counterflow flames

Here, the effects of heat release on interactions between vorticity (ω) and strain rate (s) in turbulent premixed CH 4 /O 2 /N 2 counterflow flames are investigated using simultaneous OH laser-induced fluorescence (LIF) and tomographic particle image velocimetry (TPIV) measurements. A comparison between the flames and a corresponding turbulent non-reacting variable density N 2 -vs-products counterflow reveals the impact of heat release on vorticity-strain rate alignment statistics. Vorticity and strain rate statistics in the flames and non-reacting flow are conditioned on distance from the local flame front and gas mixing layer interface (GMLI) contours, respectively. The magnitude, alignment, and spatial distribution of the vorticity and principal strain rates (s 1 , s 2 , s 3 ) are rather different when heat release is present. Density variations without heat release enhance the ω-s 2 alignment while significantly reducing the ω-s 3 alignment and modestly reducing the ω-s 1 alignment. In contrast, heat release at the flame front further reduces the ω-s 1 alignment but increases the ω-s 3 alignment and suppresses the preferential ω-s 2 alignment. Furthermore, increasing turbulence diminishes the effect of heat release on this preferential alignment. In regions with the largest vorticities, both the reacting and non-reacting counterflows show an increase in the probability of ω-s 2 alignment. All counterflow cases have a net positive vortex-stretching contribution to the enstrophy production with a peak production rate at the flame front or GMLI, but the peak values depend on the density variation, heat release, and turbulence level. Elucidation of the complex interplay between these factors contributes to the understanding of the dynamics of turbulence-flame interactions.

42 ENGINEERING↗
DOE OSTI2026

Identification and mitigation of memory block timing issue in ITk ABCStar during ASIC production

The ABCStar is a mixed-signal front-end readout ASIC for the strips sensor portion of the ATLAS ITk detector being developed as part of the High-Luminosity LHC upgrade. In pre-production testing, a subtle design flaw was uncovered in the ABCStar that was reducing wafer yields in some manufactured lots from the expected 90% to as low as 2%. The root cause was determined to be a timing issue in the logic synthesized to control previously silicon proven memory blocks re-used for this ASIC. The solutions proposed included manufacturing process changes by the wafer foundry, changes to the operating parameters for the ABCStar in the detector, and the possibility that a redesign might be required. The two mitigation efforts were undertaken in parallel, with the process modification route a less desirable solution since already manufactured wafers would need to be scrapped in favour of the new ones. Based on a knowledge of the existing process, and testing done on the worst performing wafers, it was proposed that raising the core operating voltage of the ABCStar from 1.20V to 1.25V could address the timing issue by sufficiently speeding up its transistors. An extensive testing program that included the effects of temperature and radiation expected over the lifetime of the ITk detector was conducted to validate that approach. Those tests and studies proved that even the worst performing wafers would have yields over 80% with the 1.25V core voltage, and neither the modified process nor redesign would be required for ensuring reliable operation of the ITk. Based on testing, a further timing mitigation was implemented to provide an additional margin of reliability by increasing the duty cycle of the clock to the ABCStar. Testing of all ABCStar wafers has been completed and the production of the detector modules using these ASICs is now well underway as a result of the efforts detailed herein.

FOS: Physical sciences↗
NASA NTRS1978

Fractional crystallization of iron meteorites, an experimental study

Measurements of the distribution coefficients of various elements (Ni, Co, P, Pt, Ir, Ge, Cr, and Au) between solid and liquid phases were made for FeNi alloys. A high thermal gradient crystal growing furnace was employed in order to promote plane front solidification, rapid mixing in the liquid and a minimum of diffusion in the solid. Eight measurements of the nickel distribution coefficient gave an average value of 0.91 plus or minus 0.03 and this value was independent of the Ni concentration in the range 5.5 to 8.75 wt. % Ni. Equilibrium distribution coefficients were also calculated by correcting for the diffusion boundary layer present in the liquid. These values are within 20% of the measured distribution coefficient values. Except for Co and Ni, all the measured distribution coefficient values are closer to 1.0 than those calculated by Scott (1972) for the IIIAB irons. Using the experimental values and the calculated equilibrium values for ternary additions to FeNi melts, it appears that fractional crystallization processes will not explain the observed chemical variations in the major iron meteorite groups. Alternative explanations are suggested which also recognize the fact that cooling rates vary widely in several chemical groups.

Goldstein, J. I.↗
NASA NTRS1982

Delineation of estuarine fronts in the German Bight using airborne laser-induced water Raman backscatter and fluorescence of water column constituents

The acquisition and application of airborne laser induced emission spectra from German Bight water during the 1979 MARSEN experiment is detailed for the synoptic location of estuarine fronts. The NASA Airborne Oceanographic Lidar (AOL) was operated in the fluorosensing mode. A nitrogen laser transmitter at 337.1 nm was used to stimulate the water column to obtain Gelbstoff or organic material fluorescence spectra together with water Raman backscatter. Maps showing the location and relative strength of estuarine fronts are presented. The distribution of the fronts indicates that mixing within the German Bight takes place across a relatively large area. Reasonable agreement between the patterns observed by the AOL and published results are obtained. The limitations and constraints of this technique are indicated and improvements to the AOL fluorosensor are discussed with respect to future ocean mapping applications.

Hoge, F. E.↗
DOE OSTI2023

A Cryogenic Readout IC with 100 KSPS in-Pixel ADC for Skipper CCD-in-CMOS Sensors

The Skipper CCD-in-CMOS Parallel Read-Out Circuit (SPROCKET) is a mixed-signal front-end design for the readout of Skipper CCD-in-CMOS image sensors. SPROCKET is fabricated in a 65 nm CMOS process and each pixel occupies a 50$\mu$m $\times$ 50$\mu$m footprint. SPROCKET is intended to be heterogeneously integrated with a Skipper-in-CMOS sensor array, such that one readout pixel is connected to a multiplexed array of nine Skipper-in-CMOS pixels to enable massively parallel readout. The front-end includes a variable gain preamplifier, a correlated double sampling circuit, and a 10-bit serial successive approximation register (SAR) ADC. The circuit achieves a sample rate of 100 ksps with 0.48 $\mathrm{e^-_{rms}}$ equivalent noise at the input to the ADC. SPROCKET achieves a maximum dynamic range of 9,000 $e^-$ at the lowest gain setting (or 900 $e^-$ at the lowest noise setting). The circuit operates at 100 Kelvin with a power consumption of 40 $\mu W$ per pixel. A SPROCKET test chip was submitted in September 2022, and test results will be presented at the conference.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Mixing Driven by Critical Reflection of Near-Inertial Waves over the Texas–Louisiana Shelf

Studies of internal wave-driven mixing in the coastal ocean have been focused on internal tides, while wind-driven near-inertial waves (NIWs) have received less attention in this regard. Here, this study demonstrates a scenario of NIW-driven mixing over the Texas–Louisiana shelf. Supported by a high-resolution simulation over the shelf, the NIWs driven by land–sea breeze radiate downward at a sharp front and enhance the mixing in the bottom boundary layer where the NIWs are focused because of slantwise critical reflection. The criterion for slantwise critical reflection of NIWs is ω=$\sqrt{f{^2}+s^{2}_{bot}N{^2}(1-s_{ρ}/s_{bot})}$ (where ω is the wave frequency, s bot is the bottom slope, and s ρ is the isopycnal slope), under the assumption that the mean flow is in a thermal wind balance and only varies in the slope-normal direction. The mechanism driving the enhanced mixing is explored in an idealized simulation. During slantwise critical reflection, NIWs are amplified with enhanced shear and periodically destratify a bottom boundary layer via differential buoyancy advection, leading to periodically enhanced mixing. Turbulent transport of tracers is also enhanced during slantwise critical reflection of NIWs, which has implications for bottom hypoxia over the Texas–Louisiana shelf.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Numerical Mixing Suppresses Submesoscale Baroclinic Instabilities Over Sloping Bathymetry

Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Analysis of Composite Skin-Stiffener Debond Specimens Using Volume Elements and a Shell/3D Modeling Technique

The debonding of a skin/stringer specimen subjected to tension was studied using three-dimensional volume element modeling and computational fracture mechanics. Mixed mode strain energy release rates were calculated from finite element results using the virtual crack closure technique. The simulations revealed an increase in total energy release rate in the immediate vicinity of the free edges of the specimen. Correlation of the computed mixed-mode strain energy release rates along the delamination front contour with a two-dimensional mixed-mode interlaminar fracture criterion suggested that in spite of peak total energy release rates at the free edge the delamination would not advance at the edges first. The qualitative prediction of the shape of the delamination front was confirmed by X-ray photographs of a specimen taken during testing. The good correlation between prediction based on analysis and experiment demonstrated the efficiency of a mixed-mode failure analysis for the investigation of skin/stiffener separation due to delamination in the adherents. The application of a shell/3D modeling technique for the simulation of skin/stringer debond in a specimen subjected to three-point bending is also demonstrated. The global structure was modeled with shell elements. A local three-dimensional model, extending to about three specimen thicknesses on either side of the delamination front was used to capture the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from shell/3D simulations were in good agreement with results obtained from full solid models. The good correlations of the results demonstrated the effectiveness of the shell/3D modeling technique for the investigation of skin/stiffener separation due to delamination in the adherents.

Krueger, Ronald↗