首页/文章/ 详情

数值分析(Numerical analysis)高影响力研究分析

9月前浏览449

文一:

 

浸入式边界法

摘要:

本文研究了浸没边界(IB)方法的数学结构,该方法用于流体-结构相互作用的计算机模拟,特别是在生物流体动力学中。这类问题的IB公式,从最小二乘原理推导而来,既涉及欧拉变量,也涉及拉格朗日变量,由狄拉克-德尔塔函数连接。IB方程的空间离散化基于欧拉变量的固定笛卡尔网格和拉格朗日变量的移动曲线网格。这两种类型的变量通过相互作用方程联系在一起,该方程涉及对狄拉克德尔塔函数的平滑近似。欧拉/拉格朗日恒等式控制数据从一个网格到另一个网格的传输。时间离散化采用二阶龙格-库塔方法。指出了当前和未来的研究方向,并简要讨论了IB方法的应用。

简介:

浸没边界(IB)方法被引入研究心脏瓣膜周围的流动模式,并已发展成为解决流体-结构相互作用问题的一种普遍有用的方法。IB方法既是一种数学公式,也是一种数值格式。数学公式是欧拉和拉格朗日变量的混合。这些与相互作用方程有关,其中狄拉克-德尔塔函数起着重要作用。在IB公式驱动的数值方案中,欧拉变量定义在固定的笛卡尔网格上,拉格朗日变量定义在曲线网格上,该曲线网格在固定的笛卡儿网格中自由移动,而不受任何约束以适应它。数值格式的相互作用方程涉及狄拉克-德尔塔函数的平滑近似,根据我们将讨论的某些原理构建。本文主要研究IB方法的数学结构。这既包括运动方程的IB形式,也包括IB数值格式。尽管为感兴趣的读者提供了对这些主题的参考,但省略了实现的细节,并且仅以概要的方式讨论应用程序。

文二:

 

鞍点问题的数值解法

摘要:

鞍点型大型线性系统在整个计算科学和工程中有着广泛的应用。由于它们的不确定性和通常较差的谱特性,这种线性系统对求解器开发人员来说是一个重大挑战。近年来,人们对鞍点问题的兴趣激增,并为这类系统提出了许多求解技术。本文的目的是提出并讨论鞍点形式线性系统的大量求解方法,重点是求解大型和稀疏问题的迭代方法。

简介:

近年来,大量的工作被投入到解决大型线性系统的鞍点形式的问题。产生这种兴趣的原因是,这些问题出现在各种各样的技术和科学应用中。例如,混合有限元法在流体和固体力学等工程领域的日益普及已经成为鞍点系统(Brezzi 和 Fortin 1991,Elman,Silvester 和 Wathen 2005c)的主要来源。这种兴趣激增的另一个原因是内点算法在线性和非线性优化中的非凡成功,它们的核心需要解决一系列鞍点形式的系统(Nocedal and Wright 1999,Wright 1992,Wright 1997)。

由于鞍点系统无处不在的特性,其数值解的方法和结果已经出现在各种各样的书籍、期刊和会议记录中,证明了对这一主题进行全面的调查是合理的。本文的目的是回顾许多最有前途的解决方法,重点是大型和稀疏问题的迭代方法。虽然许多这些解决方案已经开发了具体的应用(例如,斯托克斯类型的问题在流体动力学) ,它是可能的讨论他们在一个相当一般的设置使用标准的数值线性代数概念,最突出的可能是舒尔补。然而,当选择一个预处理(或开发一个新的) ,知识的起源,特定的问题在手头是必不可少的。虽然没有单一的“最佳”方法存在,但非常有效。我们已经为一些重要类别的问题开发出了解决方法因此,用一些空间来讨论出现的鞍点问题在一些选定的应用程序。希望本次调查能够对从业人员有所帮助正在寻求指导,为自己选择一种解决方法数值线性代数和科学计算的研究人员,特别是研究生,作为这个非常丰富而重要的主题。

