Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “UNIFIED MODEL”

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 829 records · Page 46

Matter inflation with A{sub 4} flavour symmetry breaking

We discuss model building in tribrid inflation, which is a framework for realising inflation in the matter sector of supersymmetric particle physics models. The inflaton is a D-flat combination of matter fields, and inflation ends by a phase transition in which some Higgs field obtains a vacuum expectation value. We first describe the general procedure for implementing tribrid inflation in realistic models of particle physics that can be applied to a wide variety of BSM particle physics models around the GUT scale. We then demonstrate how the procedure works for an explicit lepton flavour model based on an A{sub 4} family symmetry. The model is both predictive and phenomenologically viable, and illustrates how tribrid inflation connects cosmological and particle physics parameters. In particular, it predicts a relation between the neutrino Yukawa coupling and the running of the spectral index α{sub s}. We also show how topological defects from the flavour symmetry breaking can be avoided automatically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SUGRA new inflation with Heisenberg symmetry

We propose a realisation of ''new inflation'' in supergravity (SUGRA), where the flatness of the inflaton potential is protected by a Heisenberg symmetry. Inflation can be associated with a particle physics phase transition, with the inflaton being a (D-flat) direction of Higgs fields which break some symmetry at high energies, e.g. of GUT Higgs fields or of Higgs fields for flavour symmetry breaking. This is possible since compared to a shift symmetry, which is usually used to protect a flat inflaton potential, the Heisenberg symmetry is compatible with a (gauge) non-singlet inflaton field. In contrast to conventional new inflation models in SUGRA, where the predictions depend on unknown parameters of the Kaehler potential, the model with Heisenberg symmetry makes discrete predictions for the primordial perturbation parameters which depend only on the order n at which the inflaton appears in the effective superpotential. The predictions for the spectral index n{sub s} can be close to the best-fit value of the latest Planck 2013 results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

How to approach continuum physics in the lattice Weinberg-Salam model

We investigate the lattice Weinberg-Salam model without fermions numerically for the realistic choice of coupling constants correspondent to the value of the Weinberg angle {theta}{sub W{approx}}30 deg., and bare fine structure constant around {alpha}{approx}(1/150). We consider the values of the scalar self-coupling corresponding to Higgs mass M{sub H{approx}}100, 150, 270 GeV. It has been found that nonperturbative effects become important while approaching continuum physics within the lattice model. When the ultraviolet cutoff {Lambda}=({pi}/a) (where a is the lattice spacing) is increased and achieves the value around 1 TeV, one encounters the fluctuational region (on the phase diagram of the lattice model), where the fluctuations of the scalar field become strong. The classical Nambu monopole can be considered as an embryo of the unphysical symmetric phase within the physical phase. In the fluctuational region quantum Nambu monopoles are dense, and therefore, the use of the perturbation expansion around the trivial vacuum in this region is limited. Further increase of the cutoff is accompanied by a transition to the region of the phase diagram, where the scalar field is not condensed (this happens at the value of {Lambda} around 1.4 TeV for the considered lattice sizes). Within this region further increase of the cutoff is possible, although we do not observe this in detail due to the strong fluctuations of the gauge boson correlator. Both above mentioned regions look unphysical. Therefore we come to the conclusion that the maximal value of the cutoff admitted within lattice electroweak theory cannot exceed the value of the order of 1 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Operations on Graphical Models with Plates

This paper explains how graphical models, for instance Bayesian or Markov networks, can be extended to model problems in data analysis and learning. This provides a unified framework that combines lessons learned from the artificial intelligence, statistical and connectionist communities. This also offers a set of principles for developing a software generator for data analysis, whereby a learning or discovery system can be compiled from specifications. Many of the popular learning algorithms can be compiled in this way from graphical specifications. While in a sense this paper is a multidisciplinary review of learning, the main contribution here is the presentation of the material within the unifying framework of graphical models, and the observation that, as a result, the process of developing learning algorithms can be partly automated.

Buntine, Wray L.↗

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗

Atomistic Method Applied to Computational Modeling of Surface Alloys

