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 73 records · Page 4

Direct numerical simulation of a high-pressure hydrogen micromix combustor: Flame structure and stabilisation mechanism

A high-pressure hydrogen micromix combustor has been investigated using direct numerical simulation with detailed chemistry to examine the flame structure and stabilisation mechanism. The configuration of the combustor was based on the design by Schefer et al., using numerical periodicity to mimic a large square array. A precursor simulation of an opposed jet-in-crossflow was first conducted to generate appropriate partially-premixed inflow boundary conditions for the subsequent reacting simulation. The resulting flame can be described as an predominantly-lean inhomogeneously-premixed lifted jet flame. Five main zones were identified: a jet mixing region, a core flame, a peripheral flame, a recirculation zone, and combustion products. The core flame, situated over the jet mixing region, was found to burn as a thin reaction front, responsible for over 85% of the total fuel consumption. The peripheral flame shrouded the core flame, had low mean flow with high turbulence, and burned at very lean conditions (in the distributed burning regime). It was shown that turbulent premixed flame propagation was an order-of-magnitude too slow to stabilise the flame at these conditions. Stabilisation was identified to be due to ignition events resulting from turbulent mixing of fuel from the jet into mean recirculation of very lean hot products. Ignition events were found to correlate with shear-driven Kelvin-Helmholtz vortices, and increased in likelihood with streamwise distance. At the flame base, isolated events were observed, which developed into rapidly burning flame kernels that were blown downstream. Further downstream, near-simultaneous spatially-distributed ignition events were observed, which appeared more like ignition sheets. The paper concludes with a broader discussion that considers generalising from the conditions considered here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intentional Uranium Tagging for Material Provenance and Pathway Forensics (LA19-Intentional-Forensics-NDD3Bb) (Annual Report for FY21)

This report describes the outcomes of the third and final year of a project to research the feasibility of tagging uranium materials, especially nuclear fuels. The experimental focus remained on metallic uranium forms under prospective surface and bulk tagging scenarios. Overall, the results showed that the tags could be successfully imparted and characterized in both as-built and degraded conditions. This was in line with expectations coming into this project, indicating promise for both surface and bulk tagging of metallic forms of uranium. Multiple surface tagging techniques and detection strategies were explored in FY21, with an emphasis on improving tag quality, readability, and detection in the field. Non-radioactive materials were used as a testbed. Selective deposition via laser beam was determined to be successful in imparting a readable titanium deposit on a stainless steel base plate, and can be read with high resolution characterization techniques (e.g., scanning electron microscopy) and field capable tools (e.g., eddy current testing). Other deposition techniques, such as selective deposition via electron beam and photoluminescent tags, were explored in FY21, and while success for these techniques would be dependent on additional work, these techniques showed potential for surface tagging applications. To survey bulk taggant elements for bulk uranium metal, 16 tagging elements were spread among 18 depleted uranium castings (4 baseline, 3 mix, 1 dilution, and 10 recycle). Most of these were made and characterized in FY21. Taggant acceptability was based upon manufacturability, detectability, and persistence from the standpoint of two detection options: bulk chemical analysis (for “chemical taggants”) and microstructural analysis (for “second phase taggants”). Taggant detection in both up-front manufacturing and in the face of “degradations” such as dilution, mixing, and recycling was generally good. Two independent laboratories carried out chemical analysis on most of the castings, and often at several locations within a casting, and the results are discussed. Scanning electron microscopy+EDS microanalysis revealed the second phases mostly contained the expected tagging elements. The shapes and 2 spatial distributions of these micron-sized carbides, oxides, and intermetallic particles offers opportunities for further science-based investigation and tagging optimization. Overall, V and Co currently appear as the best choices for chemical taggants while Al, Ti, Mn, Co, Pd, and Tb all look good as second phase taggants. The other elements considered here – Sc, Ni, Ge, Nb, Ce, Ta, W, Ir, and Au – while not being ruled out, require more study to become viable options. It is of special note that in the recycling study only one of the 12 elements fell out of detection even after 10 meltings, demonstrating their persistence. An Appendix tabulates all chemical analysis results to enable more quantitative and statistical studies of detection opportunities and limitations, as a part of a related project.

36 MATERIALS SCIENCE↗

Impacts of Stochastic Mixing in Idealized Convection-Permitting Simulations of Squall Lines