文三:

 

序列二次规划

摘要:

自20世纪70年代末推广以来,序列二次规划(SQP)可以说是解决非线性约束优化问题最成功的方法。与大多数优化方法一样,SQP不是一种单一的算法,而是一种概念方法,从中进化出了许多特定的算法。在坚实的理论和计算基础的支持下,商业和公共领域的SQP算法都已被开发并用于解决大量重要的实际问题。最近,已经设计并测试了大规模版本,结果很有希望。

文四:

 

离散力学与变分积分

摘要:

本文综述了基于离散变分原理的有限维机械系统积分算法。变分技术给出了许多辛格式的统一处理,包括高阶格式,以及离散Noether定理的自然处理。该方法还允许我们以自然的方式包括力、耗散和约束。在作为例子的许多具体方案中,提出了Verlet、SHAKE、RATTLE、Newmark和辛划分的Runge–Kutta方案。

简介:

本文综述了离散力学和机械积分器的变分方法。这一观点不仅适用于守恒系统,也适用于强迫系统和耗散系统,因此适用于力学中的控制问题(例如)和传统的守恒问题。正如我们将要展示的那样,变分方法对许多关于离散力学和机械系统积分方法的文献给出了一个全面而统一的观点,我们认为这些都是紧密相关的主题。

该方法自然产生的一些重要主题是辛-能量-动量方法、误差分析、约束、cing和Newmark算法。除了介绍这些方法外,我们还将这些技术与其他最近令人兴奋的发展联系起来,包括用于非线性波动方程和非线性壳动力学等问题的多辛时空积分器(也称为AVI,或异步变分积分器)的PDE设置。事实上,我们的观点之一是,通过从一开始就将积分器建立在有趣的基本力学概念和方法的基础上,可以简化到其他领域,例如连续体力学和具有强迫和约束的系统。

在过去的几年里,这门学科已经发展成为一个非常庞大和活跃的研究领域,有很多观点和主题。我们将在这里集中讨论我们自己的观点,即变分观点。当然,我们必须省略一些重要的主题,但包括我们自己的几个主题。

在本文的最后部分和下面的简史中,我们确实接触了一些但不是全部的其他主题。

正如在标准力学中一样,有些事情从哈密顿的角度来看更容易,而另一些事情从拉格朗日的角度来看则更容易。从这两种观点来看辛积分器可以更深入地了解它们的性质和导数。我们试图在这篇文章中给出一个平衡的观点。

我们将假设配置流形是有限维的。这意味着,在一开始,我们将处理普通微分方程的上下文。然而,正如我们所指出的,我们的方法与变分时空多辛方法密切相关,这是一种适用于无限维PDE上下文的方法,因此对本文方法论的投资可以简化到相应PDE上下文。

文五:

 

反问题: 贝叶斯视角

摘要:

微分方程中的反问题具有巨大的实践意义,也产生了大量的数学和计算创新。通常需要某种形式的正则化来识别不适定行为。在这篇文章中,我们回顾了正则化的贝叶斯方法,发展了一个关于这个主题的函数空间观点。这种方法允许对所有可能的解决方案及其相对概率进行全面表征,同时迫使重要的建模问题以明确和精确的方式得到解决。尽管实现起来很昂贵,但在许多应用领域,这种方法开始处于可用计算资源的范围内。它还允许对不确定性和风险进行量化,而这些应用程序越来越不需要这种量化。此外,该方法在概念上对于理解更简单、计算方便的反问题方法很重要。

我们证明,当以贝叶斯方式表述时,各种反问题都有一个共同的数学框架,我们高度重视由此产生的适定性理论。适定性理论为我们描述的许多稳定性和近似结果提供了基础。我们还回顾了在采用贝叶斯方法求解逆问题时使用的一系列算法方法。其中包括MCMC方法、滤波和变分方法。

简介:

