edgeIC

Initial conditions on geometry edge

Since R2023a

Description

An edgeIC object specifies the type of initial condition on an edge of a geometry. An femodel object contains an array of edgeIC objects in its EdgeIC property.

For separate assignments to a geometric region (or the entire geometry) and the boundaries of that region, the solver uses the specified assignment on the region and chooses the assignment on the boundary as follows. The solver gives an EdgeIC assignment precedence over a FaceIC assignment, even if you specify a FaceIC assignment after an EdgeIC assignment. The assignments in order of precedence are VertexIC (highest precedence), EdgeIC, FaceIC, and CellIC (lowest precedence).

Creation

Syntax

model.EdgeIC(EdgeID) = edgeIC(Name=Value)

model.EdgeIC = edgeIC(Name=Value)

Description

model.EdgeIC(EdgeID) = edgeIC(Name=Value) creates an edgeIC object and sets properties using one or more name-value arguments. This syntax assigns the specified structural or thermal initial condition to the specified edges of the geometry stored in the femodel object model. For example, model.EdgeIC([1 2]) = edgeIC(Temperature=25) specifies the temperature on edges 1 and 2.

model.EdgeIC = edgeIC(Name=Value) assigns the specified initial condition to all edges of the geometry. For example, model.EdgeIC = edgeIC(Temperature=25) specifies the initial temperature on all edges.

Input Arguments

EdgeID — Edge IDs

vector of positive integers

Edge IDs, specified as a vector of positive integers. Find the edge IDs using pdegplot with the EdgeLabels value set to "on".

Data Types: double

Properties

Displacement — Initial displacement

numeric vector | function handle

Initial displacement, specified as a numeric vector or function handle. A numeric vector must contain two elements for a 2-D model and three elements for a 3-D model. The elements represent the components of initial displacement.

Use a function handle to specify spatially varying initial displacement. The function must return a two-row matrix for a 2-D model and a three-row matrix for a 3-D model. Each column of the matrix corresponds to the initial displacement at the coordinates provided by the solver. For details, see Nonconstant Parameters of Finite Element Model.

The default initial displacement is zero.

Velocity — Initial velocity

numeric vector | function handle

Initial velocity, specified as a numeric vector or function handle. A numeric vector must contain two elements for a 2-D model and three elements for a 3-D model. The elements represent the components of initial velocity.

Use a function handle to specify spatially varying initial velocity. The function must return a two-row matrix for a 2-D model and a three-row matrix for a 3-D model. Each column of the matrix corresponds to the initial velocity at the coordinates provided by the solver. For details, see Nonconstant Parameters of Finite Element Model.

The default initial velocity is zero.

Temperature — Initial temperature or initial guess for temperature

real number | function handle

Initial temperature or initial guess for temperature, specified as a real number or function handle. Use a function handle to specify spatially varying initial temperature. For details, see Nonconstant Parameters of Finite Element Model.

Examples

Initial Displacement

Specify initial displacement for an femodel object representing a transient structural problem.

Create an femodel object for solving a transient structural problem, and assign the unit square geometry to the model.

model = femodel(AnalysisType="structuralTransient", ...
                Geometry=@squareg);

Plot the geometry with the edge labels.

pdegplot(model.Geometry,EdgeLabels="on");
xlim([-1.1 1.1])
ylim([-1.1 1.1])

The default initial displacement on the entire geometry is zero. Specify a small nonzero initial displacement on edge 1.

model.EdgeIC(1) = edgeIC(Displacement=[0 10^(-5)]);
model.EdgeIC

ans = 
  1x4 edgeIC array

Properties for analysis type: structuralTransient

Index     Displacement     Velocity
  1      [0 1.0000e-05]       []   
  2            []             []   
  3            []             []   
  4            []             []   

  Show all properties

Initial Temperature

Specify initial temperature for an femodel object representing a transient thermal problem.

Create an femodel object for solving a transient thermal problem, and assign the unit circle geometry to the model.

model = femodel(AnalysisType="thermalTransient", ...
                Geometry=@circleg);

Plot the geometry with the edge labels.

pdegplot(model.Geometry,EdgeLabels="on");

Because there is no default initial temperature, you must first specify the initial temperature for the entire geometry.

model.FaceIC = faceIC(Temperature=0);
model.FaceIC

ans = 
  1x1 faceIC array

Properties for analysis type: thermalTransient

Index    Temperature
  1           0     

  Show all properties

Specify a small nonzero initial temperature on all edges of the circle.

model.EdgeIC = edgeIC(Temperature=0.1);
model.EdgeIC

ans = 
  1x4 edgeIC array

Properties for analysis type: thermalTransient

Index    Temperature
  1        0.1000   
  2        0.1000   
  3        0.1000   
  4        0.1000   

  Show all properties

Version History

Introduced in R2023a