This study investigates impacts of altering subgrid-scale mixing in “convection-permitting” kilometer-scale horizontal-grid-spacing (Δ h ) simulations by applying either constant or stochastic multiplicative factors to the horizontal mixing coefficients within the Weather Research and Forecasting Model. In quasi-idealized 1-km Δ h simulations of two observationally based squall-line cases, constant enhanced mixing produces larger updraft cores that are more dilute at upper levels, weakens the cold pool, rear-inflow jet, and front-to-rear flow of the squall line, and degrades the model’s effective resolution. Reducing mixing by a constant multiplicative factor has the opposite effect on all metrics. Completely turning off parameterized horizontal mixing produces bulk updraft statistics and squall-line mesoscale structure closest to an LES “benchmark” among all 1-km simulations, although the updraft cores are too undilute. The stochastic mixing scheme, which applies a multiplicative factor to the mixing coefficients that varies stochastically in time and space, is employed at 0.5-, 1-, and 2-km Δ h . It generally reduces midlevel vertical velocities and enhances upper-level vertical velocities compared to simulations using the standard mixing scheme, with more substantial impacts at 1- and 2-km Δ h compared to 0.5-km Δ h . Further, the stochastic scheme also increases updraft dilution to better agree with the LES for one case, but has less impact on the other case. Stochastic mixing acts to weaken the cold pool but without a significant impact on squall-line propagation. It also does not affect the model’s overall effective resolution unlike applying constant multiplicative factors to the mixing coefficients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Acoustic radar investigations of boundary layer phenomena

A comparison is made between acoustic radar echoes and conventional meteorological data obtained from the WKY tower, for the purpose of better understanding the relationships between acoustic radar echoes and boundary layer processes. Two thunderstorm outflow cases are presented and compared to both acoustic radar data and Charba's gust front model. The acoustic radar echoes reveal the boundary between warm and cold air and other areas of mixing and strong thermal gradient quite well. The thunderstorm outflow of 27 June 1972 is found to compare with in most respects to Charba's gust front model. The major difference is the complete separation of the head from the main body of cold air, probably caused by erosion of the area behind the head by mixing with the ambient air. Two cases of nocturnal inversions caused by advection of warmer air aloft are presented. It is found that areas of turbulent mixing or strong thermal gradient can be identified quite easily in the acoustic radar record.

Marks, J. R.↗

Pollutant Formation in Monodisperse Fuel Spray Combustion

The combustion of liquid sprays represents an extremely important class of combustion processes. In the transition region, encompassing droplet sizes in the range of 25-80 micron diameter, the mixing and evaporation processes are both incomplete at the flame front and burning occurs in a combined diffusive and premixed fashion. Under these conditions, the relative importance of heterogeneous and homogeneous effects in dominating the combustion process is switched and gives rise to a number of interesting phenomena. NO (sub x) formation in monodisperse spray combustion was investigated with the following specific objectives: (1) to quantitatively determine the effect of droplet size, number density, etc. on NO sub x formation in monodisperse fuel spray combustion; and (2) to isolate the important physical and chemical phenomena in NO sub x formation in these combustion systems.

Cernansky, N. P.↗

Variability in Ocean Color Associated with Phytoplankton and Terrigenous Matter: Time Series Measurements and Algorithm Development at the FRONT Site on the New England Continental Shelf

Fronts in the coastal ocean describe areas of strong horizontal gradients in both physical and biological properties associated with tidal mixing and freshwater estuarine output (e.g. Simpson, 1981 and O Donnell, 1993). Related gradients in optically important constituents mean that fronts can be observed from space as changes in ocean color as well as sea surface temperature (e.g., Dupouy et al., 1986). This research program is designed to determine which processes and optically important constituents must be considered to explain ocean color variations associated with coastal fronts on the New England continental shelf, in particular the National Ocean Partnership Program (NOPP) Front Resolving Observational Network with Telemetry (FRONT) site. This site is located at the mouth of Long Island sound and was selected after the analysis of 12 years of AVHRR data showed the region to be an area of strong frontal activity (Ullman and Cornillon, 1999). FRONT consists of a network of modem nodes that link bottom mounted Acoustic Doppler Current Profilers (ADCPs) and profiling arrays. At the center of the network is the Autonomous Vertically Profiling Plankton Observatory (AVPPO) (Thwaites et al. 1998). The AVPPO consists of buoyant sampling vehicle and a trawl-resistant bottom-mounted enclosure, which holds a winch, the vehicle (when not sampling), batteries, and controller. Three sampling systems are present on the vehicle, a video plankton recorder, a CTD with accessory sensors, and a suite of bio-optical sensors including Satlantic OCI-200 and OCR-200 spectral radiometers and a WetLabs ac-9 dual path absorption and attenuation meter. At preprogrammed times the vehicle is released, floats to the surface, and is then winched back into the enclosure with power and data connection maintained through the winch cable. Communication to shore is possible through a bottom cable and nearby surface telemetry buoy, equipped with a mobile modem, giving the capability for near-real time data transmission and interactive sampling control.

