Simulation of the microstructure evolution during solidifi cation is greatly benefi cial to the control of solidifi cation microstructures. A phase-fi eld method based on the full threaded tree(FTT) for the simulation...Simulation of the microstructure evolution during solidifi cation is greatly benefi cial to the control of solidifi cation microstructures. A phase-fi eld method based on the full threaded tree(FTT) for the simulation of casting solidifi cation microstructure was proposed in this paper, and the structure of the full threaded tree and the mesh refi nement method was discussed. During dendritic growth in solidifi cation, the mesh for simulation is adaptively refi ned at the liquid-solid interface, and coarsened in other areas. The numerical results of a threedimension dendrite growth indicate that the phase-fi eld method based on FTT is suitable for microstructure simulation. Most importantly, the FTT method can increase the spatial and temporal resolutions beyond the limits imposed by the available hardware compared with the conventional uniform mesh. At the simulation time of 0.03 s in this study, the computer memory used for computation is no more than 10 MB with the FTT method, while it is about 50 MB with the uniform mesh method. In addition, the proposed FTT method is more effi cient in computation time when compared with the uniform mesh method. It would take about 20 h for the uniform mesh method, while only 2 h for the FTT method for computation when the solidifi cation time is 0.17 s in this study.展开更多
A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal...A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.展开更多
Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation o...Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.展开更多
Adaptive mesh refinement (AMR) is fairly practiced in the context of high-dimensional, mesh-based computational models. However, it is in its infancy in that of low-dimensional, generalized-coordinate-based computatio...Adaptive mesh refinement (AMR) is fairly practiced in the context of high-dimensional, mesh-based computational models. However, it is in its infancy in that of low-dimensional, generalized-coordinate-based computational models such as projection-based reduced-order models. This paper presents a complete framework for projection-based model order reduction (PMOR) of nonlinear problems in the presence of AMR that builds on elements from existing methods and augments them with critical new contributions. In particular, it proposes an analytical algorithm for computing a pseudo-meshless inner product between adapted solution snapshots for the purpose of clustering and PMOR. It exploits hyperreduction—specifically, the energy-conserving sampling and weighting hyperreduction method—to deliver for nonlinear and/or parametric problems the desired computational gains. Most importantly, the proposed framework for PMOR in the presence of AMR capitalizes on the concept of state-local reduced-order bases to make the most of the notion of a supermesh, while achieving computational tractability. Its features are illustrated with CFD applications grounded in AMR and its significance is demonstrated by the reported wall-clock speedup factors.展开更多
A finite element method for analysis of pollutant dispersion in shallow water is presented. The analysis is divided into two parts ： （ 1 ） computation of the velocity flow field and water surface elevation, and （2...A finite element method for analysis of pollutant dispersion in shallow water is presented. The analysis is divided into two parts ： （ 1 ） computation of the velocity flow field and water surface elevation, and （2） computation of the pollutant concentration field from the dispersion model. The method was combined with an adaptive meshing technique to increase the solution accuracy, as well as to reduce the computational time and computer memory. The finite element formulation and the computer programs were validated by several examples that have known solutions. In addition, the capability of the combined method was demonstrated by analyzing pollutant dispersion in Chao Phraya River near the gulf of Thailand.展开更多
The wing of a butterfly consists of partially overlapping forewing and hindwing,and forewing sweeping can dynamically change the shape of the whole wing.In this work,the effect of forewing sweeping on aerodynamic perf...The wing of a butterfly consists of partially overlapping forewing and hindwing,and forewing sweeping can dynamically change the shape of the whole wing.In this work,the effect of forewing sweeping on aerodynamic performance of a butterfly like model is studied using a solver based on immersed boundary method and adaptive mesh.For aerodynamic performance,adding a"forward-backward-forward"sweeping motion to the forewing makes it more suitable for fast cruising flight,and compared to the situation without forewing sweeping,the drag is reduced by 46%and the lift to drag ratio is increased by 45%.On the contrary,adding a"backward-for ward-backward"sweeping motion to the forewing increases lift and makes it more suitable for climb flight.For downstroke and middle to late upstroke,the forewing sweeping affects the Leading-Edge Vortex(LEV)through two factors:sweeping velocity and forward sweeping angle,and their effects are coupled.A large forward sweeping velocity can enhance the strength of LEV,while a large forward sweeping angle can weaken it.For early upstroke,the forewing sweeping can affect the wake capture mechanism,sweeping backward can enhance it while sweeping forward can weaken it.The findings in this work provide insight into the design of butterfly like Micro AirVehicles(MAVs).展开更多
Semi-crystalline polymer laser powder bed fusion(L-PBF)has recently attracted increasing interest due to its potential for fabricating complex geometry.However,a more comprehensive understanding of the underlying phys...Semi-crystalline polymer laser powder bed fusion(L-PBF)has recently attracted increasing interest due to its potential for fabricating complex geometry.However,a more comprehensive understanding of the underlying physics during L-PBF is required to better control the properties of the final part.This work proposed a multi-layer numerical model to study the temperature and phase evolution during the polyamide-12(PA12)L-PBF process.The Descend and Parallel Chord methods were introduced to improve the convergence of the non-linear thermal solver.The level-set-based mesh adaptation strategy,governed by multi-physical fields,was applied to alleviate the calculation and accurately track the phase evolution.The processing simulation on the dog-bone model revealed that preheating temperature significantly influences the crystallization behavior.Finally,the multi-layer simulation demonstrated that such a developed numerical model can be used to study the phase transformation during powder layer updating and the cyclic laser sintering phenomena.Moreover,the numerical study suggested that crystallization occurs slowly during the L-PBF process.展开更多
Multi-scale problems in Computational Fluid Dynamics(CFD)often require numerous simulations across various design parameters.Using a fixed mesh for all cases may fail to capture critical physical features.Moving mesh ...Multi-scale problems in Computational Fluid Dynamics(CFD)often require numerous simulations across various design parameters.Using a fixed mesh for all cases may fail to capture critical physical features.Moving mesh adaptation provides an optimal resource allocation to obtain high-resolution flow-fields on low-resolution meshes.However,most existing methods require manual experience and the flow posteriori information poses great challenges to practical applications.In addition,generating adaptive meshes directly from design parameters is difficult due to highly nonlinear relationships.The diffusion model is currently the most popular model in generative tasks that integrates the diffusion principle into deep learning to capture the complex nonlinear correlations.A dual diffusion framework,Para2Mesh,is proposed to predict the adaptive meshes from design parameters by exploiting the robust data distribution learning ability of the diffusion model.Through iterative denoising,the proposed dual networks accurately reconstruct the flow-field to provide flow features as supervised information,and then achieve rapid and reliable mesh movement.Experiments in CFD scenarios demonstrate that Para2Mesh predicts similar meshes directly from design parameters with much higher efficiency than traditional method.It could become a real-time adaptation tool to assist engineering design and optimization,providing a promising solution for high-resolution flow-field analysis.展开更多
The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which co...The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around crack tips and large elements in the other regions. The resulting stress intensity factors and simulated crack propagation behavior are used to evaluate the effectiveness of the procedure. Three sample problems of a center cracked plate, a single edge cracked plate and a compact tension specimen, are simulated and their results assessed.展开更多
In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstructio...In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.展开更多
In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the s...In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes：a） calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b） fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c） adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.展开更多
An adaptive finite element method for high-speed flow-structure interaction is pre- sented.The cell-centered finite element method is combined with an adaptive meshing technique to solve the Navier-Stokes equations fo...An adaptive finite element method for high-speed flow-structure interaction is pre- sented.The cell-centered finite element method is combined with an adaptive meshing technique to solve the Navier-Stokes equations for high-speed compressible flow behavior.The energy equation and the quasi-static structural equations for aerodynamically heated structures are solved by applying the Galerkin finite element method.The finite element formulation and computational procedure are de- scribed.Interactions between the high-speed flow,structural heat transfer,and deformation are studied by two applications of Mach 10 flow over an inclined plate,and Mach 4 flow in a channel.展开更多
Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind sch...Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.展开更多
A self adaptive three-dimensional baroclinic model is designed. A horizontal temperature gradient is used to control the grid size, which can improve computational precision in the fronts without inordinately increasi...A self adaptive three-dimensional baroclinic model is designed. A horizontal temperature gradient is used to control the grid size, which can improve computational precision in the fronts without inordinately increasing computation in the whole area. A simulation of the development and disappearance of the front in the Huanghai Sea is conducted with this model. Simulations of temperature distribution throughout the year are also conducted. The comoutational result agrees well with the observation.展开更多
The adaptive mesh refinement （AMR） method is applied in the 2-D Euler multi-component elasticplastic hydrodynamics code （MEPH2Y）. It is applied on detonation. Firstly, the AMR method is described, including a cons...The adaptive mesh refinement （AMR） method is applied in the 2-D Euler multi-component elasticplastic hydrodynamics code （MEPH2Y）. It is applied on detonation. Firstly, the AMR method is described, including a conservative spatial interpolation, the time integration methodology with the adapitve time increment and an adaptive computational region method. The advantage of AMR technique is exhibited by numerical examples, including the 1-D C-J detonation and the 2-D implosion ignited from a single point. Results show that AMR can promote the computational efficiency, keeping the accuracy in interesting regions.展开更多
The adaptive mesh mesh technique is developed and applied in three dimensional numerical simulation of welding process on the base of the commercial software. Special user subroutine is worked out to accom- plish t...The adaptive mesh mesh technique is developed and applied in three dimensional numerical simulation of welding process on the base of the commercial software. Special user subroutine is worked out to accom- plish this function.This technique can make the dense mesh moving simultaneously with the heat source while the other area of the structure with much coarser mesh, greatly reducing the number of nodes and elements in the analysis.Temperature field,displacement and stress distributions during welding pro- cess me analyzed by FEM method with adaptive mesh and the analysis is also conducted with normal FEM method. The temperature field,displacement and stress distributions obtained with both methods are shown in contrast. The results show that the temperature fields and the displacement distributions of simulation on adaptive mesh correspond wery well with that of without adaptive mesh. Though the stress distributions have some difference,but the trends of the stress distribution are corresponding.The com- parison of the computing time of the two meshes indicates that the adaptive that the adaptive mesh can greatly reduce the calculation time when used for welding process.展开更多
In this paper,a simple and robust shock-capturing method is developed for the Flux Reconstruction(FR)framework by combining the Adaptive Mesh Refinement(AMR)technique with the positivity-preserving property.The adapti...In this paper,a simple and robust shock-capturing method is developed for the Flux Reconstruction(FR)framework by combining the Adaptive Mesh Refinement(AMR)technique with the positivity-preserving property.The adaptive technique avoids the use of redundant meshes in smooth regions,while the positivity-preserving property makes the solver capable of providing numerical solutions with physical meaning.The compatibility of these two significant features relies on a novel limiter designed for mesh refinements.It ensures the positivity of solutions on all newly created cells.Therefore,the proposed method is completely positivity-preserving and thus highly robust.It performs well in solving challenging problems on highly refined meshes and allows the transition of cells at different levels to be completed within a very short distance.The performance of the proposed method is examined in various numerical experiments.When solving Euler equations,the technique of Local Artificial Diffusivity(LAD)is additionally coupled to damp oscillations.More importantly,when solving Navier-Stokes equations,the proposed method requires no auxiliaries and can provide satisfying numerical solutions directly.The implementation of the method becomes rather simple.展开更多
High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical i...High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical issues associated with high-order schemes on an adaptively refined mesh, such as stability and accuracy are addressed. Several CAA benchmark problems are used to demonstrate the feasibility and efficiency of the approach.展开更多
The present study provides a three-dimensional volume-of-fluid method based on the adaptive mesh refinement technique.The projection method on the adaptive mesh is introduced for solving the incompressible Navier-Stok...The present study provides a three-dimensional volume-of-fluid method based on the adaptive mesh refinement technique.The projection method on the adaptive mesh is introduced for solving the incompressible Navier-Stokes equations.The octree structure mesh is employed to solve the flow velocities and the pressure.The developed solver is applied to simulate the deformation of the cubic droplet driven by the surface tension without the effect of the gravity.The numerical results well predict the shape evolution of the droplet.展开更多
An essentially conservative adaptive space time conservation element and solution element （CE/SE） method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its imple...An essentially conservative adaptive space time conservation element and solution element （CE/SE） method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its implementation is based on redefined conservation elements （CEs） and solution elements （SEs）, optimized interpolations and a Courant number insensitive CE/SE scheme. This approach is used in two applications, the Woodward double Mach reflection and a two- component Richtmyer-Meshkov instability experiment. This scheme reveals the essential features of the investigated cases, captures small unstable structures, and yields a solution that is consistent with the results from experiments or other high order methods.展开更多
基金financially supported by Program for New Century Excellent Talents in University(No.NCET-090396)the Foundation for Innovative Research Groups of the Natural Science Foundation of Hubei Province,China(2010CDA067)State Major Science and Technology Special Project Foundation for High-end Numerical Machine and Basic Manufacturing Equipment(2011ZX04014-052,2012ZX04012-011)
文摘Simulation of the microstructure evolution during solidifi cation is greatly benefi cial to the control of solidifi cation microstructures. A phase-fi eld method based on the full threaded tree(FTT) for the simulation of casting solidifi cation microstructure was proposed in this paper, and the structure of the full threaded tree and the mesh refi nement method was discussed. During dendritic growth in solidifi cation, the mesh for simulation is adaptively refi ned at the liquid-solid interface, and coarsened in other areas. The numerical results of a threedimension dendrite growth indicate that the phase-fi eld method based on FTT is suitable for microstructure simulation. Most importantly, the FTT method can increase the spatial and temporal resolutions beyond the limits imposed by the available hardware compared with the conventional uniform mesh. At the simulation time of 0.03 s in this study, the computer memory used for computation is no more than 10 MB with the FTT method, while it is about 50 MB with the uniform mesh method. In addition, the proposed FTT method is more effi cient in computation time when compared with the uniform mesh method. It would take about 20 h for the uniform mesh method, while only 2 h for the FTT method for computation when the solidifi cation time is 0.17 s in this study.
文摘A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.
基金supported by the National Natural Science Foundation of China(Grant Nos.12372200 and 12072242).
文摘Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.
基金support by the Air Force Office of Scientific Research under Grant No.FA9550-20-1-0358 and Grant No.FA9550-22-1-0004.
文摘Adaptive mesh refinement (AMR) is fairly practiced in the context of high-dimensional, mesh-based computational models. However, it is in its infancy in that of low-dimensional, generalized-coordinate-based computational models such as projection-based reduced-order models. This paper presents a complete framework for projection-based model order reduction (PMOR) of nonlinear problems in the presence of AMR that builds on elements from existing methods and augments them with critical new contributions. In particular, it proposes an analytical algorithm for computing a pseudo-meshless inner product between adapted solution snapshots for the purpose of clustering and PMOR. It exploits hyperreduction—specifically, the energy-conserving sampling and weighting hyperreduction method—to deliver for nonlinear and/or parametric problems the desired computational gains. Most importantly, the proposed framework for PMOR in the presence of AMR capitalizes on the concept of state-local reduced-order bases to make the most of the notion of a supermesh, while achieving computational tractability. Its features are illustrated with CFD applications grounded in AMR and its significance is demonstrated by the reported wall-clock speedup factors.
文摘A finite element method for analysis of pollutant dispersion in shallow water is presented. The analysis is divided into two parts ： （ 1 ） computation of the velocity flow field and water surface elevation, and （2） computation of the pollutant concentration field from the dispersion model. The method was combined with an adaptive meshing technique to increase the solution accuracy, as well as to reduce the computational time and computer memory. The finite element formulation and the computer programs were validated by several examples that have known solutions. In addition, the capability of the combined method was demonstrated by analyzing pollutant dispersion in Chao Phraya River near the gulf of Thailand.
基金supported by the National Natural Science Foundation of China(No.52175277)The computations were conducted in Computing Center in Xi'an.The authors thankfully acknowledge these institutions.
文摘The wing of a butterfly consists of partially overlapping forewing and hindwing,and forewing sweeping can dynamically change the shape of the whole wing.In this work,the effect of forewing sweeping on aerodynamic performance of a butterfly like model is studied using a solver based on immersed boundary method and adaptive mesh.For aerodynamic performance,adding a"forward-backward-forward"sweeping motion to the forewing makes it more suitable for fast cruising flight,and compared to the situation without forewing sweeping,the drag is reduced by 46%and the lift to drag ratio is increased by 45%.On the contrary,adding a"backward-for ward-backward"sweeping motion to the forewing increases lift and makes it more suitable for climb flight.For downstroke and middle to late upstroke,the forewing sweeping affects the Leading-Edge Vortex(LEV)through two factors:sweeping velocity and forward sweeping angle,and their effects are coupled.A large forward sweeping velocity can enhance the strength of LEV,while a large forward sweeping angle can weaken it.For early upstroke,the forewing sweeping can affect the wake capture mechanism,sweeping backward can enhance it while sweeping forward can weaken it.The findings in this work provide insight into the design of butterfly like Micro AirVehicles(MAVs).
文摘Semi-crystalline polymer laser powder bed fusion(L-PBF)has recently attracted increasing interest due to its potential for fabricating complex geometry.However,a more comprehensive understanding of the underlying physics during L-PBF is required to better control the properties of the final part.This work proposed a multi-layer numerical model to study the temperature and phase evolution during the polyamide-12(PA12)L-PBF process.The Descend and Parallel Chord methods were introduced to improve the convergence of the non-linear thermal solver.The level-set-based mesh adaptation strategy,governed by multi-physical fields,was applied to alleviate the calculation and accurately track the phase evolution.The processing simulation on the dog-bone model revealed that preheating temperature significantly influences the crystallization behavior.Finally,the multi-layer simulation demonstrated that such a developed numerical model can be used to study the phase transformation during powder layer updating and the cyclic laser sintering phenomena.Moreover,the numerical study suggested that crystallization occurs slowly during the L-PBF process.
基金co-supported by the Aeronautical Science Foundation of China(Nos.2018ZA52002 and 2019ZA052011)。
文摘Multi-scale problems in Computational Fluid Dynamics(CFD)often require numerous simulations across various design parameters.Using a fixed mesh for all cases may fail to capture critical physical features.Moving mesh adaptation provides an optimal resource allocation to obtain high-resolution flow-fields on low-resolution meshes.However,most existing methods require manual experience and the flow posteriori information poses great challenges to practical applications.In addition,generating adaptive meshes directly from design parameters is difficult due to highly nonlinear relationships.The diffusion model is currently the most popular model in generative tasks that integrates the diffusion principle into deep learning to capture the complex nonlinear correlations.A dual diffusion framework,Para2Mesh,is proposed to predict the adaptive meshes from design parameters by exploiting the robust data distribution learning ability of the diffusion model.Through iterative denoising,the proposed dual networks accurately reconstruct the flow-field to provide flow features as supervised information,and then achieve rapid and reliable mesh movement.Experiments in CFD scenarios demonstrate that Para2Mesh predicts similar meshes directly from design parameters with much higher efficiency than traditional method.It could become a real-time adaptation tool to assist engineering design and optimization,providing a promising solution for high-resolution flow-field analysis.
文摘The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around crack tips and large elements in the other regions. The resulting stress intensity factors and simulated crack propagation behavior are used to evaluate the effectiveness of the procedure. Three sample problems of a center cracked plate, a single edge cracked plate and a compact tension specimen, are simulated and their results assessed.
基金supported by the National Natural Science Foundation of China(Nos.91752114 and 11672160)
文摘In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.
文摘In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes：a） calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b） fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c） adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.
基金The project supported by the Thailand Research Fund(TRF)
文摘An adaptive finite element method for high-speed flow-structure interaction is pre- sented.The cell-centered finite element method is combined with an adaptive meshing technique to solve the Navier-Stokes equations for high-speed compressible flow behavior.The energy equation and the quasi-static structural equations for aerodynamically heated structures are solved by applying the Galerkin finite element method.The finite element formulation and computational procedure are de- scribed.Interactions between the high-speed flow,structural heat transfer,and deformation are studied by two applications of Mach 10 flow over an inclined plate,and Mach 4 flow in a channel.
基金supported by the Department of Science & Technology, Government of India under research grant SR/S4/MS:318/06.
文摘Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.
基金The National Science and Technology Support Program of China under contract No.2011BAC03B02the National Natural Science Foundation of China under contract No. 40976002
文摘A self adaptive three-dimensional baroclinic model is designed. A horizontal temperature gradient is used to control the grid size, which can improve computational precision in the fronts without inordinately increasing computation in the whole area. A simulation of the development and disappearance of the front in the Huanghai Sea is conducted with this model. Simulations of temperature distribution throughout the year are also conducted. The comoutational result agrees well with the observation.
基金Sponsored by the National Natural Science Foundation of China(10676120)Laboratory of Computational Physics Foundation(9140C690101070C69)
文摘The adaptive mesh refinement （AMR） method is applied in the 2-D Euler multi-component elasticplastic hydrodynamics code （MEPH2Y）. It is applied on detonation. Firstly, the AMR method is described, including a conservative spatial interpolation, the time integration methodology with the adapitve time increment and an adaptive computational region method. The advantage of AMR technique is exhibited by numerical examples, including the 1-D C-J detonation and the 2-D implosion ignited from a single point. Results show that AMR can promote the computational efficiency, keeping the accuracy in interesting regions.
文摘The adaptive mesh mesh technique is developed and applied in three dimensional numerical simulation of welding process on the base of the commercial software. Special user subroutine is worked out to accom- plish this function.This technique can make the dense mesh moving simultaneously with the heat source while the other area of the structure with much coarser mesh, greatly reducing the number of nodes and elements in the analysis.Temperature field,displacement and stress distributions during welding pro- cess me analyzed by FEM method with adaptive mesh and the analysis is also conducted with normal FEM method. The temperature field,displacement and stress distributions obtained with both methods are shown in contrast. The results show that the temperature fields and the displacement distributions of simulation on adaptive mesh correspond wery well with that of without adaptive mesh. Though the stress distributions have some difference,but the trends of the stress distribution are corresponding.The com- parison of the computing time of the two meshes indicates that the adaptive that the adaptive mesh can greatly reduce the calculation time when used for welding process.
基金supported by the National Natural Science Foundation of China(No.11721202).
文摘In this paper,a simple and robust shock-capturing method is developed for the Flux Reconstruction(FR)framework by combining the Adaptive Mesh Refinement(AMR)technique with the positivity-preserving property.The adaptive technique avoids the use of redundant meshes in smooth regions,while the positivity-preserving property makes the solver capable of providing numerical solutions with physical meaning.The compatibility of these two significant features relies on a novel limiter designed for mesh refinements.It ensures the positivity of solutions on all newly created cells.Therefore,the proposed method is completely positivity-preserving and thus highly robust.It performs well in solving challenging problems on highly refined meshes and allows the transition of cells at different levels to be completed within a very short distance.The performance of the proposed method is examined in various numerical experiments.When solving Euler equations,the technique of Local Artificial Diffusivity(LAD)is additionally coupled to damp oscillations.More importantly,when solving Navier-Stokes equations,the proposed method requires no auxiliaries and can provide satisfying numerical solutions directly.The implementation of the method becomes rather simple.
基金supported by the National Natural Science Foundation of China (11150110134)the Science Foundation of Aeronautics of China (20101271004)
文摘High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical issues associated with high-order schemes on an adaptively refined mesh, such as stability and accuracy are addressed. Several CAA benchmark problems are used to demonstrate the feasibility and efficiency of the approach.
基金This work was supported by the National Natural Science Foun-dation of China(No.41776194).
文摘The present study provides a three-dimensional volume-of-fluid method based on the adaptive mesh refinement technique.The projection method on the adaptive mesh is introduced for solving the incompressible Navier-Stokes equations.The octree structure mesh is employed to solve the flow velocities and the pressure.The developed solver is applied to simulate the deformation of the cubic droplet driven by the surface tension without the effect of the gravity.The numerical results well predict the shape evolution of the droplet.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.10732010,10972010,and 11028206)the Opening Project of State Key Laboratory of Explosion Science and Technology,China (Grant No.KFJJ13-5M)
文摘An essentially conservative adaptive space time conservation element and solution element （CE/SE） method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its implementation is based on redefined conservation elements （CEs） and solution elements （SEs）, optimized interpolations and a Courant number insensitive CE/SE scheme. This approach is used in two applications, the Woodward double Mach reflection and a two- component Richtmyer-Meshkov instability experiment. This scheme reveals the essential features of the investigated cases, captures small unstable structures, and yields a solution that is consistent with the results from experiments or other high order methods.
