BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michael Ortiz

Michael Ortiz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAphysics.comp-phcond-mat.mtrl-sciComputational Engineering, Finance, and Sciencemath.APNumerical AnalysisBiological PhysicsCell BehaviorInformation Theorymath.ITmath.PRmath.STnlin.PSphysics.app-phphysics.data-anphysics.flu-dynStatistics Theory

Catalog footprint

What is connected

14works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2021arXiv

A proof of Taylor scaling for curvature-driven dislocation motion through random arrays of obstacles

We prove Taylor scaling for dislocation lines characterized by line-tension and moving by curvature under the action of an applied shear stress in a plane containing a random array of obstacles. Specifically, we show--in the sense of optimal scaling--that the critical applied shear stress for yielding, or percolation-like unbounded motion of the dislocation, scales in proportion to the square root of the obstacle density. For sufficiently small obstacle densities, Taylor scaling dominates the linear-scaling that results from purely energetic considerations and, therefore, characterizes the dominant rate-limiting mechanism in that regime.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A spatially adaptive phase-field model of fracture

Phase-field models of fracture introduce smeared cracks of width commensurate with a regularisation length parameter $ε$ and obeying a minimum energy principle. Mesh adaptivity naturally suggests itself as a means of supplying spatial resolution were needed while simultaneously keeping the computational size of the model as small as possible. Here, a variational-based spatial adaptivity is proposed for a phase-field model of fracture. The conventional phase-field model is generalised by allowing a spatial variation of the regularisation length $ε$ in the energy functional. The optimal spatial variation of the regularisation length then follows by energy minimisation in the same manner as the displacement and phase fields. The extended phase-field model is utilised as a basis for an adaptive mesh refinement strategy, whereby the mesh size is required to resolve the optimal length parameter locally. The resulting solution procedure is implemented in the framework of the finite element library FEniCS. Selected numerical tests suggest that the spatially adaptive phase-field model exhibits the same convergence rate as the conventional phase-field model, albeit with a vastly superior constant, which results in considerable computational savings.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mollified finite element approximants of arbitrary order and smoothness

The approximation properties of the finite element method can often be substantially improved by choosing smooth high-order basis functions. It is extremely difficult to devise such basis functions for partitions consisting of arbitrarily shaped polytopes. We propose the mollified basis functions of arbitrary order and smoothness for partitions consisting of convex polytopes. On each polytope an independent local polynomial approximant of arbitrary order is assumed. The basis functions are defined as the convolutions of the local approximants with a mollifier. The mollifier is chosen to be smooth, to have a compact support and a unit volume. The approximation properties of the obtained basis functions are governed by the local polynomial approximation order and mollifier smoothness. The convolution integrals are evaluated numerically first by computing the boolean intersection between the mollifier and the polytope and then applying the divergence theorem to reduce the dimension of the integrals. The support of a basis function is given as the Minkowski sum of the respective polytope and the mollifier. The breakpoints of the basis functions, i.e. locations with non-infinite smoothness, are not necessarily aligned with polytope boundaries. Furthermore, the basis functions are not boundary interpolating so that we apply boundary conditions with the non-symmetric Nitsche method as in immersed/embedded finite elements. The presented numerical examples confirm the optimal convergence of the proposed approximation scheme for Poisson and elasticity problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data-Driven Finite Elasticity

We extend to finite elasticity the Data-Driven formulation of geometrically linear elasticity presented in Conti, Müller, Ortiz, Arch.\ Ration.\ Mech.\ Anal.\ 229, 79-123, 2018. The main focus of this paper concerns the formulation of a suitable framework in which the Data-Driven problem of finite elasticity is well-posed in the sense of existence of solutions. We confine attention to deformation gradients $F \in L^p(Ω;\mathbb{R}^{n\times n})$ and first Piola-Kirchhoff stresses $P \in L^q(Ω;\mathbb{R}^{n\times n})$, with $(p,q)\in(1,\infty)$ and $1/p+1/q=1$. We assume that the material behavior is described by means of a material data set containing all the states $(F,P)$ that can be attained by the material, and develop germane notions of coercivity and closedness of the material data set. Within this framework, we put forth conditions ensuring the existence of solutions. We exhibit specific examples of two- and three-dimensional material data sets that fit the present setting and are compatible with material frame indifference.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data-driven fracture mechanics

