In this paper,a number of ordinary differential equation(ODE)conversion techniques for trans- formation of nonstandard ODE boundary value problems into standard forms are summarised,together with their applications to...In this paper,a number of ordinary differential equation(ODE)conversion techniques for trans- formation of nonstandard ODE boundary value problems into standard forms are summarised,together with their applications to a variety of boundary value problems in computational solid mechanics,such as eigenvalue problem,geometrical and material nonlinear problem,elastic contact problem and optimal design problems through some simple and representative examples,The advantage of such approach is that various ODE bounda- ry value problems in computational mechanics can be solved effectively in a unified manner by invoking a stand- ard ODE solver.展开更多
A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the gov...A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.展开更多
In this work we provide analytic and numerical solutions for the Bjorken flow,a standard benchmark in relativistic hydrodynamics providing a simple model for the bulk evolution of matter created in collisions between ...In this work we provide analytic and numerical solutions for the Bjorken flow,a standard benchmark in relativistic hydrodynamics providing a simple model for the bulk evolution of matter created in collisions between heavy nuclei.We consider relativistic gases of both massive and massless particles,working in a(2+1)and(3+1)Minkowski space-time coordinate system.The numerical results from a recently developed lattice kinetic scheme show excellent agreement with the analytic solutions.展开更多
The finite element method is a key player in computational electromag-netics for designing RF(Radio Frequency)components such as waveguides.The frequency-domain analysis is fundamental to identify the characteristics ...The finite element method is a key player in computational electromag-netics for designing RF(Radio Frequency)components such as waveguides.The frequency-domain analysis is fundamental to identify the characteristics of the components.For the conventional frequency-domain electromagnetic analysis using FEM(Finite Element Method),the system matrix is complex-numbered as well as indefinite.The iterative solvers can be faster than the direct solver when the solver convergence is guaranteed and done in a few steps.However,such complex-numbered and indefinite systems are hard to exploit the merit of the iterative solver.It is also hard to benefit from matrix factorization techniques due to varying system matrix parts according to frequency.Overall,it is hard to adopt conventional iterative solvers even though the system matrix is sparse.A new parallel iterative FEM solver for frequency domain analysis is implemented for inhomogeneous waveguide structures in this paper.In this implementation,the previous solution of the iterative solver of Matlab(Matrix Laboratory)employ-ing the preconditioner is used for the initial guess for the next step’s solution process.The overlapped parallel stage using Matlab’s Parallel Computing Toolbox is also proposed to alleviate the cold starting,which ruins the convergence of early steps in each parallel stage.Numerical experiments based on waveguide structures have demonstrated the accuracy and efficiency of the proposed scheme.展开更多
With increasing engineering demands,there need high order accurate schemes embedded with precise physical information in order to capture delicate small scale structures and strong waves with correct“physics”.There ...With increasing engineering demands,there need high order accurate schemes embedded with precise physical information in order to capture delicate small scale structures and strong waves with correct“physics”.There are two families of high order methods:One is the method of line,relying on the Runge-Kutta(R-K)time-stepping.The building block is the Riemann solution labeled as the solution element“1”.Each step in R-K just has first order accuracy.In order to derive a fourth order accuracy scheme in time,one needs four stages labeled as“1111=4”.The other is the one-stage Lax-Wendroff(LW)type method,which is more compact but is complicated to design numerical fluxes and hard to use when applied to highly nonlinear problems.In recent years,the pair of solution element and dynamics element,labeled as“2”,are taken as the building block.The direct adoption of the dynamics implies the inherent temporal-spatial coupling.With this type of building blocks,a family of two-stage fourth order accurate schemes,labeled as“22=4”,are designed for the computation of compressible fluid flows.The resulting schemes are compact,robust and efficient.This paper contributes to elucidate how and why high order accurate schemes should be so designed.To some extent,the“22=4”algorithm extracts the advantages of the method of line and one-stage LW method.As a core part,the pair“2”is expounded and LW solver is revisited.The generalized Riemann problem(GRP)solver,as the discontinuous and nonlinear version of LW flow solver,and the gas kinetic scheme(GKS)solver,the microscopic LW solver,are all reviewed.The compact Hermite-type data reconstruction and high order approximation of boundary conditions are proposed.Besides,the computational performance and prospective discussions are presented.展开更多
In this paper, we apply arbitrary Riemann solvers, which may not satisfy the Maire＇s requirement, to the Maire＇s node-based Lagrangian scheme developed in [P. H. Maire et al., SIAM J. Sci. Comput, 29 （2007）, 1781-...In this paper, we apply arbitrary Riemann solvers, which may not satisfy the Maire＇s requirement, to the Maire＇s node-based Lagrangian scheme developed in [P. H. Maire et al., SIAM J. Sci. Comput, 29 （2007）, 1781-1824]. In particular, we apply the so-called Multi-Fluid Channel on Averaged Volume （MFCAV） Riemann solver and a Riemann solver that adaptively combines the MFCAV solver with other more dissipative Riemann solvers to the Maire＇s scheme. It is noted that neither of the two solvers satisfies the Maire＇s requirement. Numerical experiments are presented to demonstrate that the application of the two Riemann solvers is successful.展开更多
In this paper, a S-CR method with inexact solvers on the subdomains is presented, and then its convergence property is proved under very general conditions. This result is important because it guarantees the effective...In this paper, a S-CR method with inexact solvers on the subdomains is presented, and then its convergence property is proved under very general conditions. This result is important because it guarantees the effectiveness of the Schwarz alternating method when executed on message-passing distributed memory multiprocessor system.展开更多
An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction...An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction procedure of the original system of PDEs describing coupled thermomechanical behavior to a system of Differential Algebraic Equations （DAEs） where the stress-strain relationships are treated as algebraic equations. The resulting system of DAEs was then solved with a Backward Differentiation Formula （BDF） using a fully implicit algorithm. The described procedure was explained in detail, and its effectiveness was demonstrated on the solution of a transient uncoupled thermoelastic problem, for which an analytical solution is known, as well as on a fully coupled problem in the two-dimensional case.展开更多
In this work, a numerical solution of the incompressible Stokes equations is proposed. The method suggested is based on an algorithm of discretization by the unstable of Q1 – P0 velocity/pressure ?nite element approx...In this work, a numerical solution of the incompressible Stokes equations is proposed. The method suggested is based on an algorithm of discretization by the unstable of Q1 – P0 velocity/pressure ?nite element approximation. It is shown that the inf-sup stability constant is O(h) in two dimensions and O( h2) in three dimensions. The basic tool in the analysis is the method of modi?ed equations which is applied to ?nite difference representations of the underlying ?nite element equations. In order to evaluate the performance of the method, the numerical results are compared with some previously published works or with others coming from commercial code like Adina system.展开更多
Solving for detailed chemical kinetics remains one of the major bottlenecks for computational fluid dynamics simulations of reacting flows using a finite-rate-chemistry approach.This has motivated the use of neural ne...Solving for detailed chemical kinetics remains one of the major bottlenecks for computational fluid dynamics simulations of reacting flows using a finite-rate-chemistry approach.This has motivated the use of neural networks to predict stiff chemical source terms as functions of the thermochemical state of the combustion system.However,due to the nonlinearities and multi-scale nature of combustion,the predicted solution often diverges from the true solution when these machine learning models are coupled with a computational fluid dynamics solver.This is because these approaches minimize the error during training without guaranteeing successful integration with ordinary differential equation solvers.In the present work,a novel neural ordinary differential equations approach to modeling chemical kinetics,termed as ChemNODE,is developed.In this machine learning framework,the chemical source terms predicted by the neural networks are integrated during training,and by computing the required derivatives,the neural network weights are adjusted accordingly to minimize the difference between the predicted and ground-truth solution.A proof-of-concept study is performed with ChemNODE for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions.It is shown that ChemNODE accurately captures the chemical kinetic behavior and reproduces the results obtained using the detailed kinetic mechanism at a fraction of the computational cost.展开更多
The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is base...The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is based on the approach of Vila[25]and another is proposed in this article which mimics the classical 1D LaxWendroff scheme.In both approaches the classical SPH artificial viscosity term is removed preserving nevertheless the linear stability of the methods,demonstrated via the von Neumann stability analysis.Moreover,the issue of the consistency for the equations of gas dynamics is analyzed.An alternative approach is proposed that consists of using Godunov-type SPH schemes in Lagrangian coordinates.This not only provides an improvement in accuracy of the numerical solutions,but also assures that the consistency condition on the gradient of the kernel function is satisfied using an equidistant distribution of particles in Lagrangian mass coordinates.Three different Riemann solvers are implemented for the first-order Godunov type SPH schemes in Lagrangian coordinates,namely the Godunov flux based on the exact Riemann solver,the Rusanov flux and a new modified Roe flux,following the work of Munz[17].Some well-known numerical 1D shock tube test cases[22]are solved,comparing the numerical solutions of the Godunov-type SPH schemes in Lagrangian coordinates with the first-order Godunov finite volume method in Eulerian coordinates and the standard SPH scheme with Monaghan’s viscosity term.展开更多
An indirect-forcing immersed boundary method for solving the incompressible Navier-Stokes equations involving the interfaces and irregular domains is developed.The rigid boundaries and interfaces are represented by a ...An indirect-forcing immersed boundary method for solving the incompressible Navier-Stokes equations involving the interfaces and irregular domains is developed.The rigid boundaries and interfaces are represented by a number of Lagrangian control points.Stationary rigid boundaries are embedded in the Cartesian grid and singular forces at the rigid boundaries are applied to impose the prescribed velocity conditions.The singular forces at the interfaces and the rigid boundaries are then distributed to the nearby Cartesian grid points using the immersed boundary method.In the present work,the singular forces at the rigid boundaries are computed implicitly by solving a small system of equations at each time step to ensure that the prescribed velocity condition at the rigid boundary is satisfied exactly.For deformable interfaces,the forces that the interface exerts on the fluid are computed from the configuration of the elastic interface and are applied to the fluid.The Navier-Stokes equations are discretized using finite difference method on a staggered uniform Cartesian grid by a second order accurate projection method.The ability of the method to simulate viscous flows with interfaces on irregular domains is demonstrated by applying to the rotational flow problem,the relaxation of an elastic membrane and flow in a constriction with an immersed elastic membrane.展开更多
New direct spectral solvers for the 3D Helmholtz equation in a finite cylindrical region are presented.A purely variational(no collocation)formulation of the problem is adopted,based on Fourier series expansion of the...New direct spectral solvers for the 3D Helmholtz equation in a finite cylindrical region are presented.A purely variational(no collocation)formulation of the problem is adopted,based on Fourier series expansion of the angular dependence and Legendre polynomials for the axial dependence.A new Jacobi basis is proposed for the radial direction overcoming the main disadvantages of previously developed bases for the Dirichlet problem.Nonhomogeneous Dirichlet boundary conditions are enforced by a discrete lifting and the vector problem is solved by means of a classical uncoupling technique.In the considered formulation,boundary conditions on the axis of the cylindrical domain are never mentioned,by construction.The solution algorithms for the scalar equations are based on double diagonalization along the radial and axial directions.The spectral accuracy of the proposed algorithms is verified by numerical tests.展开更多
This paper summarizes suitable material models for creep and damage of concrete which are coupled with heat and moisture transfer.The fully coupled approach or the staggered coupling is assumed.Governing equations are...This paper summarizes suitable material models for creep and damage of concrete which are coupled with heat and moisture transfer.The fully coupled approach or the staggered coupling is assumed.Governing equations are spatially discretized by the finite element method and the temporal discretization is done by the generalized trapezoidal method.Systems of non-linear algebraic equations are solved by the Newton method.Development of an efficient and extensible computer code based on the C++programming language is described.Finally,successful analyses of two real engineering problems are described.展开更多
Most traditional artificial intelligence(AI)systems of the past decades are either very limited,or based on heuristics,or both.The new millennium,however,has brought substantial progress in the field of theoretically ...Most traditional artificial intelligence(AI)systems of the past decades are either very limited,or based on heuristics,or both.The new millennium,however,has brought substantial progress in the field of theoretically optimal and practically feasible algorithms for prediction,search,inductive inference based on Occam’s razor,problem solving,decision making,and reinforcement learning in environments of a very general type.Since inductive inference is at the heart of all inductive sciences,some of the results are relevant not only for AI and computer science but also for physics,provoking nontraditional predictions based on Zuse’s thesis of the computer-generated universe.We first briefly review the history of AI since Godel’s 1931 paper,then discuss recent post-2000 approaches that are currently transforming general AI research into a formal science.展开更多
基金The project is supported by National Natural Science Foundation of China
文摘In this paper,a number of ordinary differential equation(ODE)conversion techniques for trans- formation of nonstandard ODE boundary value problems into standard forms are summarised,together with their applications to a variety of boundary value problems in computational solid mechanics,such as eigenvalue problem,geometrical and material nonlinear problem,elastic contact problem and optimal design problems through some simple and representative examples,The advantage of such approach is that various ODE bounda- ry value problems in computational mechanics can be solved effectively in a unified manner by invoking a stand- ard ODE solver.
基金Supported by the National Natural Science Foundation of China(11272153)
文摘A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.
基金funding from the European Research Council under the European Union’sHorizon 2020 framework programme(No.P/2014-2020)/ERC Grant Agreement No.739964(COPMAT).
文摘In this work we provide analytic and numerical solutions for the Bjorken flow,a standard benchmark in relativistic hydrodynamics providing a simple model for the bulk evolution of matter created in collisions between heavy nuclei.We consider relativistic gases of both massive and massless particles,working in a(2+1)and(3+1)Minkowski space-time coordinate system.The numerical results from a recently developed lattice kinetic scheme show excellent agreement with the analytic solutions.
基金supported by Institute of Information&communications Technology Planning&Evaluation(ITP)grant funded by the Korea govermment(MSIT)(No.2019-0-00098,Advanced and Integrated Software Development for Electromagnetic Analysis)supported by Research Assistance Program(2021)in the Incheon National University.
文摘The finite element method is a key player in computational electromag-netics for designing RF(Radio Frequency)components such as waveguides.The frequency-domain analysis is fundamental to identify the characteristics of the components.For the conventional frequency-domain electromagnetic analysis using FEM(Finite Element Method),the system matrix is complex-numbered as well as indefinite.The iterative solvers can be faster than the direct solver when the solver convergence is guaranteed and done in a few steps.However,such complex-numbered and indefinite systems are hard to exploit the merit of the iterative solver.It is also hard to benefit from matrix factorization techniques due to varying system matrix parts according to frequency.Overall,it is hard to adopt conventional iterative solvers even though the system matrix is sparse.A new parallel iterative FEM solver for frequency domain analysis is implemented for inhomogeneous waveguide structures in this paper.In this implementation,the previous solution of the iterative solver of Matlab(Matrix Laboratory)employ-ing the preconditioner is used for the initial guess for the next step’s solution process.The overlapped parallel stage using Matlab’s Parallel Computing Toolbox is also proposed to alleviate the cold starting,which ruins the convergence of early steps in each parallel stage.Numerical experiments based on waveguide structures have demonstrated the accuracy and efficiency of the proposed scheme.
基金This work is supported by NSFC(nos.11771054,91852207)and Foundation of LCP.
文摘With increasing engineering demands,there need high order accurate schemes embedded with precise physical information in order to capture delicate small scale structures and strong waves with correct“physics”.There are two families of high order methods:One is the method of line,relying on the Runge-Kutta(R-K)time-stepping.The building block is the Riemann solution labeled as the solution element“1”.Each step in R-K just has first order accuracy.In order to derive a fourth order accuracy scheme in time,one needs four stages labeled as“1111=4”.The other is the one-stage Lax-Wendroff(LW)type method,which is more compact but is complicated to design numerical fluxes and hard to use when applied to highly nonlinear problems.In recent years,the pair of solution element and dynamics element,labeled as“2”,are taken as the building block.The direct adoption of the dynamics implies the inherent temporal-spatial coupling.With this type of building blocks,a family of two-stage fourth order accurate schemes,labeled as“22=4”,are designed for the computation of compressible fluid flows.The resulting schemes are compact,robust and efficient.This paper contributes to elucidate how and why high order accurate schemes should be so designed.To some extent,the“22=4”algorithm extracts the advantages of the method of line and one-stage LW method.As a core part,the pair“2”is expounded and LW solver is revisited.The generalized Riemann problem(GRP)solver,as the discontinuous and nonlinear version of LW flow solver,and the gas kinetic scheme(GKS)solver,the microscopic LW solver,are all reviewed.The compact Hermite-type data reconstruction and high order approximation of boundary conditions are proposed.Besides,the computational performance and prospective discussions are presented.
文摘In this paper, we apply arbitrary Riemann solvers, which may not satisfy the Maire＇s requirement, to the Maire＇s node-based Lagrangian scheme developed in [P. H. Maire et al., SIAM J. Sci. Comput, 29 （2007）, 1781-1824]. In particular, we apply the so-called Multi-Fluid Channel on Averaged Volume （MFCAV） Riemann solver and a Riemann solver that adaptively combines the MFCAV solver with other more dissipative Riemann solvers to the Maire＇s scheme. It is noted that neither of the two solvers satisfies the Maire＇s requirement. Numerical experiments are presented to demonstrate that the application of the two Riemann solvers is successful.
文摘In this paper, a S-CR method with inexact solvers on the subdomains is presented, and then its convergence property is proved under very general conditions. This result is important because it guarantees the effectiveness of the Schwarz alternating method when executed on message-passing distributed memory multiprocessor system.
文摘An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction procedure of the original system of PDEs describing coupled thermomechanical behavior to a system of Differential Algebraic Equations （DAEs） where the stress-strain relationships are treated as algebraic equations. The resulting system of DAEs was then solved with a Backward Differentiation Formula （BDF） using a fully implicit algorithm. The described procedure was explained in detail, and its effectiveness was demonstrated on the solution of a transient uncoupled thermoelastic problem, for which an analytical solution is known, as well as on a fully coupled problem in the two-dimensional case.
文摘In this work, a numerical solution of the incompressible Stokes equations is proposed. The method suggested is based on an algorithm of discretization by the unstable of Q1 – P0 velocity/pressure ?nite element approximation. It is shown that the inf-sup stability constant is O(h) in two dimensions and O( h2) in three dimensions. The basic tool in the analysis is the method of modi?ed equations which is applied to ?nite difference representations of the underlying ?nite element equations. In order to evaluate the performance of the method, the numerical results are compared with some previously published works or with others coming from commercial code like Adina system.
基金This work was supported by the U.S.Department of Energy,Office of Science under contract DE-AC02-06CH11357The research work was funded by Argonne’s Laboratory Directed Research and Development(LDRD)Innovate project#2020-0203.The authors acknowledge the computing resources available via Bebop,a high-performance computing cluster operated by the Laboratory Computing Resource Center(LCRC)at Argonne National Laboratory.
文摘Solving for detailed chemical kinetics remains one of the major bottlenecks for computational fluid dynamics simulations of reacting flows using a finite-rate-chemistry approach.This has motivated the use of neural networks to predict stiff chemical source terms as functions of the thermochemical state of the combustion system.However,due to the nonlinearities and multi-scale nature of combustion,the predicted solution often diverges from the true solution when these machine learning models are coupled with a computational fluid dynamics solver.This is because these approaches minimize the error during training without guaranteeing successful integration with ordinary differential equation solvers.In the present work,a novel neural ordinary differential equations approach to modeling chemical kinetics,termed as ChemNODE,is developed.In this machine learning framework,the chemical source terms predicted by the neural networks are integrated during training,and by computing the required derivatives,the neural network weights are adjusted accordingly to minimize the difference between the predicted and ground-truth solution.A proof-of-concept study is performed with ChemNODE for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions.It is shown that ChemNODE accurately captures the chemical kinetic behavior and reproduces the results obtained using the detailed kinetic mechanism at a fraction of the computational cost.
文摘The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach.First,we focus on the problem of the stability for two different SPH schemes,one is based on the approach of Vila[25]and another is proposed in this article which mimics the classical 1D LaxWendroff scheme.In both approaches the classical SPH artificial viscosity term is removed preserving nevertheless the linear stability of the methods,demonstrated via the von Neumann stability analysis.Moreover,the issue of the consistency for the equations of gas dynamics is analyzed.An alternative approach is proposed that consists of using Godunov-type SPH schemes in Lagrangian coordinates.This not only provides an improvement in accuracy of the numerical solutions,but also assures that the consistency condition on the gradient of the kernel function is satisfied using an equidistant distribution of particles in Lagrangian mass coordinates.Three different Riemann solvers are implemented for the first-order Godunov type SPH schemes in Lagrangian coordinates,namely the Godunov flux based on the exact Riemann solver,the Rusanov flux and a new modified Roe flux,following the work of Munz[17].Some well-known numerical 1D shock tube test cases[22]are solved,comparing the numerical solutions of the Godunov-type SPH schemes in Lagrangian coordinates with the first-order Godunov finite volume method in Eulerian coordinates and the standard SPH scheme with Monaghan’s viscosity term.
基金The last author’s research is supported by the grant AcRF RG59/08 M52110092.
文摘An indirect-forcing immersed boundary method for solving the incompressible Navier-Stokes equations involving the interfaces and irregular domains is developed.The rigid boundaries and interfaces are represented by a number of Lagrangian control points.Stationary rigid boundaries are embedded in the Cartesian grid and singular forces at the rigid boundaries are applied to impose the prescribed velocity conditions.The singular forces at the interfaces and the rigid boundaries are then distributed to the nearby Cartesian grid points using the immersed boundary method.In the present work,the singular forces at the rigid boundaries are computed implicitly by solving a small system of equations at each time step to ensure that the prescribed velocity condition at the rigid boundary is satisfied exactly.For deformable interfaces,the forces that the interface exerts on the fluid are computed from the configuration of the elastic interface and are applied to the fluid.The Navier-Stokes equations are discretized using finite difference method on a staggered uniform Cartesian grid by a second order accurate projection method.The ability of the method to simulate viscous flows with interfaces on irregular domains is demonstrated by applying to the rotational flow problem,the relaxation of an elastic membrane and flow in a constriction with an immersed elastic membrane.
文摘New direct spectral solvers for the 3D Helmholtz equation in a finite cylindrical region are presented.A purely variational(no collocation)formulation of the problem is adopted,based on Fourier series expansion of the angular dependence and Legendre polynomials for the axial dependence.A new Jacobi basis is proposed for the radial direction overcoming the main disadvantages of previously developed bases for the Dirichlet problem.Nonhomogeneous Dirichlet boundary conditions are enforced by a discrete lifting and the vector problem is solved by means of a classical uncoupling technique.In the considered formulation,boundary conditions on the axis of the cylindrical domain are never mentioned,by construction.The solution algorithms for the scalar equations are based on double diagonalization along the radial and axial directions.The spectral accuracy of the proposed algorithms is verified by numerical tests.
基金Financial support for this work was provided by the Czech Science Foundation,projects n.103/08/1119 and 105/10/1682support was provided by the Ministry of Education,Youth and Sports of the Czech Republic,project No.1M0579.
文摘This paper summarizes suitable material models for creep and damage of concrete which are coupled with heat and moisture transfer.The fully coupled approach or the staggered coupling is assumed.Governing equations are spatially discretized by the finite element method and the temporal discretization is done by the generalized trapezoidal method.Systems of non-linear algebraic equations are solved by the Newton method.Development of an efficient and extensible computer code based on the C++programming language is described.Finally,successful analyses of two real engineering problems are described.
文摘Most traditional artificial intelligence(AI)systems of the past decades are either very limited,or based on heuristics,or both.The new millennium,however,has brought substantial progress in the field of theoretically optimal and practically feasible algorithms for prediction,search,inductive inference based on Occam’s razor,problem solving,decision making,and reinforcement learning in environments of a very general type.Since inductive inference is at the heart of all inductive sciences,some of the results are relevant not only for AI and computer science but also for physics,provoking nontraditional predictions based on Zuse’s thesis of the computer-generated universe.We first briefly review the history of AI since Godel’s 1931 paper,then discuss recent post-2000 approaches that are currently transforming general AI research into a formal science.
