Questions tagged [numerical-methods]

Questions concerning problems which cannot be exactly solved, or where the exact solution may be much more difficult to acquire than a solution by numerical methods.

21 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
2
votes
2answers
54 views

How to convert the velocity field to a pressure field for obtaining the sound level?

I am studying jet mixing noise and I am given a velocity field, with $V_x$, $V_y$, $V_z$ in the Cartesian coordinates. I wanted to know if there is a method by which I could find some information ...
2
votes
0answers
34 views

Convergence problem in RANS simulation of unsteady flow

I am simulating an air flow through a two dimensional rectangular channel. The channel splits into two channels, one of them is followed by a 90 degree bend downwards and leads into a cylindrical ...
2
votes
0answers
625 views

How does Abaqus calculate reaction force?

I wrote a UEL for one element and compared to a reference UEL (FEM), where the results agreed satisfactorily, except the reaction force. The stress, strain, strain energy all were pretty accurate. I ...
2
votes
0answers
37 views

Discretise this transport equation

I am trying to discretise the PDE: $$\phi \frac{\partial c}{\partial t}+\frac{\partial j}{\partial x}=0$$ where $c$ is a function of $x$ and $t$, and $j=qc-D(q)\frac{\partial c}{\partial x}$, $q$ is ...
1
vote
0answers
35 views

Does numerical dissipation in FE simulations affect system's elastic energy? Does it depend on thickness?

I am running a finite-element simulation in Abaqus, for which I also have a theoretical model. To be precise, I am simulating wrinkling in a thin annular sheet, with shell elements. The model makes a ...
1
vote
1answer
34 views

Why is the amount of force/moment transfer occurring between bodies in FEA is independent of mesh size?

So I read somewhere that the mesh size won't effect the amount of force transfer occurring between bodies if several bodies are in contact, and a global force is applied somewhere. I wanted to know ...
1
vote
0answers
42 views

Errors obtained in the numerical solution of the 1D convection-diffusion equation

I have been trying to solve the 1D convection-diffusion equation In implementing the above, I have tried to follow closely the method described here: https://uk.mathworks.com/support/search.html/...
1
vote
1answer
52 views

Finite difference discretization of the Cauchy-Riemann PDEs

I made a forward fd-discretization of the Cauchy-Riemann PDEs but I am struggling to implement this in python. I have a quadratic mesh with heighτ = $2*\pi$. The dirichlet boundary conditions are at $...
1
vote
1answer
80 views

Analytical solution of thermal stress in a bi-metallic cylinder

As a side project I have looked to solve using Python an analytical solution for one-dimensional transient heat transfer in a thick cylinder with cladding on the inside surface, wherein the fluid in ...
1
vote
0answers
24 views

Water distribution through N x Needle orifices

Consider a water distribution setup consisting of a feeder line blocked at one end but having 100 needle outlets all the same size. The feeder line cross sectional area is much greater than a needle ...
1
vote
0answers
34 views

How to formulate a 3D version of the Navier-Stokes equations in 2D for numerical solution?

I'm considering a topology optimization problem where the design domain is a rotating annulus (in $\mathbb{R}^2$). One of the adjoint equations of the Navier-Stokes system whose solution is to be used ...
1
vote
0answers
15 views

Does excitation matter in a time-domain simulation/computation of a transfer function of an LTI system?

Let's say I am running a FDTD simulation of a wave-equation to determine a transfer function of an LTI system: \begin{equation} H(f) = \dfrac{Y(f)}{X(f)}\ \end{equation} where $Y(f)$ and $X(f)$ are ...
1
vote
0answers
2k views

How to calculate effective strain in plane stress?

What is the correct formula to calculate the effective strain under plane stress conditions? In plane stress it is assumed that the out of plane stresses are zero. The out of plane strains are not ...
1
vote
0answers
59 views

Sliding mode surface derivative exact matrix format (ds/dx)

Is there any suggestion on sliding surface state derivative matrix (ds/dx) format? ...
1
vote
0answers
43 views

Verify parameters variation in local sensitivity coefficients for microkinetic model

I’m performing a sensitivity analysis for a model of a set of elementary chemical reactions. I’m considering the variation in species concentration ($c_{k}$) with respect to the forward reaction rate ...
0
votes
2answers
53 views

What is the difference between Frictional non-linearity and Geometric non-linearity in FEA? (In terms of Stiffness matrix)

For geometric non-linearity, we know that as the geometry changes during the analysis, the global stiffness matrix changes since it already is a function of geometry of the structure. What I want to ...
0
votes
0answers
22 views

Perfectly Matched Layer (PML) for Plates in COMSOL

I am modelling infinite plate domain in COMSOL for frequency simulation. I need to apply absorbing boundary conditions all around the boundaries and I could not find PML option for plates. Could ...
0
votes
0answers
10 views

How are eigenvalues of a cell in the finite volume discretization calculated? (for euler fluid equations)

I am building a program to numerically solve the euler fluid equations based on finite volume discretization and am having trouble on a step where the eigenvalues of a cell need to be calculated ...
0
votes
0answers
23 views

Defining appropriate test function spaces for the finite element solution of Euler's fluid equations

I have the following coupled equations for the conservation of mass and momentum of a compressible fluid \begin{equation} \rho_t + (\rho u)_z = 0, \end{equation} $$ (\rho u)_t + (\rho u^2)_z + \...
0
votes
0answers
43 views

Simufact / MARC and Fortran

I am using Simufact from MSC Software to perform numerical simulations. To be able to create my own sub-routines I need to have Fortran compiler installed. Simufact manual mentions only Intel Fortran +...
0
votes
0answers
40 views

How to code this iterative scheme for solving 2-group diffusion equation?

I am trying to write a python code to solve the neutron diffusion equation to model neutron flux distribution in a one-dimensional two-group setting. The governing equations of the system are: $$-...