We present a new data-driven paradigm for variational brittle fracture mechanics. The fracture-related material modeling assumptions are removed and the governing equations stemming from variational principles are combined with a set of discrete data points, leading to a model-free data-driven method of solution. The solution at a given load step is identified as the point within the data set that best satisfies either the Kuhn-Tucker conditions stemming from the variational fracture problem or global minimization of a suitable energy functional, leading to data-driven counterparts of both the local and the global minimization approaches of variational fracture mechanics. Both formulations are tested on different test configurations with and without noise and for Griffith and R-curve type fracture behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ice penetration by a bluff-body melting probe

We analyze the operation of melting probes as a Stefan problem for the liquid/solid interface surrounding the probe. We assume that the liquid layer is thin and, therefore, amenable to analysis by lubrication theory. The resulting Stefan problem is solvable in closed form. The solution determines the dependence of the penetration speed on the temperature differential between the probe and the surrounding ice, the size, shape and weight of the probe, the viscosity of liquid water and the thermal properties of solid ice.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Oncotripsy: Targeting cancer cells selectively via resonant harmonic excitation

We investigate a method of selectively targeting cancer cells by means of ultrasound harmonic excitation at their resonance frequency, which we refer to as oncotripsy. The geometric model of the cells takes into account the cytoplasm, nucleus and nucleolus, as well as the plasma membrane and nuclear envelope. Material properties are varied within a pathophysiologically-relevant range. A first modal analysis reveals the existence of a spectral gap between the natural frequencies and, most importantly, resonant growth rates of healthy and cancerous cells. The results of the modal analysis are verified by simulating the fully-nonlinear transient response of healthy and cancerous cells at resonance. The fully nonlinear analysis confirms that cancerous cells can be selectively taken to lysis by the application of carefully tuned ultrasound harmonic excitation while simultaneously leaving healthy cells intact.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Large scale ab-initio simulations of dislocations

We present a novel methodology to compute relaxed dislocations core configurations, and their energies in crystalline metallic materials using large-scale \emph{ab-intio} simulations. The approach is based on MacroDFT, a coarse-grained density functional theory method that accurately computes the electronic structure but with sub-linear scaling resulting in a tremendous reduction in cost. Due to its implementation in \emph{real-space}, MacroDFT has the ability to harness petascale resources to study materials and alloys through accurate \emph{ab-initio} calculations. Thus, the proposed methodology can be used to investigate dislocation cores and other defects where long range elastic defects play an important role, such as in dislocation cores, grain boundaries and near precipitates in crystalline materials. We demonstrate the method by computing the relaxed dislocation cores in prismatic dislocation loops and dislocation segments in magnesium (Mg). We also study the interaction energy with a line of Aluminum (Al) solutes. Our simulations elucidate the essential coupling between the quantum mechanical aspects of the dislocation core and the long range elastic fields that they generate. In particular, our quantum mechanical simulations are able to describe the logarithmic divergence of the energy in the far field as is known from classical elastic theory. In order to reach such scaling, the number of atoms in the simulation cell has to be exceedingly large, and cannot be achieved with the state-of-the-art density functional theory implementations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Data-driven computational mechanics

We develop a new computing paradigm, which we refer to as data-driven computing, according to which calculations are carried out directly from experimental material data and pertinent constraints and conservation laws, such as compatibility and equilibrium, thus bypassing the empirical material modeling step of conventional computing altogether. Data-driven solvers seek to assign to each material point the state from a prespecified data set that is closest to satisfying the conservation laws. Equivalently, data-driven solvers aim to find the state satisfying the conservation laws that is closest to the data set. The resulting data-driven problem thus consists of the minimization of a distance function to the data set in phase space subject to constraints introduced by the conservation laws. We motivate the data-driven paradigm and investigate the performance of data-driven solvers by means of two examples of application, namely, the static equilibrium of nonlinear three-dimensional trusses and linear elasticity. In these tests, the data-driven solvers exhibit good convergence properties both with respect to the number of data points and with regard to local data assignment. The variational structure of the data-driven problem also renders it amenable to analysis. We show that, as the data set approximates increasingly closely a classical material law in phase space, the data-driven solutions converge to the classical solution. We also illustrate the robustness of data-driven solvers with respect to spatial discretization. In particular, we show that the data-driven solutions of finite-element discretizations of linear elasticity converge jointly with respect to mesh size and approximation by the data set.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coarse-graining Kohn-Sham Density Functional Theory

