Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded boundary”

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 235 records · Page 13

Determination of Elastic Coefficients for a Composite with Crossing Steel - 20086

For a composite material with an embedding of crossing reinforcement bars, the elastic coefficients are calculated by solving a boundary-value problem (elastostatic problem). The background material is concrete whereas the bars are made of steel. The bars cross each other horizontally at a right angle. The boundary-value problem is solved numerically by using 3-D quadratic basis functions in ABAQUS{sup R}. Under external shear loading conditions, the effective elastic coefficients on the macroscale are quite larger for a cross-section in which the bars are immediately below the ones in the cross-section normal to it. The shear modulus in the off-plane direction is about one and half times larger than that in the in-plane direction. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Two-Axis Direct Fluid Shear Stress Sensor

This innovation is a miniature or micro sized semiconductor sensor design that provides two axis direct non-intrusive measurement of skin friction or wall shear stress in fluid flow. The sensor is fabricated by micro-electro-mechanical system (MEMS) technology, enabling small size and low cost reproductions. The sensors have been fabricated by utilizing MEMS fabrication processes to bond a sensing element wafer to a fluid coupling wafer. This layering technique provides for an out of plane dimension that is on the same order of length as the inplane dimensions. The sensor design has the following characteristics: a shear force collecting plate with dimensions that can be tailored to various application specific requirements such as spatial resolution, temporal resolution and shear force range and resolution. This plate is located coplanar to both the sensor body and flow boundary, and is connected to a dual axis gimbal structure by a connecting column or lever arm. The dual axis gimbal structure has torsional hinges with embedded piezoresistive torsional strain gauges which provide a voltage output that is correlated to the applied shear stress (and excitation current) on force collection plate that is located on the flow boundary surface (hence the transduction method). This combination of design elements create a force concentration and resolution structure that enables the generation of a large stress on the strain gauge from the small shear stress on the flow boundary wall. This design as well as the use of back side electrical contacts establishes a non-intrusive method to quantitatively measure the shear force vector on aerodynamic bodies.

Adcock, Edward E.↗

Analysis of the Harrier forebody/inlet design using computational techniques

Under the support of this Cooperative Agreement, computations of transonic flow past the complex forebody/inlet configuration of the AV-8B Harrier II have been performed. The actual aircraft configuration was measured and its surface and surrounding domain were defined using computational structured grids. The thin-layer Navier-Stokes equations were used to model the flow along with the Chimera embedded multi-grid technique. A fully conservative, alternating direction implicit (ADI), approximately-factored, partially flux-split algorithm was employed to perform the computation. An existing code was altered to conform with the needs of the study, and some special engine face boundary conditions were developed. The algorithm incorporated the Chimera technique and an algebraic turbulence model in order to deal with the embedded multi-grids and viscous governing equations. Comparison with experimental data has yielded good agreement for the simplifications incorporated into the analysis. The aim of the present research was to provide a methodology for the numerical solution of complex, combined external/internal flows. This is the first time-dependent Navier-Stokes solution for a geometry in which the fuselage and inlet share a wall. The results indicate the methodology used here is a viable tool for transonic aircraft modeling.

Chow, Chuen-Yen↗

Characterization of Aerodynamic Performance of Boundary-Layer-Ingesting Inlet Under Crosswind