Morrison, John R.↗

A Conceptual Model for Integrating Design Thinking and Lean Startup Methods into the Innovation Process

Design thinking is a methodology that comes from the industrial design realm and is centred on culling better insight into needs from users. Another popular methodology is based on gaining insight on the potential of an opportunity through experimentation, testing, and iterating with users. These are commonly referred to as lean startup methods. However, from a research perspective, we still do not know the most effective way to implement these user-focused design methods within the innovation process within organizations, and which aspects of the design process are the most impactful in developing new opportunities. In this research, we propose a high-level conceptual process model on how user-focused design methods such as design thinking and lean startup methods can be integrated into the up-front innovation process within organizations. This is based on a mix of digital ethnography, business historiography, and ethnographic, embedded participant observation of an organization's development of an innovation process. We review the conceptual model, associated activities, and process considerations. The article concludes with thoughts on future research.

Design Innovation↗

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of Fuel Vaporization, Fuel-Air Mixing, and Combustion in Integrated Mixer-Flame Holders

Requirements to limit pollutant emissions from the gas turbine engines for the future High-Speed Civil Transport (HSCT) have led to consideration of various low-emission combustor concepts. One such concept is the Integrated Mixer-Flame Holder (IMFH). This report describes a series of IMFH analyses performed with KIVA-II, a multi-dimensional CFD code for problems involving sprays, turbulence, and combustion. To meet the needs of this study, KIVA-II's boundary condition and chemistry treatments are modified. The study itself examines the relationships between fuel vaporization, fuel-air mixing, and combustion. Parameters being considered include: mixer tube diameter, mixer tube length, mixer tube geometry (converging-diverging versus straight walls), air inlet velocity, air inlet swirl angle, secondary air injection (dilution holes), fuel injection velocity, fuel injection angle, number of fuel injection ports, fuel spray cone angle, and fuel droplet size. Cases are run with and without combustion to examine the variations in fuel-air mixing and potential for flashback due to the above parameters. The degree of fuel-air mixing is judged by comparing average, minimum, and maximum fuel/air ratios at the exit of the mixer tube, while flame stability is monitored by following the location of the flame front as the solution progresses from ignition to steady state. Results indicate that fuel-air mixing can be enhanced by a variety of means, the best being a combination of air inlet swirl and a converging-diverging mixer tube geometry. With the IMFH configuration utilized in the present study, flashback becomes more common as the mixer tube diameter is increased and is instigated by disturbances associated with the dilution hole flow.

Deur, J. M.↗

Deep Chandra Observation and Numerical Studies of the Nearest Cluster Cold Front in the Sky

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are approximately 10 per cent brighter than the surrounding gas and are separated by approximately 15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches approximately 5-10 per cent of the thermal pressure, reducing the gas density between the bright features.

hydrodynamics↗

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗

Design of an Ultra-Efficient GaN High Power Amplifier for Radar Front-Ends Using Active Harmonic Load-Pull

This work presents a new measurement technique, mixed-signal active harmonic load-pull (MSALP) developed by Anterverta-mw in partnership with Maury Microwave, that allows for wide-band ultra-high efficiency amplifiers to be designed using GaN technology. An overview of the theory behind active load-pull is presented and why load-pull is important for high-power device characterization. In addition, an example procedure is presented that outlines a methodology for amplifier design using this measurement system. Lastly, measured results of a 10W GaN amplifier are presented. This work aims to highlight the benefit of using this sophisticated measurement systems for to optimize amplifier design for real radar waveforms that in turn will simplify implementation of space-based radar systems

Active Load-pull↗

Modeling complex chemical effects in turbulent nonpremixed combustion