We present a real-space formulation for coarse-graining Kohn-Sham Density Functional Theory that significantly speeds up the analysis of material defects without appreciable loss of accuracy. The approximation scheme consists of two steps. First, we develop a linear-scaling method that enables the direct evaluation of the electron density without the need to evaluate individual orbitals. We achieve this by performing Gauss quadrature over the spectrum of the linearized Hamiltonian operator appearing in each iteration of the self-consistent field method. Building on the linear-scaling method, we introduce a spatial approximation scheme resulting in a coarse-grained Density Functional Theory. The spatial approximation is adapted so as to furnish fine resolution where necessary and to coarsen elsewhere. This coarse-graining step enables the analysis of defects at a fraction of the original computational cost, without any significant loss of accuracy. Furthermore, we show that the coarse-grained solutions are convergent with respect to the spatial approximation. We illustrate the scope, versatility, efficiency and accuracy of the scheme by means of selected examples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal Uncertainty Quantification

We propose a rigorous framework for Uncertainty Quantification (UQ) in which the UQ objectives and the assumptions/information set are brought to the forefront. This framework, which we call \emph{Optimal Uncertainty Quantification} (OUQ), is based on the observation that, given a set of assumptions and information about the problem, there exist optimal bounds on uncertainties: these are obtained as values of well-defined optimization problems corresponding to extremizing probabilities of failure, or of deviations, subject to the constraints imposed by the scenarios compatible with the assumptions and information. In particular, this framework does not implicitly impose inappropriate assumptions, nor does it repudiate relevant information. Although OUQ optimization problems are extremely large, we show that under general conditions they have finite-dimensional reductions. As an application, we develop \emph{Optimal Concentration Inequalities} (OCI) of Hoeffding and McDiarmid type. Surprisingly, these results show that uncertainties in input parameters, which propagate to output uncertainties in the classical sensitivity analysis paradigm, may fail to do so if the transfer functions (or probability distributions) are imperfectly known. We show how, for hierarchical structures, this phenomenon may lead to the non-propagation of uncertainties or information across scales. In addition, a general algorithmic framework is developed for OUQ and is tested on the Caltech surrogate model for hypervelocity impact and on the seismic safety assessment of truss structures, suggesting the feasibility of the framework for important complex systems. The introduction of this paper provides both an overview of the paper and a self-contained mini-tutorial about basic concepts and issues of UQ.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Convergence Analysis of Meshfree Approximation Schemes

This work is concerned with the formulation of a general framework for the analysis of meshfree approximation schemes and with the convergence analysis of the Local Maximum-Entropy (LME) scheme as a particular example. We provide conditions for the convergence in Sobolev spaces of schemes that are n-consistent, in the sense of exactly reproducing polynomials of degree less or equal to n, and whose basis functions are of rapid decay. The convergence of the LME in a locally Sobolev space follows as a direct application of the general theory. The analysis shows that the convergence order is linear in h, a measure of the density of the point set. The analysis also shows how to parameterize the LME scheme for optimal convergence. Because of the convex approximation property of LME, its behavior near the boundary is singular and requires additional analysis. For the particular case of polyhedral domains we show that, away from a small singular part of the boundary, any Sobolev function can be approximated by means of the LME scheme. With the aid of a capacity argument, we further obtain approximation results with truncated LME basis functions in H1 and for spatial dimension d > 2.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mesoscopic approach to granular crystal dynamics

We present a mesoscopic approach to granular crystal dynamics, which comprises a three-dimensional finite-element model and a one-dimensional regularized contact model. The approach investigates the role of vibrational-energy trapping effects in the dynamic behavior of one-dimensional chains of particles in contact (i.e., granular crystals), under small to moderate impact velocities. The only inputs of the models are the geometry and the elastic material properties of the individual particles that form the system. We present detailed verification results and validate the model comparing its predictions with experimental data. This approach provides a physically sound, first-principle description of dissipative losses in granular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

An Efficient Adaptive Procedure for Three-Dimensional Fragmentation Simulations

We present a simple set of data structures, and a collection of methods for constructing and updating the structures, designed to support the use of cohesive elements in simulations of fracture and fragmentation. Initially all interior faces in the triangulation are perfectly coherent, i.e., conforming in the usual finite element sense. Cohesive elements are inserted adaptively at interior faces when the effective traction acting on those face reaches the cohesive strength of the material. The insertion of cohesive elements changes the geometry of the boundary and, frequently, the topology of the model as well. The data structures and methods presented here are straightforward to implement and enable the efficient tracking of complex fracture and fragmentation processes. The efficiency and versatility of the approach is demonstrated with the aid of two examples of application to dynamic fracture.

0 citations0 reviewsNo code linkedOpen access