The formation of surface alloys is a growing research field that, in terms of the surface structure of multicomponent systems, defines the frontier both for experimental and theoretical techniques. Because of the impact that the formation of surface alloys has on surface properties, researchers need reliable methods to predict new surface alloys and to help interpret unknown structures. The structure of surface alloys and when, and even if, they form are largely unpredictable from the known properties of the participating elements. No unified theory or model to date can infer surface alloy structures from the constituents properties or their bulk alloy characteristics. In spite of these severe limitations, a growing catalogue of such systems has been developed during the last decade, and only recently are global theories being advanced to fully understand the phenomenon. None of the methods used in other areas of surface science can properly model even the already known cases. Aware of these limitations, the Computational Materials Group at the NASA Glenn Research Center at Lewis Field has developed a useful, computationally economical, and physically sound methodology to enable the systematic study of surface alloy formation in metals. This tool has been tested successfully on several known systems for which hard experimental evidence exists and has been used to predict ternary surface alloy formation (results to be published: Garces, J.E.; Bozzolo, G.; and Mosca, H.: Atomistic Modeling of Pd/Cu(100) Surface Alloy Formation. Surf. Sci., 2000 (in press); Mosca, H.; Garces J.E.; and Bozzolo, G.: Surface Ternary Alloys of (Cu,Au)/Ni(110). (Accepted for publication in Surf. Sci., 2000.); and Garces, J.E.; Bozzolo, G.; Mosca, H.; and Abel, P.: A New Approach for Atomistic Modeling of Pd/Cu(110) Surface Alloy Formation. (Submitted to Appl. Surf. Sci.)). Ternary alloy formation is a field yet to be fully explored experimentally. The computational tool, which is based on the BFS (Bozzolo, Ferrante, and Smith) method for the calculation of the energetics, consists of a small number of simple PCbased computer codes that deal with the different aspects of surface alloy formation. Two analysis modes are available within this package. The first mode provides an atom-by-atom description of real and virtual stages 1. during the process of surface alloying, based on the construction of catalogues of configurations where each configuration describes one possible atomic distribution. BFS analysis of this catalogue provides information on accessible states, possible ordering patterns, and details of island formation or film growth. More importantly, it provides insight into the evolution of the system. Software developed by the Computational Materials Group allows for the study of an arbitrary number of elements forming surface alloys, including an arbitrary number of surface atomic layers. The second mode involves large-scale temperature-dependent computer 2. simulations that use the BFS method for the energetics and provide information on the dynamic processes during surface alloying. These simulations require the implementation of Monte-Carlo-based codes with high efficiency within current workstation environments. This methodology capitalizes on the advantages of the BFS method: there are no restrictions on the number or type of elements or on the type of crystallographic structure considered. This removes any restrictions in the definition of the configuration catalogues used in the analytical calculations, thus allowing for the study of arbitrary ordering patterns, ultimately leading to the actual surface alloy structure. Moreover, the Monte Carlo numerical technique used for the large-scale simulations allows for a detailed visualization of the simulated process, the main advantage of this type of analysis being the ability to understand the underlying features that drive these processes. Because of the simplicity of the BFS method for e energetics used in these calculations, a detailed atom-by-atom analysis can be performed at any point in the simulation, providing necessary insight on the details of the process. The main objective of this research program is to develop a tool to guide experimenters in understanding and interpreting often unexpected results in alloy formation experiments. By reducing the computational effort without losing physical accuracy, we expect that powerful simulation tools will be developed in the immediate future, which will allow material scientists to easily visualize and analyze processes at a level not achievable experimentally.

Bozzolo, Guillermo H.↗

Dissecting jets and missing energy searches using $n$-body extended simplified models

Simplified Models are a useful way to characterize new physics scenarios for the LHC. Particle decays are often represented using non-renormalizable operators that involve the minimal number of fields required by symmetries. Generalizing to a wider class of decay operators allows one to model a variety of final states. This approach, which we dub the $n$-body extension of Simplified Models, provides a unifying treatment of the signal phase space resulting from a variety of signals. In this paper, we present the first application of this framework in the context of multijet plus missing energy searches. The main result of this work is a global performance study with the goal of identifying which set of observables yields the best discriminating power against the largest Standard Model backgrounds for a wide range of signal jet multiplicities. Our analysis compares combinations of one, two and three variables, placing emphasis on the enhanced sensitivity gain resulting from non-trivial correlations. Utilizing boosted decision trees, we compare and classify the performance of missing energy, energy scale and energy structure observables. We demonstrate that including an observable from each of these three classes is required to achieve optimal performance. In conclusion, this work additionally serves to establish the utility of $n$-body extended Simplified Models as a diagnostic for unpacking the relative merits of different search strategies, thereby motivating their application to new physics signatures beyond jets and missing energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Long-lived time-dependent remnants during cosmological symmetry breaking: From inflation to the electroweak scale