NASA has been studying future transport concepts, envisioned to be technically realizable in the timeframe of 2020-2030, to meet environmental and performance goals. One concept receiving considerable interest involves a propulsion system embedded into a hybrid wing-body aircraft. While offering significant advantages in fuel savings and noise reduction by this concept, there are several technical challenges that are not encountered in the current fleet and must be overcome so as to deliver target performance and operability. One of these challenges is associated with an inlet system that ingests a significantly thick boundary layer, developing along the wing-body surface, into a serpentine diffuser before the flow meeting fan blades. The flow is subject to considerable total pressure loss and distorted at the fan face, much more significantly than in the inlet system of conventional aircraft. In our previous studies [1, 2], we have shown that through innovative design changes on the airframe surface, it is possible to simultaneously increase total pressure recovery and decrease distortion in the flow, without resorting to conventional penalty-ridden flow control concepts, such as vortex generator or boundary layer bleeding/suction. In the current study, we are interested in understanding the following issues: how the embedded propulsion system performs under a crosswind condition by studying in detail the flow characteristics of two inlets, the baseline and another optimized previously under the cruise condition. With the insight, it is hoped that it can help in the follow-on study by devising effective strategies to minimize flow distortion arising from the integration of an embedded-engine system into an airframe to the level acceptable to the operation and fuel consumption before 2030. To achieve these demanding goals, non-conventional concepts are called for; but technology gap is too big that it requires evolutionary approach by focusing various concepts and technologies needed in the next three generations of aircraft, respectively named as N+1, N+2, and N+3. Noticeably, considerable reduction in each category of 1 is required in N+2 (relative to Boeing 777-200 and GE90 engines) and N+3 (relative to Boeing 737-800 and CFM56-7B engines). In this study, concepts for N+2 is our interest. A concept that has potential to achieve these metrics and has been under intensive study is the hybrid wing body (HWB) airframe with a tightly integrated propulsion system, see Fig. 1. The inlet is non-circular at the entrance and the entering flow, no longer uniform or free of disturbances, and is now carrying with it a boundary layer developing along the fuselage; the inlet is thus known as boundary-layer-ingesting (BLI) inlet.

Liou, Meng-Sing↗

Proton Transport

The transport of protons across membranes is an essential process for both bioenergetics of modern cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATP), synthesized from adenosine diphosphate. ATP, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this Study was: how the same process can be accomplished with the aid of similar but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC) which is a good model of the biological membranes forming cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M$_2$ protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M$_2$ protein is 97 amino acids in length, but a fragment 25 amino acids long. which contains a transmembrane domain of 19 amino acids flanked by three amino acids on each side. is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules. and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small ions, such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel with the histidine residues in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel loses its structural integrity and one helix moves away from the remaining three. Additional information is contained within the original extended abstract.

Pohorille, Andrew↗

The M2 Proton Channel of Influenza Virus: How Does It Work?

The transport of protons across membranes is an essential process for both bioenergetics of modem cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATR), synthesized from adenosine diphosphate. ATR, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this study was: how the same process can be accomplished with the aid of similar, but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC), which is a good model of the biological membranes focusing cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M2 protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M2 protein is 97 amino acids in length, but a fragment 25 amino acids long, which contains a transmembrane domain of 19 amino acids flanked by 3 amino acids on each side, is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules, and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional, important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel, in which the histidine residues were in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel looses its structural integrity and one helix moves away from the remaining three. This result indicates that such a mechanism of proton transport is unlikely. The simulations revealed that translocation along a network of water molecules in the channel and tautomerization of the histidine residues in the M2 proteins in the most likely mechanism of proton transport. The results not only explain how a remarkably simple protein system can efficiently aid in the formation of proton gradients across cell walls, but also suggest how this system can be genetically re-engineered to become a directional, reversible proton pump. Such a pump can provide energy to laboratory-built models of simple cellular systems. If they were successfully constructed it would greatly advance our understanding of the beginnings of life and find important applications in medicine and pharmacology.

Pohorille, Andrew↗

Thriving in the Carbon-Aware Market: How to Account for Emissions in the Era of Carbon-Centered Trade Policies

Emerging global policies, such as the European Union's enacted Carbon Border Adjustment Mechanism and similar policies under development in Canada, Australia and the United Kingdom, will place a premium on goods traded into their territories with higher embedded emissions than those produced domestically. U.S. manufacturers could stand to benefit from such policies; given the investments U.S. industry has made to reduce the energy and emission intensities of its operations. For example, the overall GHG intensity of U.S. steel production in 2019 was ~0.96 t CO2/t steel, less than half that of China (~1.97 t CO2/t steel), and bested only by Italy. To realize these benefits, transparent, accurate, interoperable and accepted embedded emissions accounting and calculation methods are required. Achieving this requires overcoming challenges related to data availability, boundary definitions, and product definitions among others, both at individual facilities as well as through value chains. We will present technical findings on methodology considerations and data-availability constraints for determining the emissions of traded goods, using steel as a pilot and leveraging publicly available data. Issues such as determining the appropriate scope for emissions accounting, implications of the specificity of product chosen, emissions allocation in multi-product facilities, and enumeration of emissions for products manufactured across multiple facilities will be discussed. By sharing the results of our efforts, we aim to inform the development and execution of embedded emissions accounting methods from a technical perspective such that U.S. manufacturers can thrive in emerging global markets.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗

Inlet Flow Control and Prediction Technologies for Embedded Propulsion Systems

Fail-safe, hybrid, flow control (HFC) is a promising technology for meeting high-speed cruise efficiency, low-noise signature, and reduced fuel-burn goals for future, Hybrid-Wing-Body (HWB) aircraft with embedded engines. This report details the development of HFC technology that enables improved inlet performance in HWB vehicles with highly integrated inlets and embedded engines without adversely affecting vehicle performance. In addition, new test techniques for evaluating Boundary-Layer-Ingesting (BLI)-inlet flow-control technologies developed and demonstrated through this program are documented, including the ability to generate a BLI-like inlet-entrance flow in a direct-connect, wind-tunnel facility, as well as, the use of D-optimal, statistically designed experiments to optimize test efficiency and enable interpretation of results. Validated improvements in numerical analysis tools and methods accomplished through this program are also documented, including Reynolds-Averaged Navier-Stokes CFD simulations of steady-state flow physics for baseline, BLI-inlet diffuser flow, as well as, that created by flow-control devices. Finally, numerical methods were employed in a ground-breaking attempt to directly simulate dynamic distortion. The advances in inlet technologies and prediction tools will help to meet and exceed "N+2" project goals for future HWB aircraft.

McMillan, Michelle L.↗

Aerodynamic Design and Optimization of Fan Stage for Boundary Layer Ingestion Propulsion System

The present paper addresses the process of preliminary design of a low-pressure fan and outlet guide vane (OGV) of a boundary layer ingestion (BLI) propulsion system. The tail-cone thruster systems of NASA's STARC_ABL (Single-aisle Turboelectric Aircraft with an Aft Boundary-Layer propulsor) adopts an axi-symmetric BLI type inlet as opposed to other embedded engine systems. Thus, the focus of the present work is placed on maximizing the efficiency of the fan and OGV stages under a significant radial distortion. A parameterization with B-spline function for camber line angles, metal chord, thickness distribution and stacking axis of blades is presented. The flowpath lines are also parameterized by B-spline function and aggregated in the design system of blades. The design optimization with evolutionary algorithm is performed with constraints of fan pressure ratio, OGV exit swirl angle and nozzle exit properties. The inlet conditions for the turbo-machinery CFD (Computational Fluid Dynamics) domain and the design goal of the fan stage are driven by a propulsion airframe integration (PAI) model that uses a 3-D unstructured RANS (Reynolds Average Navier Stokes) solver and actuator disk model. The expected power saving of the BLI propulsor is quantified via PAI analysis and the resulting preliminary design of the fan stages is compared with a clean-inlet flow propulsor.

Turbo-machinery↗

Computational Micromechanics Model for the Analysis of Fiber Kinking in Unidirectional FRP

A computational micromechanics (CMM) model is developed to analyze fiber kinking, which is a failure mechanism that takes place in fiber-reinforced composites when they are loaded under longitudinal compression. The CMM model consists of a single AS4 carbon fiber with an initial misalignment embedded in an 8552 polymer matrix. The deformation of the model is governed by periodic boundary conditions (PBC). The relatively simple CMM model enables the evaluation of the role played by initial misalignment of the fiber, shear yielding of the matrix and fiber-matrix debonding. A novel microscale experimental technique devoted to the characterization of the longitudinal compressive strength of the fibers, X(sup f)(sub c) , was developed. By exercising the model and comparing it with several existing models in the literature, the nonlinear shear response of the composite lamina was shown to play a fundamental role not only in the prediction of the compressive strength X(sub c), but also during the post-peak regime in terms of residual stress σ(sub r) and fiber rotation φ. Finally, the influence of the fiber-matrix interface damage (not considered in most other fiber kinking models) on the fiber kinking phenomenon was assessed through a parametric study.