数学科学家面临的一个重大挑战是开发一个连贯的数学和算法框架,使研究人员能够将复杂的数学模型与(通常是庞大的)数据集相结合,这些数据集现在在工程、科学和技术的许多领域都是常见的。在这篇文章中,我们用贝叶斯统计学的语言提出了一系列反问题,主要是由微分方程和数据的结合引起的。这样做的目的有两个:(i)突出从业者在过去几十年中取得重大进展的众多应用领域中产生的共同数学结构,从而促进不同应用领域之间的思想交流;(ii)为问题开发一个抽象的函数空间设置,以评估现有算法的效率,并开发新的算法。应用范围广泛,包括大气科学、海洋学、水文学、地球物理学、化学和生物化学、材料科学、系统生物学、交通流、计量经济学、图像处理和信号处理等领域。

本文中贝叶斯反问题子主题具体发展的指导原则是避免离散化,直到最后可能的时刻。这一原理在整个数值分析中都具有巨大的力量。例如,由于传播速度有限,一阶波动方程不能在任意小的时间内控制到给定的最终状态。然而,一阶波动方程的每一个有限差分空间离散化都会产生一个常微分方程的线性系统,该系统在任何有限时间内都可以控制到给定的最终状态;在离散化之前提出可控性问题是理解的关键(Zuazua 2005)。作为另一个例子,考虑热平衡。如果这是通过θ方法在时间上离散的(其中θ∈[0,1]和θ=0是显式欧拉,θ=1是隐式欧拉),但在空间上不离散,则函数空间上的所得算法仅定义为θ∈[1,1];因此可以推断出,如果θ∈[0,1,则必然存在Courant限制

2)(Richtmyer和Morton 1967),甚至在引入空间离散化之前。在牛顿方法的研究中可以找到另一个例子:在离散化之前,该算法在函数空间上的概念应用,当将其作为非线性微分方程边值问题的迭代方法时,可以产生相当大的见解(Deuflhard 2004)。将解密推迟到算法公式的最后是有益的问题列表,几乎是无穷无尽的。因此,同样的想法在反问题的求解中产生了洞察力,我们在贝叶斯环境中证实了这一想法,这也许并不奇怪。

来源:STEM与计算机方法
非线性化学海洋UM求解技术理论自动驾驶材料数字孪生物流控制人工智能
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-01
最近编辑:9月前
江野
博士 等春风得意,等时间嘉许。
获赞 58粉丝 96文章 153课程 0
点赞
收藏
作者推荐

固体大变形模拟研究的最新进展