Through a detailed numerical investigation in three spatial dimensions, we demonstrate that long-lived time-dependent field configurations emerge dynamically during symmetry breaking in an expanding de Sitter spacetime. We investigate two situations: a single scalar field with a double-well potential and an SU(2) non-Abelian Higgs model. For the single scalar, we show that large-amplitude oscillon configurations emerge spontaneously and persist to contribute about 1.2% of the energy density of the Universe. We also show that for a range of parameters, oscillon lifetimes are enhanced by the expansion and that this effect is a result of parametric resonance. For the SU(2) case, we see about 4% of the final energy density in oscillons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Oklahoma Center for High Energy Physics (OCHEP)

The DOE EPSCoR implementation grant, with the support from the State of Oklahoma and from the three universities, Oklahoma State University, University of Oklahoma and Langston University, resulted in establishing of the Oklahoma Center for High Energy Physics (OCHEP) in 2004. Currently, OCHEP continues to flourish as a vibrant hub for research in experimental and theoretical particle physics and an educational center in the State of Oklahoma. All goals of the original proposal were successfully accomplished. These include foun- dation of a new experimental particle physics group at OSU, the establishment of a Tier 2 computing facility for the Large Hadron Collider (LHC) and Tevatron data analysis at OU and organization of a vital particle physics research center in Oklahoma based on resources of the three universities. OSU has hired two tenure-track faculty members with initial support from the grant funds. Now both positions are supported through OSU budget. This new HEP Experimental Group at OSU has established itself as a full member of the Fermilab D0 Collaboration and LHC ATLAS Experiment and has secured external funds from the DOE and the NSF. These funds currently support 2 graduate students, 1 postdoctoral fellow, and 1 part-time engineer. The grant initiated creation of a Tier 2 computing facility at OU as part of the Southwest Tier 2 facility, and a permanent Research Scientist was hired at OU to maintain and run the facility. Permanent support for this position has now been provided through the OU university budget. OCHEP represents a successful model of cooperation of several universities, providing the establishment of critical mass of manpower, computing and hardware resources. This led to increasing Oklahoma's impact in all areas of HEP, theory, experiment, and computation. The Center personnel are involved in cutting edge research in experimental, theoretical, and computational aspects of High Energy Physics with the research areas ranging from the search for new phenomena at the Fermilab Tevatron and the CERN Large Hadron Collider to theoretical modeling, computer simulation, detector development and testing, and physics analysis. OCHEP faculty members participating on the D0 collaboration at the Fermilab Tevatron and on the ATLAS collaboration at the CERN LHC have made major impact on the Standard Model (SM) Higgs boson search, top quark studies, B physics studies, and measurements of Quantum Chromodynamics (QCD) phenomena. The OCHEP Grid computing facility consists of a large computer cluster which is playing a major role in data analysis and Monte Carlo productions for both the D0 and ATLAS experiments. Theoretical efforts are devoted to new ideas in Higgs bosons physics, extra dimensions, neutrino masses and oscillations, Grand Unified Theories, supersymmetric models, dark matter, and nonperturbative quantum field theory. Theory members are making major contributions to the understanding of phenomena being explored at the Tevatron and the LHC. They have proposed new models for Higgs bosons, and have suggested new signals for extra dimensions, and for the search of supersymmetric particles. During the seven year period when OCHEP was partially funded through the DOE EPSCoR implementation grant, OCHEP members published over 500 refereed journal articles and made over 200 invited presentations at major conferences. The Center is also involved in education and outreach activities by offering summer research programs for high school teachers and college students, and organizing summer workshops for high school teachers, sometimes coordinating with the Quarknet programs at OSU and OU. The details of the Center can be found in http://ochep.phy.okstate.edu.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Confronting generalized hidden local symmetry chiral model with the ALEPH data on the decay {tau}{sup -{yields}{pi}+{pi}-{pi}-{nu}}{sub {tau}}

The generalized hidden local symmetry (GHLS) model is the chiral model of pseudoscalar, vector, and axial vector mesons and their interactions. It contains also the couplings of strongly interacting particles with electroweak gauge bosons. Here, GHLS model is confronted with the ALEPH data on the decay {tau}{sup -{yields}{pi}-{pi}-{pi}+{nu}}{sub {tau}}. It is shown that the invariant mass spectrum of final pions in this decay calculated in GHLS framework with the single a{sub 1}(1260) resonance disagrees with the experimental data at any reasonable number of free GHLS parameters. Two modifications of GHLS model based on inclusion of two additional heavier axial vector mesons are studied. One of them, which gives a good description of the ALEPH data with all the parameters kept free, is shown to result in very large {Gamma}{sub a{sub 1}{sup {+-}}{sub {yields}{pi}}{sup {+-}}{sub {gamma}}}partial width. The other scheme with the GHLS parameters fixed in a way that the universality is preserved and the observed central value of {Gamma}{sub a{sub 1}{sup {+-}}{sub {yields}{pi}}{sup {+-}}{sub {gamma}}}is reached, which results in a good description of the three pion spectrum in {tau}{sup -{yields}{pi}+{pi}-{pi}-{nu}}{sub {tau}}decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A unified algorithm for predicting partition coefficients for PBPK modeling of drugs and environmental chemicals