Polymer-matrix composites (PMC)↗

Physical properties of Saturn's rings

Spaceborne astronomy data on Saturn's rings is summarized. The main rings are composed primarily of ice, possibly of low density and including small amounts of reddish material; particle radii are primarily centimeter-to-several meters, but significant amounts of micron-sized dust exists in certain regions including regions of resonance-induced collisions. The particle surface albedo and size distribution varies globally and locally. Evidence suggests several embedded moonlets of 10 km radius within Encke's division. The cause of the abrupt regional boundaries and associated differences in particle properties, especially the cause of the very abrupt and qualitatively similar inner A and B ring boundaries; the gap-ringlet structure and broad plateau features observed in the C and Cassini regions; the F and G rings (assuming that Enceladus is responsible for the E ring); and the cause of spokes and the SED are not understood.

Cuzzi, J. N.↗

Interactive algebraic grid-generation technique

An algebraic grid generation technique and use of an associated interactive computer program are described. The technique, called the two boundary technique, is based on Hermite cubic interpolation between two fixed, nonintersecting boundaries. The boundaries are referred to as the bottom and top, and they are defined by two ordered sets of points. Left and right side boundaries which intersect the bottom and top boundaries may also be specified by two ordered sets of points. when side boundaries are specified, linear blending functions are used to conform interior interpolation to the side boundaries. Spacing between physical grid coordinates is determined as a function of boundary data and uniformly space computational coordinates. Control functions relating computational coordinates to parametric intermediate variables that affect the distance between grid points are embedded in the interpolation formulas. A versatile control function technique with smooth-cubic-spline functions is presented. The technique works best in an interactive graphics environment where computational displays and user responses are quickly exchanged. An interactive computer program based on the technique and called TBGG (two boundary grid generation) is also described.

Smith, R. E.↗

The morphology of strike-slip faults - Examples from the San Andreas Fault, California

The dilatational strains associated with vertical faults embedded in a horizontal plate are examined in the framework of fault kinematics and simple displacement boundary conditions. Using boundary element methods, a sequence of examples of dilatational strain fields associated with commonly occurring strike-slip fault zone features (bends, offsets, finite rupture lengths, and nonuniform slip distributions) is derived. The combinations of these strain fields are then used to examine the Parkfield region of the San Andreas fault system in central California.

Bilham, Roger↗

Local grid refinement for transonic flow problems

The present use of locally refined Cartesian grids to solve transonic flow problems about three-dimensional aircraft configurations obviates surface-conforming grid generation through an embedding of surface-geometry paneling in the grid. Accurate resolution of flow close to the boundary, and in regions with strong velocity gradients, is achieved via hierarchical local refinement which subdivides a given grid cell into eight cells. Fast and reliable convergence is obtained by combining several preconditioners and damping strategies. Methods are suggested for preclusion of global convergence problems.

Melvin, Robin G.↗

Modeling the growth and interaction of fractures

The research objectives are to develop accurate, efficient numerical tool for fracture mechanics analysis with the capabilities of the method including the modeling of multiple site damage and the simulation of crack propagation. The fracture is modeled using a distribution of edge dislocations. Fundamental solutions for dislocations in an infinite medium, in the neighborhood of a circular inhomogeneity, and near a bi-material interface are incorporated into the procedure to allow exceptionally accurate modeling of fracture geometries involving these types of inhomogeneities. To model arbitrary (finite) domains under general boundary conditions, a combination of the dislocation methodology and boundary element techniques is used. The dislocation scheme is a natural means to model the crack(s), and the boundary element method is an effective way to model the region in which the crack is embedded. The hybrid scheme is relatively easy to use and is very accurate. The research accomplishments are presented.

Mear, Mark E.↗