文一: 织构和孪晶对镁合金实心杆扭转行为的影响: 与单轴拉伸/压缩比较的晶体塑性方法摘要:采用具有孪晶和去孪晶方案的弹粘塑性自洽模型,结合特定扭转的有限元方法,对镁合金AZ31实心棒在不同初始织构(包括理想基底、轧制、挤压和随机)下的变形行为进行了建模。除了具有几乎均匀变形(拉伸和压缩)的载荷外,还研究了具有非均匀变形的载荷路径,特别是在自由端和固定端扭转的情况下。模拟表明,随着初始织构和加载路径的不同,变形机制发生了显著变化,导致应力-应变响应、Swift效应和变形织构发生变化。特别是,我们的晶体平面模型很好地捕捉到了扭转条件下固体棒内部应力和孪晶体积分数的不均匀分布。这一点尤其重要,因为我们的扭转特定晶体塑性模型提供了一种数值工具,可用于通过实心棒的扭转来调节梯度孪晶结构,以提高镁合金的力学性能。 图:实心扭杆力学模型的示意图和相关模型数量的定义。 图:根据 AZ31镁合金实心棒在(a)拉伸、(b)压缩和自由端扭转下的{0001}极图,在考虑孪生和(d)没有考虑孪生的情况下,在(a)拉伸、(b)压缩的轴向应变 ε = 0.15和自由端扭转(c)的剪切应变 Γ = 0.7。注意到(c)和(d)中的极图同时在 θ-r 平面和 θ-z 平面上出现。 图:在(a)拉/压和自由/固定端扭作用下加载的实心杆试件示意图。镁合金 AZ31B 实心棒的初始(b)理想基础,(c)轧制,(d)挤压,(e)按{0001}和{1010}极图计算的随机织构。文二: 基于显微结构的大应变塑性变形多尺度建模——将全场晶体塑性谱求解器与隐式有限元求解器耦合摘要:我们在隐式有限元(FE)框架中提出了全场晶体塑性模型的全嵌入实现,该组合实现了模拟大应变塑性变形的多尺度方法。在宏观有限元模型的每个积分点,一个基于快速傅立叶变换(FFT)的光谱求解器从下面的全场多晶代表体积单元(RVE)模型中输入均匀化响应,该模型通过使用晶体塑性本构公式求解。在开源软件DAMASK中实现的现象硬化定律和基于位错密度的硬化模型都被用来提供中尺度的本构响应。FE-FFT模型的准确性已经通过FCC多晶体的几种加载场景的单元测试进行了基准测试,包括简单拉伸、简单压缩和简单剪切。应用多尺度模型模拟了四种应用情况,即FCC板的平面应变变形、FCC圆柱体的压缩、HCP棒的四点弯曲和双相钢的梁弯曲。证明了该模型在预测中尺度微观结构演化以及宏观和中尺度力学响应方面的卓越能力。 图:在FEM中嵌入频谱求解器的多尺度模型各阶段的流程图。 图:(a)根据沿 Z 方向的逆极图(IPF)着色的代表性体元素(b)分配给 RVE 的随机纹理(每个方向仅显示一个极点)。 图:单元试验(a)简单拉伸(b)简单压缩(c)简单剪切,显示单元形状和变形后嵌入的 RVE。RVE 根据沿参考 Z 方向的 IPF 着色。 图:用ABAQUS-DAMASK模拟圆形刚性模具纯铝的平面应变变形。(a)最大主应变(b)剪切应变分量LE13的分布。 图:在每个积分点分配了不同RVE的圆柱体的有限元模型。圆柱体的上半部分配有细粒度RVE,而下半部分配有粗粒度RVE。RVE根据IPF沿Z方向进行着色。 图:纯铝平面应变压缩的 ABAQUS-DAMASK 模拟: 显示在 x-z 平面上的 Von Mises 应力分布,以及嵌入 RVE 的最终形状和相应的晶体结构,显示在所指定的六个位置的 x-y 平面上的投影极图。 图:圆形刚性模具纯铝平面应变变形的 ABAQUS-DAMASK 模拟。(a)最大主应变(b)剪应变分量 LE13的分布。文三: 添加制造的316L不锈钢的速率和温度相关塑性:表征、建模和在壳晶格破碎中的应用摘要:对通过选择性激光熔炼获得的添加不锈钢316L的准静态和高应变速率响应进行了数值和实验相结合的研究。实验程序包括单轴拉伸、剪切、缺口拉伸和迷你Nakazima试样的实验,涵盖了广泛的应力状态和应变速率(从10−3到103/s)。建立了一个具有Swift-Voce硬化和Johnson-Cook速率和温度依赖性的各向异性二次塑性模型,以描述不同应力状态和应变速率下组成基体材料的行为。在弹性各向同性壳格结构上进行了低和高加载速度下的压缩实验,以进一步验证已确定的塑性模型在结构应用中的有效性。研究发现,所选择的塑性模型能够很好地描述不同速度和方向加载的光滑壳晶格的反作用力和变形模式。实验表明,无论加载速度如何,添加制造的壳晶格都能够维持60%以上的宏观压缩应变,而不会出现细胞壁的可见断裂。与相同质量的板晶格结构的结果进行比较,阐明了壳晶格的巨大能量吸收潜力。 图:(a)单元网格拓扑与有限元网格的壳格。(b)未经测试的单轴压缩样本(单元格的3 × 3 × 3堆叠)的视图,这些样本是用316L 选择性激光熔化制造的。 图:在两种加载速度下对45°取向的光滑壳晶格进行压缩实验的结果:在(a)低速和(b)高速下,选定应变水平为10%-50%时的变形过程图片。 图:(a) 每个薄板的高度为170mm、宽度为210mm、壁厚为1mm的平行薄板结构的设计。0°方向与打印方向对齐。(b) 从顶部看,打印矩形块的倾斜特写视图。 图:从(c)添加制造的块体的(a)顶部、(b)正面和(d)侧面提取样品的组成材料的EBSD图像。沿印刷方向观察到细长晶粒,平均晶粒尺寸约为20μm,晶粒尺寸至少有一个数量级的扩展。 图:中速下的实验结果文四: 大变形弹塑性的拉格朗日力学和隐式材料点法摘要:物质点法非常适合于建模涉及大变形的问题,而传统的基于网格的方法 会很困难。然而,就该方法的平衡公式而言,总拉格朗日方法和更新拉格朗日方法分别是不合适的和不理想的。这是由于材料点方法的基函数,特别是基函数的导数,通常在未成形的、有时是规则的背景网格上定义。使用材料点处的变形来映射基函数空间导数是可能的,但这引入了额外的算法复杂性和计算费用。本文提出了一种新的拉格朗日平衡表述,它是材料点法的理想表述,因为它在加载步骤开始时满足未变形背景网格上的平衡。该公式是使用准静态隐式算法实现的,该算法包括一致切线的推导,以实现全局平衡迭代的最佳收敛。该方法应用于许多大变形弹塑性问题,特别关注该方法与标准、广义插值和CPDI2材料点法的分析解的收敛性。对于广义插值方法,研究了不同的域更新方法,结果表明,在特定的简单变形场下,所有现有方法都是退化的。提出了一种新的域更新方法来克服这些问题。所提出的材料点法框架可以应用于所有现有的材料点方法,并用于隐式和显式分析,但其优点主要与前者有关。 图:变形框架的全局-局部映射和加载步骤开始时的配置 图:不同变形模式下的域更新过程,其中浅灰色和深灰色阴影区域分别为原始域和更新域。 图:弹性环: 位移场与体力分量。对于身体的力,最大的正值和负值分别是黄色和蓝色。 图:弹塑性坍塌:根据σyy着色的变形材料点位置。 图:弹性梁: 变形材料的点位置,根据 σxx 着色,h = 0.25 m。活动背景网格单元已被阴影灰色。文五: 回顾泰勒冲击试验: 高应变率下金属塑性参数的确定摘要:泰勒冲击试验是表征金属在高应变速率(高达105 s-1)下塑性行为的最简单实验方法。在此,圆柱形杆被冲击在刚性固定砧上,并且通过将实验获得的杆变形与其理论对应物进行比较来估计动态屈服强度。尽管这是一个简单而有力的概念,但到目前为止,缺乏一个全面的理论模型(将材料本构关系与杆的变形模式联系起来)限制了泰勒冲击试验在材料表征中的使用。在本研究中,重新审视了泰勒冲击试验,并试图建立一个有效确定金属本构参数的反框架。逆分析是通过扩展卡尔曼滤波器技术进行的,其中基于Chakraborty等人(2015)中给出的推导来构建正向模型。与通常的做法不同,为了更好地估计材料常数,变形时间历程也包含在本逆公式中。为了进行演示,考虑了五种不同的应力模型,即Johnson–Cook模型、Zerilli–Armstrong模型、Steinberg–Cochran–Guinan–Lund模型、机械阈值应力模型和Preston Tonks Wallace模型。估计的参数与参考值基本一致。结果表明,本次尝试最明显的特点是,一次泰勒检验就足以以合理的精度确定参数 图:确定不同应变率下材料本构参数的实验方法 图:假定泰勒试样的变形形状 图:泰勒冲击试验来源:STEM与计算机方法

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