The algorithms in the literature focusing to predict tissue:blood PC (P{sub tb}) for environmental chemicals and tissue:plasma PC based on total (K{sub p}) or unbound concentration (K{sub pu}) for drugs differ in their consideration of binding to hemoglobin, plasma proteins and charged phospholipids. The objective of the present study was to develop a unified algorithm such that P{sub tb}, K{sub p} and K{sub pu} for both drugs and environmental chemicals could be predicted. The development of the unified algorithm was accomplished by integrating all mechanistic algorithms previously published to compute the PCs. Furthermore, the algorithm was structured in such a way as to facilitate predictions of the distribution of organic compounds at the macro (i.e. whole tissue) and micro (i.e. cells and fluids) levels. The resulting unified algorithm was applied to compute the rat P{sub tb}, K{sub p} or K{sub pu} of muscle (n = 174), liver (n = 139) and adipose tissue (n = 141) for acidic, neutral, zwitterionic and basic drugs as well as ketones, acetate esters, alcohols, aliphatic hydrocarbons, aromatic hydrocarbons and ethers. The unified algorithm reproduced adequately the values predicted previously by the published algorithms for a total of 142 drugs and chemicals. The sensitivity analysis demonstrated the relative importance of the various compound properties reflective of specific mechanistic determinants relevant to prediction of PC values of drugs and environmental chemicals. Overall, the present unified algorithm uniquely facilitates the computation of macro and micro level PCs for developing organ and cellular-level PBPK models for both chemicals and drugs.

60 APPLIED LIFE SCIENCES↗

Fatigue life prediction modeling for turbine hot section materials

This paper presents a summary of the life prediction methods developed under the NASA Lewis Research Center's Hot Section Technology (HOST) program. A major objective of the fatigue and fracture efforts under the HOST program was to significantly improve the analytic life prediction tools used by the aeronautical gas turbine engine industry. This has been achieved in the areas of high-temperature thermal and mechanical fatigue of bare and coated high-temperature superalloys. Such technical improvements will eventually reduce life cycle costs. The cyclic crack initiation and propagation resistance of nominally isotropic polycrystalline alloys and highly anisotropic single crystal alloys has been addressed. A sizeable data base has been generated for three alloys [cast PWA 1455 (B-1900 + Hf), wrought Inconel 718, and cast single crystal PWA 1480] in bare and coated conditions. Two coatings systems, diffusion aluminide (PWA 273) and plasma sprayed MCrAIY overlay (PWA 286) were employed. Life prediction modeling efforts were devoted to creep-fatigue interaction, oxidation, coatings interactions, multiaxiality of stress-strain states, mean stress effects, cumulative damage, and thermomechanical fatigue. The fatigue crack initiation life models developed to date include the Cyclic Damage Accumulation (CDA) Model of Pratt & Whitney and the Total Strain Version of Strainrange Partitioning (TS-SRP) of NASA Lewis for nominally isotropic materials, and the Tensile Hysteretic Energy Model of Pratt & Whitney for anisotropic superalloys. The fatigue model being developed by the General Electric Company is based upon the concepts of Path-lndependent Integrals (PII) for describing cyclic crack growth under complex nonlinear response at the crack tip due to thermomechanical loading conditions. A micromechanistic oxidation crack extension model has been derived by researchers at Syracuse University. The models are described and discussed in the paper. Only limited verification has been achieved to-date as several of the technical programs are still in progress and the verification tasks are scheduled, quite naturally, near the conclusion of the program. To-date, efforts have concentrated on development of independent models for cyclic constitutive behavior, cyclic crack initiation, and cyclic crack propagation. The transition between crack initiation and crack propagation has not been thoroughly researched as yet, and the integration of these models into a unified life prediction method has not been addressed.

Metal fatigue↗

The V471 Tauri system: A multi-data-type probe

