We give a study to the algorithm for semi-linear parabolic PDEs in Henry-Labordere (2012) and then generalize it to the non-Markovian case for a class of Backward SDEs (BSDEs). By simulating the branching process, the...
详细信息
We give a study to the algorithm for semi-linear parabolic PDEs in Henry-Labordere (2012) and then generalize it to the non-Markovian case for a class of Backward SDEs (BSDEs). By simulating the branching process, the algorithm does not need any backward regression. To prove that the numerical algorithm converges to the solution of BSDEs, we use the notion of viscosity solution of path dependent PDEs introduced by Ekren et al. (to appear) [5] and extended in Ekren et al. (2012) [6,7]. (C) 2013 Elsevier B.V. All rights reserved.
In this paper, we proposed a novel numerical algorithm for cascaded Raman fiber laser (CRFL) using the approximate analytic results as the initial values for shooting method, which effectively reduces the calculating ...
详细信息
In this paper, we proposed a novel numerical algorithm for cascaded Raman fiber laser (CRFL) using the approximate analytic results as the initial values for shooting method, which effectively reduces the calculating time from several hours to a few minutes. With the algorithm, we obtained a numerical solution of the bilateral-pumping Ge-doped fifth-order CRFL which can avoid the so-called "end face damage" phenomenon efficiently. At the same time, we also simulated the unilateral-pumping one as a comparison, which showed both the characteristics of them are similar with each other. (c) 2006 Elsevier B.V. All rights reserved.
This letter presents an iterative estimation algorithm for modeling a class of output nonlinear systems. The basic idea is to derive an estimation model and to solve an optimization problem using the gradient search. ...
详细信息
This letter presents an iterative estimation algorithm for modeling a class of output nonlinear systems. The basic idea is to derive an estimation model and to solve an optimization problem using the gradient search. The proposed iterative numerical algorithm can estimate the parameters of a class of Wiener nonlinear systems from input-output measurement data. The proposed algorithm has faster convergence rates compared with the stochastic gradient algorithm. The numerical simulation results indicate that the proposed algorithm works well. (C) 2012 Elsevier Ltd. All rights reserved.
Boundary condition plays an important role in prediction of the effective thermal conductivity of composites. In this research, the periodic boundary condition and the representative volume element (RVE) based finite ...
详细信息
Boundary condition plays an important role in prediction of the effective thermal conductivity of composites. In this research, the periodic boundary condition and the representative volume element (RVE) based finite element (FE) homogenization method are adopted to evaluate the effective thermal conductivities of the composites reinforced by the spherical, ellipsoidal and cylindrical inclusions, and the emphases are on the numerical implementation algorithm and validation of the periodic boundary condition. The heat flux continuity of the node pairs on the opposite surfaces of the RVEs of the composites is analyzed and the effective thermal conductivity of the composites are homogenized. The results show that the heat flux continuity of the node pairs on the opposite surfaces of the RVEs of the composites can be guaranteed by the proposed numerical implementation algorithm for the periodic boundary condition, and that the predicted effective thermal conductivities of the composites agree well with those determined by the Lewis-Nielsen model and the experimental tests. Therefore, the RVE based FE homogenization method with the periodic boundary condition can accurately evaluate the effective thermal conductivity of the composites with discontinuous inclusions. (C) 2019 Elsevier Ltd. All rights reserved.
A coupled viscoelastic-viscoplastic (VE-VP) model is implemented and studied. The total strain is the sum of VE and VP parts, and the Cauchy stress is given by a linear YE model as a Boltzmann integral of the history ...
详细信息
A coupled viscoelastic-viscoplastic (VE-VP) model is implemented and studied. The total strain is the sum of VE and VP parts, and the Cauchy stress is given by a linear YE model as a Boltzmann integral of the history of VE strains. The proposed computational algorithm features fully implicit integration, return mapping based on a two-step YE predictor/VP corrector strategy, and a consistent tangent operator. The algorithm is applied to J(2) VP coupled with VE. Very compact expressions are obtained which are form-identical to classical elasto-viscoplasticity (EVP) provided that the constant linear elastic shear and bulk moduli are replaced with incremental relaxation moduli which are appropriate functions of the time increment. Two different integration methods to obtain the incremental moduli are proposed and assessed. Closed-form solutions for uniaxial tension and simple shear are developed, based on an original solution method for integro-differential equations. The analytical results enable to illustrate the constitutive model and provide unambiguous benchmarks for numerical algorithms. Model predictions are compared with experimental data and reasonable correlation is obtained. (C) 2011 Elsevier B.V. All rights reserved.
The method of enforcing sparsity during magnetic resonance imaging reconstruction has been successfully applied to partially parallel imaging (PPI) techniques to reduce noise and artifact levels and hence to achieve e...
详细信息
The method of enforcing sparsity during magnetic resonance imaging reconstruction has been successfully applied to partially parallel imaging (PPI) techniques to reduce noise and artifact levels and hence to achieve even higher acceleration factors. However, there are two major problems in the existing sparsity-constrained PPI techniques: speed and robustness. By introducing an auxiliary variable and decomposing the original minimization problem into two subproblems that are much easier to solve, a fast and robust numerical algorithm for sparsity-constrained PPI technique is developed in this work. The specific implementation for a conventional Cartesian trajectory data set is named self-feeding Sparse Sensitivity Encoding (SENSE). The computational cost for the proposed method is two conventional SENSE reconstructions plus one spatially adaptive image denoising procedure. With reconstruction time approximately doubled, images with a much lower root mean square error (RMSE) can be achieved at high acceleration factors. Using a standard eight-channel head coil, a net acceleration factor of 5 along one dimension can be achieved with low RMSE. Furthermore, the algorithm is insensitive to the choice of parameters. This work improves the clinical applicability of SENSE at high acceleration factors. Magn Reson Med 64:1078-1088, 2010. (C) 2010 WileyLiss, Inc.
In order to solve a class of linear nonlocal boundary value problems, a new reproducing kernel space satisfying nonlocal conditions is constructed carefully. This makes it easy to solve the problems. Furthermore, the ...
详细信息
In order to solve a class of linear nonlocal boundary value problems, a new reproducing kernel space satisfying nonlocal conditions is constructed carefully. This makes it easy to solve the problems. Furthermore, the exact solutions of the problems can be expressed in series form. The numerical results demonstrate that the new method is quite accurate and efficient for solving fourth-order nonlocal boundary value problems. (C) 2010 Elsevier Ltd. All rights reserved.
The ''partly coupled'' algorithm and numerical schemes are developed to solve a coupled model of air-sea-wave interaction, using the matrix and simple pivotal-condensation methods. The ''partly...
详细信息
The ''partly coupled'' algorithm and numerical schemes are developed to solve a coupled model of air-sea-wave interaction, using the matrix and simple pivotal-condensation methods. The ''partly coupled'' algorithm is used to handle boundary conditions at the air-sea interface to allow for solving systems of the atmospheric and oceanic equations, as two separate systems. This ''partly coupled'' algorithm is simpler, more effective, and requires less computing efforts to solve an air-sea-wave coupled model compared to the ''fully coupled'' algorithm developed by Ly [1] for solving a coupled model of the baroclinic atmospheric and oceanic boundary layers. In the numerical schemes, the buoyancy terms of dissipation and TKE budget equations are used as conditions for rearranging equation coefficient vectors to satisfy stability conditions. The mathematical coupled model of the air-sea-wave system is also presented. The model equations are written in the same form for both atmospheric and oceanic boundary layers, and include equations for momentum, the k-epsilon turbulence scheme, and stratification in the atmospheric and oceanic boundary layers. In the model, ocean wind-induced wave effects are modeled by using boundary conditions for the turbulent energy dissipation at the interface in a k-epsilon turbulence closure of the coupled model. numerical experiments are designed for various geostrophic winds, wave heights, and wave ages under neutral atmospheric stratification to demonstrate an example in solving a coupled model of the air-sea-wave system. The numerical simulations show that surface wind-induced waves have strong effects on the basic physical characteristics of the air-sea system.
A new numerical algorithm is presented in this study to solve the partial differential equations that govern fluid flows. This new algorithm is based on first transforming the partial differential equations by introdu...
详细信息
A new numerical algorithm is presented in this study to solve the partial differential equations that govern fluid flows. This new algorithm is based on first transforming the partial differential equations by introducing an exponential function to eliminate the convection terms. A fourth-order central differencing scheme and a second-order central differencing scheme are used to discretize the transformed equations. The algorithm is then applied to simulate fluid flows with exact solutions to validate this new algorithm. The fluid flows used in this study are a self-designed quasi-fluid flow problem, stagnation in plane flow (Hiemenz flow), and flow between two concentric cylinders. Comparisons against the exact solution are made for the results obtained using the new numerical algorithm as well as the power-law scheme. The comparisons indicate that the present fourth-order scheme performs the best and the present second-order scheme is the next most accurate. (c) 2006 Elsevier Inc. All rights reserved.
In this paper, we consider a general class of discrete-time optimal control problems subject to all-time-step constraints on the state and control variables. The derivations of the gradient formulas for the cost and c...
详细信息
In this paper, we consider a general class of discrete-time optimal control problems subject to all-time-step constraints on the state and control variables. The derivations of the gradient formulas for the cost and constraint functions for this constrained discrete-time optimal control problem are rather involved. We present a simple approach to the derivations of these gradient formulas based on reversed automatic differentiation. On this basis, a numerical algorithm is developed to solve this all-time-step constrained discrete-time optimal control problem. We then consider a class of continuous-time optimal control problems subject to continuous state inequality constraints. This constrained continuous-time optimal control problem is discretized into a discrete time optimal control problem with all-time-step constraints using the Euler discretization method. Then, the algorithm developed for constrained discrete time optimal control problem is applied to solve this discretized optimal control problem. numerical examples are presented to verify the applicability of the methods.
暂无评论