Virtually all of the energy derived from the consumption of combustibles occurs in systems which utilize turbulent fluid motion. Since combustion is largely related to the mixing of fluids and mixing processes are orders of magnitude more rapid when enhanced by turbulent motion, efficiency criteria dictate that chemically powered devices necessarily involve fluid turbulence. Where combustion occurs concurrently with mixing at an interface between two reactive fluid bodies, this mode of combustion is called nonpremixed combustion. This is distinct from premixed combustion where flame-fronts propagate into a homogeneous mixture of reactants. These two modes are limiting cases in the range of temporal lag between mixing of reactants and the onset of reaction. Nonpremixed combustion occurs where this lag tends to zero, while premixed combustion occurs where this lag tends to infinity. Many combustion processes are hybrids of these two extremes with finite non-zero lag times. Turbulent nonpremixed combustion is important from a practical standpoint because it occurs in gas fired boilers, furnaces, waste incinerators, diesel engines, gas turbine combustors, and afterburners etc. To a large extent, past development of these practical systems involved an empirical methodology. Presently, efficiency standards and emission regulations are being further tightened (Correa 1993), and empiricism has had to give way to more fundamental research in order to understand and effectively model practical combustion processes (Pope 1991). A key element in effective modeling of turbulent combustion is making use of a sufficiently detailed chemical kinetic mechanism. The prediction of pollutant emission such as oxides of nitrogen (NO(x)) and sulphur (SO(x)) unburned hydrocarbons, and particulates demands the use of detailed chemical mechanisms. It is essential that practical models for turbulent nonpremixed combustion are capable of handling large numbers of 'stiff' chemical species equations.

Smith, Nigel S. A.↗

High-Fidelity Simulations of Gravity Currents Using a High-Order Finite-Difference Spectral Vanishing Viscosity Approach

This numerical work investigates the potential of a high-order finite-difference spectral vanishing viscosity approach to simulate gravity currents at high Reynolds numbers. The method introduces targeted numerical dissipation at small scales through altering the discretisation of the second derivatives of the viscous terms in the incompressible Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of the SVV method is straightforward and comes with a negligible additional computational cost. In order to assess the ability of this high-order finite-difference spectral vanishing viscosity approach, we performed large-eddy simulations (LES) of a gravity current in a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data are compared with a direct numerical simulation (DNS) based on more than 800 million mesh nodes, and with experimental measurements. A framework for the energy budget is introduced to investigate the behaviour of the gravity current. First, it is found that the DNS is in good agreement with the experimental data for the evolution of the front location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes compared to the DNS, can reproduce the main features of the gravity currents, with the SVV model yielding slightly more accurate results. It is also found that the dynamic Smagorinsky model performs better than its static version. For the present study, the static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the SVV model, because the latter does not require the calculation of explicit SGS terms in the Navier-Stokes equations nor spatial filtering operations.

direct numerical simulation↗

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model con figurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

58 GEOSCIENCES↗

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

54 ENVIRONMENTAL SCIENCES↗

Glacial-to-interglacial Changes in NADW Fluxes?

Interglacial gradients in delta 13C between Atlantic and Pacific deep waters reflect differences between low-nutrient, 13C-enriched North Atlantic Deep Water (NADW) and high-nutrient, 13C-depleted Pacific Deep Water. Reduced Atlantic-Pacific delta 13C and cadmium gradients at the last glacial maximum have been used to suggest substantial replacement of NADW with nutrient-rich Antarctic Bottom Water (Boyle and Keigwin, 1982; Shackleton et al., 1983). We show that the Atlantic delta 13C signal is linked directly to North Atlantic polar-front migration, as reflected by planktonic foraminiferal faunas.

Mix, A. C.↗

A New Approach to Model Delamination Growth in Fatigue Using the Virtual Crack Closure Technique Without Re-Meshing

A modeling approach is proposed to simulate delamination propagation in fatigue that combines the Virtual Crack Closure Technique (VCCT) with a progressive nodal release strategy. The progressive nodal release alleviates the artificial stress concentrations found when using the VCCT with instantaneous release to model 3D delamination without re-meshing. This enables crack shapes that do not conform to the underlying mesh to be readily simulated. The progressive release is formulated such that no artificial growth or healing occurs as a result of a change in loading conditions. Furthermore, expressions for calculating energy release rates at partially released nodes are proposed and applied to arbitrarily shaped delamination fronts. The shape of the delamination front is considered explicitly in the energy release rate calculations, and is determined by considering the normalized crack positions associated with the partial nodal release. Contrary to previous implementations, the local delamination propagation direction is not assumed to follow mesh lines, but instead is computed as part of the iterative procedure proposed to determine the maximum energy release along the crack front. The results obtained suggest the approach can accurately simulate Mode I, Mixed Mode I/II and Mode II fatigue delamination growth and be used to simulate growth of arbitrarily shaped cracks that do not conform to the underlying mesh.

N V De Carvalho↗