V471 Tauri, a white dwarf–red dwarf eclipsing binary (EB) in the Hyades, is well known for stimulating development of common envelope theory, whereby novae and other cataclysmic variables form from much wider binaries by catastrophic orbit shrinkage. Our evaluation of a recent imaging search that reported negative results for a much postulated third body shows that the object could have escaped detection or may have actually been seen. The balance of evidence continues to favor a brown dwarf companion about 12 AU from the EB. A recently developed algorithm finds unified solutions from three data types. New radial velocities (RVs) of the red dwarf and BVR{sub C}I{sub C} light curves are solved simultaneously along with white dwarf and red dwarf RVs from the literature, uvby data, the Microvariability and Oscillations of Stars mission light curve, and 40 years of eclipse timings. Precision-based weighting is the key to proper information balance among the various data sets. Timewise variation of modeled starspots allows unified solution of multiple data eras. Light-curve amplitudes strongly suggest decreasing spottedness from 1976 to about 1980, followed by approximately constant spot coverage from 1981 to 2005. An explanation is proposed for lack of noticeable variation in 1981 light curves, in terms of competition between spot and tidal variations. Photometric–spectroscopic distance is estimated. The red dwarf mass comes out larger than normal for a K2 V star, and even larger than adopted in several structure and evolution papers. An identified cause for this result is that much improved red dwarf RV curves now exist.

79 ASTRONOMY AND ASTROPHYSICS↗

Electroweak symmetry breaking and proton decay in SO(10) supersymmetric GUT with TeV W{sub R}

In a recent paper, we proposed a new class of supersymmetric SO(10) models for neutrino masses where the TeV-scale electroweak symmetry is SU(2){sub L}xSU(2){sub R}xU(1){sub B-L}, making the associated gauge bosons W{sub R} and Z{sup '} accessible at the Large Hadron Collider. We showed that there exists a domain of Yukawa coupling parameters and symmetry breaking patterns which give an excellent fit to all fermion masses including neutrinos. In this sequel, we discuss an alternative Yukawa pattern which also gives a good fermion mass fit, and then study the predictions of both models for the proton lifetime. Consistency with current experimental lower limits on the proton lifetime require the squark masses of the first two generations to be larger than {approx}1.2 TeV. We also discuss how one can have simultaneous breaking of both SU(2){sub R}xU(1){sub B-L} and standard electroweak symmetries via radiative corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Unified Picture for Magnetic Correlations in Iron-Based Superconductors

The varying metallic antiferromagnetic correlations observed in iron-based superconductors are unified in a model consisting of both itinerant electrons and localized spins. The decisive factor is found to be the sensitive competition between the superexchange antiferromagnetism and the orbital-degenerate double-exchange ferromagnetism. Our results reveal the crucial role of Hund's rule coupling for the strongly correlated nature of the system and suggest that the iron-based superconductors are closer kin to manganites than cuprates in terms of their diverse magnetism and incoherent normal-state electron transport. This unified picture would be instrumental for exploring other exotic properties and the mechanism of superconductivity in this new class of superconductors.

36 MATERIALS SCIENCE↗

Multi-Domain Routing in Delay Tolerant Networks

The goal of Delay Tolerant Networking (DTN) is to provide the missing ingredient for the ever-growing collection of communicating nodes in our solar system to become a Solar System Internet (SSI). Great strides have been made in modeling particular types of DTNs, such as schedule- or discovery-based. Now, analogously to the Internet, these smaller DTNs can be considered routing domains which must be stitched together to form the overall SSI. In this paper, we propose a framework for cross-domain routing in DTNs as well as methodologies for detecting these sub-domains. Example time-varying networks are given to demonstrate the techniques proposed. A basic component is the mathematical theory of sheaves, which unifies the underlying model of DTN routing algorithms, by giving rise to routing sheaves – these can be defined for the dynamic and scheduled networks as noted above, and can also be used to define the interfaces between these domains in order to route across them. An immediate application would be routing across discovery-based networks connected by scheduled networks. These DTN subdomains remain elusive, however, and need to become well-defined and properly sized for tractable computability. In particular, a balance must be determined between areas that are too large (i.e. large matrix computations) versus areas that are too small (i.e. “many” single-noded domains). Moreover, the connections between the domains should, at least locally, be chosen to optimize data flow and connectivity: we address this in three ways. First, tools from persistent homology are given to understand underlying structures, reminiscent of hierarchies in the Internet Protocol (IP) addressing. Second, we construct a notion of temporal graph curvature based on network geometry to analyze flows induced by dynamical processes on these networks. Finally, Schrodinger Bridges, a tool arising from statistical physics, are proposed as a method of constructing flows on time-evolving networks with desirable properties such as speed, robustness, and load sensitivity. We construct an approach to temporal hypergraphs to simultaneously model unicast, multicast, and broadcast, using the language of scheme theory, and then consider DTN network coding as a way to achieve network-level computation and organization. The paper concludes with a discussion and ideas for future work.

Alan Hylton↗

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